• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.53-61
• /
• 1999

Two-dimensional, unsteady, incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. The compressible code involves a conventional upwind-differenced scheme for the convective terms and LU-SGS scheme for temporal integration. The incompressible code with pseudo-compressibility method also adopts the same schemes as the compressible code. Three two-equation turbulence models are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by predicting the flow around the RAE 2822 transonic airfoil and the NACA 4412 airfoil, respectively. In addition, both the incompressible and compressible code are used to compute the flow over the NLR 7301 airfoil with flap to study the compressible effect near the high-loaded leading edge. The grid systems are efficiently generated using Chimera overlapping grid scheme. Overall, the κ-ω SST model shows closer agreement with experiment results, especially in the prediction of adverse pressure gradient region on the suction surfaces of high-lift airfoils.

• Korea-Australia Rheology Journal
• /
• v.16 no.3
• /
• pp.117-128
• /
• 2004

The effect of molecular parameters on the steady vortex behaviors in the inertial viscoelastic flow past a cylinder has been investigated. FENE-CR model was considered as a constitutive equation. A recently developed iterative solution method (Kim et al., (in press)) was found to be successfully applicable to the computation of inertial viscoelastic flows. The high-resolution computations were carried out to understand the detailed flow behaviors based on the efficient iterative solution method armed with ILU(0) type pre-conditioner and BiCGSTAB method. The discrete elastic viscous split stress-G/streamline upwind Petrov Galerkin (DEVSS-G/SUPG) formulation was adopted as a stabilization method. The vortex size decreased as elasticity increases. However, the vortex enhancement was also observed in the case of large extensibility, which means that the vortex behavior is strongly dependent upon the material parameters. The longitudinal gradient of normal stress was found to retard the formation of vortex, whereas the extensional viscosity played a role in the vortex enhancement. The present results are expected to be helpful for understanding the inertial vortex dynamics of viscoelastic fluids in the flow past a confined cylinder.

• Journal of computational fluids engineering
• /
• v.8 no.4
• /
• pp.16-25
• /
• 2003

The design method for transonic turbine blades has been developed based on Wavier-Stokes equations. The present computing process is done on the four separate steps, i.e., determination of the blade profile, generation of the computational grids, cascade flow simulation and analysis of the computed results in the sense of the aerodynamic performance. The blade shapes are designed using the cubic polynomials under the control of the design parameters. Numerical methods for the flow equations are based on Van-Leer's FVS with an upwind TVD scheme on the finite volume. In the present study, numerical simulation has been done to investigate the effects of the design parameters on the aerodynamic peformance of the axial-flow turbine blades. Applications are made to the VKI transonic rotor blades. Computed results are analyzed with respect to four parameters and compared with the experimental data.

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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.1
• /
• pp.43-57
• /
• 1999

The fluid flow and heat transfer in a thin liquid film are investigated numerically. The flow Is assumed to be two-dimensional laminar and surface tension is considered. The most important characteristics of this flow is the existence of a hydraulic jump through which the flow undergoes very sharp and discontinuous change. Arbitrary Lagrangian-Eulerian(ALE) method is used to describe moving free boundary and a modified SIMPLE algorithm based on streamline upwind Petrov-Galerkin(SUPG) finite element method is used for time marching iterative solution. The numerical results obtained by solving unsteady full Navier-Stokes equations are presented for planar and radial flows subject to constant wall temperature or constant wall heat flux, and compared with available experimental data. It Is discussed systematically how the inlet Reynolds and Froude numbers and surface tension affect the formation of a hydraulic jump. In particular, the effect of temperature dependent fluid properties is also discussed.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.9
• /
• pp.1744-1759
• /
• 1992

A multidimensional simulation of turbulent flow and combustion with swirl in the cylinder of SI engine is implemented to clarify the effects of swirl. present simulation employs the ICED-ALE numerical technique, the skew-upwind difference scheme, a modified k_.epsilon. turbulence model, a combustion model of Arrhenius type and turbulence-mixing-control type. First, the calculations for swirling flow in an axisymmetric cylinder are carried out. The results are compared with the experimental data to validate the numerical analysis. Second, the calculations for intake, compression and combustion processes in an axisymmetric cylinder are performed. The effects of swirl on turbulent flow and combustion are examined through the parametric study of swirl number 0.0, 0.6, 1.2 and 2.4. As a result, it is numerically shown that the turbulent kinetic energy and the swirl velocity, which are produced during the intake process, affect the combustion process.

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

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2857-2861
• /
• 2007

Flow and aerodynamic characteristics were analyzed numerically for a commercial passenger airplane, Boeing 747-400, flying in the cruising condition. The model geometry with 100:1 in scale was obtained by the photo scanning measurement with the maximum error of 1.4% comparing with the real airplane dimension. The three-dimensional inviscid steady compressible governing equations were solved by the finite volume method in the unstructured grid system. The convective terms were treated by the Crank-Nicholson and first-order upwind schemes. In the computational results, the strong wing-tip vortices were clearly observed and the pressure contours on the airplane surface were suggested. The lift and drag forces in the wing with engines increase by 1.49% and 3.9%, respectively compared with the case without engines. The aerodynamic forces were estimated quantitatively for each element which consists of the airplane.

• Wind and Structures
• /
• v.9 no.6
• /
• pp.473-492
• /
• 2006

This paper presents numerical results of pedestrian-level wind environment around the base of a row of tall buildings by CFD. Four configurations of building arrangement are computed including a single square tall building. Computed results of pedestrian-level wind flow patterns and wind speeds are compared to previous wind tunnel measurement data to enable an assessment of CFD predictions. The CFD model uses the finite-volume method with RNG $k-{\varepsilon}$ model for turbulence closure. It is found that the numerical results can reproduce key features of pedestrian-level wind environment such as corner streams around corners of upwind building, sheltered zones behind buildings and channeled high-speed flow through a building gap. However, there are some differences between CFD results and wind tunnel data in the wind speed distribution and locations of highest wind speeds inside the corner streams. In locations of high ground-level wind speeds, CFD values match wind tunnel data within ${\pm}10%$.

• The KSFM Journal of Fluid Machinery
• /
• v.6 no.1 s.18
• /
• pp.23-29
• /
• 2003

Three-dimensional flow through a tubular centrifugal fan with airfoil type blades is analyzed, and the effects of turbulence model and numerical scheme on the results are investigated. Standard $k-{\epsilon}$ model and k - w model are tested as turbulence closures. The numerical schemes for convection terms, i.e., Upwind Differencing Scheme (UDS), Mass Weighted Skewed upstream differencing scheme (MWS), Linear Profile Skewed upstream differencing scheme (LPS), and Modified Linear Profile Skewed upstream differencing scheme (MLPS) are also tested, and the performances of these schemes coupled with two turbulence models are evaluated. The static pressure distributions are compared with experimental data obtained in this work, which shows that the $k-{\epsilon}$ model gives better results than the k-w model.