• Journal of ILASS-Korea
• /
• v.26 no.3
• /
• pp.127-134
• /
• 2021

Fire sprinkler initial spray was analyzed by Large eddy simulation (LES) and Volume of Fluid (VoF) integrated method. The IsoAdvector geometric VoF was used to identify the liquid-gas interface clearly even with the large Courant-Friedrichs-Lewy number. To reduce the computational costs, sector meshes and Adaptive Mesh Refinement up to level 3 were used. Base mesh size was 1 mm, which is roughly equivalent to the initial sprinkler droplet. Top surface radius of boss and deflector size were modified to investigate the effects of sprinkler head design on primary breakup process. When top surface radius of boss was increased, vertical liquid sheet was formed. This phenomenon reduced the sheet breakup radius, sheet thickness and velocity. Due to reduced liquid sheet thickness, a large amount of ligaments was created from the liquid sheet. As a result, there was a dramatic decrease in volume per surface area, indicating an increase in breakup process. Spray pattern viewed in radial direction also changed when top surface radius of boss increased. When top surface radius of boss was increased, a T-shaped pattern was observed while a V-shaped pattern was observed in all other cases. When the deflector size increases, the spray pattern remains V-shaped, even if the top surface radius of boss increased. Further studies on promoting atomization of the water supplied to the lower part of the sprinkler head in the T-shape pattern should be conducted.

• Atmosphere
• /
• v.27 no.4
• /
• pp.445-453
• /
• 2017

This study investigates turbulent flow around a square cylinder mounted on a flat surface at high Reynolds number using a large-eddy simulation (LES) model, particularly focusing on vortex streets behind the square cylinder. Total 9 simulation cases with different inflow wind directions, inflow wind speeds, and cylinder widths in the x- and y-directions are considered to examine the effects of inflow wind direction, speed, and cylinder widths on turbulent flow and vortex streets. In the control case, the inflow wind parallel to the x-direction has a maximum speed of $5m\;s^{-1}$ and the width and height of the cylinder are 50 m and 200 m, respectively. In all cases, down-drafts in front of the cylinder and updrafts, wakes, and vortex streets behind the cylinder appear. Low-speed flow below the cylinder height and high-speed flow above it are mixed behind the cylinder, resulting in strong negative vertical turbulent momentum flux at the boundary. Accordingly, the magnitude of the vertical turbulent momentum flux is the largest near the cylinder top. In the case of an inflow wind direction of $45^{\circ}$, the height of the boundary is lower than in other cases. As the inflow wind speed increases, the magnitude of the peak in the vertical profile of mean turbulent momentum flux increases due to the increase in speed difference between the low-speed and high-speed flows. As the cylinder width in the y-direction increases, the height of the boundary increases due to the enhanced updrafts near the top of the cylinder. In addition, the magnitude of the peak of the mean turbulent momentum flux increases because the low-speed flow region expands. Spectral analysis shows that the non-dimensional vortex generation frequency in the control case is 0.2 and that the cylinder width in the y-direction and the inflow wind direction affect the non-dimensional vortex generation frequency. The non-dimensional vortex generation frequency increases as the projected width of the cylinder normal to the inflow direction increases.

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

• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.98-98
• /
• 2019

중력류 또는 밀도류는 주변 유체에 비해 상대적으로 밀도가 큰 유체가 밀도차에 의한 추진력으로 흐르는 것이다. 중력류의 수치모델링에는 두 가지 어려움이 있다. 즉, 적합한 지배방정식을 구성하여 적용하는 것 그리고 난류의 영향을 합리적으로 반영하는 것이다. 기존 중력류 해석을 위한 지배방정식들은 유체의 연속방정식과 운동량 방정식 그리고 밀도 또는 농도의 이송방정식을 조합하여 구성된다. 이들 지배방정식을 이용한 연구들은 대부분 두 유체 사이의 밀도차가 충분히 작아서 밀도 변동(variations)의 영향은 오로지 부력항에서만 유지된다는 Boussinesq 근사에 근거를 둔다. 그리고 이송방정식에서 밀도 또는 농도의 확산계수을 점성계수의 함수로 표현하기 위해서 Schmidt 수를 이용한다. 수치모델링에서 Schimdt 수는 상수값을 적용하지만, 이 값은 밀도의 연직방향 경사에 근거한 부력빈도(buoyancy frequency)와 난류량의 따라 큰 차이를 보이는 것으로 알려져있다. 한편, 표준 통계학적 난류모델과 벽함수를 적용한 수치모델링은 초기 중력에 의해서 무너지는(slumping) 단계를 넘어 관성력으로 추진되는 단계와 점성 효과가 지배적인 단계에서는 정확도에 현저히 낮아지기 때문에 대부분 큰와모의(large-eddy simulation, LES) 또는 DNS(direct numerical simulation)수준의 고해상도(high-resolution) 해석기법을 적용하여 공학적인 문제에 적용하는 데는 한계가 있다. 이 연구에서는 Boussinesq 근사와 Schmidt 수를 사용하지 않으며, LES 보다 적용이 용이한 DES (detached-eddy simulation)기법을 조합한 다상흐름 수치모델을 적용하여 중력류를 해석을 시도하였다. 수치해석결과를 실험값과 함께 기존 수치모델링 기법으로 구한 수치해와 비교분석하여 이 연구에서 개발 및 적용된 수치모델링 기법의 적용성을 평가한다.

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

• Wind and Structures
• /
• v.38 no.2
• /
• pp.129-146
• /
• 2024

Aiming at the problem of non-stationary wind field simulation of downbursts, a non-stationary down-burst generation system was designed by adding a nozzle and program control valve to the inlet of the original wall jet model. The computational fluid dynamics (CFD) method was used to simulate the downburst. Firstly, the two-dimensional (2D) model was used to study the outflow situation, and the database of working conditions was formed. Then the combined superposition of working conditions was carried out to simulate the full-scale measured downburst. The three-dimensional (3D) large eddy simulation (LES) was used for further verification based on this superposition condition. Finally, the wind tunnel test is used to further verify. The results show that after the valve is opened, the wind ve-locity at low altitude increases rapidly, then stays stable, and the wind velocity at each point fluctuates. The velocity of the 2D model matches the wind velocity trend of the measured downburst well. The 3D model matches the measured downburst flow in terms of wind velocity and pulsation characteris-tics. The time-varying mean wind velocity of the wind tunnel test is in better agreement with the meas-ured time-varying mean wind velocity of the downburst. The power spectrum of fluctuating wind ve-locity at different vertical heights for the test condition also agrees well with the von Karman spectrum, and conforms to the "-5/3" law. The vertical profile of the maximum time-varying average wind veloci-ty obtained from the test shows the basic characteristics of the typical wind profile of the downburst. The effectiveness of the downburst generation system is verified.

• Fire Science and Engineering
• /
• v.32 no.3
• /
• pp.42-50
• /
• 2018

The effect of sub-grid turbulence and combustion models on the mean flame height in a heptane pool fire according to the Fire Dynamics Simulator (FDS) version (5 and 6) based on Large Eddy Simulation (LES) was examined. The heat release rate for the fire simulation was provided through experiments performed under identical conditions and the predictive performance of the mean flame height according to FDS version was evaluated by a comparison with the existing correlation. As a result, the Smagorinsky and Deardorff turbulence models applied to FDS 5 and 6, respectively, had no significant effects on the mean flow field, flame shape and flame height. On the other hand, the difference in pool fire characteristics including the mean flame height was due mainly to the difference in the mixture fraction and Eddy Dissipation Concept (EDC) combustion models applied to FDS 5 and 6, respectively. Finally, compared to FDS 6, FDS 5 provided the predictive result of a significantly longer flame height and more consistent mean flame height than the existing correlation.

• Wind and Structures
• /
• v.20 no.1
• /
• pp.37-58
• /
• 2015

As wind flows around a sharp-edged body, the resulting separated flow becomes complicated, with multiple separations and reattachments as well as vortex recirculation. This widespread and unpredictable phenomenon has long been studied academically as well as in engineering applications. In this study, the flow characteristics around rectangular prisms with five different aspect ratios were determined through wind tunnel experiments and a detached eddy simulation, that placed the objects in a simulated deep turbulent boundary layer at $Re=4.6{\times}10^4$. A series of rectangular prisms with the same height (h = 80 mm), different longitudinal lengths (l = 0.5h, h, and 2h), or different transverse widths (w = 0.5h, h, and 2h) were employed to observe the effects of the aspect ratio. Furthermore, five wind directions ($0^{\circ}$, $10^{\circ}$, $20^{\circ}$, $30^{\circ}$, and $45^{\circ}$) were selected to observe the effects of the wind direction. The simulated results of the surface pressure were compared to the wind tunnel experiment results and the existing results of previous papers. The vortex and spectrum were also analyzed to determine the detailed flow structure around the body. The paper also highlights the pressure distribution around the rectangular prisms with respect to the different aspect ratios. With an increasing transverse width, the surface suction pressure on the top and side surfaces becomes stronger. In addition, depending on the wind direction, the pressure coefficient experiences a large variation and can even change from a negative to a positive value on the side surface of the cube model.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.547-558
• /
• 2008

With the advent of super cargo ship due to the explosive increase in the amount of cargo shipped via seas, some mega ports are under construction in South Korea, to accommodate the super cargo ship, and some of them already enter their final phase. To sustain the harbor tranquility, mega ports usually comprise huge vertical type breakwaters which are intrinsically vulnerable to the attack of long waves. In this rationale, we present the chamber type breakwater with a circular curtain wall - Eco-breaker 2, to alleviate the reflection of long waves and numerically investigate the hydraulic characteristics of Eco-breaker 2. As a wave driver, we use the Navier-Stokes eq., the most robust wave driver, using SPH (Smoothed Particle Hydrodynamics) and LES (Large Eddy Simulation). For the verification of numerical results, we also carried out hydraulic model test. It is shown that Eco-breaker 2 can effectively alleviate the reflection of long waves with its inherited large organized eddies encompassing the water chamber and some region off the curtain wall of varying size. It is also shown that the scope and strength of large organized eddies strongly depends on the incident wave period, and the reflection coefficient can be lowered to 0.18 by tuning the size of water chamber such that resident time at the chamber is just short of the half period of incident waves. Based on these results, we present the specification of Eco-breaker 2 to boost its use on the development of water environment friendly harbor worldwide.