Direct Prediction of Cricondentherm and Cricondenbar Coordinates of Natural Gas Mixtures using Cubic Equation of State

摘要

A numerical algorithm is presented for direct calculation of the cricondenbar and cricondentherm coordinates of natural gas mixtures of known composition based on the Michelsen method. In the course of determination of these coordinates, the equilibrium mole fractions at these points are also calculated. In this algorithm, the property of the distance from the free energy surfaces to a tangent plane in equilibrium condition is added to saturation calculation as an additional criterion. An equation of state (EoS) was needed to calculate all required properties. Therefore, the algorithm was tested with Soave-Redlich-Kwong (SRK), Peng-Robinson (PR), and modified Nasrifar-Moshfeghian (MNM) equations of state. For different EoSs, the impact of the binary interaction coefficient (k _ij) was studied. The impact of initial guesses for temperature and pressure was also studied. The convergence speed and the accuracy of the results of this new algorithm were compared with experimental data and the results obtained from other methods and simulation softwares such as Hysys, Aspen Plus, and EzThermo.