Modelling the non-equilibrium two-phase flow during depressurisation of CO2 pipelines

Abstract The development, testing and validation of a two-fluid transient flow model for simulating outflow following the failure of high pressure CO 2 pipelines is presented. Thermal and mechanical non-equilibrium effects during depressurisation are accounted for by utilising simple constitutive relations describing inter-phase mass, heat and momentum transfer in terms of relaxation to equilibrium. Pipe wall/fluid heat exchange on the other hand is modelled by coupling the fluid model with a finite difference transient heat conduction model. The two-fluid transient flow model's performance is tested by comparison of the predicted transient pressure and temperature profiles along the pipeline against those based on the simplified homogeneous equilibrium model (HEM) as well as real data captured during the full bore rupture of a 260 m long, 233 mm internal diameter pipeline containing CO 2 at 36 bara and 273 °C. The two-fluid model is found to produce a reasonably good degree of agreement with the experimental data throughout the depressurisation process. The HEM based flow model on the other hand performs well only near the rupture plane and during the early stages of the depressurisation process.

[1]  Keun-Shik Chang,et al.  A modified HLLC-type Riemann solver for the compressible six-equation two-fluid model , 2013 .

[2]  S. T. Munkejord.  Comparison of Roe-type methods for solving the two-fluid model with and without pressure relaxation , 2007 .

[3]  F. Dittus,et al.  Heat transfer in automobile radiators of the tubular type , 1930 .

[4]  J. R. García-Cascales,et al.  Advances in the Characterization of Transient Two-Phase Flow , 2007 .

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

[6]  Svend Tollak Munkejord,et al.  Thermo- and fluid-dynamical modelling of two-phase multi-component carbon dioxide mixtures , 2010 .

[7]  Henrik Kruse,et al.  Calculating the consequences of a CO2-pipeline rupture , 1996 .

[8]  Ioannis G. Economou,et al.  Evaluation of Statistical Associating Fluid Theory (SAFT) and Perturbed Chain-SAFT Equations of State for the Calculation of Thermodynamic Derivative Properties of Fluids Related to Carbon Capture and Sequestration , 2011 .

[9]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

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

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

[12]  J. Thome,et al.  New prediction methods for CO2 evaporation inside tubes: Part I – A two-phase flow pattern map and a flow pattern based phenomenological model for two-phase flow frictional pressure drops , 2008 .

[13]  David I. Ketcheson,et al.  Highly Efficient Strong Stability-Preserving Runge-Kutta Methods with Low-Storage Implementations , 2008, SIAM J. Sci. Comput..

[14]  M. Liou.  A Sequel to AUSM , 1996 .

[15]  Haroun Mahgerefteh,et al.  A homogeneous relaxation flow model for the full bore rupture of dense phase CO2 pipelines , 2013 .

[16]  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 .

[17]  H Mahgerefteh,et al.  Modelling the impact of stream impurities on ductile fractures in CO2 pipelines , 2012 .

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

[19]  André Faaij,et al.  A state-of-the-art review of techno-economic models predicting the costs of CO2 pipeline transport , 2013 .

[20]  Simon Roussanaly,et al.  Costs benchmark of CO2 transport technologies for a group of various size industries , 2013 .

[21]  Don W. Green,et al.  Perry's Chemical Engineers' Handbook , 2007 .

[22]  Meng-Sing Liou,et al.  A sequel to AUSM, Part II: AUSM+-up for all speeds , 2006, J. Comput. Phys..

[23]  Mihaela Popescu.  Modeling of fluid dynamics interacting with ductile fraction propagation in high pressure pipeline , 2009 .

[24]  R. Reid,et al.  The Properties of Gases and Liquids , 1977 .

[25]  J. Cortes.  On the construction of upwind schemes for non-equilibrium transient two-phase flows , 2002 .

[26]  M. Mølnvik,et al.  Thermo- and fluid-dynamical modeling of two-phase multicomponent carbon dioxide mixtures , 2009 .

[27]  A. Terenzi,et al.  Influence of Real-Fluid Properties in Modeling Decompression Wave Interacting with Ductile Fracture Propagation , 2005 .

[28]  R. Pletcher,et al.  Computational Fluid Mechanics and Heat Transfer. By D. A ANDERSON, J. C. TANNEHILL and R. H. PLETCHER. Hemisphere, 1984. 599 pp. $39.95. , 1986, Journal of Fluid Mechanics.

[29]  John R. Thome,et al.  State-of-the-art of two-phase flow and flow boiling heat transfer and pressure drop of CO2 in macro- and micro-channels. , 2005 .

[30]  U. Ziegler,et al.  The NIRVANA code: Parallel computational MHD with adaptive mesh refinement , 2008, Comput. Phys. Commun..

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

[32]  John B. Bell,et al.  Parallelization of structured, hierarchical adaptive mesh refinement algorithms , 2000 .

[33]  Yang-Yao Niu,et al.  A further work on multi‐phase two‐fluid approach for compressible multi‐phase flows , 2008 .

[34]  Danièle Revel,et al.  IPCC Special Report on Renewable Energy Sources and Climate Change Mitigation , 2011 .

[35]  D. Peng,et al.  A New Two-Constant Equation of State , 1976 .

[36]  Meng-Sing Liou,et al.  How to solve compressible multifluid equations: a simple, robust and accurate method , 2007 .

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