A Hyperbolicity Analysis of the 1991 OLGA’s Model for Isothermal Flow

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

Multiphase flows are encountered in many engineering problems. Particularly in the oil and gas industry, many applications involve the transportation of a mixture of oil and natural gas in long pipelines from offshore platforms to the continent. Numerical simulations of steady and unsteady flows in pipelines are usually based on one-dimensional models, such as the two-fluid model, the drift-flux model and the homogeneous equilibrium model. The 1991’s version of the well-known and widely-used commercial software OLGA describes a system of non-linear equations of the two-fluid-model type, with an extra equation for the presence of liquid droplets. It is well known that one-dimensional formulations may be physically inconsistent due to the loss of hyperbolicity. In these cases, the associated eigenvalues become complex numbers and the model loses physical meaning locally. This paper presents a numerical study of the 1991’s version of the software OLGA, for an isothermal flow of stratified pattern, in a horizontal pipeline. For each point of interest in the stratified-pattern flow map, the eigenvalues are numerically calculated in order to verify if the eigenvalues are real and also to assess their signs. The results indicate that the model is conditionally hyperbolic and loses hyperbolicity in a vast area of the stratified region under certain flow conditions. Even though the model is not unconditionally hyperbolic, some simulations here performed for typical offshore pipeline flows are shown to be in the hyperbolic region.

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):  
Andrew Lehmann ◽  
Mark Wardle

AbstractWe characterise steady, one-dimensional fast and slow magnetohydrodynamic (MHD) shocks using a two-fluid model. Fast MHD shocks are magnetically driven, forcing ions to stream through the neutral gas ahead of the shock front. This magnetic precursor heats the gas sufficiently to create a large, warm transition zone where all fluid variables only weakly change in the shock front. In contrast, slow MHD shocks are driven by gas pressure where neutral species collide with ion species in a thin hot slab that closely resembles an ordinary gas dynamic shock.We computed observational diagnostics for fast and slow shocks at velocities vs=2–4 km/s and preshock Hydrogen nuclei densities nH = 102-4 cm−3. We followed the abundances of molecules relevant for a simple oxygen chemistry and include cooling by CO, H2 and H2O. Estimates of intensities of 12CO rotational lines show that high-J lines, above J = 6 → 5, are more strongly excited in slow MHD shocks.


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.


2019 ◽  
Vol 97 (12) ◽  
pp. 3156-3174
Author(s):  
Nannan Liu ◽  
Wei Wang ◽  
Yingying Liu ◽  
Liang Ma ◽  
Jing Gong

Processes ◽  
2020 ◽  
Vol 8 (7) ◽  
pp. 810
Author(s):  
Jade Gesare Abuga ◽  
Tiri Chinyoka

The flow of viscoelastic fluids may, under certain conditions, exhibit shear-banding characteristics that result from their susceptibility to unusual flow instabilities. In this work, we explore both the existing shear banding mechanisms in the literature, namely; constitutive instabilities and flow-induced inhomogeneities. Shear banding due to constitutive instabilities is modelled via either the Johnson–Segalman or the Giesekus constitutive models. Shear banding due to flow-induced inhomogeneities is modelled via the Rolie–Poly constitutive model. The Rolie–Poly constitutive equation is especially chosen because it expresses, precisely, the shear rheometry of polymer solutions for a large number of strain rates. For the Rolie–Poly approach, we use the two-fluid model wherein the stress dynamics are coupled with concentration equations. We follow a computational analysis approach via an efficient and versatile numerical algorithm. The numerical algorithm is based on the Finite Volume Method (FVM) and it is implemented in the open-source software package, OpenFOAM. The efficiency of our numerical algorithms is enhanced via two possible stabilization techniques, namely; the Log-Conformation Reformulation (LCR) and the Discrete Elastic Viscous Stress Splitting (DEVSS) methodologies. We demonstrate that our stabilized numerical algorithms accurately simulate these complex (shear banded) flows of complex (viscoelastic) fluids. Verification of the shear-banding results via both the Giesekus and Johnson-Segalman models show good agreement with existing literature using the DEVSS technique. A comparison of the Rolie–Poly two-fluid model results with existing literature for the concentration and velocity profiles is also in good agreement.


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.


Sign in / Sign up

Export Citation Format

Share Document