• 대한기계학회논문집B
• /
• 제31권9호
• /
• pp.814-821
• /
• 2007

A hybrid method is presented for efficient computation of turbulent flow noise at low Mach numbers. In this method, the turbulent flow field is computed by incompressible large eddy simulation (LES), while the acoustic field is computed with the linearized perturbed compressible equations (LPCE) derived in this study. Since LPCE is computed on the rather coarse acoustic grid with the flow variables and source term obtained by the incompressible LES, the computational efficiency of calculation is greatly enhanced. Furthermore, LPCE suppress the instability of perturbed vortical mode and therefore secure consistent and stable acoustic solutions. The proposed LES/LPCE hybrid method is applied to three low Mach number turbulent flow noise problems: i) circular cylinder, ii) isolated flat plate, and iii) interaction between cylinder wake and airfoil. The computed results are closely compared with the experimental measurements.

• 한국자동차공학회논문집
• /
• 제6권2호
• /
• pp.12-20
• /
• 1998

Exter ballistics of a typical high-speed projectile is studied through a flow-visualization experiment and an unstructured grid Navier-Srokes computation. Experiment produced a schlieren photograph that adequately shows the characteristic features of this complex flow, namely two kinds of oblique cone shocks and turbulent wake developing into the downstream. A hybrid scheme of finite volume-element method is used to simulate the compressible Reynolds-Averaged Navier-Stok- es solution on unstructured grids. Osher's approximate Riemann solver is used to discretize the cinvection term. Higher-order spatial accuracy is obtained by MUSCL extension and van Albada ty- pe flux limiter is used to stabilize the numerical oscillation near the solution discontinuity. Accurate Gakerkin method is used to discretize the viscous term. Explict fourth-order Runge-Kutta method is used for the time-stepping, which simplifies the application of MUSCL extension. A two-layer k-$\varepsilon$ turbulence model is used to simulate the turbulent wakes accurately. Axisymmetric folw and two-dimensional flow with an angle of attack have been computed. Grid-dependency is also checked by carrying out the computation with doubled meshes. 2-D calculation shows that effect of angle of attack on the flow field is negligible. Axi-symmetric results of the computation agrees well with the flow visualization. Primary oblique shock is represented within 2-3 meshes in numerical results, and the varicose mode of the vortex shedding is clearly captured in the turbulent wake region.

• 대한기계학회논문집B
• /
• 제22권6호
• /
• pp.771-787
• /
• 1998

A Navier-Stokes code with a modified low Reynolds number k-.epsilon. turbulence model was used to study the unsteady transitional boundary layer flow due to rotor-stator interaction. The modification, proposed by Launder, to improve prediction of stagnation flows was incorporated to the low Reynolds number k-.epsilon. turbulence model by Fan-Lakshminarayana-Barnett. Numerical solution is shown to capture well the calmed laminar flow as well as the wake induced transitional strip due to rotor-stator interaction and shows improvement, in terms of onset of transition and its length, over previous Euler/boundary layer solution. The turbulent kinetic energy shows local maximum along the upstream rotor wake in the wake induced transitional strip and this characteristics is observed untill the end of transition. The wake induced strip also shown apparent even in the laminar sublayer as the upstream rotor wake penetrates inside the boundary layer.

• 한국전산유체공학회지
• /
• 제4권2호
• /
• pp.17-30
• /
• 1999

The unsteady transitional boundary layer due to rotor-stator interaction was studied at two operation points, the design and one off design points. The off design point leads to lower blade loading and lower Reynolds number. A Navier-Stokes code developed in the previous study was parallelized to expedite computations. A low Reynolds number turbulence model was used to close the momentum equations. All computations show good agreement with experimental data. The wake induced transitional strip on the suction side of the stator is clearly captured at design point operation. There is no noticeable change in shape and phase angle of the wake induced strip even in the laminar sublayer. The wake induced transitional strip at off design point shows more complex structure. The wake induced transitional strip is observed only in the turbulent layer, and becomes obscure in the laminar sublayer and buffer layer. This behavior is probably consequent upon that the transition is governed by both wake induced strip and natural transition mechanism by Tollmien-Schlichting wave.

• 한국가시화정보학회지
• /
• 제21권2호
• /
• pp.24-33
• /
• 2023

This work is to observe the wake flow generated behind a ramp. We have conducted a large eddy simulation with two ramp models having different heights with two different inflow conditions. Reynolds number based on the height of the large ramp (LR) and small ramp (SR) are Re_h = 2.8×10⁴ and 1.4×10⁴ respectively. The wake flow visualization shows the formation of streamwise counter-rotating vortices pairs at the downstream of the obstacle. These primary vortices are stretched and lifted up when moving downstream. In order to observe the effect of the inflow condition on the wake transition, two different inlet flow conditions are given on the inlet section as an inlet boundary condition. Induced counter-rotating vortices pairs due to sharp-edged triangular ramp obstacles are developed and propagated downstream. In the result, the large ramp shows a more complicated wake structure of the boundary layer than the small ramp.

• 대한조선학회논문집
• /
• 제47권1호
• /
• pp.20-29
• /
• 2010

The complicated flow characteristics of upper propeller wake influenced by hull wake are investigated in detail in the present study. A two-frame PIV (particle image velocimetry) technique was employed to visualize the upper propeller wake region. As the upper hull wake affects strongly propeller inflow, upper propeller wake shows much unstable vortical behavior, especially in the tip vortices. Velocity field measurements were conducted in a cavitation tunnel with a simulated hull wake. Generally, the hull wake generated by the hull of a marine ship may cause different loading distributions on the propeller blade in both upper and lower propeller planes. The unstable upper propeller wake caused by the ship's hull is expressed in terms of turbulent kinetic energy (TKE) and is identified by using the proper orthogonal decomposition (POD) method to characterize the coherent flow structure in it. Instabilities appeared in the eigen functions higher than the second one, giving unsteadiness to the downstream flow characteristics. The first eigen mode would be useful to find out the tip vortex positions immersed in the unstable downstream region.

• 한국가시화정보학회지
• /
• 제19권3호
• /
• pp.92-98
• /
• 2021

Large Eddy Simulation (LES) has been popularly applied and used in the last several decades to simulate turbulent boundary layer in the numerical domain. A fully developed turbulent boundary layer has also been applied to predict the complicated wake flow behind bluff bodies. In this study we aimed to generate an artificial turbulent boundary layer, which is based on an exponential correlation function, and generates a series of realistic three-dimensional velocity data in two-dimensional inlet section which are correlated both in space and in time. The results suggest its excellent capability for high Reynolds number flows. To make an effective generation, a hexahedral mesh has been used and Cholesky decomposition was applied to possess suitable turbulent statistics such as the randomness and correlation of turbulent flow. As a result, the flow characteristics in the domain and fluctuating pressure near the wall are very close to those of fully developed turbulent boundary layers.

• Wind and Structures
• /
• 제15권1호
• /
• pp.17-26
• /
• 2012

Modeling swirling wakes is of considerable interest to wind farm designers. The present work is an attempt to develop a computational tool to understand free, far-wake development behind a single rotating wind turbine. Besides the standard momentum and continuity equations from the boundary layer theory in two dimensions, an additional equation for the conservation of angular momentum is introduced to study axisymmetric swirl effects on wake growth. Turbulence is simulated with two options: the standard ${\kappa}-{\varepsilon}$ model and the Reynolds Stress transport model. A finite volume method is used to discretize the governing equations for mean flow and turbulence quantities. A marching algorithm of expanding grids is employed to enclose the growing far-wake and to solve the equations implicitly at every axial step. Axisymmetric far-wakes with/without swirl are studied at different Reynolds numbers and swirl numbers. Wake characteristics such as wake width, half radius, velocity profiles and pressure profiles are computed. Compared with the results obtained under similar flow conditions using the computational software, FLUENT, this far-wake model shows simplicity with acceptable accuracy, covering large wake regions in far-wake study.

• 대한기계학회논문집B
• /
• 제32권7호
• /
• pp.518-528
• /
• 2008

The effects of surface roughness on a turbulent boundary layer (TBL) were investigated using particle image velocimetry (PIV). The roughness elements used were periodically arranged two-dimensional spanwise rods, and the roughness height was ${\kappa}/{\delta}$. Introduction of the roughness elements increased the wake strength and the turbulent stress not only in the roughness sublayer but also in the outer layer. This indicates the existence of interaction between inner and outer layers for 2D rod-roughened wall. Roughness effects on a turbulence structure near the wall were obtained by PIV measurements. Iso-contours of mean velocities and Reynolds stresses in the roughness sublayer showed a very good agreement with previous DNS results.

• Wind and Structures
• /
• 제5권2_3_4호
• /
• pp.317-328
• /
• 2002

A numerical investigation on the turbulent flows over a three-dimensional steep hill is presented. The numerical model developed for the present work is based on the finite volume method and the SIMPLE algorithm with a non-staggered grid system. Standard $k-{\varepsilon}$ model and Shih's non-linear model are tested for the validation of the prediction accuracy in the 3D separated flow. Comparisons of the mean velocity and turbulence profiles between the numerical predictions and the measurements show good agreement. The Shih's non-linear model is found to predict mean flow and turbulence better than the Standard $k-{\varepsilon}$. Flow patterns have also been examined to explain the difference in the cavity zone between 2D and 3D hills.