• Journal of the Korean Society of Propulsion Engineers
• /
• v.14 no.2
• /
• pp.1-9
• /
• 2010

LES(Large eddy simulation) of breakup and droplet atomization of a liquid jet into cross turbulent flow was performed. Two phase flow of gas and liquid phases were modeled by the mixed numerical scheme of both Eulerian and Lagrangian methods for gas and liquid droplet respectively. The breakup process of a liquid column and droplets was observed by implementing the blob-KH wave breakup model. The penetration depth into cross flow was comparable with experimental data for several variants of the liquid-gas momentum flux ratio by varying liquid injection velocity. SMD(Sauter Mean Diameter) distribution downstream of jet was analyzed.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.5
• /
• pp.387-393
• /
• 2019

In this study, the free motion of a riser due to vortex shedding was numerically simulated with Large Eddy Simulation (LES) and Detached Eddy Simulation (DES) turbulence models. A numerical simulation program was developed by applying the Rhie-Chow interpolation method to the pressure correction of the OpenFOAM standard solver pimpleDyMFoam. To verify the developed program, the vortex shedding around the fixed riser at Re = 3900 was calculated, and the results were compared with the existing experimental and numerical data. Moreover, the vortex-induced vibration of a riser supported by a linear spring was numerically simulated while varying the spring constant. The results are compared with published direct numerical simulation (DNS) results. The present calculation results show that the numerical method is appropriate for simulating the vortex-induced motion of a riser, including lock-in phenomena.

• International Journal of Aeronautical and Space Sciences
• /
• v.18 no.2
• /
• pp.197-205
• /
• 2017

We investigate the slat noise generation mechanism by using large-eddy simulation (LES) and simple source modeling based on linearized Euler equations. An incompressible LES of an MD 30P30N three-element airfoil in the high-lift configuration is conducted at $Re_c=1.7{\times}10^6$. Using the total derivative of the hydrodynamic pressure (DP/Dt) acquired from the incompressible LES, representative noise sources in the slat cove region are characterized in terms of simple sources such as frequency-specific monopoles and dipoles. Acoustic radiation around the 30P30N multi-element airfoil is effectively computed using the Brinkman penalization method incorporated with the linearized Euler equation. The directivity pattern of $p^{\prime}_{rms}$ at $r=20c_{slat}$ in the multiple sources is closely compared to that obtained by the application of the LES/Ffowcs-Williams and Hawking's methods to the entire flow field. The power spectrum of p' at ${\theta}=290^{\circ}$ is in good agreement with the data reported in BANC-III, especially the broadband part of the spectrum with a decaying slope ${\propto}f^{-3}$.

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

A two-phase Large Eddy Simulation(LES) has been conducted to investigate breakup and atomization of a liquid jet in a cross turbulent flow in a rectangular duct. Gas-droplet two-phase flow was solved by a coupled Eulerian-Lagrangian method which tracks every individual particles. Effects of liquid breakup models, sub-grid scale models, and a order of spatial discretization was investigated. The penetration depth in cross flow was comparable with experimental data by varying breakup model and LES scheme. SMD(Sauter Mean Diameter) distribution downstream of jet was analyzed.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.05a
• /
• pp.210-213
• /
• 2010

Recent experimental data shows that the noticeable feature of irregular roughened spots on the fuel surface occurs during the combustion test. The generation of these unexpected patterns is likely to be resulted from the disturbed boundary layer due caused by wall blowing which is intended to simulate the process of fuel vaporization. LES without chemical reaction was conducted to investigate the flow characteristics at the near-fuel surface and the behavior of turbulent structures which is evolved by the wall blowing at the Reynolds number of 23,000. Cylindrical geometry was considered to get the most reality of the calculation results because real hybrid rocket motor is circular grain configuration. It was shown that the wall blowing pushed turbulent structures upwards making them tilted and this skewed displacement, in effect, left the foot prints of the structures on the surface. This change of kinematics may explain the formation of irregular isolated spots on the fuel surface observed in the experiment.

• Journal of the Society of Naval Architects of Korea
• /
• v.59 no.3
• /
• pp.141-148
• /
• 2022

A hybrid turbulence model has developed by combining a sub-grid scale model using dynamic k equation in LES with k-𝜔 SST model of RANS equation. To ascertain potential applicability of the hybrid turbulence model, fully developed turbulent channel flows at Re_𝜏=180 have been simulated of which computational domain has a top wall with coarse cells and a bottom wall with fine cells. The streamwise mean velocity and turbulent intensity profiles showed a good agreement with DNS data when using the hybrid model rather than using a single model in k-𝜔 SST or dynamic k equation models. Computational simulations of turbulent flows around KVLCC2 with a pre-swirl duct have been mainly performed using the hybrid turbulence model. Compared to the results obtained from RANS simulation with k-𝜔 SST model as well as LES with dynamic k equation SGS model, turbulent wakes of the duct in the present simulation using the hybrid turbulence model were very similar to that of LES. Also, the resistances acting on hull, rudder and duct in hybrid turbulence model were similar to those in RANS simulation whereas the viscous forces acting on the hull in LES had a significant error due to coarse cells inappropriate to the sub-grid scale model.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.2B
• /
• pp.249-259
• /
• 2008

The present paper studies the effect of the slope gradient of a fully permeable submerged breakwater using a newly developed numerical model that is able to consider the flow through a porous midium with inertial, laminar and turbulent resistance terms, i.e. simulate directly WAve-Structure (submerged breakwater)-Sand seabed interaction and can determine the eddy viscosity with LES turbulence model in 2-Dimensional wave field (LES-WASS-2D). The developed model was validated through the comparison with an existing experimental data, and further used for various numerical experiments in oder to investigate the complicated hydrodynamics on the varying slope gradient of permeable submerged breakwater. We found an acceptable phenomenon, as we expect intuitively, that reflection and transmission coefficients decrease simultaneously as slope gradient decrease. In addition, the breaking point, the circulation flow and mean vorticity around a submerged breakwater are throughly discussed.

• 한국신재생에너지학회:학술대회논문집
• /
• 2010.11a
• /
• pp.181.1-181.1
• /
• 2010

기후변화의 주요인이 되는 온실가스 감축을 목표로 화석연료를 대체하기 위한 대체 에너지 개발을 위한 많은 노력이 진행되고 있다. 풍력 에너지와 같은 신재생에너지는 이러한 하나의 해결 수단이 될 수 있으며 풍력 에너지 사업의 활성화를 위해서는 정확한 풍력 정보 제공이 우선이다. 풍력-기상자원지도는 풍력 발전에 유용한 정보 제공을 위한 목적으로 중규모 수치 모델을 이용하여 작성된다. 본 연구에서는 중규모 수치 모델의 정확도 향상을 위한 자료동화 방법으로써 Four-Dimensional Data Assimilation (FDDA) 방법을 이용한다. 풍력-기상자원지도는 공간분해능 1 km 해상도로 개발된다. 풍력-기상자원지도는 1998-2008년까지의 평균적인 상태에 대하여 모의를 하기 위하여 통계적인 방법으로 11년 기간의 평균과 유사한 기간을 선정하였다. 풍력-기상자원지도는 연 평균, 월 평균 풍속과 주 풍향, 주풍향 발생 비율 등의 정보를 제공한다. 우리나라 풍속의 평균 분포는 내륙 산악지역, 남해안, 제주도에서 강풍이 발생하며 주 풍향은 대체로 북서풍이다. 주 풍향의 발생비율은 산악 지역과 남 동해안에서 높아 풍력 발전의 최적지 정보를 제공한다. 1 km 해상도의 모델과 관측의 오차는 서해안 등의 해안지역보다 강원 산악지역에서 오차가 더욱 증가하였다. 이러한 산악 지역의 오차는 복잡한 지형에서는 1km 미만의 수 백 m 해상도 수치모의가 필요함을 지시한다. 따라서 본 연구에서는 WRF-LES 모형을 이용하여 333m 해상도의 기상자원지도를 개발한다. 333m 해상도의 자원지도 영역은 강원도 지역에 대하여 모의되었다. 333m 해상도의 풍속 분포는 1km 해상도의 풍속 분포와 비교하였을 때 풍속의 분포가 보다 세밀하게 표현되었다. 정량적인 검증을 하였을 때 관측소에 따라 차이는 있었으나 1km 해상도에서 과대 모의된 풍속의 분포가 현저히 개선이 되었으며, 시간적인 경향도 잘 일치함을 보였다.

• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.1-10
• /
• 2016

The analysis of characteristics of turbulent flow and thermal boundary layer for natural convection caused by fire along vertical wall is performed. The 4m-high vertical copper plate is heated and kept at a uniform surface temperature of $60^{\circ}C$ and the surrounding fluid (air) is kept at $16.5^{\circ}C$. The flow and temperature is solved by large eddy simulation(LES) of FDS code(Ver.6), in which the viscous-sublayer flow is calculated by Werner-Wengle wall function. The whole analyzed domain is assumed as turbulent region to apply wall function even through the laminar flow is transient to the turbulent flow between $10^9$<$Gr_z$<$10^{10}$ in experiments. The various grids from $7{\times}7{\times}128$ to $18{\times}18{\times}128$ are applied to investigate the sensitivity of wall function to $x^+$ value in LES simulation. The mean velocity and temperature profiles in the turbulent boundary layer are compared with experimental data by Tsuji & Nagano and the results from other LES simulation in which the viscous-sublayer flow is directly solved with many grids. The relationship between heat transfer rate($Nu_z$) and $Gr_zPr$ is investigated and calculated heat transfer rates are compared with theoretical equation and experimental data.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.16 no.9
• /
• pp.6305-6314
• /
• 2015

This study presents numerical simulations of mean flow and turbulence structure of an open channel with submerged vegetation. Filtered Navier-Stokes equations are solved using large-eddy simulation (LES). The immersed boundary method (IBM) is employed based on a Cartesian grid. The numerical result is compared with experimental data of Liu et al. (2008) and shows that simulated results coincided reasonably with experimental data within the average error of 10%. Strong vortices are generated at the interface between vegetated and non-vegetated regions with spanwise extent. The generation of turbulence induced by shear at the interface is interfered with wake turbulence, resulting turbulence intensity maximum. Turbulence produced by shear affects the flow in vegetated region and the penetration depth increases with an increase in the submergence ratio. This result can be used to understand sediment transport mechanisms in the vegetated region.