Thermodynamic interpolation for the simulation of two-phase flow of non-ideal mixtures

This paper describes the development and application of a technique for the rapid interpolation of thermodynamic properties of mixtures for the purposes of simulating two-phase flow. The technique is based on adaptive inverse interpolation and can be applied to any Equation of State and multicomponent mixture. Following analysis of its accuracy, the method is coupled with a two-phase flow model, based on the homogeneous equilibrium mixture assumption, and applied to the simulation of flows of carbon dioxide (CO2) rich mixtures. This coupled flow model is used to simulate the experimental decompression of binary and quinternary mixtures. It is found that the predictions are in good agreement with the experimental data and that the interpolation approach provides a flexible, robust means of obtaining thermodynamic properties for use in flow models.

[1]  Georgios M. Kontogeorgis,et al.  Thermodynamic Models for Industrial Applications , 2010 .

[2]  Hiroshi Akima,et al.  Algorithm 761: Scattered-data surface fitting that has the accuracy of a cubic polynomial , 1996, TOMS.

[3]  M. Michelsen.  The isothermal flash problem. Part I. Stability , 1982 .

[4]  Denis V. Voskov,et al.  Fully compositional and thermal reservoir simulation , 2014, Comput. Chem. Eng..

[5]  Haroun Mahgerefteh,et al.  Efficient numerical solution for highly transient flows , 2006 .

[6]  Solomon F. Brown,et al.  Simulation of two-phase flow through ducts with discontinuous cross-section , 2015 .

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

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

[9]  Hisham N.H. Saadawi,et al.  Thermodynamic and transport property models for carbon capture and sequestration (CCS) processes with emphasis on CO2 transport , 2013 .

[10]  Saeedeh Mohebbinia.  Advanced equation of state modeling for compositional simulation of gas floods , 2013 .

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

[12]  M. Pourkashanian,et al.  The range and level of impurities in CO2 streams from different carbon capture sources , 2015 .

[13]  Hamdi A. Tchelepi,et al.  Generalized negative-flash method for multiphase multicomponent systems , 2010 .

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

[15]  M. Heikal,et al.  Numerical Simulation of Cavitation Flows Based on Their Hydrodynamic Similarity , 2006 .

[16]  M. Michelsen.  The isothermal flash problem. Part II. Phase-split calculation , 1982 .

[17]  J. Houseworth.  Comment on "Potential for environmental impact due to acid gas leakage from wellbores at EOR injection sites near Zama Lake, Alberta" by D.M. LeNeveu (2012) , 2012 .

[18]  Trond Andresen,et al.  Mathematical modeling of CO2 based heat pumping systems: New developments for simulation tools to aid the design of systems for non-residential buildings , 2009 .

[19]  Michael Dumbser,et al.  Efficient implementation of high order unstructured WENO schemes for cavitating flows , 2013 .

[20]  Ioannis G. Economou,et al.  Evaluation of Cubic, SAFT, and PC-SAFT Equations of State for the Vapor–Liquid Equilibrium Modeling of CO2 Mixtures with Other Gases , 2013 .

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

[22]  Eric S. Fraga,et al.  A geometrically based grid refinement technique for multiphase flows , 2015, Comput. Chem. Eng..

[23]  Eric S. Fraga,et al.  Global sensitivity analysis of the impact of impurities on CO2 pipeline failure , 2013, Reliab. Eng. Syst. Saf..

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

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

[26]  E. Stenby,et al.  PVT modeling of reservoir fluids using PC-SAFT EoS and Soave-BWR EoS , 2015 .

[27]  The isoenergetic-isochoric flash , 1997 .

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

[29]  Svend Tollak Munkejord,et al.  Depressurization of CO2-rich mixtures in pipes: Two-phase flow modelling and comparison with experiments , 2015 .

[30]  O. Jones,et al.  Nucleation and flashing in nozzles—2. Comparison with experiments using a five-equation model for vapor void development , 1993 .

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

[32]  Gabriele Sadowski,et al.  Perturbed-Chain SAFT: An Equation of State Based on a Perturbation Theory for Chain Molecules , 2001 .

[33]  Marcelo Castier,et al.  Solution of the isochoric–isoenergetic flash problem by direct entropy maximization , 2009 .

[34]  Clark W. Bullard,et al.  Pressure drop and flashing mechanisms in refrigerant expansion devices , 2003 .

[35]  G. Kontogeorgis,et al.  Thermodynamic Models for Industrial Applications: From Classical and Advanced Mixing Rules to Association Theories , 2010 .

[36]  N. Chen.  An Explicit Equation for Friction Factor in Pipe , 1979 .