• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.64-70
• /
• 1996

Three dimensional turbulent flow fields around ships are simulated by a numerical method. Reynolds Averaged Navier-Stokes equations are used where Reynolds stresses are approximated by Baldwin-Lomax and Sub-Grid Scale(SGS) turbulence models. Body-fitted coordinate system is introduced to conform three dimensional ship geometries. The governing equations are discretized by a finite volume method. Temporal derivatives are approximated by the forward differencing and the convection terms are approximated by the QUICK or Kawamura scheme. The 2nd-order centered differencing is used for other spatial derivatives. Pressure and velocity fields are simultaneously iterated by the Highly Simplified Marker-And-Cell method. To verify the numerical method and turbulence models, flow fields around ships are simulated and compared to the experiments.

• 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.2 no.2
• /
• pp.35-43
• /
• 1997

An efficient inverse method for 1.he supersonic/hypersonic axisymmetric body design is developed for the parabolized Navier-Stokes equations. The developed method is examined numerically for three extreme testcases in the supersonic(M/sub ∞/=3.0) and hypersonic(M/sub ∞/=6.28) speeds. The first one is a negative pressure distribution near a vacuum pressure and the second one is a positive pressure distribution over the whole region of the body. The last one is the case of abrupt change of pressure distribution to zero in the forward region of the body. These testcases show the robustness of the method. By introducing a regular-falsi method and by using a not-fully converged inverse solution, the convergence behavior was greatly improved.

• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.13-16
• /
• 2008

After the collapse of the Tacoma bay bridge at Tacoma Washington, the accurate prediction of aerodynamics became crucial to the sound design of bridges. CFD(Computational Fluid Dynamics) becomes important tool for the prediction on wind effects on the bridge due to the recent development of CFD. The usage of CFD is further prompted by the advantages in using CFD, such as low-cost and fast feed-back of design. In this paper, an unsteady compressible Reynolds averaged Navier-Stokes code is used for the computation of the flow over bridges. Coakley's ��q-${\omega}$ �� two-equation turbulence model is used for the turbulent eddy viscosity. For accurate and stable computations, the local preconditioning method is adapted to the code. Aerodynamic characteristics of a couple bridges are presented to show the validity and the accuracy of the method.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.35-42
• /
• 1998

Time marching algorithms applied to compressible Navier-Stokes equation have a convergence problem at low Mach number. It is mainly due to the eigenvalue stiffness and pressure singularity as Mach number approaches to zero. Among the several methods to overcome the shortcomings of time marching scheme, time derivative preconditioning method have been used successfully. In this numerical analysis, we adopted a preconditioner of K.H. Chen and developed a two-dimensional, axisymmetric Navier-Stokes program. The steady state driven cavity flow and backward facing step flow problems were computed to confirm the accuracy and the robustness of preconditioned algorithm for low Mach number flows. And the transonic and supersonic flows insice the JPL axisymmetric nozzle internal flow is exampled to investigate the effects of preconditioning at high Mach number flow regime. Test results showed excellent agreement with the experimental data.

• 한국전산유체공학회:학술대회논문집
• /
• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.182-187
• /
• 2003

Numerical analysis of three-dimensional vicous flow is used to compute the design speed operating line of a transonic axial-flow compressor. The Navier-Stokes equation was solved by an explicit finite-difference numerical scheme and the Baldwin-Lomax turbulence model was applied. A spatially-varying time-step and an implicit residual smoothing were used to improve convergence. Two-stage axial compressor of a turboshaft engine developed KARI was chosen for the analysis. Numerical results show reasonably good agreements with experimental measurements made by KARI. Numerical solutions indicate that there exist a strong shock-boundary layer interaction and a subsequent large flow separation. It is also observed that the shock is moved ahead of the blade passage at near-stall condition.