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

In the present work, we investigate the hydrodynamic behavior of a turbulent bore, such as tsunami bore and tidal bore, generated by the removal of a gate with water impounded on one side. The bore generation system is similar to that used in a general dam-break problem. In order to the numerical simulation of the formation and propagation of a bore, we consider the incompressible flows of two immiscible fluids, liquid and gas, governed by the Navier-Stokes equations. The interface tracking between two fluids is achieved by the volume-of-fluid (VOF) technique and the M-type cubic interpolated propagation (MCIP) scheme is used to solve the Navier-Stokes equations. The MCIP method is a low diffusive and stable scheme and is generally extended the original one-dimensional CIP to higher dimensions, using a fractional step technique. Further, large eddy simulation (LES) closure scheme, a cost-effective approach to turbulence simulation, is used to predict the evolution of quantities associated with turbulence. In order to verify the applicability of the developed numerical model to the bore simulation, laboratory experiments are performed in a wave tank. Comparisons are made between the numerical results by the present model and the experimental data and good agreement is achieved.

• Fire Science and Engineering
• /
• v.27 no.4
• /
• pp.27-34
• /
• 2013

The high prediction performance of fire detector models is essentially needed to assure the reliability of fire and evacuation modeling in the process of PBD (Performance Based fire safety Design). The main objective of the present study is to measure input information in order to predict the accurate activation time of smoke detector into a Large Eddy Simulation (LES) fire model such as FDS (Fire Dynamics Simulator). To end this, FDE (Fire Detector Evaluator) which can measure the device properties of detector was developed, and the input information of Heskestad and Cleary's models was measured for a ionization smoke detector. In addition, the activation times of smoke detectors predicted using default values into FDS and measured values in the present study were systematically compared. As a result, the device properties of smoke detector examined in the present study showed a significant difference compared to the default values used into FDS, which resulted in the considerable difference of up to 15 minutes or more in terms of the activation time of smoke detector. The database (DB) on device properties of various smoke and heat detectors will be built to improve the reliability of PBD in future studies.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.25 no.6
• /
• pp.374-385
• /
• 2013

Sediment particle motions have been numerically simulated over a sinusoidal ripple. Turbulent boundary layer flows are generated by Large Eddy Simulation, and the sediment particle motions are simulated using Lagrangian particle tracking method. Two unsteady flow conditions are used in the experiment by employing two different wave amplitudes while keeping other conditions such as wave period same. As expected, the amount of suspended sediment particles is clearly dependent on the wave amplitude as it is increasing with increasing flow intensity. However, it is also observed that the pattern of suspension may be different as well due to the only different condition caused by wave amplitude. Specially, the time of maximum sediment suspension within the wave period is not coincident between the two cases because sediment suspension is strongly affected by the existence of turbulent eddies that are formed at different times over the ripple between the two cases as well. The role of these turbulent eddies on sediment suspension is important as it is also confirmed in previous researches. However, it is also found the time of these eddies' formation may also dependent on the wave amplitude over rippled beds. Therefore, it has been proved that various flow as well as geometric conditions under waves has to be considered in order to have better understanding on the sediment suspension process over ripples. In addition, it is found that high turbulent energy and strong upward flow velocities occur during the time of eddy formation, which also supports high suspension rate at these time steps. The results indicate that the relationship between the structure of flows and bedforms has to be carefully examined in studying sediment suspension at coastal regions.

• Journal of the Korean Society of Safety
• /
• v.27 no.1
• /
• pp.20-25
• /
• 2012

The purpose of the study is to assess the extinguishing concentration of inert gases in engine nacelle fire. The experiment was performed with a two dimensional rectangular bluff body stabilized flames, where the fuel was ejected to counter flow and co-flow against an oxidizer stream. Two inert gases, $CO_2$ and $N_2$, were used for extinguishing agent in the oxidizer and methane was used for fuel. The main experimental parameters were the direction of injecting fuel, the kinds of agent and the velocity ratio between air and fuel streams, which controlled the mixing characteristic near bluff body and the strength of recirculation zone in the downstream. The result shows the flame structure and the mode were strongly dependent with fuel/air ratio and the fuel jet direction. For both flow configurations, the extinguishing concentration of $CO_2$ was smaller than the $N_2$ because of the large heat capacity of $CO_2$. However, the concentration of inert gasesat blowout was much smaller than those in the cup burner and coflow jet diffusion flames, which implies that the extinction mechanism of bluff body stabilized flames was mainly due to the aerodynamic aspect. Compared to co-flow fuel injection, the extinguishing concentration of inert gases under counter flow configuration was lower. The effect of direction might result from the mixing characteristic and strength of recirculation zonearound a bluff body. More details should be investigated for the characteristic of recirculation zone in the wake of bluff body using the LES(Large Eddy Simulation).

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1995.04a
• /
• pp.196-201
• /
• 1995

The broadband and discrete sources of sound in small cooling fans of propeller type and centrifugal type were investigated to understand the turbulent vortex structures from many bladed fans using ANSI test plenum for small air-moving devices (AMDs). The noise measurement method uses the plenum as a test apparatus to determine the acoustic source spectral density function at each operating conditions similar to real engineering applications based on acoustic similarity laws. The characteristics of fans including the head rise vs. volumetric flow rate performance were measured using a performance test facility. The sound power spectrum is decomposed into two non-dimensional functions: an acoustic source spectral distribution function F(St,.phi.) and an acoustic system response function G(He,.phi.) where St, He, and .phi. are the Strouhal number, the Helmholtz number, and the volumetric flow rate coefficient, respectively. The autospectra of radiated noise measurements for the fan operating at several volumetric flow rates,.phi., are analyzed using acoustical similarity. The rotating stall in the small propeller fan with a bell-mouth guided is mainly due to a leading edge separation. It creates a blockage in the passage and the reduction in the flow rate. The sound power levels with respect to the rotational speeds were measured to reveal the mechanisms of stall and/or surge for different loading conditions and geometries, for example, fans installed with a impinging plate. Lee and Meecham (1993) studied the effect of the large-scale motions like impinging normally on a flat plate using Large-Eddy Simulation(LES) and Lighthill's analogy.[ASME Winter Annual Meeting 1993, 93-WA/NCA-22]. The dipole and quadrupole sources in the fans tested are shown closely related to the vortex structures involved using cross-correlations of the hot-wire and microphone signals.

• Proceedings of the KSR Conference
• /
• 2008.06a
• /
• pp.1794-1801
• /
• 2008

Fire-driven flow and temperature distribution in a ventilated tunnel was analyzed by Large Eddy Simulation using FDS code. The simulated tunnel is 182m length, 5.4m wide and 2.4m height. A pool fire was located 112m from tunnel entrance and was taken as a heat source of $0.89m^2$. The heat is assumed to be released uniformly throughout the whole simulated time. The fire strength was 2.76MW and the fuel burnt was octane. The parallel computational method was employed to accelerate the computing time and manage the large grid points which is not possible to handle in the one CPU. The total grid points used were $2.4{\times}10^6$ and 7 CPUs were used to calculate the momentum and energy equations. The simulated results were well compared with the experiments.

• Journal of Power System Engineering
• /
• v.13 no.6
• /
• pp.13-21
• /
• 2009

This basic study is required to examine spray or jet behavior depending on fuel phase. In this study, analyses of diesel fuel(n-Tridecane, $C_{13}H_{28}$) spray and natural gas fuel(Methane, $CH_4$) jet under high temperature and pressure are performed by a general-purpose program, ANSYS CFX release 11.0, and the results of these are compared with experimental results of diesel fuel spray using the exciplex fluorescence method. The simulation results of diesel spray is analyzed by using the combination of Large-Eddy Simulation(LES) and Lagrangian Particle Tracking(LPT) and of a natural gas jet is analyzed by using Multi-Component Model(MCM). There are two study variables considered, that is, ambient pressure and injection pressure. In a macroscopic analysis, the higher ambient pressure is, the shorter spray or jet tip penetration is at each time after start of injection. And the higher injection pressure is, the longer spray or jet tip penetration is at each time after start of injection. When liquid fuel is injected, droplets of the fuel need some time to evaporate. However, when natural gas fuel is injected, the fuel does not need time to evaporate. Gas fuel consists of minute particles. Therefore, the gas fuel is mixed with the ambient gas more quickly at the initial time of injection than the liquid fuel is done. The experimental results also validate the usefulness of this analysis.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2011.05a
• /
• pp.37-37
• /
• 2011

태풍 루사나 매미에 의한 피해사례를 보면 하천 횡단 구조물 본체가 파괴되는 피해 뿐 만 아니라, 구조물과 제방과의 연결부가 세굴되어 붕괴되는 사례가 많이 발생하였다. 하천설계기준(2009)에는 이러한 보와 제방의 연결부 부분을 연결호안이라 하여 관련 기준을 제시하고 있으나, 설치구간의 길이를 정할 때 하천의 규모나 하도의 특성을 고려하지 못하고 일률적으로 결정하도록 하고 있다. 이에 건설기술연구원에서는 상류 Froude 수와 보 높이를 변수로 하는 연결호안 설치 길이에 대한 실험식(2006)을 제시하였다. 하지만 보 접속부 세굴에 영향을 주는 인자는 상류의 Froude 수와 보의 높이 뿐만 아니라, 상 하류 수위차, 보의 형상, 제방의 재료적 특성 등 많은 요인들이 존재한다. 이 중 본 연구에서는 상 하류 수위차와 보 접속부 세굴 범위의 관계를 파악하고자 한다. 보 접속부 세굴에 대한 수치모의에 앞서 Gill(1972)와 Dongol(1990)이 수행하였던 교대세굴 자료를 바탕으로 Flow-3D의 세굴에의 적합성 및 적용성을 검토하였다(그림 2참고). 검토 결과 수치모의의 입력 변수를 조정함으로써 실험값과 가깝게 나타나는 것을 확인하였다. 이렇게 조정된 변수를 이용하여, 보 접속부 세굴에 대한 수치모의는 상 하류의 수위차에 대해 다양한 검토를 하기 위해 상류의 수위를 1.0m로 고정한 채하류부의 수위를 0.2~0.8m까지 0.1m간격으로 변화시켰으며, 유사의 대표 입경은 0.63mm로 주문진 표준사를 가정하였다. 모의 시간은 평형세굴심이 발생할 때까지 수행하였고, 난류모델로는 LES(Large Eddy Simulation)를 채택하였다. 모의 결과 상 하류 수위차가 증가할수록 보 하류부 세굴 길이($L_d$)가 증가 하였으며, 세굴의 폭($B_d$) 역시 증가하였다. 그림 3은 수치모의 결과로 세굴이 발생된 제방의 모습을 도시한 것이다.

• Wind and Structures
• /
• v.36 no.3
• /
• pp.161-174
• /
• 2023

The analytical model of tornado vortices plays an essential role in tornado wind description and tornado-resistant design of civil structures. However, there is still a lack of guidance for the selection and application of tornado analytical models since they are different from each other. For single-cell tornado vortices, this study conducts a comparative study on the velocity characteristics of the analytical models based on numerically simulated tornado-like vortices (TLV). The single-cell stage TLV is first generated by Large-eddy simulations (LES). The spatial distribution of the three-dimensional mean velocity of the typical analytical tornado models is then investigated by comparison to the TLV with different swirl ratios. Finally, key parameters are given as functions of swirl ratio for the direct application of analytical tornado models to generate full-scale tornado wind field. Results show that the height of the maximum radial mean velocity is more appropriate to be defined as the boundary layer thickness of the TLV than the height of the maximum tangential mean velocity. The TLV velocity within the boundary layer can be well estimated by the analytical model. Simple fitted results show that the full-scale maximum radial and tangential mean velocity increase linearly with the swirl ratio, while the radius and height corresponding to the position of these two velocities decrease non-linearly with the swirl ratio.

• Wind and Structures
• /
• v.23 no.2
• /
• pp.127-142
• /
• 2016

Aerodynamic effects, such as drag force and flow-induced vibration (FIV), on civil engineering structures can be minimized by optimally modifying the structure shape. This work investigates the turbulent wake of a square prism with its faces modified into a sinusoidal wave along the spanwise direction using three-dimensional large eddy simulation (LES) and particle image velocimetry (PIV) techniques at Reynolds number $Re_{Dm}$ = 16,500-22,000, based on the nominal width ($D_m$) of the prism and free-stream velocity ($U_{\infty}$). Two arrangements are considered: (i) the top and bottom faces of the prism are shaped into the sinusoidal waves (termed as WSP-A), and (ii) the front and rear faces are modified into the sinusoidal waves (WSP-B). The sinusoidal waves have a wavelength of $6D_m$ and an amplitude of $0.15D_m$. It has been found that the wavy faces lead to more three-dimensional free shear layers in the near wake than the flat faces (smooth square prism). As a result, the roll-up of shear layers is postponed. Furthermore, the near-wake vortical structures exhibit dominant periodic variations along the spanwise direction; the minimum (i.e., saddle) and maximum (i.e., node) cross-sections of the modified prisms have narrow and wide wakes, respectively. The wake recirculation bubble of the modified prism is wider and longer, compared with its smooth counterpart, thus resulting in a significant drag reduction and fluctuating lift suppression (up to 8.7% and 78.2%, respectively, for the case of WSP-A). Multiple dominant frequencies of vortex shedding, which are distinct from that of the smooth prism, are detected in the near wake of the wavy prisms. The present study may shed light on the understanding of the underlying physical mechanisms of FIV control, in terms of passive modification of the bluff-body shape.