• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.53-59
• /
• 2004

A general purpose program NUFLEX for the analysis of 3-D heat/fluid flow in complex geometry with pre/post processor have been developed, which consists of a flow solver based on FVM and a dedicated pre/post processor. The program employs a general non-orthogonal grid system and solve laminar and turbulent (lows with standard and RNG $\kappa-\epsilon$ turbulence models. NUFLEX is capable of analysing two-phase flow with topologically complex interface, turbulent diffusion combustion, solidification problems and magnetic flow. For the purpose of verification of the program and testing the applicability, several practical problems are solved and compared with the available data. Comparison of the NUFLEX results with that by the STAR-CD program has been also made for the same flow configuration and grid structure.

• Journal of applied mathematics & informatics
• /
• v.7 no.3
• /
• pp.905-915
• /
• 2000

In this article we discuss some aspects of operational Tau Method on delay differential equations and then we apply this method on the differential delay equation defined by $\omega(u)\;=\frac{1}{u}\;for\;1\lequ\leq2$ and $(u\omega(u))'\;=\omega(u-1)\;foru\geq2$, which was introduced by Buchstab. As Khajah et al.[1] applied the Recursive Tau Method on this problem, they had to apply that Method under the Mathematica software to get reasonable accuracy. We present very good results obtained just by applying the Operational Tau Method using a Fortran code. The results show that we can obtain as much accuracy as is allowed by the Fortran compiler and the machine-limitations. The easy applications and reported results concerning the Operational Tau are again confirming the numerical capabilities of this Method to handle problems in different applications.

• International Journal of Automotive Technology
• /
• v.7 no.3
• /
• pp.283-288
• /
• 2006

A mathematical methodology is proposed for designing nozzle hole shapes producing controlled geometric cavitation. The proposed methodology uses an unstructured RANS flow solver, with the ability to compute sensitivity derivatives via an adjoint algorithm. The adjoint formulation for the N-S equations is presented while variation of the turbulence viscosity is not taken into account during the geometry modifications. The sensitivities are calculated in a mode independently of the shape parameterisation. The method is used to develop and evaluate conceptual shapes for nozzle hole cavitation reduction. The localized region at the hole inlet producing cavitation, is parameterised using its radius of curvature, while a cost function is formulated to eliminate the negative pressures present at this location. Sensitivity derivatives are used to assess the dependence of the localized region on the minimum pressure, and to drive the geometry to the targeted shape. The results show that the computer model can provide nozzle hole entry shapes that produce predefined flow characteristics, and thus can be used as an inverse design tool for nozzle hole cavitation control.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.66-71
• /
• 2013

In this paper, the supersonic flows around space launch vehicles have been numerically simulated by using a 3-D RANS flow solver. The focus of the study was made for investigating plume-induced flow separation(PIFS). For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras model was employed for the closure of turbulence. The Gauss-Seidel iteration was used for time integration. To validate the flow solver, calculation was made for the 0.04 scale model of the Saturn-5 launch vehicle at the supersonic flow condition without exhaust plume, and the predicted results were compared with the experimental data. Good agreements were obtained between the present results and the experiment for the surface pressure coefficient and the Mach number distribution inside the boundary layer. Additional calculations were made for the real scale of the Saturn-5 configuration with exhaust plume. The flow characteristics were analyzed, and the PIFS distances were validated by comparing with the flight data. The KSLV-1 is also simulated at the several altitude conditions. In case of the KSLV-1, PIFS was not observed at all conditions, and it is expected that PIFS is affected by the nozzle position.

• Journal of computational fluids engineering
• /
• v.8 no.3
• /
• pp.12-20
• /
• 2003

A numerical procedure for the calculation of solidification and melting phase change using PISO algorithm is presented. In case of phase change problem, the coupling between velocity/pressure/temperature and liquid fraction is important. The converged temperature and liquid fraction solution which satisfies the energy balance is acquired by applying enthalpy method into inner iteration in matrix solver. And a modified PISO algorithm version is introduced to properly solve the coupling between velocity/pressure/temperature and liquid fraction. A comparison of the proposed procedure with a standard iterative method shows improvement both in terms of computing speed and robustness.

• Journal of computational fluids engineering
• /
• v.5 no.3
• /
• pp.23-31
• /
• 2000

The three-dimensional unsteady compressible Full Navier-Stokes equation solver with sliding multi-block method has been applied to analyze three dimensional characteristics of the viscous flow field and compression wave around the high speed train which is entering into a tunnel. The numerical scheme of AF + ADI was used to efficiently solve Navier-Stokes equations in the curvilinear coordinate system. The vortex formation owing to the viscous interaction around the train was found and the generation of compression wave due to the blockage effects was observed ahead of the train in the form of plane wave. The three dimensional characteristics of the flow field compared to the analytic results were discussed in detail. The variation of pressure of tunnel wall surface and velocity profile of the train are identified as the train enters into a tunnel. The changes in aerodynamic forces and streamlines of each specific sections are also discussed.

• Journal of Ocean Engineering and Technology
• /
• v.7 no.2
• /
• pp.150-161
• /
• 1993

원초변수를 이용한 Navier-Stokes 방정식의 수치계산기법을 개발하고, 이를 응용하여 backward facing step의 층류 유동을 계산하였다. 직교좌표계에서의 비압축성 Navier-Stokes방정식을 풀기위해 시간과 공간항을 2차 정도의 유한 차분을 사용하여 이산화하였고 비교차격자계를 사용하여 양해법으로 수치 계산하였다. 운동량방정식과 연속방정식으로 부터 유도된 압력방정식(pressure-poisson equation)을 이용하여 무발산 조건을 만족시켰ㄲ다. Backward facing step의 층류 유동을 100.$\leq$R$_e$$\leq$1000 범위에 대해서 수치 계산하였으며 실험결과와 잘 일치하는 결과를 구할 수 있었다. 특히 step뒤에서 생기는 박리구간의 길이는 다른 계산결과들보다 실험치에 가까운 값을 얻을 수 있었으며, Re가 600보다 클때는 위쪽 벽에 또 다른 박리 유동이 발생되는 현상이 예측되었다.

• 연구논문집
• /
• s.28
• /
• pp.39-47
• /
• 1998

Industrial electrostatic precipitation is a very complex process, which involves multiple-way interaction between the electric field, the fluid flow, and the particulate motion. This paper describes a strongly coupled calculation procedure for the rigorous computation of particle dynamics during electrostatic precipitation. The turbulent gas flow and the particle motion under electrostatic forces are calculated by using the commercial computational fluid dynamics (CFD) package FLUENT linked to a finite-volume solver for the electric field and ion charge. Particle charge is determined from both local electrical conditions and the cell residence time which the particle has experienced through its path. Particle charge density and the particle velocity are averaged in a control volume to use Lagrangian information of the particle motion in calculating the gas and electric fields. The turbulent particulate transport and the effects of particulate space charge on the electrical current flow are investigated. The calculated results for poly-dispersed particles are compared with those for mono-dispersed particles, and significant differences are demonstrated.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.142-146
• /
• 2017

지진해일의 수리특성을 분석하기 위하여 유사한 파형특성을 가진 고립파를 많이 이용하고 있다. 그러나 고립파의 근사파형은 실제 지진해일에 비해 상당히 좁은 파형분포를 가지고 있다. 이에 수치모의에서는 기존의 고립파 근사식을 개량하여 고립파형의 지진해일을 수치적으로 생성하고 있다. 본 연구에서는 지진해일의 파형분포에 따른 월파특성을 수치적으로 조사하기 위하여 개량된 고립파 근사식을 2차원 N-S solver에서 적용하였다. 이것에 기초하여 수치파동수조에 직립호안과 그 배후에는 월파수조를 설치하고, 지진해일 월파량을 측정하였다. 수치해석결과로부터 직립호안 주변의 공간파형과 마루 위의 유속분포로부터 파형분포에 따른 월파현상을 분석할 수 있었다. 또한 기존 고립파 근사이론 대비 개량된 고립파의 체적비에 따른 월파량 변화를 정량적으로 조사하였다. 그 결과 지진해일의 체적비가 증가할수록 월파량이 거의 선형적으로 증가하는 경향을 나타내었다.

• Journal of computational fluids engineering
• /
• v.20 no.1
• /
• pp.47-56
• /
• 2015

In this paper, we study the effect of various turbulence models by comparing the aerodynamic characteristics and the flow patterns computed for aircraft models. An in-house CFD solver, MSAPv, that solves the three dimensional RANS equations with the turbulence model equations is used. The turbulence models used in this study are the Spalart-Allmaras model, Menter's $k-{\omega}$ SST model, Coakley's $q-{\omega}$ model, and Huang and Coakley's $k-{\varepsilon}$ model. DLR-F6 WB and WBNP configurations are selected for the study. We concentrate on the separated flow pattern variations with the turbulence models at the wing-body junction and the wing-pylon junction as well as drag polar curves.