A crossover SAFT-VR equation of state for pure fluids: preliminary results for light hydrocarbons

SAFT is perhaps the most versatile, fundamentally, based engineering equation of state in use today. However, in common with all analytic equations of state, SAFT exhibits classical behavior in the critical region rather than the non-analytical, singular behavior seen in real fluids. Recently, so-called crossover equations of state have been developed which solve this shortcoming by incorporating the scaling laws valid asymptotically close to the critical point while reducing to the original classical equation of state far from the critical point. We have combined the SAFT-VR equation of state with an analytical crossover technique to obtain the SAFT-VRX equation of state. The SAFT-VRX approach combines the accurate low temperature behavior of SAFT-VR with a precise representation of the critical region. Preliminary results are presented for hydrocarbon systems which illustrate the accuracy of the SAFT-VRX approach over the entire fluid phase region.     

Liquid mixtures of xenon with fluorinated species: xenon + sulfur hexafluoride

The vapor-liquid equilibrium of binary mixtures of xenon + SF6 has been measured at nine temperatures from 235.34 to 295.79 K and pressures up to 6.5 MPa. The mixture critical line is found to be continuous between the critical points of the pure components, and hence, the system can be classified as type I phase behavior in the scheme of van Konynenburg and Scott. The excess Gibbs free energies have been calculated, and the experimental results have been interpreted using the statistical associating fluid theory for potentials of variable range (SAFT-VR). Additionally, the SAFT-VR equation has been used to model other systems involving SF6 and alkanes, illustrating the predictability of the approach and further demonstrating the transferability of parameters between binary mixtures involving alkanes and xenon.     

Modeling the phase behavior of H2S+n-alkane binary mixtures using the SAFT-VR+D approach

A statistical associating fluid theory for potential of variable range has been recently developed to model dipolar fluids (SAFT-VR+D) [Zhao and McCabe, J. Chem. Phys. 2006, 125, 104504]. The SAFT-VR+D equation explicitly accounts for dipolar interactions and their effect on the thermodynamics and structure of a fluid by using the generalized mean spherical approximation (GMSA) to describe a reference fluid of dipolar square-well segments. In this work, we apply the SAFT-VR+D approach to real mixtures of dipolar fluids. In particular, we examine the high-pressure phase diagram of hydrogen sulfide+n-alkane binary mixtures. Hydrogen sulfide is modeled as an associating spherical molecule with four off-center sites to mimic hydrogen bonding and an embedded dipole moment (micro) to describe the polarity of H2S. The n-alkane molecules are modeled as spherical segments tangentially bonded together to form chains of length m, as in the original SAFT-VR approach. By using simple Lorentz-Berthelot combining rules, the theoretical predictions from the SAFT-VR+D equation are found to be in excellent overall agreement with experimental data. In particular, the theory is able to accurately describe the different types of phase behavior observed for these mixtures as the molecular weight of the alkane is varied: type III phase behavior, according to the scheme of classification by Scott and Konynenburg, for the H2S+methane system, type IIA (with the presence of azeotropy) for the H2S+ethane and+propane mixtures; and type I phase behavior for mixtures of H2S and longer n-alkanes up to n-decane. The theory is also able to predict in a qualitative manner the solubility of hydrogen sulfide in heavy n-alkanes.     

A modified clausius equation of state for calculation of multicomponent refrigerant vapor-liquid equilibria

Gow, A.S., 1993. A modified Clausius equation of state for calculation of multicomponent refrigerant vapor-liquid equilibria. Fluid Phase Equilibria, 90: 219-249.

A modified Clausius equation of state with a single temperature dependent energy-volume parameter a(T) in the attractive term was designed to describe the vapor pressure vs. temperature relationship of 39 pure refrigerant fluids including elementary cryogenic materials (e.g. He, Ar, N₂, CO₂, CH₄, etc.), chlorofluorocarbons (CFCs), hydrofluorocarbons (HFCs), hydrochlorofluorocarbons (HCFCs), fluorocarbons (FCs), and various other simple cryogenic compounds. The equation developed represents the vapor-liquid coexistence dome, and the superheated vapor compressibility factor and enthalpy for pure refrigerants.

The vapor-liquid equilibrium for refrigerant mixtures is calculated using a “phi-phi” method with “one fluid” van der Waals mixing and combining rules for the equation of state parameters a_M(T), b_M and c_M. A single interaction constant k₁₂ is used to describe non-ideal behavior of each binary. The binary interaction constant, which is a strong function of temperature, and the sign of which signifies the type of deviations from Raoult's law, is obtained by correlating experimental bubble point data for isothermal binary refrigerant mixtures. The proposed equation of state generally describes binary P-x,y data more accurately the higher the temperature for a given system. The method presented is extended to predict vapor-liquid equilibria for the R14-R23-R13 ternary system at 198.75 K using binary interaction constants at this temperature for the three binaries involved.     

Global fluid phase equilibria and critical phenomena of selected mixtures using the crossover soft-SAFT equation

We present here the extension of the crossover soft-statistical associating fluid theory (soft-SAFT) equation of state to mixtures, as well as some illustrative applications of the methodology to mixtures of particular scientific and technological interest. The procedure is based on White's work (White, J. A. Fluid Phase Equilib. 1992, 75, 53) from the renormalization group theory, as for the pure fluids, with the isomorphism assumption applied to the mixtures. The equation is applied to three groups of mixtures: selected mixtures of n-alkanes, the CO2/n-alkane homologous series, and the CO2/1-alkanol homologous series. The crossover equation is first applied to the pure components of the mixtures, CO2 and the 1-alkanol family, while an available correlation is used for the molecular parameters of the n-alkane series (Llovell et al. J. Chem. Phys 2004, 121, 10715). A set of transferable molecular parameters is provided for the 1-alkanols series; these are accurate for the whole range of thermodynamic conditions. The crossover soft-SAFT equation is able to accurately describe these compounds near to and far from the critical point. The theory is then used to represent the phase behavior and the critical phenomena of the selected mixtures. We use binary interaction parameters xi and eta for dissimilar mixtures. These parameters are fitted at some particular conditions (one subcritical temperature or binary critical data) and used to predict the behavior of the mixture at different conditions (other subcritical conditions and/or critical conditions). The equation is able to capture the continuous change in the critical behavior of the CO2/n-alkane and the CO2/1-alkanol homologous series as the chain length of the second compound increases. Excellent agreement with experimental data is obtained, even in the most nonideal cases. The new equation is proved to be a powerful tool to study the global phase behavior of complex systems, as well as other thermodynamic properties of very challenging mixtures.     

Solid-liquid equilibrium using the SAFT-VR equation of state: Solubility of naphthalene and acetic acid in binary mixtures and calculation of phase diagrams

A solid-liquid equilibrium (SLE) thermodynamic model based on the SAFT-VR equation of state (EOS) is presented. The model allows for the calculation of solid-liquid phase equilibria in binary mixtures at atmospheric pressure. The fluid (liquid) phase is treated with the SAFT-VR approach, where molecules are modelled as associating chains of tangentially bonded spherical segments interacting via square-well potentials of variable range. The equilibrium between the liquid and solid phase is treated following a standard thermodynamic method that requires the experimental temperature and enthalpy of fusion of the solute. The model is used to calculate the solubilities of naphthalene and acetic acid in common associating and non-associating organic solvents and to determine the solid-liquid phase behaviour of binary mixtures with simple eutectics. The SAFT-VR pure component model parameters are determined by comparison to experimental vapour pressure and saturated liquid density data with the choice of association models according to the nature of the molecule; in addition, an unlike adjustable parameter (k_ij) is used to model the solutions. The solubility data of naphthalene and acetic acid in both associating and non-associating solvents are reproduced essentially within the accuracy of the experimental measurements. The phase boundaries and the position of the eutectic points in the binary mixtures considered are, in most cases, reproduced with the accuracy commensurate with the industrial applications. Overall, the results presented show that the SAFT-VR EOS can be used with confidence for the prediction of the SLE of binary systems at atmospheric pressure.     

Augmented van der Waals equations of state: SAFT-VR versus Yukawa based van der Waals equation

Performance of the SAFT-VR equation of state developed for the hard sphere based simple fluids, namely the square-well, Sutherland and Yukawa fluids, is examined by comparing its results with simulation data and an augmented van der Waals (vdW) equation based on a Yukawa (Y) reference. Its shown that both for the equilibrium vapor-liquid data and data along selected isotherms in the liquid and supercritical fluid phases the vdW(Y) equation provides better results, particularly when going to lower temperatures.     

A semi-empirical equation of state applied to vapor-liquid equilibria calculations

A five-parameter equation of state is proposed to calculate the vapor-liquid equilibria of compounds in binary and multicomponent mixtures. This equation is closely related to a previous equation of state proposed by the author, the main modification being in the entropic term where the parameter m assumes a constant value for all compounds. It is shown that the van der Waals conditions at the critical point and the Morbidelli-Carra' algorithm enable the calculation of three other constants. Rules are given to calculate the remaining constant K which pertains to the enthalpic term. The proposed method only requires knowledge of the critical constants and of the normal boiling temperature as input parameters. A wide application of the new equation to both polar and non-polar binary systems indicates the following: the proposed method is predictive for ideal or nearly ideal mixtures; the correlation of mixtures of hydrocarbons having very different molar volumes can be obtained by optimizing only the binary interaction parameter linked to the enthalpic term; the new equation also correlates well with strongly non-ideal systems which exhibit a miscibility gap; the prediction of multicomponent vapor-liquid equilibria from the binary data alone is also reliable for both polar and non-polar mixtures.     

Pressure dependence of the vapor-liquid-liquid phase behavior in ternary mixtures consisting of n-alkanes, n-perfluoroalkanes, and carbon dioxide

Expansion of an organic solvent by an inert gas can be used to tune the solvent's liquid density, solubility strength, and transport properties. In particular, gas expansion can be used to induce miscibility at low temperatures for solvent combinations that are biphasic at standard pressure. Configurational-bias Monte Carlo simulations in the Gibbs ensemble were carried out to investigate the vapor-liquid-liquid equilibria and microscopic structures for two ternary systems: n-decane/n-perfluorohexane/CO2 and n-hexane/n-perfluorodecane/CO2. These simulations employed the united-atom version of the transferable potential for phase equilibria (TraPPE-UA) force field. Initial simulations for binary mixtures of n-alkanes and n-perfluoroalkanes showed that special mixing parameters are required for the unlike interactions of CHx and CFy pseudoatoms to yield satisfactory results. The calculated upper critical solution pressures for the ternary mixtures at a temperature of 298 K are in excellent agreement with the available experimental data and predictions using the SAFT-VR (statistical associating fluid theory of variable range) equation of state. The simulations yield asymmetric compositions for the coexisting liquid phases and different degrees of microheterogeneity as measured by local mole fraction enhancements.     

Application of a new simplified SAFT to VLE study of associating and non-associating fluids

A new equation of state based on the Statistical Associating Fluid Theory (SAFT) is presented to study the phase behavior of associating and non-associating fluids. In the new equation of state, the hard sphere contribution to compressibility factor of the simplified version of the SAFT (SSAFT) is replaced with that proposed by Ghotbi and Vera. The Ghotbi–Vera SSAFT (GV-SSAFT) was also extended to study the phase behavior of associating and non-associating mixtures. The GV-SSAFT like the SSAFT equation of state has three adjustable segment parameters for non-associating fluids and five parameters for associating fluids. The experimental data of liquid densities and vapor pressures for pure fluids studied in this work were used to obtain the best values for the parameters of the GV-SSAFT. The results obtained from the GV-SSAFT for liquid densities and vapor pressures of pure associating and non-associating fluids were compared with those obtained from the SSAFT equation of state. The results showed that the GV-SSAFT similar to the SSAFT can accurately correlate the experimental data of liquid density and vapor pressure for systems studied. On the other hand the results obtained from two SAFT-based equations of state are almost identical. In order to show capability of the GV-SSAFT and SSAFT equations of state, they were used to directly calculate heat of vaporization for a number of pure associating and non-associating fluids. Slightly better results for heat of vaporization comparing to the experimental data were obtained from the GV-SSAFT EOS than those obtained from the SSAFT. The GV-SSAFT was also used to study the VLE phase behavior for a number of binary associating and non-associating mixtures. The results also showed that the GV-SSAFT can be successfully used to study the phase behavior of mixtures studied in this work.     

The phase behavior of two-dimensional symmetrical mixtures

Using Monte Carlo simulation methods in the grand canonical and semigrand canonical ensembles, we study the phase behavior of two-dimensional symmetrical binary mixtures of Lennard-Jones particles. We discuss the interplay between the demixing transition in a liquid and the freezing in detail. Phase diagrams for several systems characterized by different parameters describing interactions in the system are presented. It is explicitly demonstrated that different scenarios involving demixing and freezing transitions, described in our earlier paper [A. Patrykiejew and S. Soko?owski, Phys. Rev. E, 81, 012501 (2010)], are possible. In one class of systems, the λ-line representing a continuous demixing transition in a liquid phase starts at the liquid side of either the vapor-liquid or liquid-solid coexistence. The second class involves the systems in which the λ-line begins at the liquid side of the vapor-liquid coexistence, in the lower critical end point, and then terminates at the liquid side of the liquid-solid coexistence, in the upper critical end point. It is also shown that in such systems the solid phase may undergo a demixing transition at the temperature above the upper critical end point.     

Vapor-liquid and vapor-liquid-liquid equilibria of carbon dioxide/n-perfluoroalkane/n-alkane ternary mixtures

Perfluoroalkanes have numerous applications (e.g., in the medical field and the chemical industry), and their high affinity for carbon dioxide makes them attractive as surfactants and cosolvents. Although research in this area has grown in the past few years, very little phase-equilibrium data is available in the open literature for these systems. In this work, we present, for the first time, predictions of vapor-liquid and vapor-liquid-liquid equilibria of binary and ternary systems of carbon dioxide/n-perfluoroalkane/n-alkane. Our results are based on the SAFT-VR EOS (statistical associating fluid theory of variable range, equation of state), and we study the influence of temperature, pressure, composition, and chain length on the phase diagram. The predicted phase diagrams are based on temperature-independent binary interaction parameters, and no ternary parameters are introduced. Comparisons to the available experimental and molecular simulation data show that the predicted diagrams should provide a good representation of the phase equilibria.     

Extending the GC-SAFT-VR approach to associating functional groups: Alcohols, aldehydes, amines and carboxylic acids

The statistical associating fluid theory is a widely used molecular-based equation of state that has been successfully applied to study a broad range of fluid systems. It provides a framework in which the effects of molecular shape and interactions on the thermodynamics and phase behavior of fluids can be separated and quantified. In the original approach, molecules were modeled as chains composed of identical segments; the heterogeneity of molecules in terms of structure and functional groups was described implicitly through effective parameters. To overcome this limitation, in recent works [Peng et al. Fluid Phase Equilib. 227(2), 131 (2009); Ind. Eng. Chem. Res. 49(3), 1378 (2010)] the GC-SAFT-VR approach has been developed to extend the theory to model chains composed of segments of different size and/or energy of interaction and enable the development of a group-contribution approach within the SAFT-VR framework in which molecular heterogeneity and connectivity is explicitly accounted for. The parameters for several key functional groups (CH₃, CH₂, CH, CH₂CH, CO, C₆H₅, esters, ethers, cis-alkenes and trans-alkenes groups) were determined by fitting to experimental vapor pressure and saturated liquid density data for a number of small molecules containing the functional groups of interest and transferability of the parameters tested by comparing the theoretical predictions with experimental data for pure fluids not included in the fitting process and binary mixtures of both simple fluids and the VLE and LLE of small molecules in polymer systems. In this work, we further extend the applicability of the GC-SAFT-VR approach through the study of the vapor-liquid phase behavior of associating systems, such as linear and branched alcohols, primary and secondary amines, aldehydes, and carboxylic acids, and their mixtures. In the study of these new molecules several new functional groups (OH (linear and branched), HCO, NH₂, NH and COOH) are defined and their molecular parameters characterized. The transferability of the parameters is again tested by comparing the theoretical predictions with experimental data for pure fluids and binary mixtures not included in the fitting process. The GC-SAFT-VR approach is found to predict the phase behavior of the systems studied in most cases in good agreement with experimental data and accurately captures the effects of changes in structure and molecular composition on phase behavior.     

Equation of state for square-well chain molecules with variable range II. Extension to mixtures

An equation of state (EOS) developed in our previous work for square-well chain molecules with variable range is further extended to the mixtures of non-associating fluids. The volumetric properties of binary mixtures for small molecules as well as polymer blends can well be predicted without using adjustable parameter. With one temperature-independent binary interaction parameter, satisfactory correlations for experimental vapor–liquid equilibria (VLE) data of binary normal fluid mixtures at low and elevated pressures are obtained. In addition, VLE of n-alkane mixtures and nitrogen + n-alkane mixtures at high pressures are well predicted using this EOS. The phase behavior calculations on polymer mixture solutions are also investigated using one-fluid mixing rule. The equilibrium pressure and solubility of gas in polymer are evaluated with a single adjustable parameter and good results are obtained. The calculated results for gas + polymer systems are compared with those from other equations of state.     

Vapor–liquid equilibrium in the n-butane + methanol system,measurement and modeling from 323.2 to 443.2 K

New experimental vapor–liquid equilibrium (VLE) data for the n-butane + methanol binary system are reported over a wide temperature range from 323.2 to 443.2 K and pressures up to 5.4 MPa. A static–analytic apparatus, taking advantage of two pneumatic capillary samplers, was used. The phase equilibrium data generated in this work are in relatively good agreement with previous data reported in the literature. Three different thermodynamic models have been used to represent the new experimental data. The first model is the cubic-based Peng–Robinson equation of state (EoS) combined with the Wong–Sandler mixing rules. The two other models are the non-cubic SAFT-VR and PC-SAFT equations of state. Temperature-dependent binary interaction parameters have been adjusted to the new data. The three models accurately represent the new experimental data, but deviations are seen to increase at low temperature. A similar evolution of the binary parameters with respect to temperature is observed for the three models. In particular a discontinuity is observed for the k_ij values at temperatures close to the critical point of butane, indicating the effects of fluctuations on the phase equilibria close to critical points.