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

A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed. A two-fluid three-field model was adopted for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. The hydrodynamic solver is for the 3D component of a nuclear system code and the component-scale analysis tools for transient two-phase flows. The finite volume method and unstructured grid are adopted, which are useful for the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been adapted to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing problems well.

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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.12
• /
• pp.1555-1569
• /
• 2000

Numerical simulation of fluid flow with moving free surface has been carried out. For the free surface flow, a VOF(Volume of Fluid)-based algorithm utilizing a fixed grid system has been investigated. In order to reduce numerical smearing at the free surface represented on a fixed grid system, a new free surface tracking algorithm based on the donor-acceptor scheme has been presented. Novel features of the proposed algorithm are characterized as two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed algorithm can be easily implemented in any irregular non-uniform grid systems that are usual in finite element method (FEM). Moreover, the proposed algorithm can be extended and applied to the 3-D free surface flow problem without additional efforts. For computation of unsteady incompressible flow, a finite element approximation based on the explicit fractional step method has been adopted. In addition, the SUPG(streamline upwind/Petrov-Galerkin) method has been implemented to deal with convection dominated flows. Combination of the proposed free surface tracking scheme and explicit fractional step formulation resulted in an efficient solution algorithm. Validity of the present solution algorithm was demonstrated from its application to the broken dam and the solitary wave propagation problems.

• Proceedings of the KSR Conference
• /
• 1998.11a
• /
• pp.137-143
• /
• 1998

Numerical simulation is conducted to clarify the heat transfer and fluid flow characteristics of traction motor for electric car SIMPLE algorithm based on finite volume method is used to make linear algebra equation. The governing equations are solved by TDMA(TriDiagonal Matrix Algorithm) with line-by-line method and block correction. From the results of simulation, the characteristics of cooling pattern is strongly affected by the size of hole in stator core. In the case of high rotational speed of rotor, temperature difference along the axial direction is more decreased than that of low rotational speed.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.48-57
• /
• 1999

This paper deals with the numerical simulation of the behavior of single bubble rising near the free surface. Volume fraction of fluid (VOF) method with continuum surface force (CSF) model, the well known method for two phase flow simulation is adopted. A bubble of spherical shape positioned beneath the free surface is assumed at the initial stage. The difference according to the fluid properties of surrounding medium is examined. Simulation results are depicted and explained with the time history of bubble shape, velocity field and vorticity distribution.

• The KSFM Journal of Fluid Machinery
• /
• v.1 no.1 s.1
• /
• pp.41-48
• /
• 1998

This study was carried out to investigate the effect of design parameters on the volume flow-rate and the noise level and to finally find the optimal design variables. Eighteen cross-flow fans were designed by the method of orthogonal array, and the flow-rate and the noise level were measured. These data were analyzed by the neural network system. The effects of eight design variables(scroll exit angle, scroll arc length et al.) on the fan performance and the noise level were valuated and discussed. This experiment shows that the design solutions suggested by neural network system may increase its volume flow-rate and reduce noise simultaneously.

• Korea-Australia Rheology Journal
• /
• v.18 no.2
• /
• pp.41-49
• /
• 2006

A modified Volume-of-Fluid (VOF) numerical method is used to predict the dynamics of a liquid drop of a low viscosity dilute polymer solution, forming in air from a circular nozzle. Viscoelastic effects are rep-resented using an Oldroyd-B model. Predicted drop shapes are compared with experimental observations. The main features, including the timing of the shape evolution and the 'bead-on-a-string' effect, are well reproduced by the simulations. The results confirm published conclusions of the third author, that the deformation is effectively Newtonian until near the time of Newtonian pinch-off and that the elastic stress becomes large in the pinch region due to the higher extensional flow there.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.4
• /
• pp.459-479
• /
• 2015

Two-dimensional steady laminar natural convection around a heat conducting and generating solid body inside a square enclosure with different thermal boundaries is performed. The mathematical model is governed by the coupled equation of mass, momentum and energy. These equations are discretized by finite volume method with power-law scheme and solved numerically by SIMPLE algorithm with under-relaxation technique. Effect of Rayleigh number, temperature difference ratio of solid-fluid, aspect ratio of solid-enclosure and the thermal conductivity ratio of solid-fluid are investigated numerically for Pr = 0.7. The flow and heat transfer aspects are demonstrated in the form of streamlines and isotherms respectively.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.12
• /
• pp.1017-1023
• /
• 2007

This study is intended to analyze the effect of thin ink-film thickness around rotating printing roll on the printing quality in the gravure printing process which is used for making electronics circuit like a RFID tag with a conductive ink. The present work numerically estimates the film thickness around rotating roller partially immersed in ink for which the volume of fluid (VOF) method was adopted to figure out the film formation process around rotating roller. Parameter studies were performed to compare the effect of ink properties (viscosity, surface tension), operating condition (roller rotating speed, initial immersed angle) on the film thickness. The result indicates that the film thickness has a strong dependency on the rotating speed, while the surface tension has negligible effect.

• Journal of the Korean Society of Combustion
• /
• v.6 no.2
• /
• pp.15-22
• /
• 2001

A pressure-based, unstructured finite volume method has been applied to couple the chemical kinetics and fluid dynamics and to capture effectively and accurately the steep gradient flame field. The pressure-velocity coupling is handled by two methodologies including the pressure-correction algorithm and the projection scheme. A stiff, operator-split projection scheme for the detailed nonequilibrium chemistry has been employed to treat the stiff reaction source terms. The conservative form of the governing equations are integrated over a cell-centered control volume with collocated storage for all transport variables. Computations using detailed chemistry and variable transport properties were performed for two laminar reacting flows: a counterflow hydrogen-air diffusion flame and a lifted methane-air triple flame. Numerical results favorably agree with measurements in terms of the detailed flame structure.