Virial mixtures

This chapter presents a general method for estimating nonidealities in a vapor mixture containing any number of components this method is based on the virial equation of state for ordinary substances and on the chemical theory for strongly associating species such as carboxylic acids. The method is limited to moderate pressures, as commonly encountered in typical chemical engineering equipment, and should only be used for conditions remote from the critical of the mixture. 

For a pure vapor the virial coefficients are functions only of temperature for a mixture they are also functions of composition. An important advantage of the virial equation is that there are theoretically valid relations between the virial coefficients of a mixture and its composition. These relations are ... 

A component in a vapor mixture exhibits nonideal behavior as a result of molecular interactions only when these interactions are very wea)c or very infrequent is ideal behavior approached. The fugacity coefficient (fi is a measure of nonideality and a departure of < ) from unity is a measure of the extent to which a molecule i interacts with its neighbors. The fugacity coefficient depends on pressure, temperature, and vapor composition this dependence, in the moderate pressure region covered by the truncated virial equation, is usually as follows ... 

To illustrate calculations for a binary system containing a supercritical, condensable component. Figure 12 shows isobaric equilibria for ethane-n-heptane. Using the virial equation for vapor-phase fugacity coefficients, and the UNIQUAC equation for liquid-phase activity coefficients, calculated results give an excellent representation of the data of Kay (1938). In this case,the total pressure is not large and therefore, the mixture is at all times remote from critical conditions. For this binary system, the particular method of calculation used here would not be successful at appreciably higher pressures. 

VPLQFT is a computer program for correlating binary vapor-liquid equilibrium (VLE) data at low to moderate pressures. For such binary mixtures, the truncated virial equation of state is used to correct for vapor-phase nonidealities, except for mixtures containing organic acids where the "chemical" theory is used. The Hayden-0 Connell (1975) correlation gives either the second virial coefficients or the dimerization equilibrium constants, as required. 

Subroutine BIJS2. This subroutine calculates the pure-component and cross second virial coefficients for binary mixtures according to the method of Hayden and O Connell (1975). 

CALCULATE EFF SECOND VIRIAL COEFFICIENT FOR COMP I IN MIXTURE, SS(I)... 

BUS calculated second virial coefficients for pure compoments and all binary pairs in a mixture of N components (N 20) at specified temperature. These coefficients are placed in common storage /VIRIAL/. 

Fender B E F and Halsey G D Jr 1962 Second virial coefficients of argon, krypton and argon-krypton mixtures at low temperatures J. Chem. Phys. 36 1881... 

Gas mixtures are subject to the same degree of non-ideality as the one-component ( pure ) gases that were discussed in the previous section. In particular, the second virial coefficient for a gas mixture can be written as a quadratic average... 

The CS pressures are close to the machine calculations in the fluid phase, and are bracketed by the pressures from the virial and compressibility equations using the PY approximation. Computer simulations show a fluid-solid phase transition tiiat is not reproduced by any of these equations of state. The theory has been extended to mixtures of hard spheres with additive diameters by Lebowitz [35], Lebowitz and Rowlinson [35], and Baxter [36]. 

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 ... 

Mixing mles 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 9 for mixtures are usually given by the following, where, as for the second virial coefficient, = 0-. 

Miscellaneous Generalized Correlations. Generalized charts and corresponding states equations have been pubhshed for many other properties in addition to those presented. Most produce accurate results over a wide range of conditions. Some of these properties include (/) transport properties (64,91) (2) second virial coefficients (80,92) (J) third virial coefficients (72) (4) Hquid mixture activity coefficients (93) (5) Henry s constant (94) and 6) diffusivity (95). 

Second Virial Coefficient. A group contribution method including polar and nonpolar contributions has been proposed for second virial coefficients (241). This method has been appHed to both pure components and mixtures, the latter through prediction of cross-second virial coefficients. 

The coefficient Bij characterizes a bimolecular interaction between molecules i and J, and therefore Bij = Bji. Two lands of second virial coefficient arise Bn and By, wherein the subscripts are the same (i =j) and Bij, wherein they are different (i j). The first is a virial coefficient for a pure species the second is a mixture property, called a cross coefficient. Similarly for the third virial coefficients Cm, Cjjj, and are for the pure species and Qy = Cyi = Cjn, and so on, are cross coefficients. 

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 ... 

When i = J, all equations reduce to the appropriate values for a pure species. When i j, these equations define a set of interaction parameters having no physical significance. For a mixture, values of By and dBjj/dT from Eqs. (4-212) and (4-213) are substituted into Eqs. (4-183) and (4-185) to provide values of the mixture second virial coefficient B and its temperature derivative. Values of and for the mixture are then given by Eqs. (4-193) and (4-194), and values of In i for the component fugacity coefficients are given by Eq. (4-196). 

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 ... 

Binary interaction parameters are determined for each pq pair p q) from experimental data. Note that = k and k = k = 0. Since the quantity on the left-hand side of Eq. (4-305) represents the second virial coefficient as predicted by Eq. (4-231), the basis for Eq. (4-305) lies in Eq. (4-183), which expresses the quadratic dependence of the mixture second virial coefficient on mole fraction. 

The MS closure results from s = 2. The HNC closure results from s = 1. In the latter two expressions, additional adjustable parameters occur, namely ( for the RY closure and for the BPGG version of the MS approximation. However, even when adjustable, these parameters cannot be chosen at will, as they should be chosen such that they eliminate the so-called thermodynamic inconsistency that plagues many approximate integral equations. We recall that a manifestation of this inconsistency is that there is a difference between the pressure as computed from the virial equation (10) and as computed from the compressibility equation (20). Note that these equations have been applied to a very asymmetric mixture of hard spheres [53,54]. Some results of the MS closure are plotted in Fig. 4. The MS result for y d) = g d) is about the same as the MV result. However, the MS result for y(0) is rather poor. Using a value between 1 and 2 improves y(0) but makes y d) worse. Overall, we believe the MS/BPGG is less satisfactory than the MV closure. 

The equilibrium between a pure solid and a gaseous mixture is one of very few classes of solution for which an exact treatment can be made by the methods of statistical mechanics. The earliest work on the theory of such solutions was based on empirical equations, such as those of van der Waals,45 of Keyes,44 and of Beattie and Bridgemann.3 However, the only equation of state of a gas mixture that can be derived rigorously is the virial expansion,46 66... 

The next level of approximation is valid to higher pressures. It assumes that the gas mixture obeys the virial equation of state, with the third, fourth and higher, virial coefficients equal to zero. Thus... 

Dymond, J. H. Smith, E. B. The Virial Coefficients of Pure Gases and Mixtures Oxford Press New York, 1980. 

Special care has to be taken if the polymer is only soluble in a solvent mixture or if a certain property, e.g., a definite value of the second virial coefficient, needs to be adjusted by adding another solvent. In this case the analysis is complicated due to the different refractive indices of the solvent components [32]. In case of a binary solvent mixture we find, that formally Equation (42) is still valid. The refractive index increment needs to be replaced by an increment accounting for a complex formation of the polymer and the solvent mixture, when one of the solvents adsorbs preferentially on the polymer. Instead of measuring the true molar mass Mw the apparent molar mass Mapp is measured. How large the difference is depends on the difference between the refractive index increments ([dn/dc) — (dn/dc)A>0. (dn/dc)fl is the increment determined in the mixed solvents in osmotic equilibrium, while (dn/dc)A0 is determined for infinite dilution of the polymer in solvent A. For clarity we omitted the fixed parameters such as temperature, T, and pressure, p. 

We could, of course, attempt more sophisticated simulations of scale formation. Since the fluid mixture is quite concentrated early in the mixing, we might use a virial model to calculate activity coefficients (see Chapter 8). The Harvie-Mpller-Weare (1984) activity model is limited to 25 °C and does not consider barium or... 

Spycher, N. F. and M. H. Reed, 1988, Fugacity coefficients of H2, CO2, CH4, H2O and of H2O-CO2-CH4 mixtures, a virial equation treatment for moderate pressures and temperatures applicable to calculations of hydrothermal boiling. Geochimica et Cosmochimica Acta 52, 739-749. 

---