• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.21-25
• /
• 2007

A general purpose program NUFLEX has been extended for two-phase flows with topologically complex interface and cavitation flows with liquid-vapor phase change caused by large pressure drop. In analysis of two-phase flow, the phase interfaces are tracked by employing a LS(Level Set) method. Compared with the VOF(Volume-of-Fluid) method based on a non-smooth volume-fraction function, the LS method can calculate an interfacial curvature more accurately by using a smooth distance function. Also, it is quite straightforward to implement for 3-D irregular meshes compared with the VOF method requiring much more complicated geometric calculations. Also, the cavitation process is computed by including the effects of evaporation and condensation for bubble formation and collapse as well as turbulence in flows. The volume-faction and continuity equations are adapted for cavitation models with phase change. The LS and cavitation formulation are implemented into a general purpose program for 3-D flows and verified through several test problems.

• Journal of computational fluids engineering
• /
• v.11 no.3 s.34
• /
• pp.14-21
• /
• 2006

A numerical simulation of the binary collision dynamics of water drops for size ratios of 1 and 0.75, for the Weber number range of 5 to 100, and for all impact parameter is reported. Two different types of separating collisions, namely reflexive and stretching separations, are identified. A numerical method is based on a fractional-step method with a finite volume formulation and the interface is tracked with Volume of Fluid(VOF) method, including surface tension. Numerical results for size ratios 1 and 0.75 are reasonablely compared with Ashgriz and Poo's experimental results.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.86-95
• /
• 2008

For analyses of multi-phase flows in a water-cooled nuclear power plant, a three-dimensional SIMPLE-algorithm based hydrodynamic solver CUPID-S has been developed. As governing equations, it adopts a two-fluid three-field model for the two-phase flows. The three fields represent a continuous liquid, a dispersed droplets, and a vapour field. The governing equations are discretized by a finite volume method on an unstructured grid to handle the geometrical complexity of the nuclear reactors. The phasic momentum equations are coupled and solved with a sparse block Gauss-Seidel matrix solver to increase a numerical stability. The pressure correction equation derived by summing the phasic volume fraction equations is applied on the unstructured mesh in the context of a cell-centered co-located scheme. This paper presents the numerical method and the preliminary results of the calculations.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.39-46
• /
• 2011

A numerical study of a turbulent natural convection in an enclosure with the lattice Boltzmann method (LBM) is presented. The primary emphasis of the present study is placed on investigation of accuracy and numerical stability of the LBM for the turbulent natural convection flow. A HYBRID method in which the thermal equation is solved by the conventional Reynolds averaged Navier-Stokes equation method while the conservation of mass and momentum equations are resolved by the LBM is employed in the present study. The elliptic-relaxation model is employed for the turbulence model and the turbulent heat fluxes are treated by the algebraic flux model. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with the deferred correction way to ensure accuracy and stability of solutions. The present LBM is applied to the prediction of a turbulent natural convection in a rectangular cavity and the computed results are compared with the experimental data commonly used for the validation of turbulence models and those by the conventional finite-volume method. It is shown that the LBM with the present HYBRID thermal model predicts the mean velocity components and turbulent quantities which are as good as those by the conventional finite-volume method. It is also found that the accuracy and stability of the solution is significantly affected by the treatment of the convection term, especially near the wall.

• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.67-72
• /
• 2005

This research is based on the iterative space-marching method for incompressible and compressible Navier-Stokes equations[1-4]. A principle of parallel computational schemes construction for steady and unsteady problems is suggested. It is analytically proven that convergence of these schemes is unconditional for incompressible case. When the parallel scheme is used the total volume of computations is the sum of a large number of independent and equal parts. Estimation of the speed-up K shows that K > 1000 in ideal case. First results of using the parallel schemes are presented.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.209-213
• /
• 2007

This paper presents the high order conservative interpolation methods. This method has been successfully implemented into a cell-centered finite volume methods code. Results using high order and maintaining conservation interpolation method are more accurate compared to bilinear interpolation.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.910-915
• /
• 2003

A Navier-Stokes equation solver for incompressible viscous flows with free surface is developed and tested. This is based upon a fractional time step method and a non-staggered finite volume formulation for unstructured meshes. For time advancement scheme, Adams -Bashforth method for convective term and Crank-Nicolson method for diffusive term are applied. The interface between two fluids with different fluid properties is tracked with Piecewise Linear Interface Calculation(PLIC) Volume-of-Fluid(VOF) methods. Computational results are presented for some test problems: the broken dam, the sloshing in a rectangular tank, the filling of a cylindrical tank.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.57-60
• /
• 2006

This paper presents an assessment of vertex-centered and cell-centered finite-volume methods on unstructured meshes. The results indicate that the vertex-centered method is more reliable than the cell-centered method.

• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.36-42
• /
• 2006

The existing numerical approximations of convection flux, especially the spatial higher-order difference schemes, in unstructured cell-centered finite volume methods are examined in detail with each other and evaluated with respect to the accuracy through their application to a 2-D benchmark problem. Six higher-order schemes are examined, which include two second-order upwind schemes, two central difference schemes and two hybrid schemes. It is found that the 2nd-order upwind scheme by Mathur and Murthy(1997) and the central difference scheme by Demirdzic and Muzaferija(1995) have more accurate prediction performance than the other higher-order schemes used in unstructured cell-centered finite volume methods.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.31-36
• /
• 2003

A field reconstruction scheme for a cell centered finite volume method on unstructured meshes is developed. Regardless of mesh quality, this method is exact within a machine accuracy if the solution is linear, which means it has full second order accuracy. It does not have any limitation on cell shape except convexity of the cells and recovers standard discretization stencils at structured orthogonal grids. Accuracy comparisons with other popular reconstruction schemes are performed on a simple example.