Hyperbolicity Analysis of a One-Dimensional Two-Fluid Two-Phase Flow Model for Stratified-Flow Pattern

Author(s):  
Carina N. Sondermann ◽  
Rodrigo A. C. Patrício ◽  
Aline B. Figueiredo ◽  
Renan M. Baptista ◽  
Felipe B. F. Rachid ◽  
...  

Two-phase flows in pipelines occur in a variety of processes in the nuclear, petroleum and gas industries. Because of the practical importance of accurately predicting steady and unsteady flows along the line, one-dimensional two-fluid flow models have been extensively employed in numerical simulations. These models are usually written as a system of non-linear hyperbolic partial-differential equations, but some of the available formulations are physically inconsistent due to a loss of the hyperbolicity property. In these cases, the associated eigenvalues become complex numbers and the model loses physical meaning locally. This paper presents a numerical study of a one-dimensional single-pressure four-equation two-fluid model for an isothermal stratified flow that occurs in a horizontal pipeline. The diameter, pressure and volume fraction are kept constant, whereas the liquid and gas velocities are varied to cover the entire range of superficial velocities in the stratified region. For each point, the eigenvalues are numerically computed to verify whether they are real numbers and to assess their signs. The results show that hyperbolicity is lost near the boundaries of the stratified pattern and in a vast area of the region itself. Moreover, the eigenvalue signs alternate, which has implications on the prescription of numerical boundary conditions.

Author(s):  
Raphael V. N. de Freitas ◽  
Carina N. Sondermann ◽  
Rodrigo A. C. Patricio ◽  
Aline B. Figueiredo ◽  
Gustavo C. R. Bodstein ◽  
...  

Numerical simulation is a very useful tool for the prediction of physical quantities in two-phase flows. One important application is the study of oil-gas flows in pipelines, which is necessary for the proper selection of the equipment connected to the line during the pipeline design stage and also during the pipeline operation stage. The understanding of the phenomena present in this type of flow is more crucial under the occurrence of undesired effects in the duct, such as hydrate formation, fluid leakage, PIG passage, and valve shutdown. An efficient manner to model two-phase flows in long pipelines regarding a compromise between numerical accuracy and cost is the use of a one-dimensional two-fluid model, discretized with an appropriate numerical method. A two-fluid model consists of a system of non-linear partial differential equations that represent the mass, momentum and energy conservation principles, written for each phase. Depending on the two-fluid model employed, the system of equations may lose hyperbolicity and render the initial-boundary-value problem illposed. This paper uses an unconditionally hyperbolic two-fluid model for solving two-phase flows in pipelines in order to guarantee that the solution presents physical consistency. The mathematical model here referred to as the 5E2P (five equations and two pressures) comprises two equations of continuity and two momentum conservation equations, one for each phase, and one equation for the transport of the volume fraction. A priori this model considers two distinct pressures, one for each phase, and correlates them through a pressure relaxation procedure. This paper presents simulation cases for stratified two-phase flows in horizontal pipelines solved with the 5E2P coupled with the flux corrected transport method. The objective is to evaluate the numerical model capacity to adequately describe the velocities, pressures and volume fraction distributions along the duct.


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):  
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.


2003 ◽  
Vol 125 (2) ◽  
pp. 387-389 ◽  
Author(s):  
Jin Ho Song

A linear stability analysis is performed for a two-phase flow in a channel to demonstrate the feasibility of using momentum flux parameters to improve the one-dimensional two-fluid model. It is shown that the proposed model is stable within a practical range of pressure and void fraction for a bubbly and a slug flow.


2014 ◽  
Vol 62 (3) ◽  
pp. 234-240 ◽  
Author(s):  
Gianandrea Vittorio Messa ◽  
Stefano Malavasi

Abstract The flow of a mixture of liquid and solid particles at medium and high volume fraction through an expansion in a rectangular duct is considered. In order to improve the modelling of the phenomenon with respect to a previous investigation (Messa and Malavasi, 2013), use is made of a two-fluid model specifically derived for dense flows that we developed and implemented in the PHOENICS code via user-defined subroutines. Due to the lack of experimental data, the two-fluid model was validated in the horizontal pipe case, reporting good agreement with measurements from different authors for fully-suspended flows. A 3D system is simulated in order to account for the effect of side walls. A wider range of the parameters characterizing the mixture (particle size, particle density, and delivered solid volume fraction) is considered. A parametric analysis is performed to investigate the role played by the key physical mechanisms on the development of the two-phase flow for different compositions of the mixture. The main focuses are the distribution of the particles in the system and the pressure recovery


2004 ◽  
Vol 45 (10) ◽  
pp. 1049-1066 ◽  
Author(s):  
Moon-Sun Chung ◽  
Seung-Kyung Pak ◽  
Keun-Shik Chang

1997 ◽  
Vol 63 (611) ◽  
pp. 2377-2385
Author(s):  
Kiyoshi MINEMURA ◽  
Tomomi UCHIYAMA ◽  
Katsuhiko KINOSHITA ◽  
Lin LYU ◽  
Shinji SYODA ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document