• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.4
• /
• pp.359-370
• /
• 2015

The present study aims to numerically investigate the behavior of liquid droplet driven by capillarity force imbalance on horizontal surfaces ranging from hydrophilic to hydrophobic, under various conditions. The droplet behavior has been simulated using an in-house solution code(PowerCFD), which employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with interface capturing method(CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The detailed droplet behavior was obtained under various conditions for droplets with different initial shapes, contact angles and surface tension forces(or Bond number). The mechanism of droplet transport was examined using the numerical results on the droplet shapes.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.474-481
• /
• 2001

This paper addresses a numerical method for predicting transient temperature distributions in the wall of a curved pipe subjected to internally thermal stratification flow. A simple and convenient numerical method of treating the unsteady conjugate heat transfer in the non-orthogonal coordinate systems is presented. The proposed method is implemented in a finite volume thermal-hydraulic computer code based on a cell-centered, non-staggered grid arrangement, the SIMPLEC algorithm, a higher-order bounded convection scheme, and the modified version of momentum interpolation method. Calculations are performed for the transient evolution of thermal stratification in two curved pipes, where the one has thick wall and the other has so thin wall that its presence can be negligible in the heat transfer analysis. The predicted results show that the thermally stratified flow and transient conjugate heat transfer in a curved pipe with a finite wall thickness can be satisfactorily analyzed by the present numerical method, and that the neglect of wall thickness in the prediction of pipe wall temperature distributions can provide unacceptably distorted results.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.7-13
• /
• 2011

A comparative analysis of thermal models in the lattice Boltzmann method(LBM) for the simulation of laminar natural convection in a square cavity is presented. A HYBRID method, in which the thermal equation is solved by the Navier-Stokes equation method while the mass and momentum conservation are resolved by the lattice Boltzmann method, is introduced and its merits are explained. 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 a deferred correction method to ensure stability of the solutions. The HYBRID method and the double-population method are applied to the simulation of natural convection in a square cavity and the predicted results are compared with the benchmark solutions given in the literatures. The predicted results are also compared with those by the conventional Navier-Stokes equation method. In general, the present HYBRID method is as accurate as the Navier-Stokes equation method and the double-population method. The HYBRID method shows better convergence and stability than the double-population method. These observations indicate that this HYBRID method is an efficient and economic method for the simulation of incompressible fluid flow and heat transfer problem with the LBM.

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

The non-staggered(collocated) grid approach in which all the solution variables are located at the centers of control volumes is very popular for incompressible flow analyses because of its numerical efficiency on the curvilinear or unstructured grids. Rhie and Chow's paper is the first in using non-staggered grid method for SIMPLE algorithm, where pressure weighted interpolation was used to prevent decoupling of pressure and velocity. But it has been known that this non-staggered grid method has stability problems when pressure fields are nonlinear like in natural convection flows. Also Rhie-Chow scheme generates large numerical diffusion near curved walls. The cause of these unwanted problems is too large pressure damping term compared to the magnitude of face velocity. In this study the magnitude of pressure damping term of Rhie-Chow's method is limited to 1∼10% of face velocity to prevent physically unreasonable solutions. The wall pressure extrapolation which is necessary for cell-centered FVM is another source of numerical errors. Some methods are applied in a unstructured FV solver and analyzed in view of numerical accuracy. Here, two natural convection problems are solved to check the effect of the Rhie-Chow's method on numerical stability. And numerical diffusion from Rhie-Chow's method is studied by solving the inviscid flow around a circular cylinder.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.38 no.3
• /
• pp.263-270
• /
• 2014

A fluid transport technique is a key issue for the development of microfluidic systems. In this paper, a new concept for transporting a droplet without external power sources is proposed and verified numerically. The proposed device is a heterogeneous surface which has both hydrophilic and hydrophobic horizontal surfaces. The numerical simulation to demonstrate the new concept is conducted by an in-house solution code (PowerCFD) which employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with interface capturing method (CICSAM) in a volume of fluid (VOF) scheme for phase interface capturing. It is found that the proposed concept for droplet transport shows superior performance for droplet transport in microfluidic systems.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.4
• /
• pp.19-30
• /
• 2000

A finite-volume method is developed to solve turbulent flows around modern commercial hull forms with bow and stern bulbs. The RANS equations are solved. The cell-centered finite-volume method employs QUICK and central difference scheme for convective and diffusive flux discretization, respectively. The SIMPLEC method is adopted for the velocity-pressure coupling. The developed numerical methods are applied to calculate turbulent flow around KRISO 3600TEU container ship. Surface meshes are generated into five blocks: bow and stern bulbs, overhang, fore and afterbody. 3-D field grid system with O-H topology is generated using elliptic grid generation method. Surface friction lines and wake distribution at propeller plane is compared with experiment. The calculated results show that the present method can be used to predict flow around a modern commercial hull forms with bulbs.

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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.7
• /
• pp.723-733
• /
• 2011

A numerical method of the CSF(Continuum Surface Force) model is presented for the calculation of the surface tension force and implemented in an in-house solution code(PowerCFD). The present method(code) employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with volume capturing method(CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The application of the present method to a 2-D liquid drop problem is illustrated by an equilibrium and nonequilibrium oscillating drop calculation. It is found that the present method simulates efficiently and accurately surface tension-dominant multiphase flows.

• KIEAE Journal
• /
• v.16 no.1
• /
• pp.89-94
• /
• 2016

This study is part of three-dimensional(3D) heat transfer analysis program developmental process. The program is being developed without it's own built in 3D-modeller. So 3D-model must be created from another 3D-modeller such as generic CAD programs and imported to the developed program. After that, according to the 3D-geometric data form imported model, 3D-mesh created for numerical calculation. But the 3D-model created from another 3D-modeller is likely to have errors in it's geometric data such as mismatch of position between vertexes or surfaces. these errors make it difficult to create 3D-mesh for calculation. These errors are must be detected and cured in the pre-process before creating 3D-mesh. So, in this study four kinds of filters and functions are developed and tested. Firstly, 'vertex error filter' is developed for detecting and curing for position data errors between vertexes. Secondly, 'normal vector error filter' is developed for errors of surface's normal vector in 3D-model. Thirdly, 'intersection filter' is developed for extracting and creating intersection surface between adjacent objects. fourthly, 'polygon-line filter' is developed for indicating outlines of object in 3D-model. the developed filters and functions were tested on several shapes of 3D-models. and confirmed applicability. these developed filters and functions will be applied to the developed program and tested and modified continuously for less errors and more accuracy.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.7
• /
• pp.721-730
• /
• 2012

The present paper presents a numerical study on multiphase flows induced by wall adhesion. The continuum surface force (CSF) model with the wall adhesion boundary condition model is used for calculating the surface tension force; this model is implemented in an in-house solution code (PowerCFD). The present method (code) employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with a volume capturing method (CICSAM) in a volume of fluid (VOF) scheme for phase interface capturing. The effects of wall adhesion are then numerically simulated by using the present method for a shallow pool of water located at the bottom of a cylindrical tank with no external forces such as gravity. Two different cases are computed, one in which the water wets the wall and one in which the water does not wet the wall. It is found that the present method efficiently simulates the surface tension-dominant multiphase flows induced by wall adhesion.