Development and extension of PSRK/UNIQUAC model to methane and nitrogen gases

Two successful procedures for matching the equation of state (EOS) and the excess Gibbs energy model at zero pressure belong to Michelsen [M.L. Michelsen, Fluid Phase Equilibria, 60 (1990) 213] and Holderbaum-Gmehling [T. Holderbaum, J. Gmehling, Fluid Phase Equilibria, 70 (1991) 251]. These procedures lead to volume independent mixing rules for the mixture parameter a. Michelsen proposed the MHV1 and MHV2 mixing rules and Holderbaum-Gmehling proposed the PSRK mixing rule. Keshtkar et al. [A. Keshtkar, F. Jalali, M. Moshfegian, Fluid Phase Equilibria, 140 (1997) 107] demonstrated the utility of these mixing rules together with UNIFAC and UNIQUAC models for the prediction of high-pressure vapor-liquid equilibrium (VLE) of CO₂-binary systems. Also, they showed that the PSRK/UNIQUAC model gives better VLE results than other models. In this paper, we develop and extend the above model to methane and nitrogen gases. For this purpose, it is required that the missing interaction parameters of model are fitted with experimental VLE data. The VLE calculations show that good agreement between calculated results and experimental data can be produced by using the obtained parameters.     

Modelling (vapour + liquid) and (vapour + liquid + liquid) equilibria of {water (H2O) + methanol (MeOH) + dimethyl ether (DME) + carbon dioxide (CO2)} quaternary system using the Peng–Robinson EoS with Wong–Sandler mixing rule

The (vapour + liquid) equilibria (VLE) and (vapour + liquid + liquid) equilibria (VLLE) binary data from literature were correlated using the Peng–Robinson (PR) equation of state (EoS) with the Wong–Sandler mixing rule (WS). Two group contribution activity models were used in the PRWS: UNIFAC–PSRK and UNIFAC–Lby. The systems were successfully extrapolated from the binary systems to ternary and quaternary systems. Results indicate that the PRWS–UNIFAC–PSRK generally displays a better performance than the PRWS–UNIFAC–Lby.     

Vapor–liquid equilibria of asymmetrical systems using UNIFAC-NRF group contribution activity coefficient model

The UNIFAC-NRF group contribution activity coefficient model is used for the calculation of vapor–liquid equilibria of binary systems of the heavy alkanes and light gases such as CH₄, C₂H₆, CO₂ and N₂. The linear combination mixing rule, LCVM, of the Huron–Vidal and Michelsen, Chen et al. modification of PSRK and Universal Group Contribution Equation of State of Ahlers and Gmehling are combined with the UNIFAC-NRF group activity coefficient model to correlation of the vapor–liquid equilibrium of both light and heavy hydrocarbons. The results show that the LCVM mixing rule combing with UNIFAC-NRF group contribution model correlate the asymmetric systems better than the LCVM-UNIFAC and the other EOS/G^E models. Also the group contribution model is used for the prediction of the phase envelope of the synthetic fluid with accurate results.     

PSRK: A Group Contribution Equation of State Based on UNIFAC

A group contribution equation of state called PSRK (Predictive Soave-Redlich-Kwong) which is based on the Soave-Redlich-Kwong equation (Soave, 1972) has been developed. It uses the UNIFAC method to calculate the mixture parameter a and includes all already existing UNIFAC parameters. This concept makes use of recent developments by Michelsen (1990b) and has the main advantage, that vapor-uquid-equilibria (VLB) can be predicted for a large number of systems without introducing new model parameters that must be fitted to experimental VLB-data. The PSRK equation of state can be used for VLB-predictions over a much larger temperature and pressure range than the UNIFAC γ--approach and is easily extended to mixtures containing supercritical compounds. Additional PSRK parameters, which allow the calculation of gas/gas and gas/alkane phase equilibria, are given in this paper. In addition to those mixtures covered by UNIFAC, phase equilibrium calculations may also include gases like CH₄ C₂H₆, C₃H₆, c₄H₁₀, CO₂, N₂, H₂ and CO.     

Improper matching of solvation energy components in G-based mixing rules

The physical significance of terms in two excess Gibbs free energy (G^ex)-based mixing rules, the modified Huron–Vidal (MHV1) and Wong–Sandler (WS) mixing rule, are examined through the use of solvation free energy. It is found that these mixing rules are in fact matching the charging contributions of solvation in an equation of state (EOS) to the complete solvation free energy in a liquid activity coefficient model (LM). The cavity contributions in the EOS are canceled as a result of the constant liquid molar volume to molecular volume ratio. The underlying idea of G^ex-based mixing rules that the EOS should behave like a LM at some limiting condition breaks down due to such an improper matching of solvation free energy components.     

Modelling (vapour + liquid) and (vapour + liquid + liquid) equilibria of {water (H2O) + methanol (MeOH) + dimethyl ether (DME) + carbon dioxide (CO2)} quaternary system using the Peng–Robinson EoS with Wong–Sandler mixing rule

The (vapour + liquid) equilibria (VLE) and (vapour + liquid + liquid) equilibria (VLLE) binary data from literature were correlated using the Peng–Robinson (PR) equation of state (EoS) with the Wong–Sandler mixing rule (WS). Two group contribution activity models were used in the PRWS: UNIFAC–PSRK and UNIFAC–Lby. The systems were successfully extrapolated from the binary systems to ternary and quaternary systems. Results indicate that the PRWS–UNIFAC–PSRK generally displays a better performance than the PRWS–UNIFAC–Lby.     

基于PSRK模型预测高分子溶液汽液平衡

为了更好地预测高分子溶液的汽液平衡,本文应用PSRK模型改进的混合规则[即将组合项(∑xiInb／bi)与原UNIFAC模型中的组合项(Flory-Huggins项)相抵消所获得的简化的混合规则]来计算方程的能量参数α。为了计入自由体积效应对高分子溶液的过量Gibbs自由能的贡献,在计算体积参bij时,将原混合规则中的指数修订为bij^1/2(bi^1/2 bj^1/2)/2。应用PSRK模型结合改进的混合规则预测二元高分子溶液体系的汽液平衡,结果表明新模型的计算精度是令人满意的。     

A local composition extension of the van der Waals mixing rule for PR cubic equation of state

New mixing rules (VWLC-I and II) capable of connecting van der Waals (VDW) to CEOS/A^E mixing rule models were developed. These models are able to incorporate the same multi-component mixture parameters obtained for the van der Waals and CEOS/A^E models simultaneously. The VWLC mixing rules directly incorporate local compositions into the cubic equations of state (CEOS). The energy parameters required for the local compositions are calculated from the CEOS parameters. The Peng–Robinson (PR) CEOS was used for this study. Binary interactions parameters were obtained by adjusting the vapor pressure of the binary mixture for several low and high-pressure systems. The predictive capabilities of the VWLC mixing rules were tested by vapor–liquid equilibria calculations for low and high-pressure multicomponent systems. The results were compared with the predictions of the VDW mixing rule and a Huron–Vidal (HV) kind of CEOS/A^E-NRTL mixing rule. The VWLC mixing rules are consistent models giving good results in a broad range of pressures and temperatures in binary and multicomponent mixtures. They compare favorably with the CEOS/A^E-NRTL mixing rule for low-pressure systems. In high-pressure ternary systems VWLC-I and II give good predictions, much better, in fact, than the CEOS/A^E-NRTL mixing rule.     

Interaction parameters for multi-component aromatic extraction with sulfolane

Aromatic extraction is an important operation in petrochemical processing. Design of an aromatic extractor requires the knowledge of multi-component liquid–liquid equilibrium (LLE) data. Such experimental LLE data are usually not available and therefore can be predicted using various activity coefficient models. These models require proper binary interaction parameters, which are not yet available for all aromatic extraction systems. Furthermore, the parameters available for most of the ternary systems are specific to that system only and cannot be used for other ternary or multi-component systems. An attempt has been made to obtain these parameters that are globally applicable. For this purpose, the parameter estimation procedure has been modified to estimate the parameters simultaneously for different systems involving common pairs. UINQUAC and UNIFAC models have been used for parameter estimation. The regressed parameters are shown to be applicable for the ternary as well as for the multi-component systems. It is observed that UNIQUAC parameters provide a better fit for ternary LLE data, whereas, as one moves towards the higher component systems (quaternary and quinary) the UNIFAC parameters, which are a measure of the group contributions, predict the LLE better. Effect of temperature on UNIQUAC binary interaction parameters has been studied and a linear dependence has been observed.     

PSRK group contribution equation of state: comprehensive revision and extension IV,including critical constants and α-function parameters for 1000 components

As a part of an ongoing process, the predictive Soave–Redlich–Kwong (PSRK) group contribution equation of state was extended by the introduction of additional structural groups (F₂, Cl₂, Br₂, HCN, NO₂, CF₄, O₃ and ClNO) and fitting of the corresponding group interaction parameters. Interaction parameters between already existing main groups were also optimized to the growing literature data base. Overall, 75 new parameter sets are given herein, and typical results are presented for various systems. For the sake of completeness, not only the group new interaction parameters but all available PSRK/UNIFAC interaction parameter sets (more than 900) are given as supplementary material. Moreover, the required pure component properties (critical properties, acentric factors, and Mathias–Copeman constants) were revised and are also included for about 1000 components.     

A comparison of mixing rules for the combination of COSMO-RS and the Peng–Robinson equation of state

In this work, the COSMO-RS model is combined with a volume-translated Peng–Robinson equation of state (EOS) via a G^E$G^{\text{E}}$-based mixing rule. The performance of several mixing rules previously published for this purpose is compared and semi-empirical modifications to one of them are introduced to improve its performance in our application. The new mixing rule contains three internal parameters that are adjusted to achieve consistency between the mixing rule and COSMO-RS. No experimental binary data is needed for our EOS. The new COSMO-RS-based, predictive EOS introduces a density dependence into COSMO-RS and extends its applicability to higher pressures and to mixtures containing supercritical components.     

Application of a volume-translated Peng-Robinson equation of state on vapor-liquid equilibrium calculations

A volume-translated Peng-Robinson (VTPR) equation of state (EOS) is developed in this study. Besides the two parameters in the original Peng-Robinson equation of state, a volume correction term is employed in the VTPR EOS. In this equation, the temperature dependence of the EOS energy parameter was regressed by an improved expression which yields better correlation of pure-fluid vapor pressures. The volume correction parameter is also correlated as a function of the reduced temperature. The VTPR EOS includes two optimally fitted parameters for each pure fluid. These parameters are reported for over 100 nonpolar and polar components. The VTPR EOS shows satisfactory results in calculating the vapor pressures and both the saturated vapor and liquid molar volumes. In comparison with other commonly used cubic EOS, the VTPR EOS presents better results, especially for the saturated liquid molar volumes of polar systems. VLE calculations on fluid mixtures were also studied in this work. Traditional van der Waals one-fluid mixing rules and other mixing models using excess free energy equations were employed in the new EOS. The VTPR EOS is comparable to other EOS in VLE calculations with various mixing rules, but yields better predictions on the molar volumes of liquid mixtures.     

Experimental measurement and modeling of the vapor–liquid equilibrium of carbon dioxide + chloroform

Isothermal bubble and dew points, saturated molar volumes, and mixture critical points for binary mixtures of carbon dioxide+chloroform (trichloromethane) (CO₂/CHCl₃) have been measured in the temperature region 303.15–333.15 K and at pressures up to 100 bar. Mixture critical points are reported at 313.15, 323.15, and 333.15 K. The data were modeled with the Peng–Robinson equation of state using both the van der Waals-1 (vdW-1) mixing rule and the Wong–Sandler (WS) mixing rule incorporating the UNIQUAC excess free energy model. The WS mixing rule provided a better representation of the data than did the vdW-1 mixing rule, though with three adjustable parameters instead of one. The extrapolating ability of both of the mixing rules was investigated. Using the parameters regressed at 323.15 K, the WS mixing rule yielded better extrapolations for the composition dependence at 303.15, 313.15, and 333.15 K than the vdW-1 mixing rule.     

Vapor-liquid equilibrium predictions of refrigerant mixtures from a cubic equation of state with a G-EoS mixing rule

A modified excess Gibbs energy model which is based on the local composition concept and assigns a single energy parameter per pair of components, is incorporated into the G^E—EoS thermodynamic formalism for vapor-liquid equilibrium (VLE) calculations of simple and complex refrigerant mixtures. One temperature set of data close to 273 K is used to obtain the model's parameters, which are used to extrapolate the VLE at other temperatures and pressures. A one-parameter form of the model based on the Wong-Sandler mixing rule is presented for several simple systems. The physical significance of the model's energy parameter is connected to the preference of the mixture for like to unlike interactions. The model is applied for VLE predictions of the ternary system R14-R23-R13, and the results are compared to calculations using the 3PWS model [H. Orbey. S.I. Sandler, Ind. Eng. Chem. Res. 34 (1995) 2520–2525] and the van der Waals mixing rule. Modelling of a few complex systems with only three data points given at each temperature is shown with a two parameter version of our model on the basis of the Huron-Vidal mixing rule.     

New mixing rule for predicting multi-component gas adsorption

The present work describes a predictive model for ascertaining the multi-component gas adsorption equilibria. The model utilizes special form of covolume-dependent (CVD) mixing which is combined with the generalized form of 2-D EOS. Four well known 2-D EOSs; van der Waals, Soave-Redlich-Kwong, Peng-Robinson, Eyring along with the modified CVD mixing rule were used to predict the total adsorption of several binary and ternary systems. Based on the concept of the CVD mixing rule, it was inspired that CVD mixing rule could be a binding bridge between the molecular size and the molecular interaction. To show this, the ratio of the classical mixing rule %AAD to the CVD mixing rule %AAD were plotted versus the difference of the collision or the Leonard-Jones diameters of the gas molecules in the mixtures. It shows that there is a criterion between the CVD and the classical mixing rules in terms of molecular size difference. It seems that, Δσ _LJ≈0.60 Å is the criterion. The CVD mixing rule is approximately predominant in the region of Δσ _LJ≥0.60 Å, whilst, region of Δσ _LJ≤0.60 Å is nearly governed by the classical mixing rule. All predictions by the new mixing rule and the classical mixing rule were compared with the experimental data from the case studies. The new form of the mixing rule is in good agreement with the experimental data even for the non-ideal systems; hence provides a powerful framework to predict multi-component gas adsorption.     

Phase equilibria of the ternary system vinyl acetate/(R,S)-1-phenylethanol/carbon dioxide at high pressure conditions

Phase equilibrium measurements, correlations and predictions are presented for the binary systems (R,S)-1-phenylethanol/CO₂ and vinyl acetate/CO₂ and for the ternary system vinyl acetate/(R,S)-1-phenylethanol/CO₂. Experiments for the ternary system were performed in the temperature range of 323–343 K and in the pressure range of 7–12 MPa, using a high pressure phase equilibrium apparatus with a high pressure visual variable volume cell. Phase compositions were determined by taking samples of each phase and analysing them by gas chromatography. Equilibrium data were correlated with the Peng–Robinson equation of state combined with the Mathias–Klotz–Prausnitz mixing rule. A good correlation of both phases behaviour was obtained with an average absolute deviation (AAD) of 6.80%. Predictions for the binary sub-systems and for the ternary system were performed using the Peng–Robinson and the Soave–Redlich–Kwong equation of state, with the predictive mixing rule MHV1.     

Mixing rules for van-der-Waals type equation of state based on thermodynamic perturbation theory

A concept based on the thermodynamic perturbation theory for a `simple fluid' has been applied to the attractive term of a van-der-Waals type equation of state (EOS) to derive a simple mixing rule for the a parameter. The new mixing rule is a small correction to the original one-fluid approximation to account for the influence of particles of j-type on the correlation function of ii-type in a mixture consisting of particles of i and j types. The importance of the correction has been shown by comparison of the calculated results for binary mixtures of Lennard–Jones fluids with the data obtained by numerical method (Monte-Carlo simulation). The new mixing rules can be considered as a flexible generalization of the conventional mixing rules and can be reduced to the original v-d-W mixing rules by defaulting the extra binary parameters to zero. In this way the binary parameters already available in the literature for many systems can be used without any additional regression work. Extension of the new mixing rules to a multicomponent system do not suffer from `Michelsen–Kistenmacher syndrome' and provide the correct limit for the composition dependence of second virial coefficients. Their applicability has been illustrated by various examples of vapor–liquid and liquid–liquid equilibria using a modified Patel–Teja EOS. The new mixing rules can be applied to any EOS of van-der-Waals type, i.e., EOS containing two terms which reflect the contributions of repulsive and attractive intermolecular forces.     

The experiments and correlations of the solubility of ethylene in toluene solvent

Polymer cyclic olefin copolymer (COC) is produced from the reaction of attaching ethyl groups to the norbornene monomer in liquid phase. The first step of process is dissolving ethylene in a liquid phase where toluene is present as a cosolvent. Thus, the solubility of ethylene in liquid toluene is the most important factor affecting the production of COC. In this study, the solubility of ethylene in toluene was measured in the temperature range from 323.15 to 423.15 K and pressure range from 5 to 25 bar. The experiments were conducted by the method of pressure decaying with a newly designed apparatus. The experimental results show that the solubility of ethylene in toluene increases with increasing pressure but decreases with increasing temperature.The experimental solubility data were expressed in the vapor–liquid equilibrium relationship and correlated fairly well by the bubble–pressure calculation with the Peng–Robinson equation of state (PR EOS) incorporated with the van der Waals one-fluid and the Zhong–Masuoka mixing rules with the consideration of binary interaction parameters. The results showed the van der Waals (vdW-1) mixing rule is slightly better than the Z–M mixing rule for pressure correlation but the Z–M mixing rule is slightly better for vapor composition correlation.A semi-empirical solubility equation with four parameters for the present binary system was proposed in this study. This proposed model estimates the solubility easier and as accurate as the PR EOS does for the present system.     

Equation of state for the NH3−H2O system

An equation of state (EOS) for the NH₃–H₂O system has been developed. This EOS incorporates a highly accurate end-member EOS and on an empirical mixing rule. The mixing rule is based on an analogy with high order contributions to the virial expansion for mixtures. Comparison with experimental data indicates that the mixed system EOS can predict both phase equilibria and volumetric properties for this binary system with accuracy close to that of the experimental data from 50°C and 1 bar to critical temperatures and pressures.