Computation of liquid-vapor critical points for multi-component mixtures

Abstract

Liquid-vapor critical point for a thermodynamic mixture is an important parameter for the design of mixture transportation and processing. The classical critical points model for the multicomponent mixture is based on the tangent plane distance in terms of the Helmholtz energy. In this work, we have applied the damped Newton-Raphson scheme to solve the model using the critical temperature, critical specific volume and change in mole fraction of each component as independent variables. According to numerical experiments, the convergence criteria, initial guesses of variables and damped factor formulation were determined. The method was tested on twenty five mixtures containing two to eleven components using the SRK and PR equation of state. Comparisons between the calculated critical parameters and experimental values illustrated the overall accuracy of the solution method is satisfactory. This method can solve all the variables simultaneously using the whole iteration, which is more reliable than the nested two-singe variable loops method proposed by Heidemann and Khalil. Thus it can be used for practical purposes of locating the liquid-vapor critical points of multi-component mixtures.

Computation of Critical Points of Mixtures Using Particle Swarm Optimization with Low-Discrepancy Sequences

“…More recently, Henderson et al (2010) have tested four different versions of the differential evolution algorithm, including a variant for calculating more than one critical point. Bonilla-Petriciolet (2006), and Jia et al (2012), also attacked the problem of the calculation of critical points using the Helmholtz free energy, but from the formulation proposed by Heidemann and Khalil (1980). Sánches-Mares and BonillaPetriciolet (2006) used the SA algorithm, and Jia et al (2012) employed a damped 4 Newton method.…”

Reliable mixture critical point computation using polynomial homotopy continuation

“…Such limitations are true of all nonlinear solvers that are local in nature. Nevertheless, a significant amount of effort has been dedicated toward improving the convergence properties and overall efficiency of local solver implementations. These techniques are typically well suited for large‐scale reservoir and multiphase flow simulations where critical point calculations may be performed or more times each time‐step.…”