Multicomponent Diffusion: Mixture-Averaged Approximation
The mixture-averaged diffusion model assumes that the relative mass flux due to molecular diffusion is governed by a Fick’s law type approximation.
Assuming that the diffusive flux, relative to the mass averaged velocity, is proportional to the mole fraction gradient, the mass flux is defined as:
(3-48)
Here ρi is the density, and xi the mole fraction of species i. Using the definition of the species density and mole fraction
,
Equation 3-48 can be expressed in terms of the mass fractions (ωi) in the manner of
Using Equation 3-48 together with the Maxwell-Stefan equations, where isobaric and isothermal conditions have been assumed, the following expression for the mixture-averaged diffusion coefficients can be derived (Ref. 3):
(3-49)
If instead the diffusive flux (relative to the mass averaged velocity) is assumed proportional to the mass fraction gradient, the mass flux is defined as:
(3-50)
For this assumption, using Equation 3-50 together with the Maxwell-Stefan equations, also assuming isobaric and isothermal conditions, the following expression for the mixture-averaged diffusion coefficients can be derived (Ref. 3):
(3-51)
The Mixture-Averaged Diffusion Coefficient
When using the mixture-averaged model, the diffusion is proportional to a single diffusion coefficient. The coefficient describes the diffusion of species i relative to the remaining mixture and is referred to as the mixture-averaged diffusion coefficient. The coefficient is explicitly given in terms of the multicomponent Maxwell-Stefan diffusivities Dik. As a consequence, no matrix inversion operation is required as for the Maxwell-Stefan diffusion model (when using four or more species). For low-density gas mixtures, the Dik components can be replaced by the binary diffusivities for the species pairs present.
Mixture Diffusion Correction
It should be noted that models for the multicomponent diffusivities based on Fick’s law are approximations. One drawback of this is that the net diffusive mass flux, the sum of all mass fluxes relative to the mass averaged velocity, is not constrained to zero. To account for this we add a flux, defined by a correction velocity uc, to enforce a zero net diffusive flux in the manner of:
(3-52)
Hence the correction velocity is:
(3-53)
Here ud,i is the diffusion velocity resulting from the flux assumption in Equation 3-48 or Equation 3-50. Note that the correction velocity is a constant correction (same for all species), but varies in space.
Using the correction velocity together with Equation 3-48, the resulting diffusive flux is
(3-54)
If instead Equation 3-50is used the resulting diffusive flux is
(3-55)
Additional Transport Mechanisms
Apart from molecular diffusion, transport due to thermal diffusion and migration of charged species in an electric field can be accounted for by adding the following terms to the diffusive flux:
where
zi (dimensionless) is the charge number
um, i the mobility of the ith species, and
(SI unit: V) is the electric potential.