Perturbed chain-statistical associating fluid theory extended to dipolar and quadrupolar molecular fluids

The perturbed chain statistical associating fluid theory (PC-SAFT) is extended to polar molecular fluids, namely dipolar and quadrupolar fluids. The extension is based on the perturbation theory for polar fluids by Stell and co-workers. Appropriate expressions are proposed for dipole-dipole, quadrupole-quadrupole, and dipole-quadrupole interactions. Furthermore, induced dipole interactions are calculated explicitly in the model. The new polar PC-SAFT model is relatively complex; for this purpose, a truncated polar PC-SAFT model is proposed using only the leading term in the polynomial expansion for polar interactions. The new model is used for the calculation of thermodynamic properties of various quadrupolar pure fluids. In all cases, the agreement between experimental data and model predictions is very good.     

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 crossover lattice fluid equation of state for pure fluids

A lattice fluid model is one of the most versatile, molecular-based engineering equations of state (EOS) but, in common with all analytic equations of state, the lattice fluid (LF) EOS exhibits classical behaviour in the critical region rather than the non-analytical, singular behaviour seen in real fluids. In this research, we use the LF EOS and develop a crossover lattice fluid (xLF) equation of state near to and far from the critical region which incorporates the scaling laws valid asymptotically close to the critical point while reducing to the original classical LF EOS far from the critical point. We show that, over a wide range of states, the xLF EOS yields the saturated vapour pressure data and the density data with much better accuracy than the classical LF EOS.     

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.     

Vapor-liquid transitions of dipolar fluids in disordered porous media: performance of angle-averaged potentials

Using replica integral equations in the reference hypernetted-chain (RHNC) approximation we calculate vapor-liquid spinodals, chemical potentials, and compressibilities of fluids with angle-averaged dipolar interactions adsorbed to various disordered porous media. Comparison with previous RHNC results for systems with true angle-dependent Stockmayer (dipolar plus Lennard-Jones) interactions indicate that, for a dilute hard sphere matrix, the angle-averaged fluid-fluid (ff) potential is a reasonable alternative for reduced fluid dipole moments m( *2)=mu(2)/(epsilon(0)sigma(3))< or =2.0. This range is comparable to that estimated in bulk fluids, for which RHNC results are presented as well. Finally, results for weakly polar matrices suggest that angle-averaged fluid-matrix (fm) interactions can reproduce main features observed for true dipolar (fm) interactions such as the shift of the vapor-liquid spinodals towards lower temperatures and higher densities. However, the effective attraction induced by dipolar (fm) interaction is underestimated rather than overestimated as in the case of angle-averaged ff interactions.     

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.     

由L-J流体的FMSA状态方程计算流体的PVT性质

运用Tang等提出的Lennard-Jones (L-J)流体两参数的一阶平均球形近似(FMSA)状态方程, 计算了流体的汽液共存相图和饱和蒸汽压曲线, 以及非饱和区的PVT性质, 并与文献数据进行比较. L-J参数由Tr＜0.95的汽液相共存数据回归得到. 计算结果表明, 对于分子较接近球形的流体, 除临界点附近外, 该方程可以在较大的温度和压力范围内计算真实流体的PVT性质, 结果满意. 对于球形分子, 该方程的精确度随分子尺寸的变大基本保持稳定. 该方程不适用于强极性物质. 在高密度区, 该方程的计算结果明显优于P-R方程. 对于分子偏离球形较远的流体, 该方程的适用性变差, 此时要考虑分子形状的影响, 可采用三参数的FMSA状态方程进行计算.     

Consideration of dipole–quadrupole interactions in molecular based equations of state

The thermodynamic behaviour of a number of real substances is determined by dipolar as well as quadrupolar interactions of the molecules. In equations of state (EOS) like, e.g. BACKONE separate contributions to the Helmholtz energy for the dipolar and the quadrupolar interactions are considered but no cross contributions. Here, the concept of effective dipole and quadrupole contributions is suggested in which the effective dipole strength μ_e is influenced by the quadrupole cross interaction. Similarily, the effective quadrupole strength Q_e takes into account the dipole cross interaction. In order to arrive at these effective dipolar and quadrupolar strengths, molecular simulations are performed. From the simulation results correlation equations are derived which are used in combination with BACKONE for the calculation of vapour–liquid equilibria (VLE) of real mixtures. By using these effective moments, the only required binary mixing rule parameter k_ij tends to small values of about 0.01 and becomes temperature-independent. Moreover, the VLE pressures are predicted now considerably better than without consideration of the cross contributions.     

Crossover SAFT equation of state for pure supercritical fluids

A crossover statistical associating fluid theory (SAFT) equation of state (EOS) is used to fit the parameters of eight common pure supercritical fluids (water, ammonia, carbon dioxide, R134a, ethane, propane, ethene and propene) and calculate their thermodynamic properties. Over a wide range including the critical region, the EOS reproduces the saturated pressure data with an average absolute deviation (AAD) of about 1% and the saturated densities with an AAD of about 2%. In the one-phase region, the EOS represents the experimental values of pressure with an AAD of about 1–3%. The results are satisfactory.     

Heat capacities of Stockmayer fluids from Monte Carlo simulations and perturbation theory

The heat capacities of dipolar fluids are investigated using a thermodynamic perturbation theory approach and the NVT and NpT Monte Carlo simulation methods. The theoretical results are compared to corresponding simulation data. The comparison shows that the applied perturbation theory is appropriate for the heat capacity calculations. As an application, the isobaric heat capacity of ammonia is also studied by the Stockmayer fluid model.     

Anomalous corresponding-states surface tension of hydrogen fluoride and of the Onsager model

In a corresponding-states analysis of the liquid-vapor surface tension originally suggested by Guggenheim, we study the behavior of different simple (i.e., nonpolar), polar and ionic fluids. The results are compared to the corresponding ones for model fluids of each of the three types. For simple and weakly polar fluids (both real and model), the data map onto a master curve, as demonstrated by Guggenheim. For strongly dipolar, associating fluids, which also exhibit hydrogen bonding, one finds deviations from the master curve at low temperatures and, thus, observes the characteristic sigmoid behavior of the reduced surface tension as a function of temperature. The same is obtained for the model ionic fluid, the restricted primitive model. Truly exceptionally low values of the reduced surface tension are found for hydrogen fluoride and for the Onsager model of dipolar fluids, the surface tension of which we evaluate using an approximate hypernetted chain relation to obtain the square-gradient term in a modified van der Waals theory. Remarkably, in the corresponding-states plot, the surface tensions of HF and of the Onsager model agree very closely, while being well separated from the values for the other fluids. We also study the gradual transition of a model fluid from a simple fluid to a strongly dipolar one by varying the relative strength of dipolar and dispersion forces.     

Application of the crossover lattice equation of state for fluid mixtures

In previous work, we developed the crossover lattice equation of state (xLF EOS) for pure fluids and the xLF EOS yielded the saturated vapour pressure and the density values with a much better accuracy than the classical LF EOS over a wide range. In this work, we extended xLF EOS to fluid mixtures. Classical composition-dependent mixing rules with only adjustable two binary interaction parameters same as the LF EOS are used. A comparison is made upon experimental data for fluids mixtures in the one- and two-phase regions. The xLF EOS shows more improved representations than the LF EOS, especially in the critical region.     

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.     

Electric polarization effects on the electronic spectral shift of centrosymmetric compounds

The classic dielectric dipolar Onsager model was extended to include quadrupolar interactions between solute molecules and solvents with different polarities. A multiparametric solvatochromic expression, based on the point quadrupole moment inside a spherical cavity embedded in a dielectric continuum, is applied to centrosymmetric sulfonamide porphyrins, zinc tetraphenyl porphyrin, squaraine and 9,10-dicyanoanthracene, in order to account for the quadrupolar polarization effect of solute molecules. The reaction field polarity functions created respectively by dipole and quadrupole moments are compared and found to be linearly correlated.     

Efficient prediction of thermodynamic properties of quadrupolar fluids from simulation of a coarse-grained model: the case of carbon dioxide

Monte Carlo simulations are presented for a coarse-grained model of real quadrupolar fluids. Molecules are represented by particles interacting with Lennard-Jones forces plus the thermally averaged quadrupole-quadrupole interaction. The properties discussed include the vapor-liquid coexistence curve, the vapor pressure along coexistence, and the surface tension. The full isotherms are also accessible over a wide range of temperatures and densities. It is shown that the critical parameters (critical temperature, density, and pressure) depend almost linearly on a quadrupolar parameter q=Q(*4)T*, where Q* is the reduced quadrupole moment of the molecule and T* the reduced temperature. The model can be applied to a variety of small quadrupolar molecules. We focus on carbon dioxide as a test case, but consider nitrogen and benzene, too. Experimental critical temperature, density, and quadrupolar moment are sufficient to fix the parameters of the model. The resulting agreement with experiments is excellent and marks a significant improvement over approaches which neglect quadrupolar effects. The same coarse-grained model was also applied in the framework of perturbation theory in the mean spherical approximation. As expected, the latter deviates from the Monte Carlo results in the critical region, but is reasonably accurate at lower temperatures.     

Description of self-diffusion coefficients of gases,liquids and fluids at high pressure based on molecular simulation data

The purpose of this work is to evaluate the potential of modeling the self-diffusion coefficient (SDC) of real fluids in all fluid states based on Lennard–Jones analytical relationships involving the SDC, the temperature, the density and the pressure. For that, we generated an equation of state (EOS) that interrelates the self-diffusion coefficient, the temperature and the density of the Lennard–Jones (LJ) fluid. We fit the parameters of such LJ–SDC–EOS using recent wide ranging molecular simulation data for the LJ fluid. We also used in this work a LJ pressure–density–temperature EOS that we combined with the LJ–SDC–EOS to make possible the calculation of LJ–SDC values from given temperature and pressure. Both EOSs are written in terms of LJ dimensionless variables, which are defined in terms of the LJ parameters ɛ and σ. These parameters are meaningful at molecular level. By combining both EOSs, we generated LJ corresponding states charts which make possible to conclude that the LJ fluid captures the observed behavioral patterns of the self-diffusion coefficient of real fluids over a wide range of conditions. In this work, we also performed predictions of the SDC of real fluids in all fluid states. For that, we assumed that a given real fluid behaves as a Lennard–Jones fluid which exactly matches the experimental critical temperature T_c and the experimental critical pressure P_c of the real fluid. Such an assumption implies average true prediction errors of the order of 10% for vapors, light supercritical fluids, some dense supercritical fluids and some liquids. These results make possible to conclude that it is worthwhile to use the LJ fluid reference as a basis to model the self-diffusion coefficient of real fluids, over a wide range of conditions, without resorting to non-LJ correlations for the density–temperature–pressure relationship. The database considered here contains more than 1000 experimental data points. The database practical reduced temperature range is from 0.53 to 2.4, and the practical reduced pressure range is from 0 to 68.4.     

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.