Decompression wave speed in CO2 mixtures: CFD modelling with the GERG-2008 equation of state

The development of CO2 pipelines for Carbon Capture and Storage (CCS) raises new questions regarding the control of ductile fracture propagation and fracture arrest toughness criteria. The decompression behaviour in the fluid must be determined accurately in order to estimate the proper pipe toughness. However, anthropogenic CO2 may contain impurities that can modify the fluid decompression characteristics quite significantly. To determine the decompression wave speed in CO2 mixtures, the thermodynamic properties of these mixtures must be determined by using an accurate equation of state. In this paper we present a new decompression model developed using the Computational Fluid Dynamics (CFD) package ANSYS Fluent. The GERG-2008 Equation of State (EOS) was implemented into this model through User Defined Functions (UDF) to predict the thermodynamic properties of CO2 mixtures. The model predictions were in good agreement with the experimental data of two ‘shock tube’ tests. A range of representative CO2 mixtures was examined in terms of the changes in fluid properties from the initial conditions, with time and distance, immediately after a sudden pipeline opening at one end. Phase changes that may occur within the fluid due to condensation of ‘impurities’ in the fluid were also investigated. Simulations were also conducted to examine how the initial temperature and impurities would affect the decompression wave speed.

[1]  Haroun Mahgerefteh,et al.  An Interpolation Technique for Rapid CFD Simulation of Turbulent Two-Phase Flows , 2007 .

[2]  J. F. Kiefner,et al.  Ductile fracture arrest in gas pipelines , 1976 .

[3]  G. Soave.  Equilibrium constants from a modified Redlich-Kwong equation of state , 1972 .

[4]  J. R. Chen,et al.  Modelling of two-phase blowdown from pipelines—II. A simplified numerical method for multi-component mixtures , 1995 .

[5]  Julia Race,et al.  Transporting the Next Generation of CO2 for Carbon, Capture and Storage: The Impact of Impurities on Supercritical CO2 Pipelines , 2008 .

[6]  Amyn S. Teja,et al.  A NEW CUBIC EQUATION OF STATE FOR FLUIDS AND FLUID MIXTURES , 1982 .

[7]  Hailong Li,et al.  PVTxy properties of CO2 mixtures relevant for CO2 capture, transport and storage: Review of available experimental data and theoretical models , 2011 .

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

[9]  M. McLinden,et al.  NIST Standard Reference Database 23: Reference Fluid Thermodynamic and Transport Properties-REFPROP, Version 8.0 , 2007 .

[10]  M. Liou,et al.  A New Flux Splitting Scheme , 1993 .

[11]  Cheng Lu,et al.  Investigation of the Effects of Pipe Wall Roughness and Pipe Diameter on the Decompression Wave Speed in Natural Gas Pipelines , 2012 .

[12]  H Mahgerefteh,et al.  Modelling outflow following rupture in pipeline networks , 2006 .

[13]  D. Feakins,et al.  The thermodynamics of solutions , 1989 .

[14]  Jinyue Yan,et al.  Impacts of equations of state (EOS) and impurities on the volume calculation of CO2 mixtures in the applications of CO2 capture and storage (CCS) processes , 2009 .

[15]  Meihong Wang,et al.  Simulation-based techno-economic evaluation for optimal design of CO 2 transport pipeline network , 2014 .

[16]  B. Metz.  IPCC special report on carbon dioxide capture and storage , 2005 .

[17]  D. Lempe,et al.  Density improvement of the SRK equation of state , 1997 .

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

[19]  K. Botros,et al.  Extension of the semi-empirical correlation for the effects of pipe diameter and internal surface roughness on the decompression wave speed to include High Heating Value Processed Gas mixtures , 2013 .

[20]  Hailong Li,et al.  Evaluating cubic equations of state for calculation of vapor–liquid equilibrium of CO2 and CO2-mixtures for CO2 capture and storage processes , 2009 .

[21]  Xiaoyan Ji,et al.  Energy consumption analysis for CO2 separation from gas mixtures , 2014 .

[22]  S. F. Brown,et al.  CFD modelling of outflow and ductile fracture propagation in pressurised pipelines , 2011 .

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

[24]  P. Bishnoi,et al.  Calculation of the thermodynamic sound velocity in two-phase multicomponent fluids , 1987 .

[25]  Andrew Cosham,et al.  GASDECOM: Carbon Dioxide and Other Components , 2010 .

[26]  Heqin Xu,et al.  Effect of Grade on Ductile Fracture Arrest Criteria for Gas Pipelines , 2006 .

[27]  Julian Barnett,et al.  The Decompression Behaviour of Carbon Dioxide in the Dense Phase , 2012 .

[28]  B. Lu,et al.  Simultaneous calculations of VLE and saturated liquid and vapor volumes by means of a 3P1T cubic EOS , 1987 .

[29]  K. K. Botros,et al.  Measurements of Speed of Sound in Lean and Rich Natural Gas Mixtures at Pressures up to 37 MPa Using a Specialized Rupture Tube , 2010 .

[30]  Zhenhao Duan,et al.  PVTx properties of the CO2–H2O and CO2–H2O–NaCl systems below 647 K: Assessment of experimental data and thermodynamic models , 2007 .

[31]  J. Pain,et al.  Gas Dynamics , 2018, Nature.

[32]  Haroun Mahgerefteh,et al.  A study of the effects of friction, heat transfer, and stream impurities on the decompression behavior in CO2 pipelines , 2012 .

[33]  Andre Peneloux,et al.  A consistent correction for Redlich-Kwong-Soave volumes , 1982 .

[34]  O. Redlich,et al.  On the thermodynamics of solutions; an equation of state; fugacities of gaseous solutions. , 1949, Chemical reviews.

[35]  Kenneth R. Hall,et al.  An algebraic method that includes Gibbs minimization for performing phase equilibrium calculations for any number of components or phases , 2003 .

[36]  S. H. Hashemi,et al.  Correction factors for safe performance of API X65 pipeline steel , 2009 .

[37]  P. Bishnoi,et al.  The importance of real‐fluid behavior and nonisentropic effects in modeling decompression characteristics of pipeline fluids for application in ductile fracture propagation analysis , 1988 .

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

[39]  Thomas J. Ahrens,et al.  Equation of State , 1993 .

[40]  J. Kiefner,et al.  Failure Stress Levels of Flaws in Pressurized Cylinders , 1973 .

[41]  Arne Dugstad,et al.  Update of DNV recommended practice RP-J202 with focus on CO2 Corrosion with Impurities , 2014 .

[42]  P. Bishnoi,et al.  The importance of real‐fluid behavior in predicting release rates resulting from high‐pressure sour‐gas pipeline ruptures , 1989 .

[43]  Xiong Liu,et al.  Source strength and dispersion of CO2 releases from high-pressure pipelines: CFD model using real gas equation of state , 2014 .

[44]  J. X. Wen,et al.  Predicting the Decompression Characteristics of Carbon Dioxide Using Computational Fluid Dynamics , 2012 .

[45]  Performance of five equations of state for the prediction of vle and densities of natural gas mixtures in the dense phase region , 2002 .

[46]  Haroun Mahgerefteh,et al.  Fast numerical simulation for full bore rupture of pressurized pipelines , 1999 .

[47]  K. E. Starling,et al.  Enthalpy of Mixtures by Modified BWR Equation , 1970 .