The van der Waals one-fluid model for viscosity in Lennard–Jones fluids: Influence of size and energy parameters

This work aims to estimate the limitations of the van der Waals one-fluid (vdW1) approximation in the prediction of the viscosity of Lennard–Jones (LJ) mixtures. To do so, non-equilibrium molecular dynamics simulations have been performed. Results on mixtures have been compared to those of their equivalent pure fluids (in the sense of the vdW1 model). Several systems (146 configurations) are studied, which are composed of binary and ternary mixtures in various thermodynamic states and for different combining rules. In a first step, deviations induced separately by LJ molecular parameters (size or energy) have been analyzed. It is shown that the vdW1 model is well designed for the energy parameter in every configuration. On the contrary, for the size parameter, deviations induced by this one-fluid approach are shown to be large in dense state and low temperature systems. In a second step, the coupling effects of the LJ size and energy parameters with the mass are studied. It appears that an accurate one-fluid approximation for viscosity should involve a coupling between the mass and the size parameters in its formulation (which is not the case with the vdW1 model) but not between the mass and the energy.     

Comparison of different mixing rules for prediction of density and residual internal energy of binary and ternary Lennard–Jones mixtures

Mixing rules are necessary when equations of state for pure fluids are used to calculate various thermodynamic properties of fluid mixtures. The well-known van der Waals one-fluid (vdW1) mixing rules are proved to be good ones and widely used in different equations of state. But vdW1 mixing rules are valid only when molecular size differences of components in a mixture are not very large. The vdW1 type density-dependent mixing rule proposed by Chen et al. [1] is superior for the prediction of pressure and vapor–liquid equilibria when components in the mixture have very different sizes. The extension of the mixing rule to chain-like molecules and heterosegment molecules was also made with good results. In this paper, the comparison of different mixing rules are carried out further for the prediction of the density and the residual internal energy for binary and ternary Lennard–Jones (LJ) mixtures with different molecular sizes and different molecular interaction energy parameters. The results show that the significant improvement for the prediction of densities is achieved with the new mixing rule [1], and that the modification of the mixing rule for the interaction energy parameter is also necessary for better prediction of the residual internal energy.     

Vapor–liquid equilibria predictions of hydrogen–hydrocarbon mixtures with the Huron–Vidal mixing rule

An excess Gibbs-equation of state (G^E-EoS) framework based on the Huron–Vidal mixing rule, has been applied to study vapor–liquid equilibria (VLE) of hydrogen–hydrocarbon mixtures. The mixing rule couples the Peng–Robinson–Stryjek–Vera (PRSV) EoS with a local composition solution model. The solution model is based on one-fluid theory treatment and assigns a single energy parameter to each binary pair. This energy parameter relates to the preference of the molecules for like to unlike interactions. The allocation of a system's number of interactions to the individual species in a binary mixture, incorporates the use of size parameters which gain significance only in the liquid phase. In a two parameter form, the framework has been used for the simultaneous data reduction of a large number of binary and several ternary hydrogen–hydrocarbon mixtures. These systems were taken over an extended range of pressures and temperatures. Results from the data reduction are reported in both tabular and graphical forms. Correlations for the model parameters have been identified with the acentric factor of the hydrocarbon in hydrogen–hydrocarbon binary mixtures. In a fully predictive mode, the model has shown to describe well VLE of binary hydrogen–linear alkane systems. Comparisons of these results with calculations from the Peng–Robinson (PR) EoS and the classical mixing rule (vdW) are included.     

A new three-parameter cubic equation of state for calculation physical properties and vapor–liquid equilibria

A new three-parameter cubic equation of state is presented by combination of a modified attractive term and van der Waals repulsive expression. Also a new alpha function for the attractive parameter of the new EOS is proposed. The new coefficients of alpha function and the other parameters of the attractive term are adjusted using the data of the saturated vapor pressure and liquid density of almost 60 pure compounds including heavy hydrocarbons. The new EOS is adopted for prediction of the various thermophysical properties of pure compounds such as saturated and supercritical volume, enthalpy of vaporization, compressibility factor, heat capacity and sound velocity. Following successful application of the new EOS for the pure components, using vdW one-fluid mixing rules, the new EOSs are applied to prediction of the bubble pressure and vapor mole fraction of the several binary and ternary mixtures. The accuracy of the new EOS for phase equilibrium calculation is demonstrated by comparison of the results of the present EOSs with the PT, PR, GPR and SRK cubic EOSs.     

Simulation studies on mixtures of dipolar with nonpolar linear molecules II. A mixing rule for the dipolar contribution to the Helmholtz energy

In recent equations of state an explicit expression for the dipolar contribution to the Helmholtz energy F_D is used. If such equations are applied to mixtures the problem of mixing rules for F_D arises. As a solution, a “one-fluid dipole moment μ_x” is introduced which is based on perturbation theory. The accuracy of the approach is tested by simulation results for mixtures of dipolar with nonpolar linear molecules. In addition to the simulation runs from Part I (Müller et al. (1994)), results from 97 new NVT simulations are presented. Comparisons of pressures, internal energies, dipole energies and excess free energies from simulations with those obtained from the equation for F_D in combination with the one-fluid dipole moment μ_x shows good to excellent agreement.     

Equations of state for strongly nonideal fluid mixtures: Application of local compositions toward density-dependent mixing rules

A local-composition, two-fluid model has been developed for equation-of-state calculations of fluid-phase equilibria for asymmetric mixtures; it is applicable to any equation of state of the van der Waals form. A modification of the quasichemical theory of Guggenheim is applied to mixtures at all fluid densities. Desirable boundary conditions are met at low densities, at high densities, and at high temperatures.In effect, the local-composition model uses density-dependent mixing rules. It contains no new adjustable binary parameters and can be extended to multicomponent mixtures without ternary (or higher) parameters. It appears that, when compared to conventional one-fluid models, significant improvement may be obtained in predictions for vapor-liquid equilibria of typical asymmetric mixtures.     

Vapor-liquid equilibria for alcohol/hydrocarbon systems using the CPA Equation of state

The recently developed Cubic-Plus-Association Equation of State (CPA EoS) is extended in this study to binary systems containing one associating compound (alcohol) and an inert one (hydrocarbon). CPA combines the Soave-Redlich-Kwong (SRK) equation of state for the physical part with an association term based on perturbation theory. The classical van der Waals one-fluid mixing rules are used for the attractive and co-volume parameters, and b, while the extension of the association term to mixtures is rigorous and does not require any mixing rules. Excellent correlation of Vapor-Liquid Equilibria (VLE) is obtained using a small value for the interaction parameter (k_ij) in the attractive term of the physical part of the equation of state even when it is temperature-independent. CPA yileds much better results than SRK and its performance is similar to that of other association models, like the Anderko EoS, and the more complex SAFT and Simplified SAFT EoS.     

三元——假二元系的过量体积

提出测量三元——假二元系(1-x)[(1-y)A+yC]+z[(1-y)B+yC]的过量热力学函数的方法, 测量了在298.15 K下两种假二元系: A=环己烷, B=苯, C=甲苯, 甲基环己烷的过量体积V~E, 用van der Waals单流体模型对实验结果进行了讨论。     

Ionic force field optimization based on single-ion and ion-pair solvation properties: going beyond standard mixing rules

Using molecular dynamics (MD) simulations in conjunction with the SPC/E water model, we optimize ionic force-field parameters for seven different halide and alkali ions, considering a total of eight ion-pairs. Our strategy is based on simultaneous optimizing single-ion and ion-pair properties, i.e., we first fix ion-water parameters based on single-ion solvation free energies, and in a second step determine the cation-anion interaction parameters (traditionally given by mixing or combination rules) based on the Kirkwood-Buff theory without modification of the ion-water interaction parameters. In doing so, we have introduced scaling factors for the cation-anion Lennard-Jones (LJ) interaction that quantify deviations from the standard mixing rules. For the rather size-symmetric salt solutions involving bromide and chloride ions, the standard mixing rules work fine. On the other hand, for the iodide and fluoride solutions, corresponding to the largest and smallest anion considered in this work, a rescaling of the mixing rules was necessary. For iodide, the experimental activities suggest more tightly bound ion pairing than given by the standard mixing rules, which is achieved in simulations by reducing the scaling factor of the cation-anion LJ energy. For fluoride, the situation is different and the simulations show too large attraction between fluoride and cations when compared with experimental data. For NaF, the situation can be rectified by increasing the cation-anion LJ energy. For KF, it proves necessary to increase the effective cation-anion Lennard-Jones diameter. The optimization strategy outlined in this work can be easily adapted to different kinds of ions.     

Molecular dynamics simulation of potassium along the liquid-vapor coexistence curve

The applicability of pair potential functions to liquid alkali metals is questionable. On the one hand, some recent reports in the literature suggest the validity of two-parameter pair-wise additive Lennard-Jones (LJ) potentials for liquid alkali metals. On the other hand, there are some reports suggesting the inaccuracy of pair potential functions for liquid metals. In this work, we have performed extensive molecular dynamics simulations of vapor-liquid phase equilibria in potassium to check the validity of the proposed LJ potentials and to improve their accuracy by changing the LJ exponents and taking into account the temperaturedependencies of the potential parameters. We have calculated the orthobaric liquid and vapor densities of potassium using LJ (12–6), LJ (8.5–4) and LJ (5–4), effective pair potential energy functions. The results show that using an LJ (8.5–4) potential energy function with temperature-independent parameters, ε and σ, is inadequate to account for the vapor-liquid coexistence properties of potassium. Taking into account the temperature-dependencies of the LJ parameters, ε(T) and σ(T), we obtained the densities of coexisting liquid and vapor potassium in a much better agreement with experimental data. Changing the magnitude of repulsive and attractive contributions to the potential energy function shows that a two-parameter LJ (5–4) potential can well reproduce the densities of liquid and vapor potassium. The results show that LJ (5–4) potential with temperature-dependent parameters produces the densities of liquid and vapor potassium more accurately, compared to the results obtained using LJ (12–6) and LJ (8.5–4) potential energy functions.     

A molecular theory for the thermodynamic behavior of polar mixtures. II. Development of an equation of state based on the local composition mixing rules

In this study, the essential features of a molecular theory developed earlier for the local composition model in solution thermodynamics is used as the basis for more applied calculations of vapor-liquid equilibria for mixtures of molecules vastly different in size, polarity, and strength of interaction. An accurate equation of state is introduced into the method by incorporating the Helmholtz free energy through the Gibbs-Helmholtz relation. In the local composition mixing rules, the interaction energy effects are represented by a multifluid model, while molecular size effects are represented by a one-fluid model, which in spirit corresponds to a mean density approximation for the molecular pair distribution functions. Calculations of the vapor-liquid equilibria of a wide variety of binary mixtures including nonpolar hydrocarbons, hydrogen-bonding alcohols, water, ammonia , and carbon dioxide show good agreement with experimental data.     

Mixing rules for hardsphere chain mixtures and their extension to heteronuclear hardsphere polyatomic fluids and mixtures

Mixing rules are very important for the calculation of fluid properties using different equations of state. In order to find the theoretical lead of the mixing rule for the size parameter, a mixing rule [1] for hardsphere mixtures has been proposed on the basis of Carnahan-Starling equation and Boublik-Mansoori equation. As its extension, mixing rules for hardsphere chain mixtures are proposed in this work. A mixing rule for the segment number (or chain length) is derived on the limitation of the equality of segment diameters, from the first order thermodynamic perturbation theories (TPT1) for pure chain fluids and for chain mixtures. Meanwhile, the mixing rule for the segment diameter is the same as the mixing rule for hardsphere mixtures on the limitation of monomer mixtures. The two mixing rules are checked together over wide ranges of conditions for hardsphere chain mixtures and compared with the first order thermodynamic perturbation theory (TPT1) and also with simulation data available in literature. An another interesting usage of new mixing rules is to describe the heteronuclear hardsphere polyatomic pure fluids, which consist of hardspheres with different segment diameters as in methane and ethane in which carbon and hydrogen atoms are looked as bonded spheres, and heteronuclear hardsphere polyatomic mixtures. The comparison with simulation data shows the validity of the mixing rules.     

Solubility of carbon dioxide in aqueous solutions of methanol. Predictions by molecular simulation and comparison with experimental data

The solubility of carbon dioxide in pure methanol, and in aqueous solutions of methanol, was computed using the Gibbs ensemble Monte Carlo (GEMC) technique for 313, 354, and 395 K at pressures up to 9 MPa. Three solvent mixtures (of methanol and water) with methanol mole fractions of 10, 50, and 75 mole percent (in the gas-free solvent mixture) were studied. The Monte Carlo simulations were conducted in an isothermal-isobaric ensemble applying effective pair potentials for the pure components from literature. Common mixing rules without any adjustable binary interaction parameters were used to describe the interactions between the mixture components. Overall, a good agreement between simulation results and recently published experimental data is achieved.     

Application of a new equation of state to liquid refrigerant mixtures

A new general equation of state recently reported for pure liquids has been developed to predict the volumetric and thermodynamic properties of six binary and two ternary liquid refrigerant mixtures (including HCs and HFCs mixtures) at different temperatures, pressures, and compositions. The results show this equation of state can be used to reproduce and predict different thermodynamic properties of liquid refrigerant mixtures within experimental errors. The composition dependence of the parameters of this equation of state has been assumed as quadratic functions of mole fraction. Using these mixing rules, the agreement between calculated and experimental densities is better than 0.6% for binary mixtures and 2.3% for ternary mixtures. To compare the performance of this new equation of state against other well-known methods such as the COSTALD method, the density of some refrigerant mixtures, for which the parameters of COSTALD were available, has been computed and compared with those of this new equation of state.     

New parametrization method for dissipative particle dynamics

We introduce an improved method of parametrizing the Groot-Warren version of dissipative particle dynamics (DPD) by exploiting a correspondence between DPD and Scatchard-Hildebrand regular solution theory. The new parametrization scheme widens the realm of applicability of DPD by first removing the restriction of equal repulsive interactions between like beads, and second, by relating all conservative interactions between beads directly to cohesive energy densities. We establish the correspondence by deriving an expression for the Helmoltz free energy of mixing, obtaining a heat of mixing which is exactly the same form as that for a regular mixture (quadratic in the volume fraction) and an entropy of mixing which reduces to the ideal entropy of mixing for equal molar volumes. We equate the conservative interaction parameters in the DPD force law to the cohesive energy densities of the pure fluids, providing an alternative method of calculating the self-interaction parameters as well as a route to the cross interaction parameter. We validate the new parametrization by modeling the binary system SnI(4)SiCl(4), which displays liquid-liquid coexistence below an upper critical solution temperature around 140 degrees C. A series of DPD simulations were conducted at a set of temperatures ranging from 0 degrees C to above the experimental upper critical solution temperature using conservative parameters based on extrapolated experimental data. These simulations can be regarded as being equivalent to a quench from a high temperature to a lower one at constant volume. Our simulations recover the expected phase behavior ranging from solid-liquid coexistence to liquid-liquid coexistence and eventually leading to a homogeneous single phase system. The results yield a binodal curve in close agreement with the one predicted using regular solution theory, but, significantly, in closer agreement with actual solubility measurements.     

Modification of PSRK mixing rules and results for vapor–liquid equilibria, enthalpy of mixing and activity coefficients at infinite dilution

The predictive Soave–Redlich–Kwong (PSRK) equation of state (EOS) is a well-established method for the prediction of thermodynamic properties required in process simulation. But there are still some problems to be solved, e.g. the reliability for strong asymmetric mixtures of components which are very different in size. The following modifications are introduced in the PSRK mixing rules: the Flory–Huggins term in the mixing rule for the EOS parameter a, and the combinatorial part in the UNIFAC model are skipped simultaneously; a nonlinear mixing rule for the EOS parameterb, instead of the linear mixing rule, is proposed. With these two modifications better results are obtained for vapor–liquid equilibria and activity coefficients at infinite dilution for alkane–alkane systems, specially for asymmetric systems. In order to obtain better results for enthalpy of mixing, temperature-dependent parameters are used. Group interaction parameters have been fitted for several groups, and the results are compared with the Modified UNIFAC (Dortmund), and the PSRK methods.     

Equation of state for square-well chain molecules with variable range II. Extension to mixtures

An equation of state (EOS) developed in our previous work for square-well chain molecules with variable range is further extended to the mixtures of non-associating fluids. The volumetric properties of binary mixtures for small molecules as well as polymer blends can well be predicted without using adjustable parameter. With one temperature-independent binary interaction parameter, satisfactory correlations for experimental vapor–liquid equilibria (VLE) data of binary normal fluid mixtures at low and elevated pressures are obtained. In addition, VLE of n-alkane mixtures and nitrogen + n-alkane mixtures at high pressures are well predicted using this EOS. The phase behavior calculations on polymer mixture solutions are also investigated using one-fluid mixing rule. The equilibrium pressure and solubility of gas in polymer are evaluated with a single adjustable parameter and good results are obtained. The calculated results for gas + polymer systems are compared with those from other equations of state.     

Acoustic,viscometric and optical properties of binary mixtures of tetrahydrofuran with 1-propanol and 2-propanol

The values estimated from various mixing rules for the ultrasonic velocity, viscosity and refractive index have been compared with the respective values measured earlier at 293, 303, and 313?K over the entire mole fraction range of two binary mixtures of tetrahydrofuran (THF) with 1-propanol (1-p) and 2-propanol (2-p). There is an excellent agreement between the experimental values of ultrasonic velocity and of refractive index with the respective values obtained from the mixing rules. The mixing rules for viscosity provide values agreeing broadly with those obtained from experimental measurements. The relative merits and interrelations of these mixing rules are discussed.     

Fundamental measure density functional theory studies on the freezing of binary hard-sphere and Lennard-Jones mixtures

Free energies and correlation functions of liquid and solid hard-sphere (HS) mixtures are calculated using the fundamental measure density functional theory. Using the thermodynamic perturbation theory the free energies of solid and liquid Lennard-Jones (LJ) mixtures are obtained from correlation functions of HS systems within a single theoretical approach. The resulting azeotrope- and spindle-type solid-liquid phase diagrams of HS and LJ binary mixtures are in good agreement with the corresponding ones from computer simulations.