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

• Journal of the Korean Society of Safety
• /
• v.15 no.3
• /
• pp.14-22
• /
• 2000

A numerical analysis has been performed on the two-dimensional rectangular gallium melting problem using the enthalpy method. The major advantage of this method is that the physical domain is discretized with fixed grids without transforming variables and the interface conditions of phase change are accounted for the definition of suitable source terms in the governing equations. But in the fixed method, there is some ambiguity in defining the porosity constant which has no physical interpretation. If the velocity correction is included in the momentum equation, for the appropriate range of porosity constant, the realistic predictions are obtained. The object of the present work is to predict the phase change within the molten steel with thin riser slab using the modified enthalpy-porosity method. The computational procedures for predicting velocity and temperature are based on the finite volume method and the non-staggered grid system. The influence of natural convection on the melting process is considered. A comparison with the experimental results shows that the modified method is better than the previous one.

• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.25-31
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP method. The result using improved CIP method shows the better result of the experiments, comparison with result of original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• Wind and Structures
• /
• v.26 no.5
• /
• pp.331-341
• /
• 2018

A numerical study on the flow over a square cylinder in the vicinity of a wall is conducted for different Couette-Poiseuille-based non-uniform flow with the non-dimensional pressure gradient P varying from 0 to 5. The non-dimensional gap ratio L (=$H^{\ast}/a^{\ast}$) is changed from 0.1 to 2, where $H^{\ast}$ is gap height between the cylinder and wall, and $a^{\ast}$ is the cylinder width. The governing equations are solved numerically through finite volume method based on SIMPLE algorithm on a staggered grid system. Both P and L have a substantial influence on the flow structure, time-mean drag coefficient ${\bar{C}}_D$, fluctuating (rms) lift coefficient ($C_L{^{\prime}}$), and Strouhal number St. The changes in P and L leads to four distinct flow regimes (I, II, III and IV). Following the flow structure change, the ${\bar{C}}_D$, $C_L{^{\prime}}$, and St all vary greatly with the change in L and/or P. The ${\bar{C}}_D$ and $C_L{^{\prime}}$ both grow with increasing P and/or L. The St increases with P for a given L, being less sensitive to L for a smaller P (< 2) and more sensitive to L for a larger P (> 2). A strong relationship is observed between the flow regimes and the values of ${\bar{C}}_D$, $C_L{^{\prime}}$ and St. An increase in P affects the pressure distribution more on the top surface than on bottom surface while an increase in L does the opposite.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.126-134
• /
• 2008

We developed a three dimensional Navier-Stokes code based on the level set method to simulate two phase flows with high density ratio. The Navier-Stokes equations with consideration of the surface tension effects are solved by using SIMPLE algorithm on a non-staggered grid. The present code is validated by simulating two test problems. First one is to simulate a rising bubble inside a cube. The thickness of the interface of the bubble is shown to affect the pressure distribution around the interface. As the thickness decreases, the pressure field around the interface becomes more oscillatory. As the bubble rises, a ring vortex is shown to form around the interface and the bubble eventually develops into an ellipsoidal shape. Merge of two bubbles inside a container is secondly tested to show the robustness of the present code for two phase flow simulation. Numerical results show stable and reliable behavior during the process of merging of two bubbles. The velocity and pressure fields around the interface of bubbles are shown oscillation free during the merging of two bubbles.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.8
• /
• pp.2611-2625
• /
• 1996

A numerical method on multi-block technique by Bi-CGSTAB(Bi-Conjugate Gradient STABilized) solver has been proposed. The present multi-block technique can reduce the numerical manipulation greatly because the common regions at the interface of each block are not necessary. In order to test the computational performance of present multi-block technique, the flow characteristics in a T type duct system and a N type duct system have been investigated by three kinds of methods such as the single-block method, the previous multi-block technique and the multi-block technique with Bi-CGSTAB solver. The results indicated that the required CPU time by present multi block technique was shorter than that of other two numerical methods and the convergency history was shown very stable at the present multi-block technique.

• Proceedings of the Korean Institute of Industrial Safety Conference
• /
• 1997.11a
• /
• pp.43-48
• /
• 1997

The smoke filling process for the atrium space containing a fire source is simulated using two types of deterministic fire modus: Zone model and Field model. The zone mode used is the CFAST(version 1.6) mode developed at the Building and Fire Research laboratories, NIST in the USA. The lied model is a self-developed fire field model based on Computational Fluid Dynamics(CFD) theories. This article is focused on finding out the smoke movement and temperature distribution in atrium space which is cubic in shape. A computational procedure for predicting velocity and temperature distribution in fro-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 Ire clear height and the smoke layer temperature.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2002.04a
• /
• pp.35-40
• /
• 2002

A numerical analysis based on two-dimensional and three-dimensional incompressible Navier-Stokes equations has been carried out for double-circular-arc compressor cascades and the results are compared with available experimental data at various incidence angles. The 2-D and 3-D computational codes based on SIMPLE algorithm adopt pressure weighted interpolation method for non-staggered grid and hybrid scheme for the convertive terms. Turbulence modeling is very important for prediction of cascade flows, which are extremely complex with separation and reattachment by adverse pressure gradient. In this paper k-$\varepsilon$ turbulence model with wall function is used to increase efficiency of computation times.

• Journal of computational fluids engineering
• /
• v.2 no.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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.8
• /
• pp.1046-1055
• /
• 2000

The flow and thermal phenomena in flat tubes of radiator are analyzed numerically. To predict the characteristics of heat transfer and pressure drop, the flow analysis program for three-dimensional complex geometry is developed, which adopted an non-staggered grid system and Cartesian velocities as dependent variables of the momentum equations. Using the developed program, the effect of tube specifications on the heat transfer characteristics is investigated for various flat tubes. From this study, the following results are obtained; (1) For the same hydraulic diameter($D_h{\doteq}5.2$mm), the Nusselt numbers of three basic modeis(D, J, and H-model) are 8.71, 8.92, and 10.58, respectively, and the pressure drops of D-, J-, and H-model are predicted as $-3.08{\times}10^{-2}\;Pa,\;-3.12{\times}10^{-2}\;Pa,\;and\; -3.98{\times}10^{-2}$ Pa, (2) In case of the same flat tube specification, the fins must be brazed at upper tube surface because the heat is more vividly transferred. Therefore, it is found that the H- model is the most effective tube as a heat exchanger and these results are used as a fundamental data for the design of tube.