Freezing of hard spheres confined in narrow cylindrical pores

Monte Carlo simulations for the equation of state and phase behavior of hard spheres confined inside very narrow hard tubes are presented. For pores whose radii are greater than 1.1 hard sphere diameters, a sudden change in the density and the microscopic structure of the fluid is neatly observed, indicating the onset of freezing. In the high-density structure the particles rearrange in such a way that groups of three particles fit in sections across the pore.     

Equation of state for hard-spheres and excess function calculations of binary liquid mixtures

Using our previously proposed matrix method, an equation of state for hard spheres is presented, which can reproduce the exact values of the first-eight virial coefficients. This equation meets both the low density and the close-packed limits and can predicts the first order fluid–solid phase transition of hard spheres. The results obtained show that the new equation of state can correlate the simulation data of compressibility factor up to high densities better than other equations of state.The new equation of state is extended to mixtures of hard spheres and excess functions of various binary liquid mixtures are calculated using the perturbation theory of Leonard–Henderson–Barker. The results are compared with existing theoretical and experimental data and with those calculated by other hard-sphere equations of state.It is seen that the results obtained by the new equation of state is quite satisfactory compared to other equations of state for the hard spheres and mixture of hard spheres.     

Alternative fundamental measure theory for additive hard sphere mixtures

The purpose of this short paper is to present an alternative fundamental measure theory (FMT) for hard sphere mixtures. Keeping the main features of the original Rosenfeld's FMT [Phys. Rev. Lett. 63, 980 (1989)] and using the dimensional and the low-density limit conditions a new functional is derived incorporating Boublik's multicomponent extension [Mol. Phys. 59, 371 (1986)] of highly accurate Kolafa's equation of state for pure hard spheres. We test the theory for pure hard spheres and hard sphere mixtures near a planar hard wall and compare the results with the original Rosenfeld's FMT and one of its modifications and with new very accurate simulation data. The test reveals an excellent agreement between the results based on the alternative FMT and simulation data for density profile near a contact and some improvement over the original Rosenfeld's FMT and its modification at the contact region.     

Thermodynamics of dipolar hard sphere mixtures

The thermodynamic properties of mixtures of hard spheres with imbedded point dipoles are investigated by an extension of the perturbation theory used by Rushbrooke et al. for pure fluids. Equations are presented for the general multicomponent mixture in which all the hard spheres have the same diameter. Numerical calculations are presented of the phase behaviour and excess thermodynamic mixing functions for the special case of the binary mixture in which only one species is polar. A brief discussion is given of the relationship of this model to experimental results for real fluid mixtures and of possible extensions of this work.     

Application of a new simplified SAFT to VLE study of associating and non-associating fluids

A new equation of state based on the Statistical Associating Fluid Theory (SAFT) is presented to study the phase behavior of associating and non-associating fluids. In the new equation of state, the hard sphere contribution to compressibility factor of the simplified version of the SAFT (SSAFT) is replaced with that proposed by Ghotbi and Vera. The Ghotbi–Vera SSAFT (GV-SSAFT) was also extended to study the phase behavior of associating and non-associating mixtures. The GV-SSAFT like the SSAFT equation of state has three adjustable segment parameters for non-associating fluids and five parameters for associating fluids. The experimental data of liquid densities and vapor pressures for pure fluids studied in this work were used to obtain the best values for the parameters of the GV-SSAFT. The results obtained from the GV-SSAFT for liquid densities and vapor pressures of pure associating and non-associating fluids were compared with those obtained from the SSAFT equation of state. The results showed that the GV-SSAFT similar to the SSAFT can accurately correlate the experimental data of liquid density and vapor pressure for systems studied. On the other hand the results obtained from two SAFT-based equations of state are almost identical. In order to show capability of the GV-SSAFT and SSAFT equations of state, they were used to directly calculate heat of vaporization for a number of pure associating and non-associating fluids. Slightly better results for heat of vaporization comparing to the experimental data were obtained from the GV-SSAFT EOS than those obtained from the SSAFT. The GV-SSAFT was also used to study the VLE phase behavior for a number of binary associating and non-associating mixtures. The results also showed that the GV-SSAFT can be successfully used to study the phase behavior of mixtures studied in this work.     

Surface tension correlation for pure polar fluids by a new molecular model and SRK equation of state

In this work, a new model based on molecular thermodynamic was presented to correlate the surface tension of pure polar liquids. This model was developed based on the Davis theory. According to this theory, the surface tension is defined as a function of radial distribution function (RDF) and potential function (PF) as well. The proposed model includes three additive terms; hard sphere, dispersion and polar interactions. The RDF of Kolafa equation of state and Dirac delta function as a PF were used for hard sphere interaction. The RDF expression of Xu and Hu was considered for both dispersion and polar interactions. The presented model has two adjustable parameters, size and energy, which were obtained by optimization of an objective function for each pure fluid. This proposed approach was used for 19 pure polar fluids divided into 6 groups; organic acids, alcohols, ketones, ethers, aldehydes, and water. The average absolute deviation percent (AAD%) obtained for 19 fluids are 0.74. Also the surface tension of these 19 fluids was calculated by the use of SRK EOS and Sugden empirical formula in two cases. In case 1, Sugden's Parachor was calculated from Hugill and van Welsenes correlation and in case 2, it was obtained by optimization of an objective function for each component. The values of AAD% are 43.544 and 2.281 for cases 1 and 2, respectively. These results show the new model, which includes two adjustable parameters, can correlate the surface tension of the pure polar liquids with a high accuracy.     

A study of the phase behavior of a simple model of chiral molecules and enantiomeric mixtures

We present a study of the solid-fluid and solid-solid phase equilibrium for molecular models representative of chiral molecules and enantiomeric mixtures. The models consist of four hard sphere interaction sites of different diameters in a tetrahedral arrangement with the fifth hard sphere interaction site at the center of the tetrahedron. The volumetric properties and free energies of the pure enantiomers and binary mixtures were calculated in both fluid and solid phases using isobaric Monte Carlo simulations. The models exhibit essentially ideal solution behavior in the fluid phase with little chiral discrimination. In the solid phase the effects of chirality are much greater. Solid-fluid phase behavior involving the pure enantiomer solids and also racemic compounds was calculated. The calculations indicate that, depending on the relative sizes of the hard sphere interaction sites, packing effects alone can be sufficient to stabilize a racemic compound with respect to the pure enantiomer solids.     

A novel equation of state for the prediction of thermodynamic properties of fluids

This work proposes a new equation of state (EOS) based on molecular theory for the prediction of thermodynamic properties of real fluids. The new EOS uses a novel repulsive term, which gives the correct hard sphere close packed limit and yields accurate values for hard sphere and hard chain virial coefficients. The pressure obtained from this repulsive term is corrected by a combination of van der Waals and Dieterici potentials. No empirical temperature functionality of the parameters has been introduced at this stage. The novel EOS predicts the experimental volumetric data of different compounds and their mixtures better than the successful EOS of Peng and Robinson. The prediction of vapor pressures is only slightly less accurate than the results obtained with the Peng-Robinson equation that is designed for these purposes. The theoretically based parameters of the new EOS make its predictions more reliable than those obtained from purely empirical forms.     

A three-parameter cubic equation of state for prediction of thermodynamic properties of fluids

A new cubic three-parameter equation of state has been proposed for PVT and VLE calculations of simple, high polar and associating fluids. The parameters are temperature dependent in sub-critical region, but temperature independent in super-critical region. The results for 42 simple and 14 associative pure compounds indicate that the calculated saturation properties and volumetric properties over the whole temperature range, up to high pressures, by the proposed equation of state (EOS), were in better agreement with the experimental data, compared with those obtained by the five well-known EOSs (P–R, P–T, Adachi et al., Yu–Lu, and M4). Two derivative properties, molar enthalpy and heat capacity of water and ammonia have been calculated, and demonstrated the thermodynamic consistency of the EOS parameters. Also VLE calculations have been performed for 41 binary mixtures of different type of fluids, including those of interest in petroleum industry. The results indicated the high capability of the proposed EOS for calculating the thermodynamic properties of pure and fluid mixtures.     

Group contribution multifluid lattice theory for simultaneous predictions of thermodynamic properties of pure fluids and mixtures

A unified group contribution (GC) lattice equation of state (EOS) was formulated based on the multifluid approximation of the nonrandom lattice fluid theory. The GC-EOS requires segment size and interaction energy parameter from functional group characteristics. The unique feature of the approach is that a single set of group parameters are used for both pure fluids and mixtures. The approach was found to be quantitatively applicable for predicting thermodynamic properties of real pure fluids and mixtures. Its potential utility was demonstrated for vapor pressures, vapor–liquid coexistence densities of pure fluids and phase equilibrium properties of mixtures including polymeric solutions.     

Thermodynamic calculation of supercritical-fluid equilibria: New mixing rules for equations of state

Simple cubic equations of state with conventional mixing rules have played an important role in the calculation of phase equilibria and other thermodynamic properties of non-polar fluid mixtures. In the application of supercritical fluids to separation processes, volumetric as well as phase equilibrium properties are very important for rational process design.

Heyen (1980) proposed a cubic equation of state which shows better accuracy in the calculation of volumetric properties, compared to the Peng-Robinson equation of state. In order to apply his equation to polar mixtures, Heyen recently proposed a density-independent mixing rule, but this does not obey the universally-observed quadratic mixing rule of the second virial coefficient in the low-density limit.

This paper proposes a new density-dependent mixing rule for the Heyen equation of state. The Heyen equation of state with our new mixing rule appears to calculate the phase equilibria and the volumetric properties of CO₂-containing non-polar as well as polar mixtures with good accuracy.     

On statistical mechanical equations of state for simple fluids : Effective hard spheres and quantum corrections

Nordholm, S., Greberg, H. and Penfold, R., 1991. On statistical mechanical equations of state for simple fluids. Effective hard spheres and quantum corrections. Fluid Phase Equilibria, 90: 307-332.

A comparison is made of the mean field generalised van der Waals theory, based on a variationally determined hard sphere diameter, with the recent equation of state proposed by Song and Mason incorporating a temperature-dependent hard sphere diameter and correlation effects through the second virial coefficient. The simple cell theory ansatz of the former is less accurate but permits a wide range of applications including the estimation of quantum effects on the bulk properties of light fluids at low temperatures. Results for the critical parameters of ³He, ⁴He, H₂, Ne, CH₄ and Ar are examined. The relevance to the corresponding theory of non-uniform fluids is noted.     

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.     

A perturbed hard-sphere-chain equation of state for nematic liquid crystals and their mixtures with polymers

A perturbed hard-sphere-chain (PHSC) equation of state is presented to compute nematicisotropic equilibria for thermotropic liquid crystals, including mixtures. The equation of state consists of an isotropic term and an anisotropic term given by the Maier-Saupe theory whose contribution disappears in the isotropic phase. The isotropic contribution is the recently presented PHSC equation of state for normal fluids and polymers which uses a reference equation of state for athermal hard-sphere chains and a perturbation theory for the squarewell fluid of variable well width. The PHSC equation of state gives excellent correlations of pure-component pressure-volume-temperature data in the isotropic region and, combined with the Maier-Saupe theory, correlates the dependence of nematic-isotropic transition temperature on the pressure. Theory also predicts a nematic-isotropic biphasic region and liquid-liquid phase separation in a temperature-composition diagram of binary mixtures containing a nematic liquid crystal and a normal fluid or polymer. Theory and experiment show good agreement for pure fluids as well as for mixtures.     

Primary electroviscous effect in a suspension of charged porous spheres

Primary electroviscous effect for a dilute suspension of porous spheres with fixed volumetric charge density is investigated theoretically. In the absence of flow, the electrical potential and solution charge density are assumed to satisfy the linearized Poisson-Boltzmann equation. With incorporation of the electrical body force, the Brinkman equation and the Stokes equation are used to govern the fluid flow inside and outside a sphere. The theory is formulated by assuming weak deviation of the charge cloud from its equilibrium state. However, the electrical body force is not restricted to be small compared to the viscous force in the fluid momentum equation. The results show that the double layer distortion is increased with increasing particle permeability, thereby enhancing the relative importance of its stress contribution. Nonetheless, the intrinsic viscosity remains a decreasing function of permeability, similar to the case of uncharged particles.     

The effect of rotational and translational energy exchange on tracer diffusion in rough hard sphere fluids

A study is presented of tracer diffusion in a rough hard sphere fluid. Unlike smooth hard spheres, collisions between rough hard spheres can exchange rotational and translational energy and momentum. It is expected that as tracer particles become larger, their diffusion constants will tend toward the Stokes-Einstein hydrodynamic result. It has already been shown that in this limit, smooth hard spheres adopt "slip" boundary conditions. The current results show that rough hard spheres adopt boundary conditions proportional to the degree of translational-rotational energy exchange. Spheres for which this exchange is the largest adopt "stick" boundary conditions while those with more intermediate exchange adopt values between the "slip" and "stick" limits. This dependence is found to be almost linear. As well, changes in the diffusion constants as a function of this exchange are examined and it is found that the dependence is stronger than that suggested by the low-density, Boltzmann result. Compared with smooth hard spheres, real molecules undergo inelastic collisions and have attractive wells. Rough hard spheres model the effect of inelasticity and show that even without the presence of attractive forces, the boundary conditions for large particles can deviate from "slip" and approach "stick."     

Macromolecule-Induced Clustering of Hard Spheres

The connectivity Ornstein-Zernike formalism, together with the polymer reference interaction site model (PRISM), is employed to describe connectivity and network formation in mixtures of spheres and polymers. Results are presented for the percolation of spheres induced by both flexible coil-like and rigid rod-like linear polymers; the Percus-Yevick (PY) approximation is used throughout. Our results are compared with predictions based on the adhesive hard sphere (AHS) model, and correlations with the polymer-mediated second virial coefficient between spheres are discussed. Copyright 2001 Academic Press.     

Theory of excluded volume equation of state: higher approximations and new generation of equations of state for entire density range

A novel theory of an equation of state based on excluded volume and formulated in two preceding papers for gases and gaseous mixtures is extended to the entire density range by considering higher (beginning from the third) approximations of the theory. The algorithm of constructing higher approximations is elaborated. Equations of state are deduced using the requirement of maximum simplicity and contain a single free parameter to be chosen by reason of convenience or simplicity or to be used as a fitting parameter with respect to the computer simulation database. In this way, precise equations of state are derived for the hard-sphere fluid in the entire density range. On the side, the theory reproduces most known earlier equations of state for hard spheres and determines their place in the hierarchy of approximations. Equations of state for van der Waals fluids are also presented, and their critical parameters are estimated.     

A consistent correction for Redlich-Kwong-Soave volumes

If the volumetric and phase behaviour of a fluid mixture is calculated by means of an equation of state, certain translations along the volume axis may be effected that leave the predicted phase equilibrium conditions unchanged. This property may be exploited in the form at a consistent correction to improve volume estimations by the Redlich-Kwong-Soave method. Applications of this improved method to pure liquids, mixtures of liquids or gases, and petroleum fluids show that markedly superior volume estimations are obtained, except in the neighbourhood of the pure-component critical points; nonetheless, critical volumes for mixtures can be estimated correctly.