• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.90-95
• /
• 1998

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.11
• /
• pp.1371-1378
• /
• 1999

Flow through turbomachinery has a very complex structure and Is Intrinsically unsteady. In addition, trend to highly loaded turbomachinery makes the flow extremely complex due to the interaction between rotor and stator. In this study, flows through UTRC LSRR turbine are numerically analyzed using 2 dimensional Navier-Stokes equations. The convective terms of the governing equations are discretized using the Van-Leer's FVS(Flux vector splitting) with an upwind TVD scheme. The conventional central differencing is used to discretize the diffusion terms on the finite volume. The accurate unsteady motion is achieved by using a 2nd order accurate, 3-point Euler implicit scheme. The quasi-conservative zonal scheme is used for calculating the flow variables on the zonal interface between the rotor and stator. The axial gap between stator and rotor has been configured in two variations, 15% and 65% of average chord length. The analysis program is validated using experimental results and the effect of axial gap is examined. The numerical analysis results are presented by time averaged pressure coefficient and pressure magnitude coefficient and compared with experimental results.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2005.06a
• /
• pp.210-215
• /
• 2005

A three-dimensional numerical simulation is performed to investigate on a low Reynolds number mixed convection in a horizontal rectangular channel with the upper part cooled and the lower part heated uniformly. The three-dimensional governing equations are solved using a finite volume method. For convective term, the central differencing scheme is used and for the pressure correction, the PISO algorithm is used. Solutions are obtained for A=4, Pr=0.72, 10, 909, the Reynolds number ranging from $2.1{\times}10^{-2}$ to $1.2{\times}10^{-1}$, the Rayleigh number is $3.5{\times}10^4$. It is found that vortex roll structures of mixed convection in horizontal rectangular channel can be classified into three roll structures which affected by Prandtl number and Reynolds number.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.7
• /
• pp.873-881
• /
• 1997

A numerical simulation has been carried out for three-dimensional turbulent flows around an Ahmed body. The Reynolds-averaged Navier-Stokes equation is solved with the SIMPLE method in general curvilinear coordinates system. Several k-.epsilon. turbulence models with two convective difference schemes are evaluated for the performance such as drag coefficient, velocity and pressure fields. The drag coefficient, the velocity and pressure fields are found to be changed considerably with the adopted k-.epsilon. turbulence models as well as the finite difference schemes. The results of simulation prove that the RNG k-.epsilon. model with the QUICK scheme predicts fairly well the tendency of velocity and pressure fields and gives more reliable drag coefficient. It is also demonstrated that the large difference between simulations and experiment in the drag coefficient is due to relatively high predicted values of pressure drag from vertical rear end base.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.1
• /
• pp.1-14
• /
• 1998

A two-dimensional numerical model for tidal currents based on the depth averaged equation is developed. The mode1 employs a rectangular grid system for its simplicity in the application of complicate coastal shore lines. To raise computing efficiency, implicit approximate factorization scheme is implemented in solving governing equations. An upwind-differencing is used to discretize convective terms, which provides a numerical dissipation automatically and suppresses any oscillations caused by nonlinear instabilities. Some numerical tests are made against the analytic solutions of a linearized shallow water equation to validate the developed numerical scheme, and comparisons of the model prediction with the analytic solution are satisfactory. As a real application, the tidal currents are computed on the Inchon area where the tidal currents are important for the design of new canal which is under construction.

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

The three-dimensional turbulent cascade flows with and without endwall fences are numerically investigated by solving the incompressible Navier-Stokes equations with a high-Reynolds number $k-{\varepsilon}$ turbulence closure model. A projection method based algorithm is used in the finite-volume formulation, with the second order upwind-differencing scheme for the convective terms. First, assessments on accuracy of the present method are made by comparing the static pressure distributions at the mid-span of the cascade with measured data, and also by confirming the experimental observations on the choice of an optimal fence height for the secondary flow control. In understanding the three-dimensional nature of the secondary flow in turbine cascade, the limiting streamline patterns and the static pressure contours at the suction surface of the blade as well as on the cascade endwall are employed to visualize the effectiveness of the endwall fence for the secondary flow control. Analysis on the streamwise vorticity contour maps along the cascade with the three-dimensional representation of their iso-surfaces reveals the strucuture of the complicated vortical flow in the turbine cascade with endwall fence, and also leads to an understanding on formation of the counter-rotating streamwise vortex over the endwall fence, in explaining the mechanisms of controlling the secondary flow and also for the proper selection of an optimal fence height.