A new cubic simplified perturbed hard-body equation of state

A new cubic equation of state (EOS) was developed in this study for vapor-liquid equilibrium (VLE) calculations of nonpolar fluids. The repulsive term of this EOS reexpressed the results of Walsh and Gubbins (1990) from their modified thermodynamic perturbation theory of polymerization into a simple form in which a non-spherical parameter was employed to account for the different shapes of molecules. The repulsive compressibility factors calculated from this EOS agree well with the molecular simulation data for various kinds of hard bodies ranging from a single hard sphere to tangent or fused long chain molecules. A simple attractive term was then coupled with the repulsive to complete the EOS in a cubic form. Equation parameters were determined for a diversity of nonpolar real fluids. These parameters were expressed in generalized forms for engineering computations. Satisfactory results from this EOS on the saturated properties of pure nonpolar fluids were obtained. This EOS was also extended to calculate the VLE of nonpolar fluid mixtures. The results are again satisfactory over wide ranges of temperature and pressure.     

A Thermodynamic Model for Square-well Chain Fluid: Theory and Monte Carlo Simulation

A thermodynamic model for the freely jointed square-well chain fluids was developed based on the thermodynamic perturbation theory of Barker-Henderson, Zhang and Wertheim. In this derivation Zhang's expressions for square-well monomers improved from Barker-Henderson compressibility approximation were adopted as the reference fluid, and Wertheim's polymerization method was used to obtain the free energy term due to the bond connectivity. An analytic expression for the Helmholtz free energy of the square-well chain fluids was obtained. The expression without adjustable parameters leads to the thermodynamic consistent predictions of the compressibility factors, residual internal energy and constant-volume heat capacity for dimer, 4-mer, 8-mer and 16-mer square-well fluids. The results are in good agreement with the Monte Carlo simulation. To obtain the MC data of residual internal energy and the constant-volume heat capacity needed, NVT MC simulations were performed for these square-well chain fluids.     

The extension of coordination number model by using the local composition theory: mixtures

The extension of a new coordination number model to mixture is presented in this work. Extended model agrees well with the Monte Carlo (MC) simulation results for square-well (SW) mixture fluids and shows better results compared with other models. To test our model, we compare the compressibility factors from various models for SW fluids at different λ values and for SW fluid mixtures at λ=1.5. Although our model is obtained by fitting simulation data at λ=1.5, it shows better results for the different λ values than other coordination number model. Compared with the compressibility factors of various binary mixtures of SW fluids calculated from other models, this model presents better results. Because our model considers the temperature dependency importantly by using the total site number, it predicts coordination number and compressibility factor well in the wide temperature range and enables one to derive an equation of state (EOS) through integration of the coordination number equation.     

From amorphous solid to defective crystal. A study of structural peculiarities in close packings of hard spheres

The structure of models of close packings of hard spheres is examined, with density of the packings exceeding the maximum value for uniform disordered packings (η = 0.64). High densities are achieved due to regions of the closest packing (η = 0.74) emerging in the models. Spatial geometry of good tetrahedral atomic configurations and of simplest elements of the crystal structure identified by Delaunay simplices was studied using the Voronoi— Delaunay method. Models with a packing coefficient η varying from 0.639 to 0.706 were considered. At smaller densities, a well-known disordered close “Bernal packing” is realized. At η = 0.706 (the greatest density achieved), a unified crystal structure with numerous defects is formed. At intermediate densities, stochastically oriented crystalline nuclei are observed. Specific atomic aggregates — stacks of five-membered rings in good tetrahedral configurations of spheres — are revealed in models having a substantial fraction of crystalline phase (η = 0.664). Such non-trivial structures can occur only in packings that are intermediate between amorphous and crystalline phases.     

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.     

On the Calculation of Surface Tensions of <Emphasis Type="Italic">n</Emphasis>-Alkanes Using the Modified SAFT-BACK-DFT Approach

Density functional theory is combined with the modified SAFT-BACK EOS to investigate liquid–vapor interfaces of n-alkanes. We evaluate the temperature dependence of the interfacial width and the surface tension. Differences in chain length of the alkanes lead to differences in the thermodynamic properties of the fluids. A single curve for the reduced width of the interface as a function of reduced temperature serves to correlate interfacial properties of a wide variety of linear chain fluids (excluding methane and ethane) with sufficient accuracy for our purposes.     

Calculating thermodynamic properties from perturbation theory: II. An analytic representation for the square-well chain fluid

An analytic representation of thermodynamic properties of the freely jointed square-well chain fluid is developed based on the thermodynamic perturbation theory of Barker–Henderson, Zhang and Weitheim. By using a real function expression for the radial distribution function and incorporating structural information for square-well monomer of TPT1 model, an analytic expression for the Helmholtz energy of square-well chain fluid is expanded from Zhang’s analytic expressions for thermodynamic properties of square-well monomer. The expression leads to good predictions of the compressibility factor, residual internal energy and constant-volume heat capacity for 4-mer, 8-mer and 16-mer square-well fluids when compared with the Monte Carlo (MC) simulation results. The incorporating structural information for square-well dimer of TPT-D model is also calculated. To obtain the constant-volume heat capacity needed, NVT MC simulations were performed.     

A new coordination number model by using the local composition theory

A new coordination number model for square-well (SW) fluid is proposed in this work. It is based on the local composition theory and starts from the point that total site numbers around one molecule can vary with reduced density and reduced temperature. The total site numbers are correlated with reduced density and reduced temperature by using local composition theory and Monte Carlo (MC) simulation results. This new model agrees well with the MC simulation results for the pure SW fluids at the low-density limit and in the high-density region. In addition, it can predict the transition to the face-centered cubic (fcc) fluid property which is shown in Young and Alder’s [J. Chem. Phys. 73 (1980) 2430] phase diagram. An equation of state (EOS) is derived from this new coordination number model and gives satisfactory compressibility factors for SW fluids.     

A completely analytic equation of state for the square-well chain fluid of variable well width

A completely analytic perturbation theory equation of state for the freely-jointed square-well chain fluid of variable well width (1 ≤ λ ≤ 2) is developed and tested against Monte Carlo simulation data. The equation of state is based on second-order Barker and Henderson perturbation theory to calculate the thermodynamic properties of the reference monomer fluid, and on first-order Wertheim thermodynamic perturbation theory to account for the connectivity of monomers to form chains. By using a recently developed real function expression for the radial distribution function of hard spheres in perturbation theory, we obtain analytic, closed form expressions for the Helmholtz free energy and the radial distribution function of square-well monomers of any well width. This information is used as the reference fluid in the perturbation theory of Wertheim to obtain an analytic equation of state, without adjustable parameters, that leads to good predictions of the compressibility factors and residual internal energies for 4-mer, 8-mer and 16-mer square-well fluids when compared with the simulation results. Further, very good results are obtained when this equation of state with temperature-independent parameters is used to correlate the vapor pressures and critical points of the linear alkanes from methane to n-decane.     

Study on vapor-liquid equilibria and surface tensions for nonpolar fluids by renormalization group theory and density gradient theory

An equation of state (EOS) applicable for both the uniform and nonuniform fluids is established by using the density-gradient theory (DGT). In the bulk phases, the EOS reduces to statistical associating fluid theory (SAFT). By combining the EOS with the renormalization group theory (RGT), the vapor-liquid-phase equilibria and surface tensions for 10 nonpolar chainlike fluids are investigated from low temperature up to the critical point. The obtained results agree well with the experimental data.     

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.     

Confinement, entropy, and single-particle dynamics of equilibrium hard-sphere mixtures

We use discontinuous molecular dynamics and grand-canonical transition-matrix Monte Carlo simulations to explore how confinement between parallel hard walls modifies the relationships between packing fraction, self-diffusivity, partial molar excess entropy, and total excess entropy for binary hard-sphere mixtures. To accomplish this, we introduce an efficient algorithm to calculate partial molar excess entropies from the transition-matrix Monte Carlo simulation data. We find that the species-dependent self-diffusivities of confined fluids are very similar to those of the bulk mixture if compared at the same, appropriately defined, packing fraction up to intermediate values, but then deviate negatively from the bulk behavior at higher packing fractions. On the other hand, the relationships between self-diffusivity and partial molar excess entropy (or total excess entropy) observed in the bulk fluid are preserved under confinement even at relatively high packing fractions and for different mixture compositions. This suggests that the excess entropy, calculable from classical density functional theories of inhomogeneous fluids, can be used to predict some of the nontrivial dynamical behaviors of fluid mixtures in confined environments.     

Prediction of gas condensate properties by Esmaeilzadeh–Roshanfekr equation of state

The Esmaeilzadeh–Roshanfekr (ER) equation of state (EOS) is used to predict the PVT properties of gas condensate reservoir fluids. Three gas condensate fluid samples taken from three wells in a real field in Iran, referred here as SA1, SA4 and SA8, as well as five samples from literature have been used to check the validity of the ER EOS in calculating the PVT properties of gas condensate mixtures. Some experiments such as constant composition expansion (CCE), constant volume depletion (CVD) and dew point pressures are carried out on these samples. In order to have an unbiased comparison between the ER and the Peng–Robinson (PR) equation of state, van der Waals mixing rules are used without using any adjustable parameters (k_ij = 0). Also, no pure component parameters are adjusted. The critical properties and acentric factor for plus-fraction are estimated by the Kesler–Lee, Pedersen et al. and Riazi–Daubert characterization methods. The results of dew point pressure calculations show that the ER EOS has smaller error than the PR EOS. For some mixtures, relative volume, gas compressibility factor and condensate drop-out in CVD and CCE test were also predicted. Comparison results between experimental and calculated data indicate that the ER EOS has smaller error than the PR EOS. The total average absolute deviation was found to be 0.82% and 2.97% for calculating gas compressibility factor and gas specific gravity in CVD test. Also, the total average absolute deviation was found to be 2.06% and 3.42% for calculating gas compressibility factor and relative volume in CCE test.     

Equation of state for square-well chain molecules with variable range. I: Application for pure substances

An equation of state (EOS) for square-well chain molecules with variable range developed on the basis of statistical mechanics for chemical association in our previous work is employed for the calculations of pVT properties and vapor–liquid equilibria (VLE) of pure non-associating fluids. The molecular parameters for 73 normal substances and 46 polymers are obtained from saturated vapor pressure and liquid molar volume data for normal fluids or pVT data for polymers. Linear relations are found for the molecular parameters of normal fluids with their molecular weight of homologous compounds. This indicates that the model parameters of homologous series, subsequently pVT and VLE, can be predicted when experimental data are not available. The predicted saturated vapor pressures and/or liquid volumes are satisfactory through the generalized model parameters. The calculated VLE and pVT for normal fluids and polymers by this EOS are compared with those from other engineering models, respectively.     

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.     

基于化学缔合统计理论的链状流体状态方程

基于化学缔合统计理论的链状流体状态方程(EOS)能够反映实际分子的形状、链节成链、缔合等具体信息,在实际流体热力学性质计算中有着广泛应用.一般的链状流体EOS仅考虑相邻链节间的相关性,我们则借助统计力学和计算机模拟结果在模型中纳入了相间链节间的相关性,获得的硬球链流体(HSCF)模型能够更好地预测模型流体的压缩因子和第二维里系数.以HSCF为参考,引入方阱色散微扰项获得了实际方阱链流体(SWCF)EOS;结合根据黏滞球模型导得的缔合项,进一步构建了缔合流体EOS.最近,我们根据微扰理论和积分方程方法又开发了一新的变阱宽方阱链流体(SWCF-VR)模型.SWCF和SWCF-VREOSs可很好地用于计算小分子、聚合物、离子液体等纯流体及混合物的相行为、热焓、表面张力、黏度等热力学及传递性质,显示了模型良好的工程应用价值.本文就本课题组多年来在自由空间范畴内基于化学缔合统计理论开发链状流体EOS及其实际应用作系统的总结.     

Modeling of Osmotic Pressure of Aqueous Poly(Ethylene Glycol) Solutions Using the Artificial Neural Network and Free Volume Flory Huggins Model

In this work, the modified Flory-Huggins coupled with the free-volume concept and the artificial neural network models were used to obtain the osmotic pressure of aqueous poly(ethylene glycol) solutions. In the artificial neural network, the osmotic pressure of aqueous poly(ethylene glycol) solutions depends on temperature, molecular weight and the mole fractions of poly(ethylene glycol) in aqueous solution. The network topology is optimized and the (3-1-1) architecture is found using optimization of an objective function with batch back propagation (BBP) method for 134 experimental data points. The results obtained from the neural network in obtaining of the osmotic pressure of aqueous poly(ethylene glycol) were compared with those obtained from the free volume Flory-Huggins model (FV-FH). The results showed that the modified Flory-Huggins model and also the artificial neural network can accurately predict the osmotic pressure of aqueous poly(ethylene glycol) solutions but the accuracy of ANN is much better than the modified Flory-Huggins model.     

Analysis of experimental data within the statistical theory of critical phenomena

An approach devised earlier is generalized for a binodal line and used to analyze experiments on a broad set of fluids from He⁴ to organic compounds. Experimental binodal lines and isothermal compressibility are described with precision for the above substances. The error characterizing the theory’s deviation from the experimental data is several times lower than in a similar analysis using the Ising model.     

Semi-empirical insertion probability for hard-spheres and hard-sphere chains

In spite of its simplicity and a well-defined theoretical basis, the Flory–Guggenheim approach is conventionally regarded as inapplicable to off-lattice system since the insertion probability of the approach does not account for the excluded region, existing in the off-lattice system. In this work, we propose the insertion probability accounting for the excluded region of off-lattice fluids and derive a new version of equation of state (EOS) for hard-sphere chains basing on the Flory–Guggenheim approach. To investigate the behavior of the excluded regions, a Monte Carlo sampling was performed for hard disks and the various excluded regions were found to have different density dependence. On the basis of the simulation result, we formulated the insertion probability for hard-sphere and that of hard-sphere chain which accounts for the effect of chain-connectivity on the monomer insertion. The proposed insertion probability was found to correctly predict the simulation data for monomer and correctly correlate the simulation data for chain fluids. The resulting EOS was found to meet closed-packed limit and predict the simulation data of compressibility factor for monomer and chains with a reasonable degree of accuracy. When compared with other off-lattice based EOS, it shows a comparable or better result. For second virial coefficient of chain molecules, the model was found to reasonably predict the simulation data.