Application of a global optimization algorithm to phase stability and liquid–liquid equilibrium calculations

Abstract A global optimization algorithm was developed for finding all singular points (minima, maxima and saddles) of an objective function by exploring the natural connectedness that exists between their singular points. The idea of following ridges and valleys using information gathered along the way was significantly enhanced by applying the arc length continuation method. The algorithm was applied to the Gibbs tangent plane stability test for multiphase liquid mixtures. The global optimization algorithm gives an efficient and robust scheme for locating all stationary points of the tangent plane distance function predicted by any thermodynamic model. Since it provides very good initial estimates for the liquid–liquid equilibrium calculations, it became an integral part of a combined phase equilibrium and stability algorithm. The combined algorithm is self-starting and significantly improves reliability and robustness of multiphase equilibrium calculations. It was successfully tested on a variety of problems and is applicable to any component's mixtures and any number of liquid phases. Solutions have been found for the entire phase diagram.

[1]  Gordon S. G. Beveridge,et al.  Optimization: theory and practice , 1970 .

[2]  Michael F. Malone,et al.  Computing all homogeneous and heterogeneous azeotropes in multicomponent mixtures , 1999 .

[4]  Computation of phase equilibrium: Optimization with thermodynamic inconsistency , 1985 .

[5]  P. I. Barton,et al.  Bilevel optimization formulation for parameter estimation in vapor–liquid(–liquid) phase equilibrium problems , 2009 .

[6]  Vassilis Gaganis,et al.  A general framework of model functions for rapid and robust solution of Rachford–Rice type of equations , 2012 .

[7]  Global Stability Analysis and Calculation of Liquid−Liquid Equilibrium in Multicomponent Mixtures† , 1996 .

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

[9]  G. P. Rangaiah,et al.  Evaluation of stochastic global optimization methods for modeling vapor-liquid equilibrium data , 2010 .

[10]  Christodoulos A. Floudas,et al.  Global Optimization and Analysis for the Gibbs Free Energy Function Using the UNIFAC, Wilson, and ASOG Equations , 1995 .

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

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

[13]  Yushan Zhu,et al.  Global stability analysis and phase equilibrium calculations at high pressures using the enhanced simulated annealing algorithm , 2000 .

[14]  A. W. Francis.  Ternary Systems with Three Separate Binodal Curves. , 1956 .

[15]  C. Floudas,et al.  Global optimization for the phase and chemical equilibrium problem: Application to the NRTL equation , 1995 .

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

[17]  Hans Hasse,et al.  Convex envelope method for the determination of fluid phase diagrams , 2012 .

[18]  Rainer Storn,et al.  Differential Evolution – A Simple and Efficient Heuristic for global Optimization over Continuous Spaces , 1997, J. Glob. Optim..

[19]  Abbas Firoozabadi,et al.  Initialization of phase fractions in Rachford-Rice equations for robust and efficient three-phase split calculation , 2012 .

[20]  S. Gómez,et al.  Multiphase equilibria calculation by direct minimization of Gibbs free energy with a global optimization method , 2002 .

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

[22]  Michael F. Malone,et al.  Conceptual design of distillation systems , 2001 .

[23]  Rosendo Monroy-Loperena.  A block-algebra approach to the solution of Σi(ai / λ-ci) equations , 2012, Comput. Chem. Eng..

[24]  C. Floudas,et al.  GLOPEQ: A new computational tool for the phase and chemical equilibrium problem , 1997 .

[25]  M. Stadtherr,et al.  Reliable Phase Stability Analysis for Excess Gibbs Energy Models , 2000 .

[26]  W. Seider,et al.  Homotopy-continuation method for stability analysis in the global minimization of the Gibbs free energy , 1995 .

[27]  Seif-Eddeen K. Fateen,et al.  Cuckoo Search: A new nature-inspired optimization method for phase equilibrium calculations , 2013 .

[28]  W. Arlt,et al.  Liquid-liquid equilibrium data collection , 1979 .

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

[30]  The existence and uniqueness of trajectories joining critical points for differential equations in R3 , 2001 .

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

[32]  G. Cholakov,et al.  A powerful algorithm for liquid–liquid–liquid equilibria predictions and calculations , 2000 .

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

[34]  Jani Kangas,et al.  Modified Bounded Homotopies in the Solving of Phase Stability Problems for Liquid–Liquid Phase-Splitting Calculations , 2011 .

[35]  M. Doherty,et al.  On the dynamics of distillation processes-V: The topology of the boiling temperature surface and its relation to azeotropic distillation , 1984 .

[36]  A. Kabouche,et al.  Interaction parameter estimation in liquid–liquid phase equilibrium modeling using stochastic and hybrid algorithms , 2012 .

[37]  Stanislaw K. Wasylkiewicz,et al.  Synthesis and design of heterogeneous separation systems with recycle streams , 2003 .

[38]  A. Bonilla-Petriciolet.  On the capabilities and limitations of harmony search for parameter estimation in vapor-liquid equilibrium modeling , 2012 .

[39]  Angelo Lucia,et al.  Multivariable terrain methods , 2003 .