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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1229-1239
• /
• 1999

A computational study on unsteady compressible flows has been performed by adopting a low Reynolds number $k-{\omega}$ turbulence model in conjunction with dual time stepping scheme. An explicit four-stage Runge-Kutta scheme for the Navier-Stokes equations and an approximate factorization scheme for the $k-{\omega}$ turbulence model equations are used. Computational results obtained for blade surface pressure distributions in the process of rotor-stator interaction in a turbine stage are in good agreement with extant experimental data. The effects of the wake from the stator on the boundary-layer transition over the rotor blade surface are discussed by showing that high intensity turbulence of the stator wake induces an early transition.

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

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

• International Journal of Aeronautical and Space Sciences
• /
• v.1 no.1
• /
• pp.36-47
• /
• 2000

Three-dimensional compressible turbulent flow fields within the passage of a diffusing S-duct have been simulated by solving the Navier-Stokes equations with SIMPLE scheme. The average inlet Mach number is 0.6 and the Reynolds number based on the inlet diameter is $1.76{\times}10^6$ The extended $k-{\varepsilon}$ turbulence model is applied to modeling the Reynolds stresses. Computed results of the flow in a circular diffusing S-duct provide an understanding of the flow structure within a typical engine inlet system. These are compared with experimental wall static-pressure, total-pressure fields, and secondary velocity profiles. Additionally, boundary layer thickness, skin friction values, and streamlines in the symmetric plane are presented. The computed results depict the interaction between the low energy flow by the flow separation and the high energy flow by the reversed duct curvature. The computed results obtained using the extended $k-{\varepsilon}$ turbulence model.

• Journal of Mechanical Science and Technology
• /
• v.18 no.1
• /
• pp.153-166
• /
• 2004

This paper presents a comparative study of a fully coupled, upwind, compressible Navier-Stokes code with three two-equation models and the Baldwin-Lomax algebraic model in predicting transonic/supersonic flow. The k-$\varepsilon$ turbulence model of Abe performed well in predicting the pressure distributions and the velocity profiles near the flow separation over the axisymmetric bump, even though there were some discrepancies with the experimental data in the shear-stress distributions. Additionally, it is noted that this model has y$\^$*/ in damping functions instead of y$\^$＋/. The turbulence model of Abe and Wilcox showed better agreements in skin friction coefficient distribution with the experimental data than the other models did for a supersonic compression ramp problem. Wilcox's model seems to be more reliable than the other models in terms of numerical stability. The two-equation models revealed that the redevelopment of the boundary layer was somewhat slow downstream of the reattachment portion.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2006.11a
• /
• pp.323-326
• /
• 2006

This paper presents the computational results for the two-dimensional compressible non-reacted flow in a converging-diverging micro thrust nozzle of which the ratio of exit to throat width (0.541 in.) is 1.8. The RNG model is applied to calculate the turbulence by loading the standard coefficients. The results agreed very well with the experiments in the view of the shock structure and the pressure distribution at the various pressure ratios between the stagnation and the environmental states. The plume structures are also discussed on the view of the shock-cell structure.

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

Numerical Analysis has been done for the supersonic off-design jet flow due to the pressure difference between the jet and the ambient fluid. The difference of pressure generates an oblique shock or an expansion wave at the nozzle exit, The waves reflect repeatedly at the center axis and on the sonic surface in the shear layer, and the pressure difference is resolved across these waves interacted with the turbulence mixing layer. In this paper, the axi-symmetric Navier-Stokes equation has been used with two equation $k-{\varepsilon}$ turbulence closure model. The second order TVD scheme with flux limiters, based on the flux vector split by the smooth eigenvalue split, has been used to capture internal shocks and other discontinuities. The correction term for the compressible flow and the damping function are used in the turbulence model. Numerical calculations have been done to analyze the off-design jet flow due to the pressure difference. The variation of pressure along the flow axis is compared with an experimental result and other numerical result. The characteristics of the interaction between the shock cell and the turbulence mixing layer have been analyzed.

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3367-3382
• /
• 2023

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian computational fluid dynamics method that has been widely used in the analysis of physical phenomena characterized by large deformation or multi-phase flow analysis, including free surface. Despite the recent implementation of eddy-viscosity models in SPH methodology, sophisticated turbulent analysis using Lagrangian methodology has been limited due to the lack of computational performance and numerical consistency. In this study, we implement the standard and dynamic Smagorinsky model and dynamic Vreman model as sub-particle scale models based on a weakly compressible SPH solver. The large eddy simulation method is numerically identical to the spatial discretization method of smoothed particle dynamics, enabling the intuitive implementation of the turbulence model. Furthermore, there is no additional filtering process required for physical variables since the sub-grid scale filtering is inherently processed in the kernel interpolation. We simulate lid-driven flow under transition and turbulent conditions as a benchmark. The simulation results show that the dynamic Vreman model produces consistent results with experimental and numerical research regarding Reynolds averaged physical quantities and flow structure. Spectral analysis also confirms that it is possible to analyze turbulent eddies with a smaller length scale using the dynamic Vreman model with the same particle size.

• Journal of computational fluids engineering
• /
• v.8 no.2
• /
• pp.1-7
• /
• 2003

The V2F turbulence model, which has shown very good performance in several test cases at low speeds, has been applied to supersonic ramp flow with 20. corner angle at the free stream Mach number of 2.79. The flow is known to manifest strong shock wave/turbulent boundary layer interactions. As a comparative study, low-Reynolds k-ε models are also considered. While the V2F model predicts wall-pressure distribution well, it relatively predicts larger separation bubble and higher skin-friction after the reattachment than the experimental data. Although the ellpticity of f equation is the characteristics of incompressible flows, the converged solutions are acquired in the compressible flow with shock waves. The effect of the realizability constraints used in the model is also examined. In contrast to the result of impinging jet flows, the realizability bounds proposed by Durbin deterioate the overall solutions of the supersonic ramp flow.