• Journal of computational fluids engineering
• /
• v.3 no.2
• /
• pp.39-51
• /
• 1998

This paper describes the analysis of flow field around an automobile. The governing equations of the 3-D unsteady incompressible Navier-Stokes equations are solved by the iterative time marching scheme. The Chimera grid technique has been applied to efficiently simulate the flow around the side-view mirror. To validate the capability of simulating the flow around a ground vehicle, the flows around the Ahmed body with 12.5$^{\circ}$ and 30$^{\circ}$ of slant angles are simulated and good agreements with experiment and other numerical results are achieved. To validate Chimera grid technique, the flow field around a cylinder was also calculated. The computed results are also well agreed with other numerical results and experiment. After code validations, the flow phenomena around the ground vehicle are evidently shown. The flow around the side-view mirror is also well simulated using the Chimera grid technique.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.4 no.4
• /
• pp.59-69
• /
• 2000

A numerical analysis based on two-dimensional, incompressible and compressible Navier-Stokes equations was carried out for double circular arc compressor cascade and the results are compared with available experimental data. The incompressible code based on SIMPLE algorithm adopts pressure weighted method and hybrid scheme for the convective terms. The compressible code with preconditioning method involves a upwind-biased scheme for the convective terms and LU-SGS scheme for temporal integration. Several turbulence models are evaluated by computing the turbulent viscous flows; Baldwin-Lomax, standard $\kappa$ -$\varepsilon$, $\kappa$ -$\varepsilon$ Lam. Bremhorst, standard $\kappa$-$\omega$, $\kappa$ -$\omega$ SST model.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.1
• /
• pp.1-8
• /
• 2001

The computations of the turbulent flow around the ship models with the free-surface effects were carried out. Incompressible Reynolds-Averaged Navier-Stokes equations were solved by using an explicit finite-difference method with the nonstaggered grid system. The method employed second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration. For the turbulence closure, a modified Baldwin-Lomax model was exploited. The location of the free surface was determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and a free-surface conforming grid was generated at each time step so that one of the grid boundary surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition was applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method developed in the present study, the computations were carried out for beth Wigley and Series 60 $C_B=0.6$ ship model and the computational results showed good agreements with the experimental data.

• 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 the Korean Society for Industrial and Applied Mathematics
• /
• v.3 no.2
• /
• pp.161-171
• /
• 1999

Aerodynamic sensitivity analysis codes are developed via the hand-differentiation using a direct differentiation method and an adjoint method respectively from discrete two-dimensional compressible Navier-Stokes equations. Unlike previous other researches, Baldwin-Lomax algebraic turbulence model is also differentiated by hand to obtain design sensitivities with respect to design variables of interest in turbulent flows. Discrete direct sensitivity equations and adjoint equations are efficiently solved by the same time integration scheme adopted in the flow solver routine. The required memory for the adjoint sensitivity code is greatly reduced at the cost of the computational time by allowing the large banded flux jacobian matrix unassembled. Direct sensitivity code results are found to be exactly coincident with sensitivity derivatives obtained by the finite difference. Adjoint code results of a turbulent flow case show slight deviations from the exact results due to the limitation of the algebraic turbulence model in implementing the adjoint formulation. However, current adjoint sensitivity code yields much more accurate sensitivity derivatives than the adjoint code with the turbulence eddy viscosity being kept constant, which is a usual assumption for the prior researches.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.28-40
• /
• 2002

Turbulent flows around two-dimensional wing sections in ground effect are analysed by incompressible RANS equations and a finite difference method. The Baldwin-Lomax algebraic turbulence model is used to simulate high Reynolds number flows. The main purpose of this study is to clarify the two-dimensional ground effect and its flow characteristics due to different ground boundary conditions, i.e., moving and fixed bottom boundary. As a first step, to validate the present numerical code, the computational result of Clark-Y(t/C 11.7%) is compared with published numerical results and experimental data. Then, NACA4412 section in ground effect is calculated for various ground clearances with two bottom boundary conditions. According to the computational results, the difference in the lift and moment simulated with the two bottom boundary conditions is negligible, but the drag force simulated by the fixed bottom is to some extent smaller than that by the moving bottom. Therefore, it can be concluded that the drag force measured in a wind tunnel with the fixed bottom could be smaller than that with the moving bottom.

• Journal of Ship and Ocean Technology
• /
• v.4 no.1
• /
• pp.11-20
• /
• 2000

Reynolds-averaged Navier-Stokes equations are numerically solved using a secondorder finite difference method for the analysis of turbulent flows around single and tandem hydrofoils advancing under the free surface. The location of the free surface, not known a priori, is computed from the kinematic free surface condition and the computational grid is conformed at each iteration to the free surface deformation. The eddy viscosity model of Baldwin-Lomax is employed for the turbulence closure. The method is validated through the comparision of the numerical results with the experimental data for a single hydrofoil of a Joukowski foil section. A computational study is also carried out to investigate the effect of the submergence depth and the Froude number on the lift and the drag of the hydrofoil. For tandem hydrofoils, computations are performed for several separation distances between the forward and aft foils to see the interference effect. The result shows clearly how the lift and drag change with the separation distance.

• 한국전산유체공학회:학술대회논문집
• /
• 1996.05a
• /
• pp.46-51
• /
• 1996

The flow field around a three dimensional minivan-like body has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system using second-order accurate schemes for the time derivatives, and third/second-order scheme for the spatial derivatives. The Marker-and-Cell concept is applied to efficiently solve continuity equation. The fourth -order artificial damping is added to the continuity equation for numerical stability. A H-H type multi-block grid system is generated around a three dimensional minivan-like body. Turbulent flows have been modeled by the Baldwin-Lomax turbulent model. The simulation shows three dimensional vortex-pair just behind body. And the flow separation is also observed the rear of the body. It has concluded that the results of present study properly agree with physical flow phenomena.

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

The numerical procedure has been developed for simulating incompressible viscous flow around a turbine stage with rotor-stator interaction. This study solves 2-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system. The Marker-and-Cell concept is applied to efficiently solve continuity equation. To impose an accurate boundary condition, O-H multiblocked grid system is generated. O-type grid and H-type grid is generated near and outer rotor-stator The cubic-spline interpolation is applied to handle a relative motion of a rotor to the stator. Turbulent flows have been modeled by the Baldwin- Lomax turbulent model. To validate present procedure, the time averaged pressure coefficients around the rotor and stator are compared with experiment and a good agreement obtained.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.10a
• /
• pp.15-24
• /
• 2000

The transient response of an airfoil to a rapidly deploying spoiler is numerically investigated using the turbulent compressible Navier-Stokes equations in two dimensions. Algebraic Baldwin-Lomax model, Wilcox $\kappa-\omega$ model, and SST $\kappa-\omega$ turbulence model are used to calculate the unsteady separated flow due to the rapid spoiler deployment. The spoiler motion relative to a stationary airfoil is treated by an overset grid hounded by a Dynamic Domain-Dividing Line which has been devised by the authors. The adverse effects of the spoiler influenced by the spoiler location and the hinge gap are expounded. The numerical results are in reasonably good agreement with the existing experimental data.