• Journal of Environmental Science International
• /
• v.11 no.7
• /
• pp.679-686
• /
• 2002

Turbulence greatly influence on atmospheric flow field. In the atmosphere, turbulence is represented as turbulent diffusion coefficients. To estimate turbulent diffusion coefficients in previous studies, it has been used constants or 2-level method which divides surface layer and Ekman layer. In this study, it was introduced Smagorinsky method which estimates turbulent diffusion coefficient not to divide the layer but to continue in vertical direction. We simulated 3-D flow model and TKE equation applied turbulent diffusion coefficients using two methods, respectively. Then we showed the values of TKE and the condition of each term to TKE. The results of Smagorinsky method were reasonable. But the results of 2-level method were not reasonable. Therefor, it had better use Smagorinsky method to estimate turbulent diffusion coefficients. We are expected that if it is developed better TKE equation and model with study of computational method in several turbulent diffusion coefficients for reasonably turbulent diffusion, we will able to predict precise wind field and movements of air pollutants.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.4
• /
• pp.1-8
• /
• 2004

As a computer has been continuously progressed to reduce R&D time and cost, the study of the flow physics has been significantly relied on the numerical method. Recently, Large Eddy Simulation(LES) has been widely used in CFD community to accurately capture the turbulent flows. The LES code requires high accuracy in time, as well as in space. Also, it should have strong robustness to ensure the convergence in various complicated flows. The objective of the present work is to develop a base code for LES simulation, having 2$^{nd}$ order accuracy in time and 4$^{th}$ order accuracy in space. To achieve the present objective, the four-step fractional step method was enhanced by adopting compact Pade'scheme. The standard Smagorinsky model was implemented for the first stage of the present code development. The flows over a cylinder and an airfoil were successfully simulated. and an airfoil were successfully simulated.

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

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.3
• /
• pp.193-204
• /
• 2011

To develop more robust oil boom which is vulnerable to various failure mode under severe weather condition, highly accurate wave model is developed using Spatially filtered Navier-Stokes Eq., LDS (Lagrangian Dynamic Smagorinsky model) for residual stresses, SPH (Smoothed Particle Hydrodynamics). To clarify the hydraulic characteristics of floating type oil boom, we numerically simulate the behavior of oil spill around oil boom under very energetic progressive waves. At the first stage, we firmly anchored the oil boom, and then, allowed the excursion of the oil boom. It turns out that oil boom with skirt of enough length (longer than 30% of depth) effectively confines the oil spill even against very energetic waves. We can also observe obliquely descending vertical eddies between y = 1~2 m as horizontal vortices shedding at the interface of oil spill and water are diffused toward the bottom, which is believed to be the birth, growing and break-down of Kelvin-Helmholz wave.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.11
• /
• pp.1523-1534
• /
• 2000

The performance of a number of existing dynamic subgrid-scale(SGS) models is evaluated in large-eddy simulations(LES) of two prototype transitional and turbulent shear flows, a planar jet and a channel flow. The dynamic SGS models applied include the dynamic Smagorinsky model(DSM);Germano et al. 1991, Lully 1992), the dynamic tow-component model(DTM; Akhavan et al. 2000), the dynamic mixed model(DMM;Zang et al, 1993). and the dynamic two-parameter model(DTPM; Salvetti & Banerjee 1995). The results are compared with those for DNS for their evaluation. The LES results demonstrate the superior performance of DTM with use of a sharp cutoff filter and DMM with use of a box filter, as compared to their respect counterpart DSM, in predicting the mean statistics, spectra and large-scale structure of the flow, Such features of DTM and DMM derive from the construction of the models in which tow separate terms are included to represent the SGS interactions; a Smagorinsky edd-viscosity term to account for the non-local interactions, and a local-interaction term to account for the nonlinear dynamics between the resolved and subgrid scales in the vicinity of the LES cutoff. As well, overall the SGS models using a sharp cutoff filter are more successful than those using a box filter in capturing the statistics and structure of the flow. Finally, DTPM is found to be compatible or inferior to DMM.

• Journal of Power System Engineering
• /
• v.4 no.3
• /
• pp.40-47
• /
• 2000

The turbulent flow with wake, reattachment and recirculation flow is very important from the viewpoint of engineering. But that is still difficult because of especially the unsteady problems which are related with the vehicle dynamics and the aerodynamics noise. This paper evaluate LES that can analyze about all fluid flow region including the laminar, transition and turbulent. So we compare the results of LES with those of PIV measurement and Reynolds averaging models. In conclusion, LES predicts flow behavior better than Reynolds averaging models.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.407-413
• /
• 2011

VIC (Vortex-In-Cell) method for viscous incompressible flow is presented to simulate the wake behind a modified NACA16 foil. With uniform rectangular grid, the velocity in field is calculated using streamfunction from vorticity field by solving the Poisson equation in which FFT(Fast Fourier Transform) is combined with 2nd order finite difference scheme. Here, LES(Large Eddy Simulation) with Smagorinsky model is applied for turbulence calculation. Effective viscosity is formulated using magnitude of strain tensor(or vorticity). Then the turbulent diffusion as well as viscous diffusion becomes particle strength exchange(PSE) with averaged eddy viscosity. The well-established panel method is combined to obtain the irrotational velocity and to apply the no-penetration boundary condition on the body panel. And wall diffusion is used for no-slip condition numerical results of turbulent stresses are compared with experimental results (Bourgoyne, 2003). Before comparing process, LES(Large Eddy Simulation) SGS(Subgrid scale) stress is transformed Reynolds averaged stress (Winckelmans, 2001).

• Journal of Korea Water Resources Association
• /
• v.44 no.6
• /
• pp.461-470
• /
• 2011

Abrupt and gradual levee breach analyses on the flat domain were implemented by laboratory experiments and numerical simulations. To avoid the reflective wave from the side wall the experiment was performed in a large domain surrounded by waterway. A numerical model was developed for solving the two-dimensional gradual levee breach flow. The results of the numerical simulation developed in this study showed good agreement with those of the experimental data. However, even if the numerical schemes effectively replicated the trends of the observed water depth for the first shock, there were little differences for the second shock. In addition, even though the model considered the Smagorinsky horizontal eddy viscosity, the location and height of the hydraulic jump in the numerical simulation were not fairly well agree with experimental measurements. This shows the shallow water equation solver has a limitation which does not exactly reproduce the energy dissipation from the hydraulic jump. Further study might be required, considering the energy dissipation due to the hydraulic jump or transition flow from reflective wave.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.4
• /
• pp.343-351
• /
• 2016

Cryogenic spray characteristics of a nitrogen swirl injector operating in supercritical environment have been numerically investigated. By comparing the equation of states(EOS) used for supercritical condition, SRK EOS was applied to predict the nitrogen thermodynamic property under supercritical environment. A Chung's method was implemented for the calculation of viscosity and conductivity and Takahashi's correlation based on Fuller's Theorem was implemented for the calculation of diffusion coefficient. By injecting the nitrogen with 5 bar differential pressure into 50 bar chamber filled with nitrogen, numerical simulation has been conducted. The dynamic Smagorinsky sub-grid scale (SGS) model has been compared with the algebraic Smagorinsky SGS model using FFT frequency analysis. The instability at the liquid film and gas core inside injector and the propagation of pressure oscillation into the injector has been investigated. The spreading angle of swirl injector obtained by numerical calculation has been validated with experimental result.

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

In order to investigate the hydraulic characteristics of a boundary layer streaming over the beach cusps appeared in swells prevailing mild seas, we numerically simulated the shoaling process of Edge waves over the beach cusp. Synchronous Edge waves known to sustain the beach cusps could successfully be duplicated by generating two obliquely colliding Edge waves in front of beach cusps. The amplitude A_B and length L_B of Beach Cusp were elected to be 1.25 m and 18 m, respectively based on the measured data along the Mang-Bang beach. Numerical results show that boundary layer streaming was formed at every phase of shoaling process without exception, and the maximum boundary layer streaming was observed to occur at the crest of sand bar. In RUN 1 where the shortest waves were deployed, the maximum boundary layer streaming was observed to be around 0.32 m/s, which far exceeds the amplitude of free stream by two times. It is also noted that the maximum boundary layer streaming mentioned above greatly differs from the analytical solution by Longuet-Higgins (1957) based on wave Reynolds stress. In doing so, we also identify the recovery procedure of natural beaches in swells prevailing mild seas, which can be summarized such as: as the infra-gravity waves formed in swells by the resonance wave-wave interaction arrives near the breaking line, the sediments ascending near the free surface by the Phase II waves orbital motion were carried toward the pinnacle of foreshore by the shoreward flow commenced at the steep front of breaking waves, and were deposited near the pinnacle of foreshore due to the infiltration.