Lecture Notes on the Estimation of Binary Diffusivities
Lecture Notes on the Estimation of Binary Diffusivities
R. Shankar Subramanian
Department of Chemical and Biomolecular Engineering
Clarkson University
1
If σ AB is not known for a given binary system composed of a non-polar molecular
pair, construct it using
1
σ AB
= (σ A + σ B ) Equation (24-39), page 444
2
Values of the collision integral for molecular diffusion are given in Table K-1 in
Appendix K as a function of a dimensionless temperature kT / ε AB . Table K-2 in the
same Appendix lists values of ε / k for several molecules. For a binary system
composed of a non-polar molecular pair, construct ε AB / k using
ε AB εA εB
=
k k k
Fuller-Schettler-Giddings Equation
1/2
−3 1 1
M + M
1.75
10 T
DAB = A B
P ( ∑ v ) A + ( ∑ v ) B
1/3 1/3 2
2
Predictions can be made for polar-non-polar pairs or for polar-polar pairs, but we
won’t be able to discuss them here. See pages 446-447 in the book for information
about one approach called Brokaw’s method.
1 y2′ y′ y′
= + 3 + + n
D1− mixture D1− 2 D1−3 D1− n
y2 y2
=y2′ =
y2 + y3 + + yn 1 − y1
3
Binary diffusivity in liquids
T 7.4 ×10 ( Φ B M B )
−8 1/2
DAB =
µB VA0.6
DAB : Binary diffusivity of A at infinite dilution in B cm 2 / s .
T: Absolute temperature in Kelvin
µ B : Viscosity of solvent B in cp (this means centipoise)
M B : Molecular weight of solvent B
Φ B : Association parameter of solvent B
VA : Molal volume of solute A at normal boiling point in cm3 / g − mol
Use Tables 24.4 and 24.5 on pages 449-450 to estimate atomic and molal volumes
at the normal boiling point. The association parameter for common solvents can
be found in a table that appears below Table 24.5 on page 450.
A simpler result that gives nearly the same prediction is the correlation of Hayduk
and Laudie, which is to be used for evaluating the diffusion coefficients of non-
electrolytes in water in dilute solution (in principle, at infinite dilution). It is given
as Equation (24-53) on page 450 in the book.
The notation is the same as that used in the Wilke-Chang correlation, and the
diffusivity of A at infinite dilution in B is obtained in cm 2 / s .
Note that the above estimates for a solute at infinite dilution. In practice, they can
be used as adequate approximations so long as the solution is reasonably dilute;
however, unlike binary gas phase diffusivities, liquid phase diffusivities depend
strongly on concentration. In pages 451-452 in Chapter 24 of the textbook, the
authors briefly mention a common method used to obtain binary diffusivities in
non-dilute solutions.
4
Remarks added in April 2020
In a recent review article on separation column efficiency, Taylor and Duss (2019)
discuss a general approach for estimating the diffusion coefficient of a key species
in a multicomponent mixture in both gas phase mixtures and liquid mixtures. For
example, when using the equations proposed by Hirschfelder, Bird, and Spotz or
Fuller,-Schechter-Giddings for gas phase diffusion coefficients, if A is the key
species in the mixture, we lump all the other species into B, and obtain M B from
n
xi
MB = ∑ x ′M
i = 1≠ A
i i where xi′ =
1 − xA
In the summation, the case i = A is omitted, and the result for xi′ is used only for
i ≠ A . A similar weighted result is used for the “molecular volume” of B, if needed
in the correlation.
n
VB = ∑ x ′V
i = 1≠ A
i i
These ideas also can be extended to the liquid phase in the correlations of Wilke-
Chang or Hayduk and Laudie. Wherever information for species B is needed, a
weighted average on an A-free basis obtained in the above manner is substituted.
The viscosity of the solvent B is estimated using any method applicable to
mixtures, but on an A-free basis.
For more details, please consult the article by Taylor and Duss (2019).
References
J. Welty, C.E. Wicks, G.L. Rorrer, and R.E. Wilson, Fundamentals of Momentum,
Heat, and Mass Transfer, 5th edition, Wiley, Hoboken, NJ (2007).