Our paper entitled: The spinodal of single- and multi-component fluids and its role in the development of modern equations of state has been accepted for publication in Fluid Phase Equilibria. The spinodal represents the limit of thermodynamic stability of a homogeneous fluid. In this work, we present a robust methodology to obtain the spinodal of multicomponent fluids described even with the most sophisticated equations of state (EoS) available. We elaborate how information about the spinodal and its uncertainty can contribute both in the development of modern EoS and to estimate their uncertainty in the metastable regions. Inequality constraints are presented that can be exploited in the fitting of modern EoS of single-component fluids to avoid inadmissible pseudo-stable states between the vapor and liquid spinodals.
We find that even cubic EoS violate some of these constraints. With the use of a selection of EoS representative of modern applications, we compare vapor and liquid spinodal curves, superheat and supersaturation limits from classic nucleation theory (CNT), and available experimental data for the superheat limit. Computations are performed with pure species found in natural gas, binary mixtures, as well as a multi-component natural gas mixture in order to demonstrate the scalability of the approach. We demonstrate that there are large inconsistencies in predicted spinodals from a wide range of EoS such as cubic EoS, extended corresponding state EoS, SAFT and multiparameter EoS. We demonstrate that a large part of the metastable domain of the phase diagram is currently unavailable to experiments (the area between the solid red (spinodals) and the green dashed lines (CNT) in the figure below, which shows the pressure of methane as a function of density and temperature). Novel techniques, experimental or with computational simulations, should be developed to characterize the thermodynamic properties in these regions, and to identify the thermodynamic states that define the spinodal. The work was performed in in collaboration with Peder Aursand, Magnus Aa. Gjennestad, Eskil Aursand and Morten Hammer from SINTEF Energy Research.