Numerical Modeling of Two-Phase Flows Using Advanced Two Fluid Systems

Author(s):  
Anela Kumbaro ◽  
Imad Toumi ◽  
Vincent Seignole

The purpose of this paper is to report on the development and assessment of approximate Riemann solver methods for the discretization of non-linear non-conservative systems arising in the simulation of two-phase flows. These methods are able to treat general two-phase flow systems with realistic state equations and are flexible enough to be applied on any mesh type, structured as well as unstructured. We will detail models that go from the basic 6 equation two-fluid model to the coupling of this system with one or more transport equations, for instance on volumetric interfacial area concentration, or on partial void fractions of groups of bubbles (MUlti-Size-Group model). This kind of transport equation is useful to predict at a finer level the interfacial patterns or bubble size distribution and takes account of coalescence or breakup rates of inclusions. We make a glimpse at the choices made regarding this aspect. Different physico-numerical benchmarks are provided in order to illustrate the numerical and physical modeling. Confrontation with experimental or analytical reference data are performed whenever possible. Computer simulations are performed using OVAP, a new multidimensional CFD code.

Author(s):  
Mamoru Ishii ◽  
Seungjin Kim ◽  
Xiaodong Sun ◽  
Takashi Hibiki

A dynamic treatment of interfacial area concentration has been studied over the last decade by employing the interfacial area transport equation. When coupled with the two-fluid model, the interfacial area transport equation replaces the flow regime dependent correlations for interfacial area concentration and eliminates potential artificial bifurcation or numerical oscillations stemming from these static correlations. An extensive database has been established to evaluate the model under various two-phase flow conditions. These include adiabatic and heated conditions, vertical and horizontal flow orientations, round, rectangular, annulus, and 8×8 rod-bundle channel geometries, and normal-gravity and reduced-gravity conditions. Currently, a two-group interfacial area transport equation is available and applicable to comprehensive two-phase flow conditions spanning from bubbly to churn-turbulent flow regimes. A framework to couple the two-group interfacial area transport equation with the modified two-fluid model is established in view of multiphase computational fluid dynamics code applications as well as reactor system analysis code applications. The present study reviews the current state-of-the-art in the development of the interfacial area transport equation, available experimental databases, and the analytical methods to incorporate the interfacial area transport equation into the two-fluid model.


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.


2016 ◽  
Vol 138 (8) ◽  
Author(s):  
Avinash Vaidheeswaran ◽  
Martin Lopez de Bertodano

Analytical expressions for interfacial pressure coefficients are obtained based on the geometry of the bubbles occurring in two-phase flows. It is known that the shape of the bubbles affects the virtual mass and interfacial pressure coefficients, which in turn determines the cutoff void fraction for the well-posedness of two-fluid model (TFM). The coefficient used in the interfacial pressure difference correlation is derived assuming potential flow around a perfect sphere. In reality, the bubbles seen in two-phase flows get deformed, and hence, it is required to estimate the coefficients for nonspherical geometries. Oblate and prolate ellipsoids are considered, and their respective coefficients are determined. It is seen that the well-posedness limit of the TFM is determined by the combination of virtual mass and interfacial pressure coefficient used. The effect of flow separation on the coefficient values is also analyzed.


Author(s):  
Xiaodong Sun ◽  
Yang Liu ◽  
Basar Ozar ◽  
Mamoru Ishii ◽  
Joseph M. Kelly

To apply the two-fluid model to a wide range of flow regimes in gas-liquid two-phase flows, the gas phase is categorized into two groups: small spherical/distorted bubbles as Group 1 and large cap/slug/churn-turbulent bubbles as Group 2 in the modeling of interfacial area transport. The interfacial transfer terms of momentum and energy for the gas phase are then divided into two groups accordingly in the implementation of the two-group interfacial area transport equation to the two-fluid model. Thus, the drag coefficients and the interfacial heat transfer for each group bubbles need to be developed. An approach has been sought for evaluating the drag coefficients of each bubble group based on a comprehensive experimental data base obtained in air-water upward flows in various size round pipes. Comparisons have been made with the theory of the drag coefficients and it was found that the agreement is not very satisfactory although the general trends can be predicted by the current approach.


Sign in / Sign up

Export Citation Format

Share Document