Extension of modified RAND to multiphase flash specifications based on state functions other than (T,P)

Abstract The recently proposed modified RAND formulation is extended from isothermal multiphase flash to several other state function based flash specifications. The obtained general formulation is applicable to chemical equilibrium although this study is focused on flash with only phase equilibrium. It is demonstrated that a common symmetric Jacobian matrix can be formulated for all of these flash specifications. Newton iteration with the common Jacobian is used to converge for the majority of cases and a Q-function maximisation with nested isothermal flash in the inner loop is used for the non-convergent exceptions. For isothermal flash with modified RAND, it can happen in rare occasions that the modified RAND step is ascending in the Gibbs energy. A correction of the step is proposed for such cases to obtain a descent direction without violating the condition used in the derivation of modified RAND. A two-phase example is used to demonstrate that the described method is suitable for ( H , P ) , ( T , V ) , ( S , V ) and ( U , V ) flash specifications and a four-phase case is examined in more detail for the difficult ( U , V ) case. Two- and three-phase examples close to critical regions are used to demonstrate the effectiveness of the correction procedure for the modified RAND step and to show that satisfactory rates of convergence are obtained.

[1]  New Formulations for Isothermal Multiphase Flash , 2017 .

[2]  V. N. Huff,et al.  General method and thermodynamic tables for computation of equilibrium composition and temperature of chemical reactions , 1951 .

[3]  R. Fletcher.  Practical Methods of Optimization , 1988 .

[4]  Optimization Based Isenthalpic Flash for Thermal Reservoir Simulations , 2017 .

[5]  M. Michelsen,et al.  Calculation of tri-critical points , 1988 .

[6]  Selmer M. Johnson,et al.  Chemical Equilibrium in Complex Mixtures , 1958 .

[7]  Michael Locht Michelsen,et al.  Thermodynamic Modelling: Fundamentals and Computational Aspects , 2004 .

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

[9]  S. Brinkley,et al.  Calculation of the Equilibrium Composition of Systems of Many Constituents , 1947 .

[10]  Dan Vladimir Nichita,et al.  Robust and efficient Trust-Region based stability analysis and multiphase flash calculations , 2014 .

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

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

[13]  Jiří Mikyška,et al.  Phase stability testing and phase equilibrium calculation at specified internal energy, volume, and moles , 2017 .

[14]  William R. Smith,et al.  Chemical Reaction Equilibrium Analysis: Theory and Algorithms , 1982 .

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

[16]  H. Greiner.  AN EFFICIENT IMPLEMENTATION OF NEWTON'S METHOD FOR COMPLEX NONIDEAL CHEMICAL EQUILIBRIA , 1991 .