Rapid multi-component phase-split calculations using volume functions and reduction methods

Abstract

We present a new family of fast and robust methods for the calculation of the vapor–liquid equilibrium at isobaric-isothermal (PT-flash), isochoric-isothermal (VT-flash), isenthalpic-isobaric (HP-flash), and isoenergetic-isochoric (UV-flash) conditions. The framework is provided by formulating phase-equilibrium conditions for multi-component mixtures in an effectively reduced space based on the molar specific value of the recently introduced volume function derived from the Helmholtz free energy. The proposed algorithmic implementation can fully exploit the optimum quadratic convergence of a Newton method with the analytical Jacobian matrix. This article provides all required exact analytic expressions for the general cubic equation of state. Computational results demonstrate the effectivity and efficiency of the new methods. Compared to conventional methods, the proposed reduced-space iteration leads to a considerable speed-up as well as to improved robustness and better convergence behavior near the spinodal and coexistence curves of multi-component mixtures, where the preconditioning by the reduction method is most effective.