• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.4 s.247
• /
• pp.364-372
• /
• 2006

The objective of this study is to evaluate the efficiency and the prediction accuracy of developed large eddy simulation (LES) program for complex turbulent flows, such as recirculating and swirling flows. To save the computational cost, a Beowulf cluster system consisting 16 processors was constructed. The flows in backward-facing step and dump combustor were examined as representative recirculating and swirling flows. Firstly, a direct numerical simulation (DNS) for laminar backward-facing step flows was previously conducted to validate the overall performance of program. Then LES was carried out for turbulent backward-facing step flows. The results of laminar flow showed a qualitative and quantitative agreement between simulations and experiments. The simulations of the turbulent flow also showed reasonable results. Secondly, LES results for non-swirling and swirling flows in a dump combustor were compared with the results of Reynolds-averaged Navier-Stokes (RANS) using standard $k-{\varepsilon}$ model. The results show that LES has a better performance in predicting the mean axial and azimuthal velocities, comer recirculation zone (CRZ) and center toroidal recirculation zone (CTRZ) than those of RANS. Finally, it was examined the capability of LES for the description of unsteady phenomena.

• Journal of the Society of Naval Architects of Korea
• /
• v.55 no.3
• /
• pp.265-273
• /
• 2018

Direct numerical and large eddy simulations of transitional flows around studs installed on flat plate and bulbous bow have been performed to investigate an effectiveness of turbulence stimulators on laminar-to-turbulence transition at a very low speed. The flow velocity was determined to be 0.366m/s corresponding to 4 knots of full-scale ship speed when the objective ship was Kriso container ship. The spatial evolution of skin friction coefficient disclosed that a fully development of turbulence was observed behind the second stud installed on flat plate while a rapid transition from laminar to turbulence gave rise to the fully development of turbulence behind the first stud installed on bulbous bow. A comparison of streamwise mean velocity profiles showed that the viscous sublayer and log-layer were in good agreement with previous results although the friction velocity of Smagrosinsky sub-grid scale model was about 10% larger than that of direct numerical simulation. While the turbulence intensities of bulbous bow was similar to those of flat plate in inner region, larger intensities of turbulence were observed in outer region of bulbous bow than those of flat plate.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.9
• /
• pp.1282-1289
• /
• 2003

Transport of a scalar quantity, such as chemical concentration or temperature, is important in many engineering applications and environmental flows. Here we report on results obtained from the large eddy simulations of flow and concentration fields inside the tank performed using a spectral multi-domain technique. The computations were driven by specifying the impeller-induced flow at the blade tip radius (Yoon et al.). This study focused on the concentration development at different molecular diffusivities in a stirred tank operated under turbulent conditions. The main objective of the work presented here is to study the large-scale mixing structure at different molecular diffusivities in a stirred tank by using the large eddy simulation. The time sequence of concentration and flow fields shows the flow dependency of the concentration development. The presence of spatial inhomogenieties is detailed by observing the time variation oflocal concentration at different positions.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.4
• /
• pp.346-352
• /
• 2009

The objective of the present work is to improve numerical predictions of unsteady turbulent swirling flows in the draft tubes of hydraulic power plants. We present Large Eddy Simulation (LES) results on a simplified draft tube consisting of a straight conical diffuser. The basis of LES is to solve the large scales of motion, which contain most of the energy, while the small scales are modeled. LES strategy is here preferred to the average equations strategies (RANS models) because it resolves directly the most energetic part of the turbulent flow. LES is now recognized as a powerful tool to simulate real applications in several engineering fields which are more and more frequently found. However, the cost of large-eddy simulations of wall bounded flows is still expensive. Bypass methods are investigated to perform high-Reynolds-number LES at a reasonable cost. In this study, computations at a Reynolds number about 2 $10^5$ are presented. This study presents the result of a new near-wall model for turbulent boundary layer taking into account the streamwise pressure gradient (adverse or favorable). Validations are made based on simple channel flow, without any pressure gradient and on the data base ERCOFTAC. The experiments carried out by Clausen et al. [1] reproduce the essential features of the complex flow and are used to develop and test closure models for such flows.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1806-1813
• /
• 2003

The stirred tank reactor is one of the most commonly used devices in industry for achieving mixing and reaction. Here we report on results obtained from the large eddy simulations of flow inside the tank performed using a spectral multi-domain technique. The computations were driven by specifying the impeller-induced flow at the blade tip radius. Stereoscopic PIV measurements (Hill et $al.^{(1)}$) along with the theoretical model of the impeller-induced flow (Yoon et $al.^{(2)}$) were used in defining the impeller-induced flow as superposition of circumferential, jet and tip vortex pair components. Large eddy simulation of flow in a stirred tank was carried out for the three different Reynolds numbers of 4000, 16000 and 64000. The effect of different Reynolds numbers is well observed in both instantaneous and time averaged flow fields. The instantaneous and mean vortex structures are identified by plotting an isosurfaces of swirling strength for all Reynolds numbers. The Reynolds number dependency of the nondimeansional eddy viscosity, resolve scale and subgrid scale dissipations is clearly shown in this study.

• International Journal of Automotive Technology
• /
• v.7 no.4
• /
• pp.415-422
• /
• 2006

Swirl and tumble flows in an engine in-cylinder have been simulated by using a three-dimensional computational fluid dynamics code, and the results are validated in comparison with experimental data. The large eddy simulation based on the Smagorinsky model and the fractional step method is adopted to describe the turbulence of in-cylinder flows and to save computing time, respectively. The main purpose of this study is connected with the effect of various conditions of intake flows on formation and development of in-cylinder tumble and swirl motions. The engine speeds considered are 1000 rpm and 3000 rpm for intake flows with inclination angles between $-10^{\circ}$ and $20^{\circ}$ at deflection angles of $0^{\circ}$, $22.5^{\circ}$, and $30^{\circ}$. The results are discussed by visualizing flow fields and by evaluating parameters in relation to vortex intensity such as swirl and tumble ratios.

• The KSFM Journal of Fluid Machinery
• /
• v.13 no.3
• /
• pp.5-11
• /
• 2010

As a basic study of the aero-acoustic noise, Large eddy simulations were carried out for a fixed circular cylinder at Renolds number (Re=$9.0\times10^4$) using commercial CFD code, FLUENT. The subgrid-scale turbulent viscosity was modeled by Smagorinsky-Lilly model adapted to structured meshes. The results of analysis showed that time-averaged value, $\bar{C}_D$ is approximately 1.47 which is considerably adjacent with the experimentally measured value of 1.32 in comparison to the values performed by previous researchers. It is observed that there are the very small acoustic pressure fluctuation with the same frequency of the Karman vortex street.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.807-810
• /
• 2002

The suitability of high-order accurate, central and upwind-biased compact difference schemes is evaluated for the large-eddy simulations of flows in complex geometry. Two flow geometries are considered: channel and circular cylinder. The effects of numerical dissipation and aliasing error on the evaluation of subgrid scale stress are investigated by extending the analysis by Ghosal (1) to centered and upwind compact schemes. It is shown that the failure of upwind schemes mainly comes from the aliasing error.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.973-983
• /
• 2003

The suitability of high-order accurate, centered and upwind-biased compact difference schemes is evaluated for large eddy simulation of turbulent flow. Two turbulent flows are considered: turbulent channel flow at Re = 23000 and flow over a circular cylinder at Re = 3900. The effects of numerical dissipation on the finite differencing and aliasing errors and the subgrid-scale stress are investigated. It is shown through the simulations that compact upwind schemes are not suitable for LES, whereas the fourth order-compact centered scheme is a good candidate for LES provided that proper dealiasing of nonlinear terms is performed. The classical issue on the aliasing error and the treatment of nonlinear terms is revisited with compact difference schemes.

• Wind and Structures
• /
• v.1 no.2
• /
• pp.127-144
• /
• 1998

In this paper, the wake structures behind a square cylinder at the Reynolds number of 22,000 are simulated using the large eddy simulation, and the main features of the wake structure associated with unsteady vortex-shedding are investigated. The Smagorinsky model is used for parametrization of the subgrid scales. The finite element method with isoparametric linear elements is employed in the computations. Unsteady computations are performed using the explicit method with streamline upwind scheme for the advection term. The time integration incorporates a subcycling strategy. No-slip condition is enforced on the wall surface. A comparative study between two-and three-dimensional computations puts a stress on the three-dimensional effects in turbulent flow simulations. Simulated three-dimensional wake structures are compared with numerical and experimental results reported by other researchers. The results include time-averaged, phase-averaged flow fields and numerically visualized vortex-shedding pattern using streaklines. The results show that dynamics of the vortex-shedding phenomenon are numerically well reproduced using the present method of finite element implementation of large eddy simulation.