Analysis of time integration methods for the compressible two-fluid model for pipe flow simulations

Abstract In this paper we analyse different time integration methods for the two-fluid model and propose the BDF2 method as the preferred choice to simulate transient compressible multiphase flow in pipelines. Compared to the prevailing Backward Euler method, the BDF2 scheme has a significantly better accuracy (second order) while retaining the important property of unconditional linear stability ( A -stability). In addition, it is capable of damping unresolved frequencies such as acoustic waves present in the compressible model ( L -stability), opposite to the commonly used Crank–Nicolson method. The stability properties of the two-fluid model and of several discretizations in space and time have been investigated by eigenvalue analysis of the continuous equations, of the semi-discrete equations, and of the fully discrete equations. A method for performing an automatic von Neumann stability analysis is proposed that obtains the growth rate of the discretization methods without requiring symbolic manipulations and that can be applied without detailed knowledge of the source code. The strong performance of BDF2 is illustrated via several test cases related to the Kelvin–Helmholtz instability. A novel concept called Discrete Flow Pattern Map (DFPM) is introduced which describes the effective well-posed unstable flow regime as determined by the discretization method. Backward Euler introduces so much numerical diffusion that the theoretically well-posed unstable regime becomes numerically stable (at practical grid and timestep resolution). BDF2 accurately identifies the stability boundary, and reveals that in the nonlinear regime ill-posedness can occur when starting from well-posed unstable solutions. The well-posed unstable regime obtained in nonlinear simulations is therefore in practice much smaller than the theoretical one, which might severely limit the application of the two-fluid model for simulating the transition from stratified flow to slug flow. This should be taken very seriously into account when interpreting results from any slug-capturing simulations.

[1]  Q. Tran,et al.  Transient simulation of two-phase flows in pipes , 1998 .

[2]  V. H. Ransom,et al.  Linear and nonlinear analysis of an unstable, but well-posed, one-dimensional two-fluid model for two-phase flow based on the inviscid Kelvin–Helmholtz instability , 2014 .

[3]  Eugene Kazantsev,et al.  Parameterizing subgrid scale eddy effects in a shallow water model , 2017 .

[4]  H. Holmås,et al.  Numerical simulation of transient roll-waves in two-phase pipe flow , 2010 .

[5]  Andreas Holm Akselsen,et al.  Efficient Numerical Methods for Waves in One-Dimensional Two-Phase Pipe Flows , 2016 .

[6]  Dag Biberg,et al.  An explicit approximation for the wetted angle in two‐Phase stratified pipe flow , 1999 .

[7]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

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

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

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

[11]  M. Montini.  Closure relations of the one-dimensional two-fluid model for the simulation of slug flows , 2011 .

[12]  Renan Martins Baptista,et al.  Numerical simulation of stratified-pattern two-phase flow in gas pipelines using a two-fluid model , 2017 .

[13]  W. Fullmer,et al.  ONE-DIMENSIONAL TWO-EQUATION TWO-FLUID MODEL STABILITY , 2013 .

[14]  Jin Ho Song,et al.  A remedy for the ill-posedness of the one-dimensional two-fluid model , 2003 .

[15]  M. Ishii,et al.  The one-dimensional two-fluid model with momentum flux parameters , 2001 .

[16]  J. Butcher.  Numerical Methods for Ordinary Differential Equations: Butcher/Numerical Methods , 2005 .

[17]  R. Mei,et al.  A study on the numerical stability of the two-fluid model near ill-posedness , 2008 .

[18]  John A. Trapp,et al.  Characteristics, Stability, and Short-Wavelength Phenomena in Two-Phase Flow Equation Systems , 1978 .

[19]  R. I. Issa,et al.  Simulation of slug flow in horizontal and nearly horizontal pipes with the two-fluid model , 2003 .

[20]  Victor H. Ransom,et al.  Hyperbolic two-pressure models for two-phase flow , 1984 .

[21]  Yehuda Taitel,et al.  Interfacial and structural stability of separated flow , 1994 .

[22]  James H. VanZwieten,et al.  Comparison of numerical methods for slug capturing with the two-fluid model , 2016 .

[23]  Haihua Zhao,et al.  RELAP-7 Theory Manual , 2015 .

[24]  E. D. Hughes,et al.  Characteristics and Stability Analyses of Transient One-Dimensional Two-Phase Flow Equations and Their Finite Difference Approximations , 1978 .

[25]  Victor H. Ransom,et al.  Hyperbolic two-pressure models for two-phase flow revisited , 1988 .

[26]  Steinar Evje,et al.  Hybrid flux-splitting schemes for a common two-fluid model , 2003 .

[27]  Haihua Zhao,et al.  Applications of high-resolution spatial discretization scheme and Jacobian-free Newton–Krylov method in two-phase flow problems , 2015 .

[28]  Alistair Fitt,et al.  The numerical and analytical solution of ill-posed systems of conservation laws , 1989 .

[29]  Yehuda Taitel,et al.  Kelvin-Helmholtz stability criteria for stratified flow: viscous versus non-viscous (inviscid) approaches , 1993 .

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