Virial coefficients mixing rules

Although PVT equations of state are based on data for pure fluids, they are frequently appHed to mixtures. 7h.e virial equations are unique in that rigorous expressions are known for the composition dependence of the virial coefficients. Statistical mechanics provide exact mixing rules which show that the nxh. virial coefficient of a mixture is nxh. degree in the mole fractions ... 

Although developed for pure materials, this correlation can be extended to gas or vapor mixtures. Basic to this extension is the mixing rule for second virial coefficients and its temperature derivative ... 

The application of cubic equations of state to mixtures requires expression of the equation-of-state parameters as func tions of composition. No exact theory like that for the virial coefficients prescribes this composition dependence, and empirical mixing rules provide approximate relationships. The mixing rules that have found general favor for the Redhch/Kwong equation are ... 

Mixing rules for the parameters in an empirical equation of state, eg, a cubic equation, are necessarily empirical. With cubic equations, linear or quadratic expressions are normally used, and in equations 34—36, parameters b and 0 for mixtures are usually given by the following, where, as for the second virial coefficient, 0 = 0ji. 

However, these new mixing rules (based both to infinite- or zero pressure limit) give, for the composition dependence of the second virial coefficient, results that are inconsistent with those obtained from statistical mechanics. 

Here CiU and C222 are the third virial coefficients for pure species 1 and 2 whereas C,12 and C122 are cross-coefficients. Published generalized correlation for third virial coefficients ) are based on a very limited supply of experimen data. Consistent with the mixing rules of Eq. (11.44) and (14.1), the temperatu ( derivatives of B and C are given exactly by... 

The classical expressions for ai, pf , and up are similar to the pure-gas expressions with the slight added complication of having two sets of isolated-molecule parameters o oo, fio,. Fo, The determination of the mixed-pair parameters involved in uq is much more difficult. The determination of these parameters from mixed viscosities and pressure virial coefficients is hampered by a shortage of experimental data, and the usual procedure is to employ a set of empirical combining rules which relate the mixed parameters to those of the pure gases. A number of such rules have been proposed for the 6-12 potential, the most widely used for dielectric calculations being... 

Two types of virial coefficients have appeared Bi 1 and. 622 for which the successive subscripts are tlie same, and B12, for which the two subscripts are different. The first type is a pure-species virial coefficient tlie second is a mixture property, known as a cross coefficient. Both are functions of temperature only. Expressions such as Eqs. (11.57) and (11.58) relate mixture coefficients to pure-species and cross coefficients. They are called mixing rules. 

The density series virial coefficients B, C, D,.. . depend on temperature and composition only. In practice, truncation is to two or three terms. The composition dependencies of B and C are given by the exact mixing rules... 

Although developed for pure materials, these correlations can be extended to gas or vapor mixtures. Basic to this extension are the mixing rules for the second virial coefficient and its temperature derivative as given by Eqs. (4-60) and (4-62). Values for the cross coefficients Bp with i and their derivatives are provided by Eq. (4-72) written in extended form ... 

Partial Molar Equation-of-State Parameters The parameters in equations of state as applied to mixtures are related to composition by mixing rules. Eor the second virial coefficient... 

The overall conclusion is that even though the conventional van der Waals mixing rules are simple to use and conform to the second virial coefficient boundary condition, they are very limited in their application and are not useful for either the correlation or the prediction of the VLE of complex mixtures. 

The WS mixing rule satisfies the low-density boundary condition that the second virial coefficient be quadratic in composition and the high-density condition that excess free energy be produced like that of currently used activity coefficient models, whereas the mixing rule itself is independent of density. This model provides a correct alternative to the earlier ad hoc density-dependent mixing rules (Copeman and... 

As indicated above, eqn. (3.3.8) is not the only option that can be used with eqn. (4.4.1). Indeed Tochigi et al. (1994) have proposed a MH V mixing rule consistent with the second virial coefficient boundary condition by combining eqns. (4.3.1 and 4.3.2 and 4.4.5). In their implementation they have eliminated the binary interaction parameter in eqn. (4.3.2), leading to the mixing rule... 

Tochigi, K., Kolar, R, Izumi, T., and Kojima, K., 1994. A note on a modified Huron-Vidal mixing rule consistent with the second virial coefficient condition. Fluid Phase Eq., 96 215-221. 

As first example we will illustrate the application of mixing rules with the Virial EOS. If only the second term is kept, then we may write Z = 1 + B ,(7 )/ E, exactly as for a pure component. However, this time the coefficient Bm is function of temperature as well as of composition. It can be demonstrated by statistical mechanics that the following mixing rule applies ... 

The factor 1.443 VT changes when other model equations of state are used. Note that (6.5.4) involves only pure component parameters, fl and fo,- so, no combining rules are needed. However, these mixing rules do not reproduce the known composition dependence of the second virial coefficient (4.5.18). 

Wong and Sandler has shown that the following mixing rule does satisfy the second virial coefficient equation ... 

The same equations are used for mixtures of gases. In this case the coefficients are combinations of the coefficients of the pure components. They are combined according to mixing rules, which are also specified by statistical mechanics for any number of components. The second virial coefficient is, for a mixture of m components,... 

In order to obtain interaction second virial coefficients for mixtures, some method is required for determining the acentric factor a>y and the pseudo-critical constants T ij and pertaining to the unlike interactions. In the present case, extended van der Waals one-fluid mixing rules are applied in terms of which... 

Wong and Sandler followed a different approach by matching the Helmholtz function at infinite pressure from the cubic equation of state and from an activity coefficient model. This approach ensures consistency with statistical mechanics requirements that the second virial coefficient of a mixture has a quadratic dependence on composition. For the case of the Peng-Robinson equation of state, the Wong-Sandler mixing rules are ... 

Any error in the virial coefficient of the pure component will also introduce an error into the calculated fugacity independent of the choice of mixing and combining rules. 

---