Two-Group Relative Velocity in Boiling Two-Phase Flow

Author(s):  
Caleb S. Brooks ◽  
Basar Ozar ◽  
Takashi Hibiki ◽  
Mamoru Ishii

In an effort to improve the prediction of void fraction and heat transfer characteristics in two-phase systems, the two-group interfacial area transport equation has been developed for use with the two-group two-fluid model. In the one-dimensional formulation, a closure relation is required for the group-1 and group-2 area-average local relative velocity. Furthermore, in the case of the modified two-fluid model with the gas-mixture momentum equation, the group-1 and group-2 void weighted gas velocities must be calculated with additional closure relations. Therefore, the drift-flux general expression is extended to two bubble groups in order to describe the group-1 and group-2 void weighted gas velocities and area-averaged local relative velocities. Correlations for group-1 and group-2 distribution parameters and drift velocities are proposed and evaluated with a dataset containing 57 boiling conditions taken in an internally heated annulus. The proposed distribution parameters show an agreement within +/−5%. The overall estimation of group-1 and group-2 void weighted gas velocities calculated with the newly proposed two-group drift-flux general expression shows an agreement within +/−15% of the measured value.

Author(s):  
Marco Germano Conte ◽  
Cristiane Cozin ◽  
Fausto Arinos Barbuto ◽  
Rigoberto E. M. Morales

Two-phase slug flow is present in many industrial processes, such as the exploitation and transportation of hydrocarbon mixtures from oil wells. This kind of flow is characterized by two distinct structures which repeat intermittently: a liquid slug with a large amount of momentum followed by a compressible gas bubble. In recent decades, a few models for simulating such complex flows were developed, as the eulerian two-fluid model and drift flux, and the lagrangian slug tracking. The aim of this work is to present a detailed study on the numerical implementation of the hybrid model proposed by Fabien Renault and Nydal which is able to track down waves that arise in the gas-liquid interface and possible slugs generated by them. This model was developed from the two-fluid model equations in which the motion generated by the dynamic pressure of the gas on the slugs is decoupled from the slow movement of the liquid below the gas. The movement of the bubbles in the liquid is then modeled similarly to shallow-water equations. The solution of the equation set is achieved in two steps. The first step provides the pressure field and the gas flow through the numerical solution of the equations for the gas, using the finite difference method. The second step solves the adapted shallow-water equations analytically. The model was coded in object-oriented Intel Visual Fortran95. Simulations to analyze the ability of the code to generate slugs for some pairs of water-air superficial velocities at atmospheric pressure were carried out. The results, as the distribution of the slug length, frequency and average values were compared to experimental results reported in the literature.


2019 ◽  
Vol 23 (5 Part B) ◽  
pp. 3003-3013
Author(s):  
Vahid Shokri ◽  
Kazem Esmaeili

In this paper, a numerical study is performed in order to investigate the effect of the liquid phase compressibility two-fluid model. The two-fluid model is solved by using conservative shock capturing method. At the first, the two-fluid model is applied by assuming that the liquid phase is incompressible, then it is assumed that in three cases called water faucet case, large relative velocity shock pipe case, and Toumi?s shock pipe case, the liquid phase is compressible. Numerical results indicate that, if an intense pressure gradient is governed on the fluid-flow, single-pressure two-fluid model by assuming liquid phase incompressibility predicts the flow variables in the solution field more accurate than single-pressure two-fluid model by assuming liquid phase compressibility.


2021 ◽  
Vol 33 (3) ◽  
pp. 033324
Author(s):  
Alejandro Clausse ◽  
Martín López de Bertodano

2013 ◽  
Author(s):  
Sung Chan Cho ◽  
Yun Wang

In this paper, two-phase flow dynamics in a micro channel with various wall conditions are both experimentally and theoretically investigated. Annulus, wavy and slug flow patterns are observed and location of liquid phase on different wall condition is visualized. The impact of flow structure on two-phase pressure drop is explained. Two-phase pressure drop is compared to a two-fluid model with relative permeability correlation. Optimization of correlation is conducted for each experimental case and theoretical solution for the flows in a circular channel is developed for annulus flow pattern showing a good match with experimental data in homogeneous channel case.


Author(s):  
David Heinze ◽  
Thomas Schulenberg ◽  
Lars Behnke

A simulation model for the direct contact condensation of steam in subcooled water is presented that allows determination of major parameters of the process, such as the jet penetration length. Entrainment of water by the steam jet is modeled based on the Kelvin–Helmholtz and Rayleigh–Taylor instability theories. Primary atomization due to acceleration of interfacial waves and secondary atomization due to aerodynamic forces account for the initial size of entrained droplets. The resulting steam-water two-phase flow is simulated based on a one-dimensional two-fluid model. An interfacial area transport equation is used to track changes of the interfacial area density due to droplet entrainment and steam condensation. Interfacial heat and mass transfer rates during condensation are calculated using the two-resistance model. The resulting two-phase flow equations constitute a system of ordinary differential equations, which is solved by means of the explicit Runge–Kutta–Fehlberg algorithm. The simulation results are in good qualitative agreement with published experimental data over a wide range of pool temperatures and mass flow rates.


Author(s):  
Hiroyuki Yoshida ◽  
Takeharu Misawa ◽  
Kazuyuki Takase

Two-fluid model can simulate two phase flow less computational cost than inter-face tracking method and particle interaction method. Therefore, two-fluid model is useful for thermal hydraulic analysis in large-scale domain such as a rod bundle. Japan Atomic Energy Agency (JAEA) develops three dimensional two-fluid model analysis code ACE-3D, which adopts boundary fitted coordinate system in order to simulate complex shape channel flow. In this paper, boiling two-phase flow analysis in a tight lattice rod bundle is performed by ACE-3D code. The parallel computation using 126CPUs is applied to this analysis. In the results, the void fraction, which distributes in outermost region of rod bundle, is lower than that in center region of rod bundle. At height z = 0.5 m, void fraction in the gap region is higher in comparison with that in center region of the subchannel. However, at height of z = 1.1m, higher void fraction distribution exists in center region of the subchannel in comparison with the gap region. The tendency of void fraction to concentrate in the gap region at vicinity of boiling starting point, and to move into subchannel as water goes through rod bundle, is qualitatively agreement with the measurement results by neutron radiography. To evaluate effects of two-phase flow model used in ACE-3D code, numerical simulation of boiling two-phase in tight lattice rod bundle with no lift force model (neglecting lift force acting on bubbles) is also performed. From the comparison of numerical results, it is concluded that the effects of lift force model are not so large on overall void fraction distribution in tight lattice rod bundle. However, higher void fraction distribution in center region of the subchannel was not observed in this simulation. It is concluded that the lift force model is important for local void fraction distribution in rod bundles.


Author(s):  
Wei Yao ◽  
Christophe Morel

In this paper, a multidimensional two-fluid model with additional turbulence k–ε equations is used to predict the two-phase parameters distribution in freon R12 boiling flow. The 3D module of the CATHARE code is used for numerical calculation. The DEBORA experiment has been chosen to evaluate our models. The radial profiles of the outlet parameters were measured by means of an optical probe. The comparison of the radial profiles of void fraction, liquid temperature, gas velocity and volumetric interfacial area at the end of the heated section shows that the multidimensional two-fluid model with proper constitutive relations can yield reasonably predicted results in boiling conditions. Sensitivity tests show that the turbulent dispersion force, which involves the void fraction gradient, plays an important role in determining the void fraction distribution; and the turbulence eddy viscosity is a significant factor to influence the liquid temperature distribution.


Author(s):  
Youn-Gyu Jung ◽  
Moon-Sun Chung ◽  
Sung-Jae Yi

This study discusses on the implementation of an upwind method for a one-dimensional two-fluid model including the surface tension effect in the momentum equations. This model consists of a complete set of six equations including two-mass, two-momentum, and two-internal energy conservation equations having all real eigenvalues. Based on this equation system with upwind numerical method, the present authors first make a pilot code and then solve some benchmark problems to verify whether this model and numerical method is able to properly solve some fundamental one-dimensional two-phase flow problems or not.


Sign in / Sign up

Export Citation Format

Share Document