A Study of Equation-Solving and Gibbs Free Energy Minimization Methods for Phase Equilibrium Calculations

• Published in: Chemical engineering research & design Vol. 80; no. 7; pp. 745 - 759
• Main Authors: Teh, Y.S., Rangaiah, G.P.
• Format: Journal Article
• Language: English
• Published: Rugby Elsevier B.V 01.10.2002; Institution of Chemical Engineers
• Subjects: Chemical thermodynamics; Chemistry; equation-solving approach; Exact sciences and technology; General and physical chemistry; General. Theory; Gibbs free energy minimization; liquid–liquid equilibrium; Phase equilibria; vapour–liquid equilibrium; vapour–liquid-liquid equilibrium; vapour–liquid equilibrium; liquid–liquid equilibrium; equation-solving approach; Gibbs free energy minimization; vapour–liquid-liquid equilibrium; Gibbs free energy; Thermodynamic model; Liquid liquid vapor equilibrium; Liquid vapor equilibrium; Multiphase equilibrium; Minimization; Equation resolution; Phase equilibrium; Liquid liquid equilibrium
• ISSN: 0263-8762; 1876-4800; 1744-3563; 1744-3598; 0957-5820
• DOI: 10.1205/026387602320776821

Phase equilibrium calculations are often involved in the design, simulation, and optimization of chemical processes. Reported methods for these calculations are based on either equation-solving or Gibbs free energy minimization approaches. The main objective of this work is to compare selected methods for these two approaches, in terms of reliability to find the correct solution, computational time, and number of K-value evaluations. For this, four equation-solving and three free minimization methods have been selected and applied to commonly encountered vapour–liquid equilibrium (VLE), liquid–liquid equilibrium (LLE), and vapour–liquid–liquid equilibrium (VLLE) examples involving multiple components and popular thermodynamic models. Detailed results show that the equation-solving method based on the Rachford–Rice formulation accompanied by mean value theorem and Wegstein's projection is reliable and efficient for two-phase equilibrium calculations not having local minima. When there are multiple minima and for three-phase equilibrium, the stochastic method, genetic algorithm (GA) followed by modified simplex method of Nelder and Mead (NM) is more reliable and desirable. Generic programs for numerical methods are ineffective for phase equilibrium calculations. These findings are of interest and value to researchers and engineers working on phase equilibrium calculations and/or developing thermodynamic models for phase behaviour.