• 한국전산유체공학회:학술대회논문집
• /
• 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.

• Journal of the Korean Society of Visualization
• /
• v.16 no.1
• /
• pp.37-41
• /
• 2018

The two-fluid model is widely used for practical applications involving multi-phase flows in chemical reactor, nuclear reactor, desalination systems, boilers, and internal combustion engine. There are several modeling terms in the two-fluid model, which must be determined properly. This study suggests the best models for turbulent vertical bubbly flow.

• Journal of Korean Port Research
• /
• v.13 no.1
• /
• pp.167-174
• /
• 1999

Numerical analysis of two-fluid flows for both water and air is carried out. Free-Surface flows with an arbitrary deformation have been simulated around two dimensional submerged hydrofoil. The computation is performed using a finite volume method with unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell-wise local mesh refinement. the integration in space is of second order based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels The linear equation systems are solved by conjugate gradient type solvers and the non-linearity of equations is accounted for through picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations the continuity equation the conservation equation of one species and the equations or two turbulence quantities.

• Journal of Environmental Science International
• /
• v.23 no.7
• /
• pp.1223-1232
• /
• 2014

Density currents have been easily observed in environmental flows, for instance turbidity currents and pollutant plumes in the oceans and rivers. In this study, we explored the propagation dynamics of density currents using the FLOW-3D computational fluid dynamics code. The renormalization group (RNG) $k-{\varepsilon}$ scheme, a turbulence numerical technique, is employed in a Reynold-averaged Navier-Stokes framework (RANS). The numerical simulations focused on two different types of intrusive density flows: (1) propagating into a two-layer ambient fluid; (2) propagating into a linearly stratified fluid. In the study of intrusive density flows into a two-layer ambient fluid, intrusive speeds were compared with laboratory experiments and analytical solutions. The numerical model shows good quantitative agreement for predicting propagation speed of the density currents. We also numerically reproduced the effect of the ratio of current depth to the overall depth of fluid. The numerical model provided excellent agreement with the analytical values. It was also clearly demonstrated that RNG $k-{\varepsilon}$ scheme within RANS framework is able to accurately simulate the dynamics of density currents. Simulations intruding into a continuously stratified fluid with the various buoyancy frequencies are carried out. These simulations demonstrate that three different propagation patterns can be developed according to the value of $h_n/H$ : (1) underflows developed with $h_n/H=0$ ; (2) overflows developed when $h_n/H=1$ ; (3) intrusive interflow occurred with the condition of 0 < $h_n/H$ < 1.

• Journal of computational fluids engineering
• /
• v.12 no.4
• /
• pp.44-53
• /
• 2007

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

• Korea-Australia Rheology Journal
• /
• v.21 no.3
• /
• pp.161-173
• /
• 2009

This study analyzes especially drag and lift models recently developed for fluid-solid, fluid-fluid or liquid-liquid two-phase flows to understand their applicability on the computational fluid dynamics, CFD modeling of pulsatile blood flow. Virtual mass effect and the effect of red blood cells, RBCs aggregation on CFD modeling of blood flow are also shortly reviewed to recognize future tendencies in this field. Recent studies on two-phase flows are found as very useful to develop more powerful drag-lift models that reflect the effects of blood cell's shape, deformation, concentration, and aggregation.

• 한국가시화정보학회:학술대회논문집
• /
• 2004.11a
• /
• pp.90-93
• /
• 2004

Micro PIV was applied to measure velocity profiles of two-fluid flows in a microchannel. In this work, the two-fluid flow of two glycerol-water mixtures was measured for three cases $(\phi=0\;and\;\phi=0.2;\;\phi=0.1\;and\;\phi=0.5;\;\phi=0\;and\;\phi=0.6)$. The flow rates of two fluids were the same. The experimental velocity profiles agreed well with numerical simulations. However, a slight deviation was found in the fluid with low concentration. Rather than the effects of the varying refractive indices inside the channel, the high velocity gradient effect was thought as the main source of the deviation.

• Nuclear Engineering and Technology
• /
• v.54 no.11
• /
• pp.4322-4328
• /
• 2022

The two-fluid model is widely used to describe two-phase flows in complex systems such as nuclear reactors. Although the two-phase flow was successfully simulated, the standard two-fluid model suffers from an ill-posed nature. There are several remedies for the ill-posedness of the one-dimensional (1D) two-fluid model; among those, artificial viscosity is the focus of this study. Some previous works added artificial diffusion terms to both mass and momentum equations to render the two-fluid model well-posed and demonstrated that this method provided a numerically converging model. However, they did not consider mass conservation, which is crucial for analyzing a closed reactor system. In fact, the total mass is not conserved in the previous models. This study improves the artificial viscosity model such that the 1D incompressible two-fluid model is well-posed, and the total mass is conserved. The water faucet and Kelvin-Helmholtz instability flows were simulated to test the effect of the proposed artificial viscosity model. The results indicate that the proposed artificial viscosity model effectively remedies the ill-posedness of the two-fluid model while maintaining a negligible total mass error.

• Nuclear Engineering and Technology
• /
• v.51 no.6
• /
• pp.1504-1513
• /
• 2019

Recently, ocean nuclear reactors have received attention due to enhanced safety features. The movable and transportable characteristics distinguish ocean nuclear reactors from land-based nuclear reactors. Therefore, for safety/design analysis of the ocean reactor, the thermos-hydraulics must be investigated in the moving system. However, there are no studies reporting the general two-fluid equations that can be used for multi-dimensional simulations of two-phase flows in moving systems. This study is to systematically formulate the multi-dimensional two-fluid equations in the non-inertial frame of reference. To demonstrate the applicability of the formulated equations, we perform a total of six different simulations in 2D tanks with translational and/or rotational motions.

• Nuclear Engineering and Technology
• /
• v.44 no.6
• /
• pp.623-638
• /
• 2012

The volume of fluid (VOF) model of FLUENT and the lattice Boltzmann method (LBM) are used to simulate two-phase flows. Both methods are validated for static and dynamic bubble test cases and then compared to experimental results. The VOF method does not reduce the spurious currents of the static droplet test and does not satisfy the Laplace law for small droplets at the acceptable level, as compared with the LBM. For single bubble flows, simulations are executed for various Eotvos numbers, Morton numbers and Reynolds numbers, and the results of both methods agree well with the experiments in the case of low Eotvos numbers. For high Eotvos numbers, the VOF results deviated from the experiments. For multiple bubbles, the bubble flow characteristics are related by the wake of the leading bubble. The coaxial and oblique coalescence of the bubbles are simulated successfully and the subsequent results are presented. In conclusion, the LBM performs better than the VOF method.