A method for simulating two-phase pipe flow with real equations of state

Abstract Common two-fluid models for pipe flow assume local non-equilibrium regarding phase transfer. To solve the two-fluid models together with accurate equations of state for real fluids will in most cases require mechanical, thermal and chemical equilibrium between the phases. The reason is that reference equations of state for real substances typically describe full thermodynamic equilibrium. In this paper, we present a method for numerically solving an equilibrium model analysed by Morin and Flatten in the paper A two-fluid four-equation model with instantaneous thermodynamical equilibrium, 2013. The four-equation two-fluid model with instantaneous thermodynamical equilibrium is derived from a five-equation two-fluid model with instantaneous thermal equilibrium. The four-equation model has one mass equation common for both phases, but allows for separate phasic velocities. For comparison, the five-equation two-fluid model is numerically solved, using source terms to impose thermodynamical equilibrium. These source terms are solved using a fractional-step method. We employ the highly accurate Span–Wagner equation of state for CO 2 , and use the simple and robust FORCE scheme with MUSCL slope limiting. We demonstrate that second-order accuracy may be achieved for smooth solutions, whereas the first-order version of the scheme even allows for a robust transition to single-phase flow, also in the presence of instantaneous phase equilibrium.

[1]  D. Bestion,et al.  The physical closure laws in the CATHARE code , 1990 .

[2]  Stanley Osher,et al.  Convergence of Generalized MUSCL Schemes , 1985 .

[3]  Svend Tollak Munkejord,et al.  Solution of the Span–Wagner Equation of State Using a Density–Energy State Function for Fluid-Dynamic Simulation of Carbon Dioxide , 2012 .

[4]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[5]  Frédéric Coquel,et al.  A Numerical Method Using Upwind Schemes for the Resolution of Two-Phase Flows , 1997 .

[6]  J. H. Stuhmiller,et al.  The influence of interfacial pressure forces on the character of two-phase flow model equations , 1977 .

[7]  Christophe Eric Corre,et al.  On the extension of the AUSM+ scheme to compressible two-fluid models , 2003 .

[8]  H. Bruce Stewart,et al.  Two-phase flow: Models and methods , 1984 .

[9]  Steinar Evje,et al.  A MUSTA Scheme for a Nonconservative Two-Fluid Model , 2009, SIAM J. Sci. Comput..

[10]  Richard Saurel,et al.  Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixtures , 2009, J. Comput. Phys..

[11]  Richard Saurel,et al.  A multiphase model for compressible flows with interfaces, shocks, detonation waves and cavitation , 2001, Journal of Fluid Mechanics.

[12]  Svend Tollak Munkejord,et al.  A new coupled fluid-structure modeling methodology for running ductile fracture , 2012 .

[13]  Raymond F. Bishop,et al.  Thermo-fluid Dynamic Theory of Two-Phase Flow , 1975 .

[14]  Tore Flåtten,et al.  Relaxation two-phase flow models and the subcharacteristic condition , 2011 .

[15]  R. Abgrall,et al.  A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows , 1999 .

[16]  Svend Tollak Munkejord,et al.  On the effect of temperature and velocity relaxation in two-phase flow models , 2012 .

[17]  W. Wagner,et al.  A New Equation of State for Carbon Dioxide Covering the Fluid Region from the Triple‐Point Temperature to 1100 K at Pressures up to 800 MPa , 1996 .

[18]  Rémi Abgrall,et al.  A comment on the computation of non-conservative products , 2010, J. Comput. Phys..

[19]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[20]  Peder Aursand,et al.  ON THE DISPERSIVE WAVE DYNAMICS OF 2 × 2 RELAXATION SYSTEMS , 2012 .

[21]  A. Robinson,et al.  On the practical importance of the SSP property for Runge–Kutta time integrators for some common Godunov‐type schemes , 2005 .

[22]  Svend Tollak Munkejord,et al.  Method Using a Density–Energy State Function with a Reference Equation of State for Fluid-Dynamics Simulation of Vapor–Liquid–Solid Carbon Dioxide , 2013 .

[23]  O. Metayer,et al.  Élaboration des lois d'état d'un liquide et de sa vapeur pour les modèles d'écoulements diphasiques Elaborating equations of state of a liquid and its vapor for two-phase flow models , 2004 .

[24]  A. Kumbaro,et al.  An Approximate Linearized Riemann Solver for a Two-Fluid Model , 1996 .

[25]  Randi Moe,et al.  The dynamic two-fluid model OLGA; Theory and application , 1991 .

[26]  Svend Tollak Munkejord,et al.  A Roe scheme for a compressible six‐equation two‐fluid model , 2013 .

[27]  Manuel Jesús Castro Díaz,et al.  Why many theories of shock waves are necessary: Convergence error in formally path-consistent schemes , 2008, J. Comput. Phys..

[28]  A. Morin,et al.  A Two-Fluid Four-Equation Model with Instantaneous Thermodynamical Equilibrium , 2016 .

[29]  Halvor Lund,et al.  A Hierarchy of Relaxation Models for Two-Phase Flow , 2012, SIAM J. Appl. Math..

[30]  R. Abgrall.  How to Prevent Pressure Oscillations in Multicomponent Flow Calculations , 1996 .

[31]  Michael L. Michelsen,et al.  State function based flash specifications , 1999 .

[32]  Rémi Abgrall,et al.  Modelling phase transition in metastable liquids: application to cavitating and flashing flows , 2008, Journal of Fluid Mechanics.

[33]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .