On the impact of using volume as an independent variable for the solution of P-T fluid-phase equilibrium with equations of state

a b s t r a c t The constant pressure–temperature (P–T) flash plays an important role in the modelling of fluid-phase behaviour, and its solution is especially challenging for equations of state in which the volume is expressed as an implicit function of the pressure. We explore the relative merits of solving the P–T flash in two ensembles: mole numbers, pressure and temperature, in which each free-energy evaluation requires the use of a numerical solver; and mole numbers, volume and temperature, in which a direct evaluation of the free-energy is possible. We examine the performance of two algorithms, HELD (Helmholtz free energy Lagrangian dual), introduced in Pereira et al. (2012), and GILD (Gibbs free energy Lagrangian dual), introduced here, for the fluid-phase equilibria of 8 mixtures comprising up to 10 components, using two equations of state. While the reliability of both algorithms is comparable, the computational cost of HELD is consistently lower; this difference becomes increasingly pronounced as the number of components is increased. © 2014 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license

[1]  H. H. Rachford,et al.  Procedure for Use of Electronic Digital Computers in Calculating Flash Vaporization Hydrocarbon Equilibrium , 1952 .

[2]  S. Gómez,et al.  Efficient location of multiple global minima for the phase stability problem , 2009 .

[3]  George Jackson,et al.  New reference equation of state for associating liquids , 1990 .

[4]  M. Stadtherr,et al.  Reliable Computation of Phase Stability and Equilibrium from the SAFT Equation of State , 2002 .

[5]  Michael A. Saunders,et al.  User''s guide for NPSOL (Ver-sion 4.0): A FORTRAN package for nonlinear programming , 1984 .

[6]  Amparo Galindo Lowri A. Davies Alej.  The thermodynamics of mixtures and the corresponding mixing rules in the SAFT-VR approach for potentials of variable range , 1998 .

[7]  G. Jackson,et al.  Modelling the fluid phase behaviour of aqueous mixtures of multifunctional alkanolamines and carbon dioxide using transferable parameters with the SAFT-VR approach , 2012 .

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

[9]  Claire S. Adjiman,et al.  The HELD algorithm for multicomponent, multiphase equilibrium calculations with generic equations of state , 2012, Comput. Chem. Eng..

[10]  M. Assael,et al.  Thermophysical Properties of Fluids: An Introduction to Their Prediction , 1996 .

[11]  P. Gill,et al.  Some theoretical properties of an augmented lagrangian merit function , 1986 .

[12]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[13]  Ignacio E. Grossmann,et al.  Computers and Chemical Engineering , 2014 .

[14]  M. Donohue,et al.  Extension of the associated perturbed anisotropic chain theory to mixtures with more than one associating component , 1988 .

[15]  Will Tribbey,et al.  Numerical Recipes: The Art of Scientific Computing (3rd Edition) is written by William H. Press, Saul A. Teukolsky, William T. Vetterling, and Brian P. Flannery, and published by Cambridge University Press, © 2007, hardback, ISBN 978-0-521-88068-8, 1235 pp. , 1987, SOEN.

[16]  Ioannis G. Economou,et al.  Mean field calculations of thermodynamic properties of supercritical fluids , 1990 .

[17]  A. Bonilla-Petriciolet,et al.  Performance of Stochastic Global Optimization Methods in the Calculation of Phase Stability Analyses for Nonreactive and Reactive Mixtures , 2006 .

[18]  George Jackson,et al.  Statistical associating fluid theory for chain molecules with attractive potentials of variable range , 1997 .

[19]  M. Stadtherr,et al.  Reliable Nonlinear Parameter Estimation in VLE Modeling , 2000 .

[20]  Claire S. Adjiman,et al.  Optimal Solvent Design for Batch Separation Based on Economic Performance , 2003 .

[21]  Angelo Lucia,et al.  Simple Process Equations, Fixed-Point Methods, and Chaos , 1990 .

[22]  K. Kobe.  The properties of gases and liquids , 1959 .

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

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

[25]  P. I. Barton,et al.  A dual extremum principle in thermodynamics , 2007 .

[26]  B. Ivanov,et al.  Phase Stability Analysis with Equations of State—A Fresh Look from a Different Perspective , 2013 .

[27]  Philip E. Gill,et al.  Numerically stable methods for quadratic programming , 1978, Math. Program..

[28]  Rafiqul Gani,et al.  Are safe results obtained when the PC-SAFT equation of state is applied to ordinary pure chemicals? , 2010 .

[29]  A. S. Cullick,et al.  New strategy for phase equilibrium and critical point calculations by thermodynamic energy analysis. Part II. Critical point calculations , 1991 .

[30]  George Jackson,et al.  SAFT: Equation-of-state solution model for associating fluids , 1989 .

[31]  Amparo Galindo,et al.  Experimental and molecular modelling study of the three-phase behaviour of (propane + carbon dioxide + water) at reservoir conditions , 2013 .

[32]  M. Donohue,et al.  A simple equation of state for non-spherical and associating molecules , 1990 .

[33]  Claire S. Adjiman,et al.  Fluid phase stability and equilibrium calculations in binary mixtures: Part I: Theoretical development for non-azeotropic mixtures , 2009 .

[34]  Michael A. Saunders,et al.  SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization , 2002, SIAM J. Optim..

[35]  Amparo Galindo,et al.  Experimental and molecular modeling study of the three-phase behavior of (n-decane + carbon dioxide + water) at reservoir conditions. , 2011, The journal of physical chemistry. B.

[36]  Frances E. Pereira,et al.  A duality-based optimisation approach for the reliable solution of (P, T) phase equilibrium in volume-composition space , 2010 .

[37]  R. J. Topliss,et al.  Computational aspects of a non-cubic equation of state for phase-equilibrium calculations. Effect of density-dependent mixing rules , 1988 .

[38]  R. P. Marques,et al.  Modeling and analysis of the isothermal flash problem and its calculation with the simulated annealing algorithm , 2001 .

[39]  Mark A. Stadtherr,et al.  Reliable prediction of phase stability using an interval Newton method , 1996 .

[40]  M. Donohue,et al.  Thermodynamics of hydrogen‐bonded molecules: The associated perturbed anisotropic chain theory , 1986 .

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

[42]  G. P. Rangaiah,et al.  A Study of Equation-Solving and Gibbs Free Energy Minimization Methods for Phase Equilibrium Calculations , 2002 .

[43]  Michael A. Saunders,et al.  Inertia-Controlling Methods for General Quadratic Programming , 1991, SIAM Rev..

[44]  A. Lucia,et al.  MULTIPHASE EQUILIBRIUM FLASH CALCULATIONS , 2000 .

[45]  Ignacio E. Grossmann,et al.  An equation-oriented approach for handling thermodynamics based on cubic equation of state in process optimization , 2010, Comput. Chem. Eng..

[46]  George Jackson,et al.  THE THERMODYNAMICS OF MIXTURES AND THE CORRESPONDING MIXING RULES IN THE SAFT-VR APPROACH FOR POTENTIALS OF VARIABLE RANGE , 1998 .

[47]  S. Gómez,et al.  Phase stability analysis with cubic equations of state by using a global optimization method , 2002 .

[48]  A. V. Levy,et al.  The Tunneling Algorithm for the Global Minimization of Functions , 1985 .

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

[50]  Christodoulos A. Floudas,et al.  Global optimization for the phase stability problem , 1995 .

[51]  Kurt Binder,et al.  Artificial multiple criticality and phase equilibria: an investigation of the PC-SAFT approach. , 2005, Physical chemistry chemical physics : PCCP.

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

[53]  Claire S. Adjiman,et al.  Fluid phase stability and equilibrium calculations in binary mixtures: Part II: Application to single-point calculations and the construction of phase diagrams , 2009 .

[54]  Dimitrios P. Tassios,et al.  An Equation of State for Associating Fluids , 1996 .

[55]  Claire S. Adjiman,et al.  Understanding the fluid phase behaviour of crude oil: Asphaltene precipitation , 2011 .

[56]  Saeed Khaleghi Rahimian,et al.  Global solution approaches in equilibrium and stability analysis using homotopy continuation in the complex domain , 2008, Comput. Chem. Eng..

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

[58]  S. Gómez,et al.  VOLUME-BASED THERMODYNAMICS GLOBAL PHASE STABILITY ANALYSIS , 2006 .

[59]  Gade Pandu Rangaiah,et al.  Implementation and evaluation of random tunneling algorithm for chemical engineering applications , 2006, Comput. Chem. Eng..

[60]  Frances E. Pereira,et al.  Transferable SAFT-VR models for the calculation of the fluid phase equilibria in reactive mixtures of carbon dioxide, water, and n-alkylamines in the context of carbon capture. , 2011, The journal of physical chemistry. B.

[61]  A. S. Cullick,et al.  New strategy for phase equilibrium and critical point calculations by thermodynamic energy analysis. Part I. Stability analysis and flash , 1991 .

[62]  J. E. Glynn,et al.  Numerical Recipes: The Art of Scientific Computing , 1989 .

[63]  G. P. Rangaiah.  Evaluation of genetic algorithms and simulated annealing for phase equilibrium and stability problems , 2001 .

[64]  K. I. M. McKinnon,et al.  A Generic Global Optimization Algorithm for the Chemical and Phase Equilibrium Problem , 1998, J. Glob. Optim..

[65]  谢安俊,et al.  模拟软件—ASPEN PLUS , 1995 .

[66]  Constantinos C. Pantelides,et al.  Efficient Solution of the Association Term Equations in the Statistical Associating Fluid Theory Equation of State , 2006 .

[67]  Imran Rahman,et al.  Evaluation of repulsive particle swarm method for phase equilibrium and phase stability problems , 2009 .

[68]  Huanquan Pan,et al.  Complex Multiphase Equilibrium Calculations by Direct Minimization of Gibbs Free Energy by Use of Simulated Annealing , 1998 .

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

[70]  Michael L. Michelsen,et al.  Instability of Successive Substitution , 1995 .

[71]  Paul M. Mathias,et al.  Computational aspects of equations of state: Fact and fiction , 1986 .

[72]  L. E. Baker,et al.  Gibbs energy analysis of phase equilibria , 1982 .