Viscosity and diffusion coefficients for the gaseous binary mixtures of benzene toluene, benzenephenol and benzenepxylene over a wide range of temperature and composition have been predicted using the semiempirical inversion method. The accuracies are within 3% and 4% for viscosities and diffusion coefficients, respectively.
INTRODUCTION
Transport properties of gaseous mixtures are of great importance in process design, mass and energy transfer because the mixtures are often encountered more than pure gases. From the theoretical point of view, due to the absent of data over a wide range of temperature, pressure and composition, the semiempirical methods are interested. The semiempirical inversion method has been used extensively to extract the pair potential function and then predict transport properties of pure and mixtures of gases.
1

10
The results of kinetic and statisticalmechanical theories provide theoretical expressions for various equilibrium and nonequilibrium properties in terms of the potential energy of interaction between a pair of molecules.
11
Thus the evaluation of such quantities from a known pair potential function is a reliable procedure and the problem is how to gain the intermolecular potential. Due to the fact of availability and measurability, transport properties of gases at low density are an important source of the information about the interactions at the molecular level. The viscosity is probably easiest to measure and most accurate than others. The general method which has been used to infer intermolecular potential from transport properties was adopting a model potential with a number of adjustable parameters which were varied until a good fit was obtained to a given set of experimental data. The main disadvantage of this method is that the parameters determined from one property, for example viscosity, are slightly different from those obtained from other ones, say diffusion.
A desirable and reliable manner of solving the problem would be to obtain potential energy u(r), directly from the experimental data without explicit assumption of a mathematical model for u(r), is the semiempirical inversion method which is developed by Smith and coworkers.
12

15
The present work while describes an iterative method for generating the effective pair potential energy for gaseous binary mixtures of benzenetoluene, benzenephenol and benzenepxylene systems, then use them to predict the low density transport properties over a wide range of temperature and composition. The accuracies for viscosity and diffusion coefficients are within 3% and 4%, respectively.
 Kinetic Theory of Gases
According to the kinetic theory of gases and the ChapmanEnskog solution
16
,
17
of the Boltzmann transport equation,
18
the transport properties of gases at low density can be related to a series of collision integrals Ω
^{(l,s)}
which depend on the intermolecular potential energy, u(r). The superscripts l and s appearing in Ω denote weighting factors that account for the mechanism of transport by molecular collision. They are l=2 and s=2 for viscosity and l=1 and s=1 for diffusion, respectively.
11
The collision integrals are defined as:
Where,
θ
is the scattering angle,
Q
^{(l)}
(
E
) is the transport collision integral,
b
the impact parameter, and
E
the relative kinetic energy of colliding partners. The deflection angle is given by:
The definition of collision integrals as dimensionless reduced quantities makes the calculations of transport properties more convenient. The reduced collision integral is defined as:
where range parameter σ denotes the intermolecular separation for which the potential is zero. Numerical differentiation of the mentioned collision integrals and usage of recursion relation can generate collision integrals higher than the ones mentioned. That is:
where reduced temperature
T
^{*}
is defined as:
in which the energy parameter ε represent the depth of the potential energy well.
We have shown that three successive numerical integrations (equations 11 to 13) should be done to obtain collision integrals. The potential energy would serve as the input information required in calculating the collision integrals and consequently the transport properties. Kinetic theory expressions for the transport properties in term of the collision integrals for pure gases and gaseous mixture are as following:
Pure gas,
Viscosity,
Selfdiffusion coefficient,
The ratios of the calculated collision integrals have been determined by the following equations:
These ratios are weak functions of
T
^{*}
and their magnitudes are approximately unity and exactly one for rigid spheres. Collision integrals and their ratios are functions of temperature and the parameters of the selected model for intermolecular forces.
For gaseous mixtures the equations for viscosity are:
Where,
A
^{*}
_{ij}
is the ratio of the collision integrals defined by equation (1.9),
x
is the mole fraction of components, and
η
_{ij}
is the interaction viscosity. Subscripts
i
and
j
represent the heavier and lighter component of the
i

j
pair, respectively. For diffusion we have:
For diffusion we have:
where
P
is the pressure and Δ
_{ij}
is a higher order correction term of the binary diffusion coefficient which can be defined as:
where
 The Inversion Method
The Inversion method uses a given set of experimental reduced viscosity collision integrals, Ω
^{*(2,2)}
, over a wide range of temperatures to invert a pair of values (Ω
^{*(2,2)}
,
T
) to
V
/
ε
as a function of
r
/
σ
using the relations:
G
is the inversion function which is a function of reduced temperature (
T
*) and can be calculated by using a model system whose intermolecular potential is known. The values of
G
for the LennardJones 126 potential as the initial model have been determined from Viehland
et al.
.
19
The new potential energies are closer approximations to the true potential than the potential of the initial model. This then can be used for calculations of improved collision integrals by using the integrals of equations (1.1) to (1.3) and the computer program developed by O’Hara and Smith.
20
,
21
This process is repeated by an iterative manner until convergence occurs. The convergence condition is the degree to which the calculated collision integrals for a given iteration are close to the experimental correlation within experimental accuracy. Our result obtained converged after two iterations.
 The extended principle of corresponding states
The law of corresponding states is based on this idea that the equilibrium and transport properties of substances whose molecules obey a potential function with two scaling parameters, energy scaling factor
ε
, and length scaling factor
σ
, can be expressed in terms of functions that are nearly universal when written in terms of scaling factors. The works on viscosity of noble gases
22
and also binary gas mixtures
23
,
24
showed that the function Ω
^{*(2,2)}
is universal. The extended law of corresponding state is a revised principle of corresponding states that characterizes each binary interaction with the aid of five scaling parameters instead of two. The five scaling parameters are:
σ
,
ε
,
ρ
^{*}
,
C
_{6}
^{*}
,
V
_{0}
^{*}
; also a quantum mechanical parameter Λ
^{*}
maybe come into play. The dimensionless parameter
C
_{6}
^{*}
, characterizes the longrange region of the potential energy curve and dimensionless parameters
ρ
^{*}
and
V
_{0}
^{*}
characterize the short range repulsive wall. The Λ
^{*}
becomes significant when the molecular partners are small and/or the temperature is low.
In our work the correlation function proposed by Najafi and coworkers
25
was applied in order to calculate the unlike pair potential energy. The proposed viscosity correlation function is as following:
where

a0=4.369×101±7.8×104

a1=4.505×101±1.3×103a2=5.326×102±8.1×104

a3=3.519×102±9.2×104a4=1.751×102±4.0×104

a5=2.773×103±7.0×105a6=1.529×104±4.3×106
RESULTS AND DISCUSSION
The extraction of information about the forces by analyzing the bulk properties will be done if the interaction potential energy be available. The degree of success then depends on the accuracy of both the measurements and theory connecting the forces to macroscopic properties. The transport properties of dilute gases, in special viscosity, which depends on the binary interactions, satisfy the above requirements. In this relation also inversion procedure plays an important role for generating an unlikemolecule potential from the viscosity data and their corresponding states correlations. This, in turn permits us to calculate collision integrals and consequently the transport properties more accurately than is possible by correlations of the corresponding states.
The reduced collision integrals and their ratios for BenzenePhenol system
The reduced collision integrals and their ratios for BenzenePhenol system
In the case of the corresponding states principle it should be mentioned that since the functional equation obtained from this principle apply in a more limited form to molecular gases than to noble gases, the arithmetic mean of the function should be used for gases 1 and 2. Also we are confined ourselves to
T
^{*}
≥1, consequently in order to proceed to the next iteration it is necessary to extrapolate
V
(
r
) at the long range region (low temperature). The extrapolation function that we used is
V
^{*}
=
C
_{6}
u
^{6}
, where
C
_{6}
is the dispersion coefficient and
u
is the reciprocal of the intermolecular distance, both in atomic units. It is remarkable to mention that the inversion procedure is insensitive to the nature of the extrapolating function.
26
In this work we obtained accurate reduced potential energies for three binary gas mixtures at low density by direct inversion of the viscosity collision integral equations. These accurate potential energies then have been used to produce improved collision integrals which in turn serve to calculate the transport properties.
1

3
contain the results of collision integrals and their ratios for three studied gaseous mixtures (benzenephenol, benzenetoluene and benzenepxylene).
The reduced collision integrals and their ratios for BenzeneToluene system
The reduced collision integrals and their ratios for BenzeneToluene system
The reduced collision integrals and their ratios for Benzenep Xylene system
The reduced collision integrals and their ratios for Benzenep Xylene system
According to these collision integrals and their ratios we calculated the viscosities of binary mixtures using equations 115 to 118 and then compared to experimental ones.
27
,
28
4

6
contain the calculated and experimental values of viscosity coefficients for studied mixtures. The absolute deviation for quantity X is calculated according to:
Experimental data for viscosity and diffusion coefficients for mixtures of benzenetoluene and benzenepxylene are taken from
Ref. 27
and for benzenephenol mixtures from
Ref. 28
.
Comparison of calculated and experimental* values of viscosity coefficients for benzenetoluene system
*Experimental data are taken from [27].
Also according to equations 119 to 123 we calculated diffusion coefficients for equimolar binary gaseous mixtures at different temperatures in 1 atmosphere pressure. Then we compared these values with experimental reported data. Results are depicted in
7

9
Comparison of calculated and experimental* values of viscosity coefficients for benzenephenol system
*Experimental data are taken from [28].
Comparison of calculated and experimental* values of viscosity coefficients for benzenepxylene system.
*Experimental data are taken from [27].
Also deviation plots based on deviation values are demonstrated in
1

6
.
As can be seen the most deviations points are within 03% for viscosity and to maximum 4% for diffusion coefficient. This is of great importance to note that the benzenephenol mixture is a nonpolarpolar system which includes large and nonspherical molecules and the results show the ability of inversion method for describing this system in addition to nonpolarspherical (noble) gas mixtures where the method is proposed at first.
22
The other two studied systems are also including large and nonspherical molecules but the result represents the advantage of the method used for predicting transport properties without introducing complicated potential functions and related mathematical difficulties.
Calculated and experimental* values of diffusion coefficients for equimolar binary gaseous mixtures of benzene toluene at different temperatures.
*Experimental data are taken from [27].
Calculated and experimental* values of diffusion coefficients for equimolar binary gaseous mixtures of benzenephenol at different temperatures.
*Experimental data are taken from [28].
Calculated and experimental* values of diffusion coefficients for equimolar binary gaseous mixtures of benzenepxylene at different temperatures.
*Experimental data are taken from [27].
Deviation plot for viscosities of benzenetoluene mixtures at different compositions and temperatures.
Deviation plot for viscosities of benzenephenol mixtures at different compositions and temperatures.
Deviation plot for viscosities of benzenepxylene mixtures at different compositions and temperatures.
Deviation plot for diffusion coefficients of equimolar mixtures of benzenetoluene at different temperatures.
Deviation plot for diffusion coefficients of equimolar mixtures of benzenephenol at different temperatures.
Deviation plot for diffusion coefficients of equimolar mixtures of benzenepxylene at different temperatures.
CONCLUSION
According to kinetic theory of gases and using the semiempirical inversion method we extract the pair potential function in the binary gaseous mixtures of benzenetoluene, benzenepxylene and benzenephenol systems. Then using collision integrals and the extended correspondingstates law we calculated viscosity and diffusion coefficients in mixtures over a wide range of temperature and composition. It is worthwhile to note that we selected mixtures having large and nonspherical molecules and also test the polarnonpolar system to examine the ability of the inversion treatment. The results showed that accuracy of method in mixtures is acceptable meanwhile the advantage of using a property (viscosity) to predict another one (diffusion coefficient) even for mixtures including polar molecules, is demonstrated as well.
Papari M. M.
,
Haghighi B.
,
Maftoon L.
,
Boushehri A.
2002
Int. J. Trans. Phenom.
4
17 
Maghari A.
,
Najafi M.
1995
J. Sci. I. R. Iran
6
99 
Hirschfelder J. O.
,
Curtiss C. F.
,
Bird R. B.
1964
Molecular Theory of Gases and Liquids
John Wiley
New York
Maitland G. C.
,
Smith E. B.
1977
Amer. Soc. Mech. Eng.
412
234 
Chapman S.
,
Cowling T. G.
1970
The Mathematical Theory of Nonuniform Gases
Cambridge University Press
London
Boltzman L.
1872
Wein. Sitz.
66
275 
Assael M. J.
,
Bich E.
,
Vogel E.
1994
J. Nonequilib. Thermodynamic
19
47 