• Fire Science and Engineering
• /
• v.15 no.1
• /
• pp.7-15
• /
• 2001

This paper describes the smoke filling process of a fire field model based on a self-developed SMEP (Smoke Movement Estimating Program) codo to the simulation of fire induced flows in the atrium space (SIVANS atrium at Japan) containing smoke radiation effect. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k-$\varepsilon$ turbulence model with buoyancy term. Also it solves the radiation equation using the discrete ordinates method. The result of the calculated smoke temperature containing radiation effect has shown a better prediction than the result calculated by only convection effect in comparison with the experimental data. This seems to come from the radiation effect of $H_2$O and $CO_2$ gas under smoke productions. Thus, the consideration of the radiation effect under smoke in fire should be necessary in order to get more realistic result. Also the numerical results indicated that the smoke layer is developing at a rate of about 0.1 m/s. It would take about 450 seconds after starting the ultra fast fire of 560 kW that the smoke layer move down to 1.5m above the escape level.

• Fire Science and Engineering
• /
• v.14 no.3
• /
• pp.6-12
• /
• 2000

This paper describes the smoke movement of a fire field model based on a self-developed SMEP(Smoke Movement Estimating Program) code to the simulation of fire induced flows in the two types of compartment space containing the radiation effect under smoke movement in room fire. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k-$\varepsilon $ turbulence model with buoyancy term. Also it solves the radiation equation using the discrete ordinates method. The result of the calculated smoke temperature containing radiation effect has shown reasonable agreement compared with the experimental data. On the other hand, a difference of a lot was found between the temperature predicted by the SMEP with only convection effect and obtained by the experimental result. This seems to come from the radiation effect of $H_2$O and $CO_2$ gas under smoke productions. Thus, the consideration of the radiation effect under smoke in fire may be necessary in order to produce more realistic result.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.13 no.1
• /
• pp.12-17
• /
• 2010

Oscillating Water Column (OWC) device has been widely employed in the wave energy conversion. Wave Focusing Device (WFD) is proposed to be helpful for improving the operating performance of OWC chamber. In the present paper, a Numerical Wave Tank (NWT) using two-phase VOF model is utilized to simulate the generation and propagation of incident regular waves, water column oscillation inside the chamber. The NWT consists of the continuity equation, Reynolds-averaged Navier-Stokes equations and two-phase VOF functions. The standard k- turbulence model, the finite volume method, NITA-PISO algorithm and dynamic mesh technique are employed. Effects of WFD on the operating performance of OWC chamber are investigated numerically.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.4 no.5
• /
• pp.130-141
• /
• 1996

Numerical study of a confined, swirling, isothermal and spray-combusting flows has been presented. The pressure-velocity coupling in the Eulerian gas-phase equation is handled by the improved PISO algorithm. The droplet dispersion by turbulence is introduced by a Stochastic Separated Flow(SSF) model. The k-$\varepsilon$ turbulence model and the eddy dissipation model are employed to account for turbulence-combustion interaction. The detailed comparison with experimental data has been made for the isothermal jet swirling flows and the nearly monodisperse spray-combusting flow in the swirl combustor.

• 한국연소학회:학술대회논문집
• /
• 2002.06a
• /
• pp.95-98
• /
• 2002

A numerical investigation on the flame propagation and extinction in a micro combustor is described. Previous measurements of $H_2-air$ flame propagation in a submilimeter scale combustor exhibited significance of wall effects on burning velocity and extinction. The heat transfer to wall becomes important not only in the cooling of burnt gases but also during the flame ropagation, which has be by and large ignored in macro scale combustor calculations. In order to take the heat loss into account the combustion calculation, we developed a numerical code with a heat transfer model that was determined empirically from measured data. PISO algorithm was used for differencing of conservation equations. $H_2-air$ reaction was modeled with 10 species - 16 steps. Comparison with measured data showed good agreement in flame propagation speed. Also the pressure decrease after flame extinction was accurately predicted by the model. A further study is desirable for a better quenching model that can predict the quenching location.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.169-174
• /
• 1995

A finite element scheme using the concept of PISO method has been developed to solve the viscous flow problems in all speed range. In this study, new pressure equation is proposed such that both the hyperbolic term related with the density variations and elliptic term reflecting the incompressibility constraint are included. Present method has been applied to incompressible flow in two-dimensional driven cavity(Re=100, 400 and 1,000), and its computed results are compared with other's. Also, Carter plate problem(M=3 and Re=1,000) is computed and the comparison is made with Carter's results. Finally, we simulate a shock-boundary layer interaction problem(M=2 and $Re=2.96{\times}10^5$) to illustrate the shock capturing capability of the present solution algorithm.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.524-529
• /
• 2010

In the present study, a CFD program is developed for the Fluid-Structure Interaction(FSI) analysis. The non-staggered, non-orthogonal, and unstructured grid system was also used to handle the complicated geometries in the program. In order to validate the capabilities of the developed CFD program, various models are investigated by using unstructured and nonorthogonal meshes. The predicted results are a good agreement with analytic solution, experimental data and commercial software. And also PISO algorithm is applied for transient flow analysis. The cyclic boundary condition and baffle cell are developed in order to improve the effectiveness of the calculation for complex geometry.

• Journal of Power System Engineering
• /
• v.13 no.5
• /
• pp.11-17
• /
• 2009

A numerical study of unsteady mixed convection in a cavity with high viscous fluid is presented. Finite volume method was employed for the discretization and PISO algorithm was used for calculating pressure term. The parameters governing the problem are the Rayleigh number ($10^3\;{\leq}\;Ra\;{\leq}\;10^5$), the Reynolds number (0 < Re $\leq$ 1), and the aspect ratio (0.5 $\leq$ AR $\leq$ 2). The fluid used is silicon oil, a high prandtl number fluid, Pr = 909.1. The results show velocity vectors and temperature distributions. It is found that the periodic flows in a cavity are observed at very low Reynolds numbers, and the period of periodic flow decreases with increasing Reynolds and Rayleigh numbers, and increases with increasing aspect ratio. Also, the Reynolds number range of periodic flow increases with increasing Rayleigh numbers and aspect ratio.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2005.06a
• /
• pp.210-215
• /
• 2005

A three-dimensional numerical simulation is performed to investigate on a low Reynolds number mixed convection in a horizontal rectangular channel with the upper part cooled and the lower part heated uniformly. The three-dimensional governing equations are solved using a finite volume method. For convective term, the central differencing scheme is used and for the pressure correction, the PISO algorithm is used. Solutions are obtained for A=4, Pr=0.72, 10, 909, the Reynolds number ranging from $2.1{\times}10^{-2}$ to $1.2{\times}10^{-1}$, the Rayleigh number is $3.5{\times}10^4$. It is found that vortex roll structures of mixed convection in horizontal rectangular channel can be classified into three roll structures which affected by Prandtl number and Reynolds number.

• Journal of Power System Engineering
• /
• v.11 no.4
• /
• pp.56-64
• /
• 2007

A numerical simulation is carried out mixed convection in horizontal channel with a heat source from below of rectangular cavity. Finite volume method was employed for the discretization and PISO algorithm was used for calculating pressure term. The parameters governing the problem are the Reynolds number ($10^{-2}{\leq}Re{\leq}50$), the Rayleigh number ($10^3{\leq}Ra{\leq}2.06{\times}10^5$), the Prandtl number ($0.72{\leq}Pr{\leq}909$), the aspect ratio ($0.5{\leq}AR=W/H{\leq}2$) and the angle of inclination ($0^{\circ}{\theta}60^{\circ}$). Mean Nusselt number distributions were obtained and effect of Reynolds number, Rayleigh number and Prandtl number on mixed convection in the horizontal channel with rectangular cavity were investigated.