• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.4
• /
• pp.387-397
• /
• 2013

An assessment of two-equation turbulence models, the low Reynolds k-${\varepsilon}$ and k-${\omega}$ SST models, with the compressibility corrections proposed by Sarkar and Wilcox, has been performed. The compressibility models are evaluated by investigating transonic or supersonic flows, including the arc-bump, transonic diffuser, supersonic jet impingement, and unsteady supersonic diffuser. A unified implicit finite volume scheme, consisting of mass, momentum, and energy conservation equations, is used, and the results are compared with experimental data. The model accuracy is found to depend strongly on the flow separation behavior. An MPI (Message Passing Interface) parallel computing scheme is implemented.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.5 no.4
• /
• pp.357-368
• /
• 1993

With the assumption of the diffusion dominated flow, a numerical model has been developed for undertow and turbulence structure under the breaking wave by using the $textsc{k}$-$\varepsilon$ turbulence model. Undertow is a strong mean current which moves seqwards below the level of wave trough in the surf zone. The turbulence, generated by wave breaking in the roller, spreads and dissipates downwards. The governing equations are composed of the equation of motion with the period-averaged shear stress due to waves; $textsc{k}$- and $\varepsilon$-equations with the turbulence energy Production due to wave breaking. They are discretised by the three-level fully implicit scheme, which can be solved by using Thomas algorithm. The model gives good agreements with measurements except for the station that is closest to the breaking point.

• International Journal of Aeronautical and Space Sciences
• /
• v.13 no.3
• /
• pp.307-316
• /
• 2012

Numerical simulations of 3D aircraft configurations are performed in order to understand the effects of turbulence models on the prediction of aircraft's aerodynamic characteristics. An in-house CFD code that solves 3D RANS equations and two-equation turbulence model equations are used. The code applies Roe's approximated Riemann solver and an AF-ADI scheme. Van Leer's MUSCL extrapolation with van Albada's limiter is also adopted. Various versions of Menter's $k-{\omega}$ SST turbulence models as well as Coakley's $q-{\omega}$ model are incorporated into the CFD code. Menter's $k-{\omega}$ SST models include the standard model, the 2003 model, the model incorporating the vorticity source term, and the model containing controlled decay. Turbulent flows over a wing are simulated in order to validate the turbulence models contained in the CFD code. The results from these simulations are then compared with computational results from the $3^{rd}$ AIAA CFD Drag Prediction Workshop. Numerical simulations of the DLR-F6 wing-body and wing-body-nacelle-pylon configurations are conducted and compared with computational results of the $2^{nd}$ AIAA CFD Drag Prediction Workshop. Aerodynamic characteristics as well as flow features are scrutinized with respect to the turbulence models. The results obtained from each simulation incorporating Menter's $k-{\omega}$ SST turbulence model variations are compared with one another.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.5
• /
• pp.1158-1174
• /
• 1988

A numerical technique is developed for the solution of fully developed turbulent recirculating flow in the passage of variable area using the non-orthogonal coordinate transformation. In the numerical analysis, primitive pressure-velocity finite difference equations were solved by SIMPLER algorithm with 2-equation turbulence model and algebraic stress model (ASM). QUICK scheme on the differencing of convective terms which is free from the inaccuracies of numerical diffusion has been applied to the variable grids and the results compared with those from HYBRID scheme. In order to test the effect of streamline curvatures on turbulent diffusion Lee and Choi streamline curvature correction model which has been obtained by modifying the Leschziner and Rodi's model is testes. The ASM was also employed and the results are compared to those from another turbulence model. The results show that difference of convective differencing schemes and turbulence models give significant differences in the prediction of velocity fields in the expansion region and outlet region of the combustor, however show little differences in the parallel flow region.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.235-238
• /
• 2005

The unsteady supersonic flow over two- and three-Dimensional cavities has been analyzed by the integration of unsteady Reynolds-Averaged Navier-Stokes(RANS) with the k - w turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in cavities. Numerical method is upwind TVD scheme based on the flux vector split with the Van Leer limiters, and time accuracy is used explicit 4th stage Runge-Kutta scheme. Cavity flows are Comparison of two- and three-dimensional. The cavity has a L/D ratio of 3 for two-dimensional case. and same L/D and W/D ratio is 1 for three-dimensional case. The Mach and Reynolds numbers are held constant at 1.5 and 450000 respectively. For the three-dimensional case, the flow field is observed to oscillate in the 'shear layer mode' with a feedback mechanism that follow Rossiter's formula. On the other hand, the self-sustained oscillating flow transitions to a 'wake mode' for the two-dimensional simulation, with more violent fluctuations inside the cavity.

• Journal of Korean Society of Water and Wastewater
• /
• v.25 no.5
• /
• pp.635-642
• /
• 2011

Hydraulic efficiency was a vital component in evaluating the disinfection capability of clearwell. Current practice evaluates these system based on the tracer test only. In this paper, CFD(Computational Fluid Dynamics) was applied on the clearwell for alternating or supplementing the tracer test. The baffle factor derived from the CFD modeling closely matched the values obtained from full scale tracer testing. And, for suggesting proper numerical model in clearwell; the turbulence model, discretization scheme, convergence criteria were investigated through separate simulation runs. The model validation was conducted by comparing the simulated data with experimental data. In the turbulence model, the realizable ${\kappa}-{\varepsilon}$ model and the standard ${\kappa}-{\varepsilon}$ model were found to be more appropriate than RNG ${\kappa}-{\varepsilon}$ model. The residuals of convergence criteria should be used as not $10^{-3}$ but $10^{-4}$ or $10^{-5}$. In discretization scheme, the difference of simulated values in 1st, 2nd, 3rd upwind scheme was found to be insignificant. Moreover, the result of this study suggest that CFD modeling can be a reliable alternative to tracer testing for evaluating the hydraulic efficiency.

• 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.31 no.3
• /
• pp.87-99
• /
• 1994

Reynolds Averaged Navier-Stokes equations are solved numerically for the computation of turbulent flow around a Wigley double model. A second order finite difference method is applied for the spatial discretization on the nonstaggered grid system and 4-stage Runge-Kutta scheme for the numerical integration in time. In order to increase the time step, residual averaging scheme of Jameson is adopted. Pressure field is obtained by solving the pressure-Poisson equation with the appropriate Neumann boundary condition. For the turbulence closure, 0-equation turbulence model of Baldwin-Lomax is used. Numerical computation is carried out for the Reynolds number of 4.5 million. Comparisons of the computed results with the available experimental data show good agreements for the velocity and pressure distributions.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.165-176
• /
• 2015

A numerical analysis was conducted to investigate the inlet buzz and combustion oscillation in an axisymmetric ramjet engine with wedge-type flame holders. The physical model of concern includes the entire engine flow path, extending from the leading edge of the inlet center-body through the exhaust nozzle. The theoretical formulation is based on the Farve-averaged conservation equations of mass, momentum, energy, and species concentration, and accommodates finite-rate chemical kinetics and variable thermo-physical properties. Turbulence closure is achieved using a combined scheme comprising of a low-Reynolds number k-${\varepsilon}$ two-equation model and Sarkar's compressible turbulence model. Detailed flow phenomena such as inlet flow aerodynamics, flame evolution, and acoustic excitation as well as their interactions, are investigated. Mechanisms responsible for driving the inlet buzz are identified and quantified for the engine operating at subcritical conditions.

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