• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.26-34
• /
• 1996

The performance of the SIMPLE, SIMPLE-C and PISO algorithms for the treatment of the pressure-velocity coupling in fluid flow problems were examined by comparing the computational effort required to obtain the same level of the convergence. Example problems are circular duct and 90-degree bent square-duct. For circular duct case, laminar and turbulent flow were computed. For 90-degree bent square-duct case, laminar flow was simulated by the time-marching method as well as the iterative method. The convergence speed of the other two algorithms are not always superior to SIMPLE algorithm. SIMPLE algorithm is faster than SIMPLE-C algorithm in the simple laminar flow calculations. The application of the PISO algorithm in three dimensional general coordinates is not so effective as in two-dimensional ones. Since computational time of PISO algorithm is increased at each time step(or iterative step) in three dimension, the total convergence speed is not decreased. But PISO algorithm is stable for large time step by using time marching method,.

• Journal of computational fluids engineering
• /
• v.8 no.3
• /
• pp.12-20
• /
• 2003

A numerical procedure for the calculation of solidification and melting phase change using PISO algorithm is presented. In case of phase change problem, the coupling between velocity/pressure/temperature and liquid fraction is important. The converged temperature and liquid fraction solution which satisfies the energy balance is acquired by applying enthalpy method into inner iteration in matrix solver. And a modified PISO algorithm version is introduced to properly solve the coupling between velocity/pressure/temperature and liquid fraction. A comparison of the proposed procedure with a standard iterative method shows improvement both in terms of computing speed and robustness.

• Journal of computational fluids engineering
• /
• v.7 no.1
• /
• pp.1-9
• /
• 2002

In this paper, a modification of PISO algorithm based on standard k-ε turbulence model was proposed. The numerical technique used in this research is finite volume method, hybrid scheme for discretizing convection term, Euler implicit scheme for discretizing time term, and non-staggered grid. The basic idea of the modification of PISO algorithm is to perform an additional corrector stage for turbulence kinetic energy and dissipation rate to correct the inconsistence of flow and turbulence. In order to validate this algorithm, simulation of flow around a square cylinder (Re=3000) was performed in two-dimensional case. The results obtained from the proposed scheme show better agreement with those from the experiment than using original PISO algorithm in coherent velocity field.

• Fire Science and Engineering
• /
• v.13 no.1
• /
• pp.21-30
• /
• 1999

In this paper, the smoke filling process for the three types of atrium spaces are simulated u using the two types of deterministic fire model; zone models and field models. The zone models u used are the FffiST, CFAST, and CCFM.VENTS m떠els develo야퍼 at the Building and Fire R Research Laboratories, NIST, USA and the NBTC one-room model of FIR.ECAIι delveloped at C CSffiO, Austr;외ia. The field models used are the fire field model developed by W. K Chow and a a self-developed Sl\1EP(Smoke Movement Estimating Program) based on computational fluid d dynamics the$\alpha$1es. The results pn려icted by the two approaches are very similar. The field model u using SIl\1PLE algorithm or SIl\1PLER algorithm requires much more computing time compared w with the use of Sl\1EP using PISO algorithm.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.2 no.1
• /
• pp.63-73
• /
• 1990

A numerical analysis of initial unsteady state flow and heat transfer in an enclosed cavity has been performed by the Modified QUICK Scheme. The stable QUICK Scheme which modified the coefficient always to be positive is included in this numerical analysis. The implicit method is applied to solve the unsteady state flow; between iterations the PISO (Pressure - Implicit with Splitting of Operators) algorithm is employed to correct and update the velocity and pressure fields on a staggered grid. The accuracy of the Modified QUICK Scheme is proved by applying fewer grid systems than those which Ghia et al. and Davis applied. The initial unsteady mixed convection in an enclosed cavity is analyzed using the above numerical procedure. This study focuses on the development of the large main vortex and secondary vortex in forced convection, the effects of the Rayleigh Number in natural convection and the relative direction of the forced and natural convection.

• Fire Science and Engineering
• /
• v.11 no.4
• /
• pp.3-14
• /
• 1997

The smoke filling process for the atrium space containing a fire source is simulated using two types of deterministic fire model : Zone model and Field model. The zone model used is the CFAST(version 1.6) model developed at the Building and Fire Research Laboratories, NIST in the USA. The field model is a self-developed frie field model based on Computational Fluid Dynamic (CFD) theories. This article is focused on finding out the smoke movement and temperature distribution in atrium space which is cubic in shape. For solving the liked set of velocity and pressure equation, the PISO algorithm, which strengthened the velocity-pressure coupling, was used. Since PISO algorithm is a time-marching procedure, computing time si very fast. A computational procedure for predicting velocity and temperature distribution in fire-induced flow is based on the solution, in finite volume method and non-staggered grid system, of 3-dimensional equations for the conservation of mass, momentum, energy, species and so forth. The fire model i.e Zone model and Field model predicted similar results for clear heights and the smoke layer temperature.

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

This study is concerned with the critical evaluation of predicative capability of a k-$\varepsilon$ turbulence model using the Renormalization Group(RNG) theory. The present numerical model for solution of the Navier-Stokes System is based on the modified PISO algorithms. Computations have been performed with the RNG-based K-$\varepsilon$ model for the two-dimensional flow over a backward-facing step, a confined coaxial jet, and a swirling flow in a swirl combustor. Numerical results are compared with experimental data in terms of mean flow velocities, turbulent kinetic energy, and turbulent stresses. Numerical results clearly indicate that the RNG-based K-$\varepsilon$ turbulence model shows a significant improvement over a standard K-$\varepsilon$ model in predicting the turbulent flows with flow separation and swirl.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.112-122
• /
• 1996

A finite element scheme using the concept of PISO method has been developed to solve the Navier-Stokes viscous flows in all speed range. This scheme includes development of new pressure equation that retains both the hyperbolic term related with the density variation and the elliptic term reflecting the incompressibility constraint. The present method is applied to the incompressible two-dimensional driven cavity flow problems(Re=100, 400 and 1,000). For compressible flows, the Carter plate problem(M=3 and Re=1,000) is computed. Finally, we have simulated the shock-boundary layer interaction(M=2 and Re=2.96×10/sup 5/), a more difficult problem, and compared its results with the experiment to demonstrate the shock capturing capability of the present solution algorithm.

• Journal of computational fluids engineering
• /
• v.6 no.2
• /
• pp.22-31
• /
• 2001

The two-dimensional unsteady incompressible viscous flow behind rectangular bluff bodies between two parallel plates was analyzed numerically. The steady state flow and the vortex flow behind rectangular bluff bodies in the channel were investigated for two regimes i.e., the laminar(Re = 100, 300, 500) and the turbulent flows(Re = 10⁴∼10/sup 6/). The vortex shedding was generated by a physical disturbance(6%) numerically imposed at the rear of the bluff bodies for a short time. It was observed that the perturbed flow became periodic after a transient period. And in the case of unsteady inflow, the sinusoidal pulsatile flow was applied as the inlet condition in the turbulent flow of Reynolds number of 1.0×10/sup 5/. FLUENT code was employed to solve the problems. The power-law scheme was used to get stable linearized equations and the PISO algorithm was applied to finding the solution of them.

• Fire Science and Engineering
• /
• v.14 no.2
• /
• pp.21-32
• /
• 2000

In this study, reduced-scale experiments as the alternative to a real-scale fire test were conducted to understand fire properties in tunnel, and their results were compared with those of numerical simulation. The 1/20 scale experiments were conducted under the Froude scaling since smoke movement in tunnel is governed by buoyancy farce. A numerical simulations were on performed 3D unstructured meshes with PISO algorithm and buoyant plume models. Results showed that data was in reasonable agreement with the numerical data of smoke velocity, temperature distribution, and clear height.