• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.104-110
• /
• 2008

Detached Eddy Simulation (DES) is applied to an axisymmetric base flow at supersonic mainstream. DES is a hybrid approach to modeling turbulence that combines the best features of the Reynolds-averaged Navier-Stokes RANS) and large-eddy simulation (LES) approaches. In the Reynolds-averaged mode, the model is currently based on either the Spalart-Allmaras (S-A) turbulence model. In the large eddy simulation mode, it is based on the Smagorinski subgrid scale model. Accurate predictions of the base flowfield and base pressure are successfully achieved by using the DES methodology with less computational cost than that of pure LES and monotone integrated large-eddy simulation (MILES) approaches. The DES accurately resolves the physics of unsteady turbulent motions, such as shear layer rollup, large-eddy motions in the downstream region, small-eddy motions inside the recirculating region. Comparison of the results shows that it is necessary to resolve approaching boundary layers and free shear-layer velocity profiles from the base edge correctly for the accurate prediction of base flows. The consideration of an empirical constant CDES for a compressible flow analysis may suggest that the optimal value of empirical constant CDES may be larger in the flows with strong compressibility than in incompressible flows.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.104-110
• /
• 2008

Detached Eddy Simulation (DES) is applied to an axisymmetric base flow at supersonic mainstream. DES is a hybrid approach to modeling turbulence that combines the best features of the Reynolds-averaged Navier-Stokes (RANS) and large-eddy simulation (LES) approaches. In the Reynolds-averaged mode, the model is currently based on either the Spalart-Allmaras (S-A) turbulence model. In the large eddy simulation mode, it is based on the Smagorinski subgrid scale model. Accurate predictions of the base flowfield and base pressure are successfully achieved by using the DES methodology with less computational cost than that of pure LES and monotone integrated large-eddy simulation (MILES) approaches. The DES accurately resolves the physics of unsteady turbulent motions, such as shear layer rollup, large-eddy motions in the downstream region, small-eddy motions inside the recirculating region. Comparison of the results shows that it is necessary to resolve approaching boundary layers and free shear-layer velocity profiles from the base edge correctly for the accurate prediction of base flows. The consideration of an empirical constant CDES for a compressible flow analysis may suggest that the optimal value of empirical constant CDES may be larger in the flows with strong compressibility than in incompressible flows.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.59-72
• /
• 2022

In this work the numerical results of the flow around a 5:1 rectangular cylinder at Reynolds numbers 3 000 and 40 000, zero angle of attack and smooth incoming flow condition are presented. Implicit Large Eddy Simulations (ILES) have been performed with a high-order accurate spatial scheme and an implicit high-order accurate time integration method. The spatial approximation is based on a discontinuous Galerkin (dG) method, while the time integration exploits a linearly-implicit Rosenbrock-type Runge-Kutta scheme. The aim of this work is to show the feasibility of high-fidelity flow simulations with a moderate number of DOFs and large time step sizes. Moreover, the effect of different parameters, i.e., dimension of the computational domain, mesh type, grid resolution, boundary conditions, time step size and polynomial approximation, on the results accuracy is investigated. Our best dG result at Re=3 000 perfectly agrees with a reference DNS obtained using Nek5000 and about 40 times more degrees of freedom. The Re=40 000 computations, which are strongly under-resolved, show a reasonable correspondence with the experimental data of Mannini et al. (2017) and the LES of Zhang and Xu (2020).

• Wind and Structures
• /
• v.8 no.1
• /
• pp.49-60
• /
• 2005

A tornado changes its wind speed and direction rapidly; therefore, it is difficult to study the effects of a tornado on buildings in a wind tunnel. The status of the tornado-structure interaction and various models of the tornado wind field found in literature are surveyed. Three dimensional computer modeling work using the turbulence model based on large eddy simulation is presented. The effect of a tornado on a cubic building is considered for this study. The Navier-Stokes (NS) equations are approximated by finite difference method, and solved by an semi-implicit procedure. The force coefficients are plotted in time to study the effect of the Rankine combined vortex model. The tornado is made to translate at a $0^{\circ}$ and $45^{\circ}$ angle, and the grid resolution is refined. Some flow visualizations are also reported to understand the flow behavior around the cube.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.426-432
• /
• 2004

Large-Eddy Simulation (LES) is applied for the simulation of compressible flat plate boundary with Reynolds number up to 5 X 10$^{5}$ . Numerical examples include shock/boundary layer interaction and boundary layer transition, aiming future application to the analysis of transonic fan/compressor cascades. The present LES code uses hybrid com-pact/WENO scheme for the spatial discretization and compact diagonalized implicit scheme for the time integration. The present code successfully predicted the bypass transition of subsonic boundary layer. As for supersonic turbulent boundary layer, mean and fluctuation velocity of the attached boundary, as well as the evolution of the friction coefficient and the displacement thickness both upstream and downstream of the separation region are all in good agreement with experiment. The separation point also agreed with the experiment. In the simulation of the shock/laminar boundary layer interaction, the dependence of the transition upon the shock strength is reproduced qualitatively, but the extent of the separation region is overpredicted. These numerical examples show that LES can predict the behavior of boundary layer including transition and shock interaction, which are hardly managed by the conventional Reynolds-averaged Navier-Stokes approach, although there needs to be more effort before achieving quantitative agreement.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2006.05a
• /
• pp.310-313
• /
• 2006

The present study aims at investigating the effectiveness of a new passive cavity flow control technique, sub-cavity. The characteristics of cavity flow oscillation with the device are compared with those with other control techniques tested previously, including a triangular bump and blowing jet. In the computation, the three-dimensional, unsteady Navier-Stokes equations governing the supersonic cavity flow are solved based on an implicit finite volume scheme spatially and multi-stage Runge-Kutta scheme temporally. Large eddy simulation (LES) is carried out to properly predict the turbulent features of cavity flow. The present results show that the pressure oscillation near the downstream edge dominates overall time-dependent cavity pressure variations, and the amplitude of the pressure oscillation can be reduced in the presence of a sub-cavity.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2005.11a
• /
• pp.117-120
• /
• 2005

The present study describes unsteady flow phenomena generated in a supersonic flow passing over a rectangular cavity and suggests a way of control of pressure oscillation, doing harm to overall performance and stable operation of aerodynamic and industrial applications. The three-dimensional, unsteady, compressible Navier-stokes equations are numerically solved based on a fully implicit finite volume scheme and large eddy simulation. The cavity flow are simulated with and without control methods, including a triangular bump and blowing jet installed near the leading edge of the cavity. The results show that the pressure oscillation is attenuated by both control techniques, especially near the trailing edge of cavity.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.392-397
• /
• 2001

A numerical study is performed to investigate the interaction between subsonic axial turbine blade boundary layer and periodically oncoming rotor induced wakes. An implicit scheme for solving the compressible Navier-Stokes equation is developed, which adopts a 4th-order compact difference for spatial discretiztion, a 2nd order Crank-Nicolson scheme for temporal discretization and the dynamic eddy viscosity model as the subgrid scale model. The efficiency and the accuracy of the proposed method are verified by applying to some benchmark problems such as laminar cylinder flow, laminar airfoil cascade flow and a transitional flat plate boundary layer flow. Computational results show good agreements with previous experimental and numerical results. Finally, flow through a stator cascade is simulated at $Re = 7.5{\times}10^5$ without free-stream turbulence intensity. The velocity fields and skin friction coefficients in the transitional region show similar trends with previous boundary layer natural transition.

• Wind and Structures
• /
• v.6 no.3
• /
• pp.209-220
• /
• 2003

A tornado changes its wind speed and direction rapidly; therefore, it is difficult to study the effects of a tornado on buildings in a wind tunnel. In this work, the status of the tornado-structure interaction is surveyed by numerical simulation. Various models of the tornado wind field found in literature are surveyed. Three-dimensional computer modeling work using the turbulence model based on large eddy simulation is presented. The effect of tornado on a cubic building is considered for this study. The Navier-Stokes (NS) equations are approximated by finite difference method, and solved by a semi-implicit procedure. The force coefficients are plotted in time to study the effect of the Rankine-Combined Vortex Model. Some flow visualizations are also reported to understand the flow behavior around the cube.

• Coupled systems mechanics
• /
• v.7 no.4
• /
• pp.421-440
• /
• 2018

In this study, a computational method for wall shear stress combined with an implicit direct-forcing immersed boundary method is presented. Near the immersed boundaries, the sub-grid stress is determined by a wall model in which the wall shear stress is directly calculated from the Lagrangian force on the immersed boundary. A coupling mathematical model of the transition process for a model Francis turbine comprising turbulent flow and rotating rigid guide vanes is established. The spatiotemporal distributions of pressure, velocity, vorticity and turbulent quantity are gained with the transient process; the drag and lift coefficients as well as other forces (moments) are also obtained as functions of the attack angle. At the same time, analysis is conducted of the characteristics of pressure pulsation, velocity stripes and vortex structure at some key parts of flowing passage. The coupling relations among the turbulent flow, the dynamical force (moment) response of blade and the rotating of guide vane are also obtained.