• 한국전산유체공학회지
• /
• 제8권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.

• 한국전산유체공학회지
• /
• 제1권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,.

• 한국전산유체공학회지
• /
• 제7권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.

• 한국전산유체공학회지
• /
• 제1권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.

• 설비공학논문집
• /
• 제2권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.

• 한국화재소방학회논문지
• /
• 제11권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.

• 한국해양공학회지
• /
• 제22권1호
• /
• pp.1-5
• /
• 2008

A 2-D Fluent-based numerical wave tank(FBNWT) capable of simulating wave propagating and overtopping is presented. The FBNWT model is based on the Reynolds averaged Naiver-Stokes equations and VOF free surface tracking method. The piston wave maker system is realized by dynamic mesh technology(DMT) and user defined function(UDF). The non-iteration time advancement(NITA) PISO algorithm is employed for the velocity and pressure coupling. The FBNWT numerical solutions of linear wave propagation have been validated by analytical solutions. Several overtopping problems are simulated and the prediction results show good agreements with the experimental data, which demonstrates that the present model can be utilized in the corresponding analysis.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.58-59
• /
• 2003

• 한국전산유체공학회지
• /
• 제2권1호
• /
• pp.1-7
• /
• 1997

A FA-FD hybrid method, developed for solving three-dimensional incompressible Navier-Stokes equations, is applied to calculate three-dimensional laminar flows through a square duct with a 90° bend. The method discretizes the convective terms in the primary flow direction with 3rd-order upwind finite-differences and the convective and diffusive terms in the transverse directions with the two-dimensional finite analytic method. The non-staggered grid system is used and the pressure-velocity coupling is achieved by a global iteration procedure based on the PISO algorithm. Detailed comparisons between the computed solutions and the available experimental data are given mainly for the velocity distributions at cross-sections in a 90° bend of a square duct with both fully developed and developing entry flows. Although the computational result shows generally a good agreement with the experimental data, there are some significant discrepancies underlining the necessity of more accurate numerical methods as well as reliable experimental data for their validation.

• International Journal of Air-Conditioning and Refrigeration
• /
• 제10권3호
• /
• pp.162-175
• /
• 2002

To investigate smoke movement with radiation in a room fires, a numerical and experi-mental analysis were performed. In this paper, results from a field model based on a self-developed SMEP (Smoke Movement Estimating Program) were compared with Stockier's ex-periment and the experiments on various sized pool fires in a room with door 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 S-N discrete ordinates method (DOM). The result of the cal-culated smoke temperature considering radiation effect has shown good agreement compared with the experimental data, although there are large discrepancy in the hot smoke layer be-tween the temperature predicted by the SMEP with only convection effect and obtained by the experimental result. This large discrepancy is caused from the radiation effect of $H_2O$ and $CO_2$ gas under smoke productions. Hence the radiation effect under smoke in fire is the point to be specially considered in order to produce more realistic result.