• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.307-310
• /
• 2002

A two-dimensional numerical method for inviscid two-fluid flows with evolution of density interface is developed, and an initially stationary two-dimensional fluid sheet surrounded by another fluid is studied. The Interface between two fluids is modeled as a vertex sheet, and the flow field u÷th the evolution of interface is solved by using vortex-in-cell/front-tracking method. The edge of the sheet Is pulled back into the sheet due to surface tension and a blob is formed at the edge. This blob and fluid sheet are connected by a thin neck. In the inviscid limit, such process of the blob and neck formation is examined in detail and their kinematic characteristics are summarized with dimensionless parameters. The edge recedes at $V=1.06({\sigma}/{\rho}h)^{0.5}$ and the capillary wave Propagating into the fluid sheet must be considered for bettor understanding of the edge receding.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.9
• /
• pp.1227-1238
• /
• 2000

This study examines the effects of the discretization schemes on numerical solutions of viscoelastic fluid flows. For this purpose, a temporally evolving mixing layer, a two-dimensional vortex pair interacting with a wall, and a turbulent channel flow are selected as the test cases. We adopt a fourth-order compact scheme (COM4) for polymeric stress derivatives in the momentum equations. For convective derivatives in the constitutive equations, the first-order upwind difference scheme (UD) and artificial diffusion scheme (AD), which are commonly used in the literature, show most stable and smooth solutions even for highly extensional flows. However, the stress fields are smeared too much and the flow fields are quite different from those obtained by higher-order upwind difference schemes for the same flow parameters. Among higher-order upwind difference schemes, a third-order compact upwind difference scheme (CUD3) shows most stable and accurate solutions. Therefore, a combination of CUD3 for the convective derivatives in the constitutive equations and COM4 for the polymeric stress derivatives in the momentum equations is recommended to be used for numerical simulation of highly extensional flows.

• Nuclear Engineering and Technology
• /
• v.47 no.4
• /
• pp.416-423
• /
• 2015

It had been disputed how to apply wall drag to the dispersed phase in the framework of the conventional two-fluid model for two-phase flows. Recently, Kim et al. [1] introduced the volume-averaged momentum equation based on the equation of a solid/fluid particle motion. They showed theoretically that for dispersed two-phase flows, the overall two-phase pressure drop by wall friction must be apportioned to each phase, in proportion to each phase fraction. In this study, the validity of the proposed wall drag model is demonstrated though one-dimensional (1D) simulations. In addition, it is shown that the existing form loss model incorrectly predicts the motion of the dispersed phase. A new form loss model is proposed to overcome that problem. The newly proposed form loss model is tested in the region covering the lower plenum and the core in a nuclear power plant. As a result, it is shown that the new models can correctly predict the relative velocity of the dispersed phase to the surrounding fluid velocity in the core with spacer grids.

• Journal of computational fluids engineering
• /
• v.20 no.3
• /
• pp.1-7
• /
• 2015

When the steam is used as working fluid in fluid machinery, different from other gases as air, phase transition (steam condensation) can occur and it affects not only the flow fields, but also machine performance & efficiency. Therefore, considering phase transition phenomena in CFD calculation is required to achieve accurate prediction of steam flow and non-equilibrium wet-steam model is needed to simulate realistic steam condensing flow. In this research, non-equilibrium wet-steam model is implemented on in-house code(T-Flow), the flow fields including phase transition phenomena in convergent-divergent nozzle are studied and compared to results of advance researches.

• Journal of Korean Society for Atmospheric Environment
• /
• v.12 no.3
• /
• pp.231-241
• /
• 1996

The effect of radiation and buoyancy on the thermophoresis phenomenon owing to the presence of highly absorbing, emitting particles (such as soot or pulverized coal) suspended in a two phase flow system was investigated numerically for a turbulent mixed convection flow. The analysis of conservation equations for a gas-particle flow system was performed on the basis of a two-fluid model from a continuum Eulerian viewpoint. The modified van Driest and Cebeci mixing length turbulence model was adopted in the anaylsis of turbulent flow. In addition, the P-1 approximation was used to evaluate the radiation heat transfer. As expected from the particle concentration and drift velocity distribution, the cumulative collection efficiency E (x) becomes larger when the buoyancy effect increases (i.e. higher Grashof number), while smaller as the radiation effect increases (i.e. higher optical thickness).

• The KSFM Journal of Fluid Machinery
• /
• v.14 no.3
• /
• pp.10-17
• /
• 2011

Free energy based lattice Boltzmann method (LBM) has been used to simulate the contact angle and the bubble necking with large density ratio. LBM with the proper contact angle model is able to reduce the spurious currents and eliminate the singularity in the contact lines. The numerical results of the contact angles are satisfied with the Youngs law. For bubble necking flows, simulations are executed for various viscosities and contact angles. The phenomena of the bubble necking are simulated successfully and the subsequent results are presented. The present method is also applicable to the nucleate boiling flows.

• The KSFM Journal of Fluid Machinery
• /
• v.11 no.3
• /
• pp.7-13
• /
• 2008

The cavitating flow simulation is of practical importance for many engineering systems, such as pump, turbine, nozzle, Infector, etc. In the present work, a solver for two-phase flows has been developed and applied to simulate the cavitating flows past hydrofoils. The governing equation is the two-phase Navier-Stokes equation, comprised of the continuity equation of liquid and vapor phase. The momentum and energy equation is in the mixture phase. The solver employs an implicit, dual time, preconditioned algorithm using finite difference scheme in curvilinear coordinates. An experimental data and other numerical data were compared with the present results to validate the present solver. It is concluded that the present numerical code has successfully accounted for two-phase Navier-Stokes model of cavitation flow.

• International Journal of Fluid Machinery and Systems
• /
• v.1 no.1
• /
• pp.1-9
• /
• 2008

In the first part of the paper, Computational Fluid Dynamics analysis of the combusting flow within a high-swirl lean premixed gas turbine combustor and over the $1^{st}$ row nozzle guide vanes is presented. In this analysis, the focus of the investigation is the fluid dynamics at the combustor/turbine interface and its impact on the turbine. The predictions show the existence of a highly-rotating vortex core in the combustor, which is in strong interaction with the turbine nozzle guide vanes. This has been observed to be in agreement with the temperature indicated by thermal paint observations. The results suggest that swirling flow vortex core transition phenomena play a very important role in gas turbine combustors with modern lean-premixed dry low emissions technology. As the predictability of vortex core transition phenomena has not yet been investigated sufficiently, a fundamental validation study has been initiated, with the aim of validating the predictive capability of currently-available modelling procedures for turbulent swirling flows near the sub/supercritical vortex core transition. In the second part of the paper, results are presented which analyse such transitional turbulent swirling flows in two different laboratory water test rigs. It has been observed that turbulent swirling flows of interest are dominated by low-frequency transient motion of coherent structures, which cannot be adequately simulated within the framework of steady-state RANS turbulence modelling approaches. It has been found that useful results can be obtained only by modelling strategies which resolve the three-dimensional, transient motion of coherent structures, and do not assume a scalar turbulent viscosity at all scales. These models include RSM based URANS procedures as well as LES and DES approaches.

• The KSFM Journal of Fluid Machinery
• /
• v.13 no.4
• /
• pp.31-36
• /
• 2010

Free energy based lattice Boltzmann method (LBM) has been used to simulate the rising bubble flows with large density ratio. LBM with compact discretization is able to reduce the spurious current of the static bubble test and be satisfied with the Laplace law. The terminal rise velocity and shape of the bubbles are dependent on Eotvos number, Morton number and Reynolds number. For single bubble flows, simulations are executed for various Eotvos number, Morton number and Reynolds number, and the results are agreed well with the experiments. For multiple bubbles, the bubble flow characteristics are related by the vortex pattern of the leading bubble. The coalescence of the bubbles are simulated successfully and the subsequent results are presented. The present method is validated for static, dynamic bubble test cases and compared to the numerical, experimental results.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.1-9
• /
• 2001

In this paper, some numerical methods recently developed for gas-liquid two-phase flows are reviewed. And then, a preconditioning method to solve cavitating flow by the author is introduced. This method employs a finite-difference Runge-Kutta method combined with MUSCL TVD scheme, and a homogeneous equilibrium cavitation model. So that it permits to treat simply the whole gas-liquid two-phase flow field including wave propagation, large density changes and incompressible flow characteristic at low Mach number. Finally, numerical results such as detailed observations of the unsteady cavity flows, a sheet cavitation break-off phenomena and some data related to performance characteristics of hydrofoils are shown.