PSRK模型的修正混合规则及其应用于聚合物溶液汽液平衡的预测

To extend the PSRK(predictive Soave-Redlich-Kwong equation of state)model to vapor-liquid equilibria of polymer solutions,a new EOS-gE mixing rule is applied in which the term ∑ xi ln(b/bi)in the PSRK mixing rule for the parameter a,and the combinatorial part in the original universal functional activity coefficient(UNIFAC)model are cancelled.To take into account the free volume contribution to the excess Gibbs energy in polymer solution,a quadratic mixing rule for the cross co-volume bij with an exponent equals to 1/2 is applied [b1/2 ij= 1/2(b1/2 i b1/2 j)].The literature reported Soave-Redlich-Kwong equation of state(SRK EOS)parameters of pure polymer are employed.The PSRK model with the modified mixing rule is used to predict the vapor-liquid equilibrium(VLE)of 37 solvent-polymer systems over a large range of temperature and pressure with satisfactory results.     

汽液相平衡计算混合法则中协体积项b的改进

混合物汽液相平衡通常需要引入合适的混合法则才能更好地进行关联。针对目前混合法则中协体积项b的改进提出了一种新的修正协体积项b的方法。该方法以Mie势能理论、London色散力理论为基础,同时引入Leach等分子形状系数进行修正,这样协体积项b中的交互作用参数l_ij形式上是纯组分分子形状系数θ和φ及临界参数的函数,无需实验数据拟合得到,从而真实反映了二元混合物的实际混合情况。将修正的协体积项b运用于MHV1、LCVM及HV混合法则并结合PR NRTL模型对不同种类的16个体系汽液相平衡进行了计算,并与采用b ∑x_ib_i的MHV1、LCVM及HV混合法则结合PR NRTL模型计算的结果进行比较,结果表明运用修正的协体积项b的混合法则+PR+NRTL模型所得的计算结果精度优于采用b ∑x_ib_i的MHV1、LCVM及HV混合法则 PR NRTL模型的计算结果。     

状态方程模拟醇胺系统的密度和汽液相平衡

通过考虑醇胺分子间的缔合作用,结合先前开发的非缔合变阱宽链流体状态方程(SWCF-VREOS)建立了一个缔合方阱链流体状态方程,并利用方程模拟了醇胺系统的密度和汽液相平衡。通过关联不同温度下醇胺的饱和蒸气压和液体体积得到了18种醇胺流体的分子参数,新方程计算的饱和蒸气压和液体密度总的平均误差分别为0.94%和0.88%。结合简单的混合规则,将此方程扩展到混合系统。研究发现,建立的方程可预测二元和三元醇胺混合物的密度。当引入一个与温度无关的可调参数时,方程能满意关联二元系统的汽液相平衡数据,并可进一步预测多元混合系统的汽液相平衡,预示着新方程可模拟醇胺系统的相行为。     

A modification of Wong-Sandler mixing rule for the prediction of vapor-liquid equilibria in binary asymmetric systems

Systems consisting of light components and heavy hydrocarbons are highly asymmetric and industrially important. Design and control of facilities for separation and purification of such mixtures require vapor-liquid equilibrium data. Coupling of the cubic equation of state (EOS) with excess Gibbs energy models (EOS/G ^ex models) failed to represent the vapor-liquid equilibria (VLE) of such systems accurately. The main purpose of this work is to present a modification of Wong-Sandler mixing rule with using the composition dependent binary interaction parameter. Vaporliquid equilibria for 30 binary systems are calculated using the SRK equation of state with proposed model and Wong-Sandler mixing rule. Calculated pressures and mole fractions of vapor phase are compared with experimental data. The average absolute percentage deviation indicates that error involved in the application of modified Wong-Sandler model is less than Wong-Sandler model in most cases.     

A novel equation of state (EOS) for prediction of solute solubility in supercritical carbon dioxide: Experimental determination and correlation

Solubility data of organophosphorous metal extractants in supercritical fluids (SCF) are crucial for designing metal extraction processes. We have developed a new equation of state (EOS) based on virial equation including an untypical parameter as BP/RT, reduced temperature and pressure for prediction of solute solubility in supercritical carbon dioxide (SC CO₂). Solubility experimental data (solubility of tributylphosphate in SC CO₂) were correlated with the two cubic equations of state (EOS) models, namely the Peng–Robinson EOS (PR‐EOS) and the Soave–Redlich–Kwong EOS (SRK‐EOS), together with two adjustable parameter van der Waals mixing and combining rules and our proposed EOS. The AARD of our EOS is significantly lower than that obtained from the other EOS models. The proposed EOS presented more accurate correlation for solubility data in SC CO₂. It can be employed to speed up the process of SCF applications in industry.     

Thermodynamic modeling of vapor-liquid equilibria and excess properties of the binary systems containing diethers and n-alkanes by cubic equation of state

A comparison of the performances of two different approaches of cubic equations of state models, based on a classical van der Waals and mixing rules incorporating theG ^E equation, was carried out for correlation of Vapor-Liquid Equilibria (VLE), H^E and C _P ^E data alone, and simultaneous correlation of VLE+H^E, VLE+C _P ^E , H^E +C _P ^E and VLE+H^E +C _P ^E data for the diethers (1,4-dioxane or 1,3-dioxolane) with n-alkane systems. For all calculations the Peng-Robinson-Stryjek-Vera cubic equation of state (PRSV CEOS) was used. A family of mixing rules for the PRSV CEOS based on the Modified van der Waals one-fluid mixing rule (MvdW1) and two well-known CEOS/G^E mixing rules (MHV1 and MHV2), was considered. The NRTL equation, as the G^E model with linear or reciprocal temperature dependent parameters, was incorporated in the CEOS/G^E models. The results obtained by the CEOS/G^E models exhibit significant improvement in comparison to the MvdW1 models.     

Fluid‐phase coexistence for the oxidation of CO2 expanded cyclohexane: Experiment,molecular simulation,and COSMO‐SAC

The gas solubility of pure oxygen and of pure carbon dioxide as well as of their gaseous mixture are measured in the ternary liquid mixture cyclohexane + cyclohexanone + cyclohexanol at 313.6 K with a high‐pressure view‐cell technique using the synthetic method. The new experimental data are used to assess the capability of molecular simulation and conductor‐like screening model (COSMO)‐SAC to predict multicomponent fluid‐phase coexistence behavior. These methods are also compared systematically on the basis of experimental binary fluid‐phase coexistence data. In that comparison also the Peng–Robinson (PR) equation of state is included as a reference. Molecular simulation and COSMO‐SAC yield good results and are found to be far superior to the PR equation of state both in predictive and in adjusted mode. © 2013 American Institute of Chemical Engineers AIChE J, 59: 2236–2250, 2013     

基于无限稀释活度系数模型的状态方程

发展了一种从纯物质的性质参数预测混合物汽液平衡的方法,它是以立方状态方程的过量吉布斯自由能 g~E结合型混合规则和无限稀释活度系数 MOSCED模型为基础,由MOSCED模型预测得到活度系数方程的相互作用参数,从而得到混合规则中的g~E,使得状态方程能在较大温度和压力范围内预测汽液平衡.用此法预测了含极性组分体系常压到较高压力的汽液平衡,结果良好,与UNIFAC基团贡献状态方程的结果相当.     

立方型状态方程含过量自由焓模型的新混合规则

提出了在参考压力下将过量自由焓(g~E)模型引入状态方程的修正Huron-Vidal混合规则,由g~E模型确定状态方程中混合物的参数,以SRK方程和Wilson模型为例,取大气压为参考压力,计算了16个体系50组汽液平衡,结果表明本文建议的混合规则可直接使用现有文献报道的常压g~E模型参数由立方型状态方程预测常压汽液平衡,并对直接外推预测高压汽液平衡作了尝试.     

AN EMPIRICAL VIRIAL-LIKE CUBIC EQUATION OF STATE FOR PHASE EQUILIBRIUM CALCULATIONS

An empirical cubic equation of state(EOS) was obtained by truncating the virial expansion in reciprocal of molar volume after the third term. The constants of the EOS was generalized in terms of critical temperature, critical pressure and Pitzer's acentric factor.

In pure component applications the EOS exhibited a performance comparable to Peng-Robinson (1976) EOS in the reduced temperature range of 0.5 to 1. The present EOS tends to predict better saturation liquid volumes at reduced temperatures below 0.8, and better estimations for second virial coefficient at high reduced temperatures.

The EOS was successfully employed for vapor liquid equilibrium calculations for some mixtures of normal or slightly polar fluids with traditional one binary parameter mixing rule at moderately high pressures. At low reduced temperatures, where conventional one adjustable parameter applications of the cubic equations compare unfavorably with dual methods based on excess Gibbs energy functions for the liquid phase, a new two constant mixing rule introduced by Stryjek and Vera (1986) was employed for the present equation of state.     

一种新的WS型混合法则

立方形状态方程被广泛应用于相平衡计算中．绝大多数状态方程是结合简单的van der Waals混合法则来进行VLE数据关联的．然而对于应用这些方程关联高度非理想混合物的复杂的相行为,仅用van der waals单维流体混合法则是不够的．甚至在van der Waals混合法则引入另外一个参数都不能够精确地关联含有极性物质的混合物．     

Modeling of solid–liquid equilibria for polyethylene and polypropylene solutions with equations of state

We modeled solid–liquid equilibria (SLEs) in polyethylene and polypropylene solutions with a Soave–Redlich–Kwong (SRK) cubic equation of state (EOS) and a perturbed‐chain statistical associating fluid theory (PC‐SAFT) EOS. Two types of mixing rules were used with SRK EOS: The Wong–Sandler mixing rule and the linear combination of the Vidal and Michelsen mixing rules (LCVM), both of which incorporated the Bogdanic and Vidal activity coefficient model. The performance of these models was evaluated with atmospheric‐pressure and high‐pressure experimental SLE data obtained from literature. The basic SLE equation was solved for the equilibrium melting temperature instead of for the composition. The binary interaction parameters of SRK and PC‐SAFT EOS were estimated to best describe the experimental equilibrium behavior of 20 different polymer–solvent systems at atmospheric pressure and 31 other polymer–solvent systems at high pressure. A comparison with experimental data showed that SRK–LCVM agreed very well with the atmospheric SLE data and that PC‐SAFT EOS was more efficient in high‐pressure conditions. © 2011 Wiley Periodicals, Inc. J Appl Polym Sci, 2011     

Modeling of Nitrogen and Carbon Dioxide Solubility in Alternative Fuels at High Pressures Using the Soave‐Redlich‐Kwong Equation of State

In this paper, the Soave‐Redlich‐Kwong equation of state with quadratic mixing rule has been tested for correlation of vapor‐liquid equilibria (VLE) at high pressures in the binary nitrogen + dimethyl ether, dimethyl ether + methanol, nitrogen + methanol, carbon dioxide + dimethyl ether and carbon dioxide + methanol systems. The interaction parameters k_ij were evaluated for each binary pair and used for prediction of VLE in the ternary nitrogen + dimethyl ether + methanol and carbon dioxide + dimethyl ether + methanol systems at high pressures. The results of correlation and prediction are discussed.     

Modeling high pressure vapor–liquid equilibrium of ternary systems containing supercritical CO2 and mixed organic solvents using Peng–Robinson equation of state

High pressure vapor–liquid equilibrium (VLE) of CO₂-expanded organic solvents was investigated using Peng–Robinson-LCVM-UNIFAC equation of state. Bubble pressure of several ternary mixtures was predicted using this model and correlations were developed based only on binary experimental data. A sensitivity study of the LCVM parameter numerical value was done by considering the coherence between the mathematical features of the mixing rule and the quality of the simulation. The results provided by PR-LCVM-UNIFAC were compared with those ones given by Peng–Robinson equation of state using the classical quadratic mixing rules (PR-CMR). Despite the use of two adjustable parameters for each binary system, PR-CMR is not able to provide good results when applied to ternary systems. The capability of PR-LCVM-UNIFAC model to predict liquid mixture density for ternary systems using parameters regressed only from bubble pressure experimental data was also investigated. Due to the lack of liquid density experimental data, it was possible to perform only a qualitative assessment of the density curves calculated by this equation of state.     

Thermodynamic modeling of vapor-liquid equilibrium of binary systems ionic liquid + supercritical {CO2 or CHF3} and ionic liquid + hydrocarbons using Peng-Robinson equation of state

Vapor-liquid equilibrium (VLE) data from literature for binary systems involving several ionic liquids were correlated. The Peng-Robinson equation of state, coupled with the van der Waals and Wong-Sandler mixing rules, was used as the thermodynamic model to evaluate the fugacity coefficients. The UNIQUAC and NRTL models were used to calculate the excess Gibbs free energy in the Wong-Sandler mixing rule. A molecular modeling strategy using the software ChemOffice was used to calculate the volume and surface area parameters of ionic liquids for UNIQUAC, while the binary interaction energy parameters for UNIQUAC and NRTL models, as well as the binary interaction parameter of the van der Waals and Wong-Sandler mixing rules were estimated through a method based on the genetic algorithm. The results show that, as expected, the Wong-Sandler mixing rules represented better the data, with both activity coefficient models showing high accuracy. However, in one case, NRTL predicted an erroneous azeotropic condition, while UNIQUAC was able to correlate the data without this error.     

新的氨水体系状态方程

本文提出一个以三参数立方型方程为基础的氨水体系状态方程。由同时拟合蒸汽压、汽相体积和液相体积确定的三个参数都是温度的函数。混合物的混合法则考虑了三分子相互作用对引力的贡献。四个交互作用参数处理成压力的线性关系后,在压力小于21atm下预测汽液平衡和热力学性质达到满意的结果。     

A practical equation of state for non‐spherical and asymmetric systems for application at high pressures. Part 2: Extension to mixtures

In part I of this series the pure component PHCT‐DNSK equation of state (EOS) was presented. In this paper the EOS is extended to describe mixtures, particularly asymmetric mixtures containing one or more low molecular weight spherical compound together with one or more high molecular weight chain‐like compound. The EOS utilises theoretically correct mixing rules and is generally able to predict the correct trends quantitatively for binary mixtures, and in most cases outperform other EOSs. With the use of a small, temperature independent, interaction parameter the EOS is able to predict the phase behaviour of the investigated systems qualitatively. The EOS is able to predict the phase behaviour of a multi‐component system containing one or more light components and a range of heavy hydrocarbons with improved accuracy compared to other EOSs at reduced computational times. © 2011 Canadian Society for Chemical Engineering     

EMPIRICAL TWO-PARAMETER MIXING RULES FOR A CUBIC EQUATION-OF-STATE

Four two-parameter mixing rules are presented along with expressions for the fugacity coefficient and compressibility factor obtained when they are used in conjunction with the Peng-Robinson equation of state. A comparison of the forms of these expressions and results obtained in the prediction of VLE for nine systems with these temperature-, pressure-, density-, and composition-dependent rules provides a guideline for the applicability, advantages and disadvantages of each. Excellent results were obtained, especially in the near critical region, when these simple pressure- and density-dependent rules were used to describe P-x data. The density-dependent mixing rule consistently yielded better predictions of VLE over a range of temperatures (using parameters obtained at a single temperature) than the results obtained by optimizing the interaction parameter of conventional mixing rule at each temperature. The composition-dependent rule was slightly less effective for most systems, but for highly polar, asymmetric systems, this rule was superior     

Calculation of the solid solubilities in supercritical carbon dioxide using a new G mixing rule

The aim of this study is to develop a new EOS/G^ex-type mixing rule with special attention to calculating the solid solubilities of aromatic hydrocarbons, aliphatic carboxylic acids, aromatic acids, and heavy aliphatic and aromatic alcohols in supercritical carbon dioxide. A volume correction term is applied with a combination of second and third virial coefficients which the equation for the third virial coefficient is quadratic, according to the suggestion by Hall and Iglesias-Silva. In this study, the cubic Peng-Robinson (PR) and Soave-Redlich-Kwong (SRK) equations of state have been used to calculate the solid solubilities of 23 solutes in supercritical CO₂, by using six mixing rules, namely, the Wong-Sandler (WS) rule, the Orbey-Sandler (OS) rule, the van der Waals one fluid rule with one (VDW1) and two (VDW2) adjustable parameters, the covolume dependent (CVD) rule and the new mixing rule. In all cases, the NRTL model was chosen as the excess Gibbs free energy model. The coefficients of the NRTL model and the binary interaction parameters of six mixing rules with two EOSs (PR and SRK EOSs) have been determined for 100 data sets of 23 binary systems over a wide range of temperatures and pressures covering more than 970 experimental data points which are reported in the literature. The results show that the PR EOS with the new mixing rule model is more accurate than the PR and SRK EOSs with the other mixing rules for solid solubility calculations in supercritical carbon dioxide.The regressed interaction parameters of the binary system, without any further modification, were then extended to four ternary mixtures, giving satisfactory results of the solid solubilities in supercritical CO₂.     

Generalized binary interaction parameters for hydrogen-heavy-n-alkane systems using Peng–Robinson equation of state

The Peng–Robinson equation of state (PR EOS) is used to model the vapor–liquid equilibria (VLE) of binary systems of hydrogen with five heavy-n-alkanes: n-decane (n-C₁₀H₂₂), n-hexadecane (n-C₁₆H₃₄), n-octacosane (n-C₂₈H₅₈), n-hexatriacontane (n-C₃₆H₇₄), and n-hexatetracontane (n-C₄₆H₉₄). Using literature experimental data for these systems, binary interaction parameters (BIPs) were calculated using the PR EOS coupled with three different alpha (α) functions (Soave, Twu, and Gasem). The calculated BIPs have been fitted to a generalized correlation that can be used to estimate the BIPs and model the VLE within the temperature range of 283.2–449.6?K, pressure range of 1.151–15.970?MPa, and hydrogen solubility range of 0.016–0.257?mole fraction. It is found that PR EOS combined with one specific form of α function is capable of reproducing the experimental VLE data with an overall %AARD of 1.1%.