An algebraic method that includes Gibbs minimization for performing phase equilibrium calculations for any number of components or phases

The most widely used technique for performing phase equilibria calculations is the K-value method (equality of chemical potentials). This paper proposes a more efficient algorithm to achieve the results that includes Gibbs minimization when we know the number of phases. Using the orthogonal derivatives, the tangent plane equation and mass balances, it is possible to reduce the Gibbs minimization procedure to the task of finding the solution of a system of non-linear equations. Such an operation is easier and faster than finding tangents or areas, and appears to converge as fast as the K-value method. Examples illustrate application of the new technique to two and three phases in equilibrium for binary and ternary mixtures.

[1]  Sandro Ridella,et al.  Minimizing multimodal functions of continuous variables with the “simulated annealing” algorithmCorrigenda for this article is available here , 1987, TOMS.

[2]  Equal Area Rule Methods for Ternary Systems , 1995 .

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

[4]  P. T. Eubank,et al.  Equal area rule and algorithm for determining phase compositions , 1995 .

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

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

[7]  P. T. Eubank,et al.  Maximum partial area rule for phase equilibrium calculations , 1996 .

[8]  Solution of the multiphase equilibrium problem for pure component, binary and ternary systems using the area method , 1996 .

[9]  B. Cairns,et al.  Multicomponent three-phase azeotropic distillation. 2. Phase-stability and phase-splitting algorithms , 1990 .

[10]  Calculation of multi-phase equilibria using the equal area rule with application to mixtures , 1996 .

[11]  R. Heidemann.  Three‐phase equilibria using equations of state , 1974 .

[12]  Wallace B. Whiting,et al.  Area method for prediction of fluid-phase equilibria , 1992 .

[13]  B. M. Fulk.  MATH , 1992 .

[14]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[15]  R. P. Stateva,et al.  A Rigorous Mathematical Proof of the Area Method for Phase Stability , 1998 .

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

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

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

[19]  W. Seider,et al.  Computation of phase and chemical equilibrium: Part I. Local and constrained minima in Gibbs free energy , 1979 .

[20]  P. T. Eubank,et al.  Areas of J. C. Maxwell for Pure-Component Fluid-Phase Equilibria from Equations of State , 1996 .

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

[22]  Yushan Zhu,et al.  A reliable prediction of the global phase stability for liquid-liquid equilibrium through the simulated annealing algorithm: Application to NRTL and UNIQUAC equations , 1999 .