• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.112-122
• /
• 1996

A finite element scheme using the concept of PISO method has been developed to solve the Navier-Stokes viscous flows in all speed range. This scheme includes development of new pressure equation that retains both the hyperbolic term related with the density variation and the elliptic term reflecting the incompressibility constraint. The present method is applied to the incompressible two-dimensional driven cavity flow problems(Re=100, 400 and 1,000). For compressible flows, the Carter plate problem(M=3 and Re=1,000) is computed. Finally, we have simulated the shock-boundary layer interaction(M=2 and Re=2.96×10/sup 5/), a more difficult problem, and compared its results with the experiment to demonstrate the shock capturing capability of the present solution algorithm.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.123-141
• /
• 1996

In this paper the upwind flux difference splitting Navier-Stokes method has been applied to study the perfect gas and the equilibrium chemically reacting hypersonic flow over an axisymmetric sphere-cone(5°) geometry. The effective gamma(γ), enthalpy to internal energy ratio was used to couple chemistry with the fluid mechanics for equilibrium chemically reacting air. The test case condition was at altitude(30km) and Mach number(15). The equilibrium shock thickness over the blunt body region was much thinner than that of perfect gas shock. The pressure difference between perfect gas and equilibrium gas was about 3 ∼ 5 percent. The heat transfer coefficient were also calculated. The results were compared with VSL results in order to validate the current numerical analysis. The results from current method were compared well VSL results ; however, not well at near nose. The proper boundary condition and grid system will be studied to improve the solution quality.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.203-212
• /
• 1995

In this paper the upwind flux difference splitting Navier-Stokes method has been applied to study the perfect gas and the equilibrium chemically reacting hypersonic flow over an axisymmetric sphere-cone($5^{\circ}$) geometry. The effective gamma($\bar{r}$), enthalpy to internal energy ratio was used to couple chemistry with the fluid mechanics for equilibrium chemically reacting air. The test case condition was at altitude(30Km) and Mach number(15). The equilibrium shock thickness over the blunt body region was much thinner than that of perfect gas shock. The pressure difference between perfect gas and equilibrium gas was about $3\sim5$ percent. The skin friction coefficient and heat transfer coefficient were also calculated.

• Journal of computational fluids engineering
• /
• v.4 no.2
• /
• pp.67-73
• /
• 1999

In the present paper, three types of the flow solvers were used to investigate the influence on the airfoil shape optimization. The adopted equations, i.e., Euler, thin layer Navier-Stokes and full Navier-Stokes ones. are solved using implicit LU-ADI decomposition scheme. The gradient projection method with the sinusoidal function was used as an optimization algorithm. The present numerical method was applied to the drag minimization problems under the initial shape of NACA0012 airfoils.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.171-176
• /
• 1999

In the present paper, three types of the flow solvers were used to investigate the influence on the airfoil shape optimization. The adopted equations, i.e., Euler , thin layer Navier- Stokes and full Navier-Stokes ones, are solved using implicit LU-ADI decomposition scheme. The feasible direction algorithm with the sinusoidal function was used as an optimization algorithm. The present numerical method was applied to the drag minimization problems under the initial shape of NACA0012 airfoils.

• 유체기계공업학회:학술대회논문집
• /
• 1999.12a
• /
• pp.27-32
• /
• 1999

A numerical study based on the three-dimensional thin-layer Navier-Stokes solver is carried out to analyze the flowfield through a single stage transonic compressor. Explicit four-step Runge-Kutta scheme with spatially variable time step and implicit residual smoothing is used. The governing equations are discretized with exploit finite difference method. Mixed-out average method is used at the interface between rotor and stator. And, an artificial dissipation model is used to assure the stability of solution. The results with k-$\omega$ turbulence model were compared to the results with Baldwin-Lomax model, and physical phenomena of transonic compressor are presented. The two turbulence models give the results that show reasonably good agreements with experimental data.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.6
• /
• pp.1-8
• /
• 2004

The objective of this study is convergence acceleration of multigrid Navier-Stokes solvers. This study has been performed to enhance the performance of preconditioned multi-stage time stepping method which is a popular smoother for the multigrid Navier-Stokes solvers. Comparative study on the convergence characteristics of the ADI and DDADI preconditioners has been conducted. It is shown that the DDADI preconditioner has better performance than the ADI by numerical tests on the 2-D compressible turbulent flows past airfoils. The Spalart-Allmaras turbulent model and the Baldwin-Lomax turbulent model have been compared with the multigrid calculations.

• 유체기계공업학회:학술대회논문집
• /
• 2005.12a
• /
• pp.169-172
• /
• 2005

Numerical analysis of turbulent flow in three-dimensional channel with V-shaped ribs extruded on both walls has been carried out. Reynolds-averaged Navier-Stokes are calculated for analysis of fluid flow and heat transfer. Shear stress transport (SST) turbulence model is used as a turbulence closure. Computational results for heat transfer rate show good agreements with experimental data.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2002.10a
• /
• pp.212-216
• /
• 2002

평행한 벽 내부의 흐르는 규제 속에 원형실린더를 놓고 운동장 해석을 수행하였다. 비압축성 Navier- Stokes 방정식을 풀었고 3차 풍상미분의 수치해법을 이용하였다. 채널 내부에서 실린더의 위치를 이동하면서 벽면의 효과 전단력, 와류의 현장을 규명하였다. 수치처리는 Marker & Cell 기법에 의한 유한차분법을 사용하였다. 본 연구를 통하여 실린더와 벽 경계 사이에서의 생성된 와가 박리 전단에 영향을 미치는 것을 알 수 있었다.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.88-97
• /
• 1996

This paper treats an adaptive finite-element method for the viscous compressible flow governed by Navier-Stokes equations in two dimensions. The numerical algorithm is the two-step Taylor-Galerkin mettled using unstructured triangular grids. To increase accuracy and stability, combined moving node method and grid refinement method have been used for grid adaption. Validation of the present algorithm has been made by comparing the present computational results with the existing experimental data and other numerical solutions. Four benchmark problems are solved for demonstration of the present numerical approach. They include a subsonic flow over a flat plate, the Carter flat plate problem, a laminar shock-boundary layer interaction. and finally a laminar flow around NACA0012 airfoil at zero angle of attack and free stream Mach number of 0.85. The results indicates that the present adaptive triangular grid method is accurate and useful for laminar viscous flow calculations.