• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.20 no.2
• /
• pp.137-143
• /
• 2008

In the present study, a forest fire spread was simulated with a three-dimensional, fully-transient, physics-based, computer simulation program. Physics-based fire simulation is based on the governing equations of fluid dynamics, combustion and heat transfer. The focus of the present study is to perform parametric study to simulate fire spread through flat and inclined wildland with vegetative fuels like trees or grass. The fire simulation was performed in the range of the wind speeds and degrees of inclination. From the results, the effect of the various parameters of the forest fire on the fire spread behavior was analyzed for the future use of the simulation in the prediction of fire behavior in the complex terrain.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.09b
• /
• pp.842-845
• /
• 2005

A method for the numerical simulation of two-dimensional free-surface flow resulting from the propagation of regular gravity waves over topography with arbitrary bottom shape is presented. The method is based on the numerical solution of the Euler equations subject to the fully nonlinear free-surface boundary conditions and the appropriate bottom, inflow and outflow conditions using a hybrid finite-differences and spectral-method scheme. The formulation includes a boundary-fitted transformation, and is suitable for extension to incorporate large-eddy simulation (LES) and large-wave simulation (LWS) terms for turbulence and breaking wave modeling, respectively. Results are presented for the simulation of the free-surface flow over two different bottom topographies, with constant slope values of 1:10 and 1:20, two different inflow wave lengths and two different inflow wave heights. An absorption outflow zone is utilized and the results indicate minimum wave reflection from the outflow boundary. Over the bottom slope, lengths of waves in the linear regime are modified according to linear theory dispersion, while wave heights remain more or less unchanged. For waves in the nonlinear regime, wave lengths are becoming shorter, while the free surface elevation deviates from its initial sinusoidal shape.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.392-397
• /
• 2001

A numerical study is performed to investigate the interaction between subsonic axial turbine blade boundary layer and periodically oncoming rotor induced wakes. An implicit scheme for solving the compressible Navier-Stokes equation is developed, which adopts a 4th-order compact difference for spatial discretiztion, a 2nd order Crank-Nicolson scheme for temporal discretization and the dynamic eddy viscosity model as the subgrid scale model. The efficiency and the accuracy of the proposed method are verified by applying to some benchmark problems such as laminar cylinder flow, laminar airfoil cascade flow and a transitional flat plate boundary layer flow. Computational results show good agreements with previous experimental and numerical results. Finally, flow through a stator cascade is simulated at $Re = 7.5{\times}10^5$ without free-stream turbulence intensity. The velocity fields and skin friction coefficients in the transitional region show similar trends with previous boundary layer natural transition.

• Fire Science and Engineering
• /
• v.17 no.2
• /
• pp.56-61
• /
• 2003

The smoke removal rate from a room with an opening was investigated for different smoke control systems by using the large eddy simulation turbulence model of the Fire Dynamics Simulator. The decreasing rate of the particles randomly distributed in the 2m X 2m X 2.4m room was com-pared for the ventilation system, pressurization system and extraction system, and for the air flowrate of the ventilation system. Difference in the smoke removal rate among the three smoke control systems was small when the opening was closed. The pressurization system showed less smoke removal rate than the other two systems when the opening existed, and hence is not recommended for subway stations with large openings. It was also shown that a less flowrate in the ventilation system leads to a much longer smoke removal time.

• Journal of Ocean Engineering and Technology
• /
• v.29 no.1
• /
• pp.9-15
• /
• 2015

We investigated the Reynolds number effects on the flow over a twisted offshore structure in the range of 3×10³≤ Re ≤ 1 × 10⁴. To analyze the effect of the twisted surface treatment, a large eddy simulation (LES) with a dynamic subgrid model was employed. A simulation of the cylindrical structure was also carried out to compare the results with those of the twisted offshore structure. As Re increased, the mean drag and lift coefficient of the twisted offshore structure increased with the same tendency as those of the cylindrical structure. However, the increases in the mean drag and lift coefficient of the twisted offshore structure were much smaller than those of the cylindrical structure. Furthermore, elongated shear layer and suppressed vortex shedding from the twisted offshore structure occurred compared to those of the cylindrical cylinder, resulting in a drag reduction and suppression of the vortex-induced vibration (VIV). In particular, the twisted offshore structure achieved a significant reduction of over 96% in VIV compared with that of the cylindrical structure, regardless of increasing Re. As a result, we concluded that the twisted offshore structure effectively controlled the flow structures with reductions in the drag and VIV compared with the cylindrical structure, irrespective of increasing Re.

• Wind and Structures
• /
• v.26 no.4
• /
• pp.215-229
• /
• 2018

The appearance of a rivulet at the upper surface of a stay cable is responsible for rain-wind-induced vibration (RWIV) of cables of cable-stayed bridges. However, the formation mechanism of the upper rivulet and its aerodynamic effects on the stay cable has not been fully understood. Large eddy simulation (LES) method is used to investigate flow around and aerodynamics of a circular cylinder with an upper rivulet at a Reynolds number of 140,000. Results show that the mean lift coefficients of the circular cylinder experience three distinct stages, zero-lift stage, positive-lift stage and negative-lift stage as the rivulet located at various positions. Both pressure-induced and friction-induced aerodynamic forces on the upper rivulet are helpful for its appearance on the upside of the stay cable. The friction-induced aerodynamic forces, which have not been considered in the previous theoretical models, may not be neglected in modeling the RWIV. In positive-lift stage, the shear layer separated from the upper rivulet can reattach on the surface of the cylinder and form separation bubbles, which result in a high non-zero mean lift of the cylinder and potentially induces the occurrence of RWIV. The separation bubbles are intrinsically unsteady flow phenomena. A serial of small eddies first appears in the laminar shear layer separated from the upper rivulet, which then coalesces and reattaches on the side surface of the cylinder and eventually sheds into the wake.

• Wind and Structures
• /
• v.24 no.1
• /
• pp.1-24
• /
• 2017

Wind environment in urban residential areas is an important index to consider when evaluating the living environment. However, due to the complexity of the flow field in residential areas, it is difficult to specify the correct inflow boundary conditions in the large eddy simulation (LES). In this paper, the weighted amplitude wave superposition (WAWS) is adopted to simulate the fluctuating velocity data, which satisfies the desired target wind field. The fluctuating velocity data are given to the inlet boundary of the LES by developing an UDF script, which is implemented into the FLUENT. Then, two numerical models - the empty numerical wind tunnel model and the numerical wind tunnel model with spires and roughness elements are established based on the wind tunnel experiment to verify the present method. Finally, the turbulence generation approach presented in this paper is used to carry out a numerical simulation on the wind environment in an urban residential area in Lisbon. The computational results are compared with the wind tunnel experimental data, showing that the numerical results in the LES have a good agreement with the experimental results, and the simulated flow field with the inlet fluctuations can generate a reasonable turbulent wind field. It also shows that strong wind velocities and turbulent kinetic energy occur at the passageways, which may affect the comfort of people in the residential neighborhood, and the small wind velocities and vortexes appear at the leeward corners of buildings, which may affect the spreading of the pollutants.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.9 no.3
• /
• pp.49-59
• /
• 2005

Turbulent heat transfer over a fully-developed wavy channel is investigated by a turbulence model. The nonlinear k- f - f$_{ model of Park et at.[1] is slightly modified and their explicit algebraic heat flux model is employed. The Reynolds number is fixed at Re$_{b}$=6760 and the wave configuration is varied in the range of 0 $\leq$ $\alpha$/$\lambda$$\leq$0.15 and 0.25 $\leq$A/H$\leq$4.0. In order to verify model performances, a large eddy simulation is performed for the selected cases. The model performance is shown to be generally satisfactory. By using k- $\varepsilon$ - f$_{ model, the enhancement of heat transfer and the characteristics of turbulent flow in wavy wall are investigated. Finally, the influence of wavy configuration on heat transfer is scrutinized.

• Journal of Korean Society for Atmospheric Environment
• /
• v.21 no.5
• /
• pp.525-535
• /
• 2005

Using a three-dimensional computational fluid dynamics (CFD) model with the $k-{\varepsilon}$ turbulence closure scheme based on the renormalization group theory, flow regimes in urban street canyons are classified according to the building and street aspect ratios. The transition between skimming flow (SF) and wake interference flow (WIF) is determined with the size of double-eddy circulation generated behind the upwind building. The transition between WIF and isolated roughness flow (IRF) is determined with the flow reattachment distance from the upwind building. The critical aspect ratios at which the flow transition occurs are found and compared with those in previous studies. The results show that the flow-regime classification method used in this study is quite reasonable and that the values of the critical aspect ratios are generally consistent with those in fluid experiments or large-eddy simulation. The regression equation describing a relation between the building and street aspect ratios at the flow-regime transition is presented.

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