A simple non-classical equation of state for fluids and fluid mixtures

A method for improving the behavior of classical equations of state (EOS) in the critical region, originally proposed by Fox [J.R. Fox, Fluid Phase Equilibria 14 (1983) 45–53], has been modified in this work for the Patel–Teja (PT) EOS [N.C. Patel, A.S. Teja, Chem. Eng. Sci. 37, 463–473]. The application of the new equation (NPT) for predicting PVT and vapor pressure behavior of pure substances, as well as vapor–liquid equilibrium behavior of binary mixtures, is demonstrated. The NPT equation is simple to use and requires the same input information as the original PT equation. However, it reproduces the correct PVT behavior in the critical region. Limitations of both the PT and NPT equations in calculating the isochoric heat capacity are discussed.     

Prediction of phase equilibrium in mixtures containing ionic liquids using UNIFAC

The UNIFAC model is extended to mixtures of ionic liquids consisting of the imidazolium cation and the hexafluorophosphate anion with alkanes, cycloalkanes, alcohols and water. Two new main groups, the imidazolium and the hexafluorophosphate groups, are introduced in UNIFAC. The required group interaction parameters between these groups and the existing UNIFAC main groups, CH₂, OH and H₂O, are determined by fitting binary liquid–liquid equilibrium and infinite dilution activity coefficient experimental data. The predictive capability of the extended UNIFAC model is examined against experimental data for vapour–liquid equilibrium, liquid–liquid equilibrium and activity coefficients at infinite dilution of binary and ternary systems containing 1-alkyl-3-alkyl′-imidazolium hexafluorophosphate ionic liquids, alkanes, cycloalkanes, alcohols and water. The results indicate that UNIFAC is a reliable model for phase equilibrium predictions in mixtures containing this type of ionic liquids.     

An analytical perturbed equation of state for hard chain fluids: Application to n-alkanes and n-perfluoroalkanes

A completely analytical equation of state for pure hard chain fluids, derived on the basis of perturbation theory and reported in our previous work, is applied for the calculation of pVT properties and the prediction of vapour–liquid equilibria of n-alkanes and n-perfluoroalkanes. The molecules are treated as a chain formed from freely joined spheres which interact via an extended site-site square-well potential. The molecular parameters of compounds are obtained from the experimental compressibility factor data above the critical temperature. These parameters are capable of relatively satisfactory prediction of the vapour–liquid equilibrium coexistence curves of compounds. Linear relationships have been found between the potential parameters of fluids and their molecular weight, which make it possible to predict the pVT data and vapour–liquid phase equilibria of heavier compounds.     

Interfacial properties of water + alcohol mixtures

Mixtures of water with alcohol are important in numerous engineering applications. Caused by the polarity of water and alcohol self-association of water and alcohol as well cross-association between water and alcohol appear in such complex mixtures. These features show significant impact on physical and chemical properties, especially phase equilibrium behaviour and hence interfacial properties. The Cahn–Hilliard theory was combined with original statistical associated fluid theory equation of states (SAFT EOS) in order to describe both the phase behaviour and interfacial properties with respect of association. The paper focuses on theoretical investigations of surface tension, density profiles, surface thickness in vapour–liquid or vapour–liquid–liquid equilibrium of mixtures of water with ethanol or 1-butanol. Results of vapour–liquid equilibrium surface tension calculations were compared with experimental data taken from the literature.     

Low pressure vapor–liquid equilibrium in ethanol + congener mixtures using the Wong–Sandler mixing rule

Phase equilibrium in binary ethanol mixtures found in alcoholic beverage production has been analyzed using a cubic equation of state (EoS) and suitable mixing and combining rules. The main objective of the study is the accurate modeling of the congener concentration in the vapor phase (substances different from ethanol), considered to be an important enological parameter in the alcohol industry. The Peng–Robinson (PR) equation of state has been used and the Wong–Sandler (WS) mixing rules, that include a model for the excess Gibbs free energy, have been incorporated into the equation of state constants. In the Wong–Sandler mixing rules the van Laar (VL) model for the excess Gibbs energy has been used. This combination of equations of state, mixing rules and combining rules are commonly applied to high pressure phase equilibrium and have not yet been treated in a systematic way to complex low pressure ethanol mixtures as done in this work. Nine binary ethanol + congener mixtures have been considered for analysis. Comparison with available literature data is done and the accuracy of the calculations is discussed, concluding that the model used is accurate enough for engineering applications.     

A crossover cubic equation of state near to and far from the critical region

In this research, we use the Patel–Teja (PT) cubic equation of state [N.C. Patel, A.S. Teja, Chem. Eng. Sci. 37 (1982) 463–473.] and develop a crossover cubic model near to and far from the critical region, which incorporates the scaling laws asymptotically close to the critical point and it transformed into original classical cubic equations of state far away from the critical point. This equation of state is used to calculate thermodynamic properties of pure systems (carbon dioxide, normal alkanes from methane to heptane). We show that, over a wide range of states, the equation of state yields the saturated vapour pressure data and the saturated density data with a much better accuracy than the original PT equation of state.     

Empirical correction to the Peng–Robinson equation of state for the saturated region

This work presents an empirical correction to improve the Peng–Robinson equation of state (PR EOS) for representing the densities of pure liquids and liquid mixtures in the saturated region using the volume translation method. A temperature-dependent volume correction is employed to improve the original PR EOS so that it can match the true critical point of pure fluids. The volume correction is generalized as a function of the critical parameters and the reduced temperature. The volume translation PR (VTPR) EOS with the generalized volume correction accurately represents the saturated liquid densities for different polar and non-polar fluids, including alkanes, cycloparaffins, halogenated hydrocarbons, olefins, cyclic olefins, aromatics and inorganic molecules. The average relative deviations for 91 pure compounds was 1.37%. The generalized VTPR EOS was also used to predict the saturated liquid density of 53 binary mixtures with a relative deviation of 0.98%. The generalized VTPR EOS can also be extended to other materials. The accuracy of the generalized VTPR EOS compares well with other methods and equations of state.     

High-pressure phase equilibria in the carbon dioxide + pyrrole system

Liquid–vapor (LV) and liquid–liquid (LL) phase equilibria in the carbon dioxide + pyrrole system were measured at temperatures between 313 K and 333 K, and pressures between 8.4 MPa and 15.1 MPa. The data were used to predict the overall phase behavior of the system using the Patel–Teja equation of state and the Mathias–Klotz–Prausnitz mixing rules with two temperature-independent parameters. The calculations suggest that the carbon dioxide + pyrrole system may exhibit type IV phase behavior according to the classification of Scott and van Konynenburg.     

Estimation of solubility parameter by the modified ER equation of state

The applications of the solubility parameter in chemical, petroleum and polymer engineering industries have been cleared up along the past 50 years. In this article, the Hildebrand solubility parameter of over 250 substances were calculated by the modified ER (Esmaeilzadeh–Roshanfekr) equation of state and some others (the Peng–Robinson, Soave–Redlich–Kwong, Patel–Teja and Schmidt–Wenzel) and compared with the experimental data. Once the less average errors of the mER method predictions were satisfied in subcritical and some supercritical fluids region, a correlation based on this EOS was presented in order to calculate the total HSP (Hansen solubility parameter) of various types of organic components categorized in 13 distinct groups including paraffins, olefins, aromatics, naphthenes, alcohols, aldehydes, ketones, ethers, esters, amines, carboxylic acids and two petroleum sub-fractions (resins and asphaltenes). The optimal values of the model parameters were obtained applying the DE (differential evolution) optimization method. The absolute average deviations of the proposed correlations results from the experimental ones lied between 0.09 and 6%.     

A modified perturbed hard-sphere-chain equation of state: consideration of attractive contribution

We adapted a new perturbation term reported by Bokis et al. for the modified perturbed hard-sphere-chain (PHSC) equation of state. The new model better agrees with the simulation data for illustrating both repulsive and attractive contributions than those of the PHSC equation of state. Three characteristic parameters, r, σ, and , are obtained from thermodynamic properties of pure saturated species and directly used to calculate coexistence curves of given binary systems. Vapor–liquid equilibria and liquid–liquid equilibria of fluids are remarkably well predicted by the model even in the case of polymer solution that shows non-ideal phase behavior.     

The phase envelopes of alternative solvents (ionic liquid,CO2) and building blocks of biomass origin (lactic acid,propionic acid)

Solid–liquid, liquid–liquid and vapour–liquid equilibrium measurements for binary and ternary systems containing building blocks of biomass origin such as propionic acid, lactic acid and alternative solvents like carbon dioxide and 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ionic liquid have been carried out at 313.15 K. The binary solid–liquid and liquid–liquid equilibrium measurements were performed at ambient pressure. The vapour–liquid equilibrium was studied in the range of pressure from 3.54 to 12 MPa while ternary systems were examined at 9, 10 and 12 MPa. The samples from the coexisting phases were taken and the compositions of both liquid and vapour phases were determined experimentally. The three-phase system was observed for lactic acid + ionic liquid + CO₂ as well. The achieved results were correlated using the Peng–Robinson equation of state with the Mathias–Klotz–Prausnitz mixing rule. The set of interaction parameters for the employed equations of state and the mixing rule for the investigated systems were obtained.     

Modelling the VLE properties using van der Waals-type equation of state

The vapour–liquid equilibrium (VLE) properties of polar and non-polar fluids have been modelled by the use of two modified van der Waals (vdW)-type equations of state (EOSs). In this article, a revised method is applied to the above-mentioned EOSs to improve the representation of VLE properties of different class of fluids. In this respect, the repulsion parameter b is considered to be temperature dependent and also a temperature-dependent revision factor α(T) is introduced to the liquid fugacity coefficient expression derived from traditional isothermal integration to reproduce the vapour pressure (P^s) of pure liquids. The present method is also extended to represent the VLE properties of binary mixtures containing noble gases, refrigerants and hydrocarbons. This method outperforms the original vdW-type EOSs in predicting the VLE and pressure-volume-temperature (PVT) properties of 22 pure substances and 7 binary mixtures.     

Crossover Peng-Robinson equation of state with introduction of a new expression for the crossover function

In this research, we use the original Peng-Robinson (PR) equation of state (EOS) for pure fluids and develop a crossover cubic equation of state which incorporates the scaling laws asymptotically close to the critical point and it is transformed into the original cubic equation of state far away from the critical point. The modified EOS is transformed to ideal gas EOS in the limit of zero density. A new formulation for the crossover function is introduced in this work. The new crossover function ensures more accurate change from the singular behavior of fluids inside the regular classical behavior outside the critical region. The crossover PR (CPR) EOS is applied to describe thermodynamic properties of pure fluids (normal alkanes from methane to n-hexane, carbon dioxide, hydrogen sulfide and R125). It is shown that over wide ranges of state, the CPR EOS yields the thermodynamic properties of fluids with much more accuracy than the original PR EOS. The CPR EOS is then used for mixtures by introducing mixing rules for the pure component parameters. Higher accuracy is observed in comparison with the classical PR EOS in the mixture critical region.     

A new model and extension of Wong–Sandler mixing rule for prediction of (vapour + liquid) equilibrium of polymer solutions using EOS/GE

The cubic equation of state (CEOS) is a powerful method for calculation of (vapour + liquid) equilibrium (VLE) in polymer solutions. Using CEOS for both the vapour and liquid phases allows one to calculate the non-ideality of polymer solutions based on a single EOS approach. However, the traditional mixing rules are not appropriate to extend the CEOS to non-ideal mixtures such as polymer solutions. Several authors have applied the EOS/G^E approach to predict (vapour + liquid) equilibria in polymer solutions, however, incorporating an appropriate excess Gibbs free energy for the new mixing rule is a major step. In this research, the NRTL-NRF model was extended in terms of volume fraction of polymer and solvent (instead of mole fraction), then equilibrium calculations were carried out using PRSV EOS and Wong–Sandler mixing rules. Using the adjustable parameters as a function of solution temperature, the NRTL-NRF model can be used as a predictive model. In comparison with NRTL model, the results of the new NRTL-NRF model show better accuracy.     

A new cubic equation of state for prediction of VLE of polymer solutions

A new cubic equation of state (CEOS) is proposed based on temperature–pressure superposition principle. A generic CEOS form, with the Peng–Robinson parameters, is used and a temperature-dependent attractive term a(T)$a(T)$ is developed, allowing an easy calculation of thermodynamic properties and vapor–liquid equilibrium. The new equation is applied to pure polymer and polymer solutions and its results are compared with those of two others equations of state. For polymer solutions, two mixing rules without binary interaction parameters were used. The vapor–liquid equilibrium (VLE) predictions showed good agreement with experimental data as well as pressure–volume–temperature (PVT) behavior of polymer liquids, attesting the appropriate form of the new equation proposed.     

Modeling high-pressure densities at wide temperature range with volume scaling: Cyclohexane + n-hexadecane mixtures

High-pressure density data for cyclohexane + n-hexadecane mixtures at a wide temperature range was modeled with several classical equations of state (EOS) and correlative models. A modification for softening the co-volume and another for a volume scaling of the Peng–Robinson EOS (VS-PR) were proposed. The VS-PR model is able to correlate the pure component experimental data employing only five adjustable parameters, with root-mean-square deviation (RMSD) between calculated and experimental densities essentially within the experimental error. This result is superior to widely used approaches, i.e., a six parameter Tait model and six parameter volume translations (temperature and pressure dependent) for Peng–Robinson and Patel–Teja EOS. The VS-PR model also represents well the isobaric thermal expansion and the isothermal compressibility coefficients of the pure cyclohexane, a small naphthenic substance as well as a long chain n-alkane hydrocarbon, n-hexadecane. When modeling the mixture data, the use of VS-PR model of pure components along with the Redlich–Kister expansion, truncated at the first term, the density was correlated within a RMSD only 60% greater than the experimental error. The proposed model is able to accurately represent all the tested mixture data with a relatively small number of parameters.     

Totally inclusive cubic equation of state with five parameters for pure fluids

A totally inclusive cubic equation of state (cubic EOS) is proposed. Although, its form is fairly simple as compared with the present cubic equations, it can include all of them as special cases. The EOS has five parameters. By fitting the experimental critical isothermal for six typical substances combining the critical conditions, the generalized expressions for the five parameters at critical temperature are established. The temperature coefficients of the five parameters for 43 substances are determined by fitting the experimental data of vapor pressure and saturated liquid density. These coefficients are correlated with the critical compressibility factor and acentric factor to obtain the generalized expressions. The predicted saturated vapor pressure, saturated liquid density, critical isothermal and coexistence curve near the critical point show that the equation gives the best results when compared with the Redlich–Kwong–Soave (RKS) and Peng–Robinson (PR) EOS.