• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.8
• /
• pp.680-687
• /
• 2007

The paper presents a numerical simulation of flow of a turbulent boundary layer, representing a typical wind environment and matching a series of wind tunnel observations. The simulations are carried out at a Reynolds number of 20,000, based on the velocity U at a pseudo-height h, and large enough that the flow be effectively Reynolds number independent. Some wall models are proposed for the LES(Large Eddy Simulation) of the turbulent boundary layer over a rough surface. The Jenson number, $J=h/z_0$, based on the roughness length $z_0$, is 600 to match the wind tunnel data. The computational mesh is uniform with a spacing of h/32, as this aids rapid convergence of the multigrid solver, and the governing equations are discretised using second order finite differences within a parallel multiblock environment. The results presented include the comparison between wind tunnel measurements and LES computations of the turbulent boundary layer over rough surface.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.1
• /
• pp.69-84
• /
• 2020

Natural shoreline repeats its re-treatment and advance in response to the endlessly varying sea-conditions, and once severely eroded under stormy weather conditions, natural beaches are gradually recovered via a boundary layer streaming when swells are prevailing after storms cease. Our understanding of the boundary layer streaming over surf-zone often falls short despite its great engineering value, and here it should be noted that the most sediments available along the shore are supplied over the surf-zone. In this rationale, numerical simulation was implemented to investigate the hydraulic characteristics of boundary layer streaming over the surf zone in this study. In doing so, comprehensive numerical models made of Spatially filtered Navier-Stokes Eq., LES (Large Eddy Simulation), Dynamic Smagorinsky turbulence closure were used, and the effects of turbulence closure such as Dynamic Smagorinsky in LES and k-ε on the numerically simulated flow field were also investigated. Numerical results show that due to the intrinsic limits of k-ε turbulence model, numerically simulated flow velocity near the bottom based on k-ε model and wall function are over-predicted than the one using Dynamic Smagorinsky in LES. It is also shown that flow velocities near the bottom are faster than the one above the bottom which are relatively free from the presence of the bottom, complying the typical boundary layer streaming by Longuet-Higgins (1957), the spatial scope where boundary layer streaming are occurring is extended well into the surf zone as incoming waves are getting longer. These tendencies are plausible considering that it is the bottom friction that triggers a boundary layer streaming, and longer waves start to feel the bottom much faster than shorter waves.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.11a
• /
• pp.31-37
• /
• 2011

Detail flow structure in a large solid rocket motor with two inhibitors has been investigated using 3D Large Eddy Simulation and Proper Orthogonal Decomposition(POD) analysis. Vortex shedding frequencies periodically occurred by inhibitors are coupled with flow acoustics induced by the impinging of vorticity on nozzle head. As a result of 3D analysis, it was observed that the nozzle exit flow causes roll-torques from the vortex being decomposed in unbalanced shape for the impinging of vorticity on the nozzle head.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.09a
• /
• pp.231-232
• /
• 2005

• Proceedings of the Korean Nuclear Society Conference
• /
• 2001.05a
• /
• pp.160.2-160
• /
• 2001

• Journal of computational fluids engineering
• /
• v.19 no.4
• /
• pp.29-36
• /
• 2014

A computational study of a thermal striping in the upper plenum of PGSFR(Prototype Generation-IV Sodium-cooled Fast Reactor) being developed at the KAERI(Korea Atomic Energy Research Institute) is presented. The LES(Large Eddy Simulation) approach is employed for the simulation of thermal striping in the upper plenum of the PGSFR. The LES is performed using the WALE (Wall-Adapting Local Eddy-viscosity) model. More than 19.7 million unstructured elements are generated in upper plenum region of the PGSFR using the CFX-Mesh commercial code. The time-averaged velocity components and temperature field in the complicated upper plenum of the PGSFR are presented. The time history of temperature fluctuation at the eight locations of solid walls of UIS(Upper Internal Structure) and IHX(Intermediate Heat eXchanger) are additionally stored. It has been confirmed that the most vulnerable regions to thermal striping are the first plate of UIS. From the temporal variation of temperature at the solid walls, it was possible to find the locations where the thermal stress is large and need to assess whether the solid structures can endure the thermal stress during the reactor life time.

• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.60-66
• /
• 2010

Unsteady Reynolds Averaged Navier-Stokes(URANS) and Large Eddy Simulation(LES) simulation of an axial flow fan are calculated upon same conditions and computational grids in order to study aeroacoustic noise of an axial flow fan numerically. Results of computed performance and predicted noise are compared with those of measurement. Both performances show accurate results with a significant difference of less than 5%. However, noise of LES result is more close to measured noise qualitatively than URANS. Levels of tonal noises of both LES and URANS are quite similar with those of measured at BPF(Blade Passing Frequency) in sound spectrum. However, as leading edge separation and tip vortex shedding phenomena of LES are showed more clearly than those of URANS, sound level of broadband noise of LES corresponds better than that of URANS, especially.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.369-378
• /
• 2002

The application of Large Eddy Simulation (LES) in a curvilinear coordinate system to the flow around a square cylinder is presented. In order to obtain sufficient resolution near the side of the cylinder, we use an O-type grid. Even with a curvilinear coordinate system, it is difficult to avoid the numerical oscillation arising in high-Reynolds-number flows past a bluff body, without using an extremely fine grid used. An upwind scheme has the effect of removing the numerical oscillations, but, it is accompanied by numerical dissipation that is a kind of an additional sub-grid scale effect. Firstly, we investigate the effect of numerical dissipation on the computational results in a case where turbulent dissipation is removed in order to clarify the differences between the effect of numerical dissipation. Next, the applicability and the limitations of the present method, which combine the dynamic SGS model with acceptable numerical dissipation, are discussed.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2000.04a
• /
• pp.75-80
• /
• 2000

We revisit the arrested Ekman boundary layer problem, using a fully non-linear numerical model with the subgrid dissipation modeled by the large eddy simulation method (LES). The main objective of this study is to find out whether the dynamic balance of the arrested Ekman boundary layer explained by MacCready and Rhines (1991) is valid for high Reynolds number. The model solution indicates that for high Reynolds number and low Richardson number flows, the density anomaly diffusion by near-wall turbulent action may become intense enough to homogenize completely the density structure within the boundary layer, in the direction perpendicular to the sloping wall. Then the buoyancy effect becomes negligible allowing a near-equilibrium Ekman boundary layer flow to persist for a long period.

• Proceedings of the SAREK Conference
• /
• 2008.11a
• /
• pp.55-60
• /
• 2008

Turbulent flow characteristics on the gap of two parallel channels are investigated using LES(large eddy simulation) approach. Two parallel channels have the same cross-section area and are connected by the narrow channel named the gap. Turbulent flow near the gap makes the flow pulsation along the streamwise direction of two channels. The flow condition is the Reynolds number of $2.5{\times}10^{-5}$. We compared the predicted results with the previous experimental results and presented the axial mean velocity, turbulent intensities, Reynolds shear stresses and turbulent kinetic energy.