MASS TRANSFER

Assoc. Prof. Dr. Emel AKYOL

6.MOLECULAR DIFFUSION IN SOLIDS
This type of diffusion in solids does not depend on the actual structure of the solid.

Generally, simplified equations are used. Using the general equation (2.5) for binary
d y A CA
N A  C DAB  (N A  N B ) (2.5)
dz C

CA
The bulk flow term, ( N A  N B ) , even if present, is usually small, since XA is quite small. Hence
C
it is neglected. Also, C is assumed constant giving for diffusion in solids,
𝐷𝐴𝐵 𝑑𝐶𝐴
𝑁𝐴 = − (6.1)
𝑑𝑧
Where DAB is diffusivity in m2/s A through B and usually is assumed constant independent of
pressure for solids. Note that DABDBA in solids.
Integration of equation (6.1) for a solid slab at steady state gives
D AB
N  (C A1  C A2) (6.2)
z 2z1
A
 MOLECULAR DIFFUSION IN SOLIDS
For the case of diffusion radially through a cylinder wall of inner radius r1 and outer r2 and length L,

𝑁𝐴 𝑑𝐶𝐴
= −𝐷𝐴𝐵
2𝜋𝑟𝐿 𝑑𝑟
(6.3)

2𝜋𝐿
𝑁𝐴 = 𝐷𝐴𝐵 (𝐶𝐴1 − 𝐶𝐴2 ) 𝑟 (6.4)
ln( 2 )
𝑟1
The diffusion coefficient DAB in the solid as stated above is not dependent upon the pressure of the gas or
liquid on the outside of the solid.
The solubility of a solute gas (A) in a solid is usually expressed as S in m3 solute (at STP of 0 oC and 1 atm)
per m3 solid per atm partial pressure of A.

𝑆𝑚3 (𝑆𝑇𝑃)/(𝑚3 𝑠𝑜𝑙𝑖𝑑.𝑎𝑡𝑚) 𝑆𝑃𝐴 𝑘𝑔 𝑚𝑜𝑙 𝐴

𝐶𝐴 = 𝑚3 𝑆𝑇𝑃
𝑃𝐴 𝑎𝑡𝑚 = (6.5)
22.414(𝑘𝑔 𝑚𝑜𝑙 𝐴) 22.414 𝑚3 𝑠𝑜𝑙𝑖𝑑
 EXAMPLE

The gas hydrogen at 17°C and 0.010 atm partial pressure is diffusing through a membrane of vulcanized neoprene
rubber 0.5 mm thick. The pressure of H2 on the other side of the neoprene is zero. Calculate the steady-state flux,
assuming that the only resistance to diffusion is in the membrane. The solubility S of H2 gas in neoprene at 17°C is
0.051 m3 (at . STP of Dec and 1 atm)/m3 solid· atm and the diffusivity DAB is 1.03 X 10 -10 m 2 /sat 17°C.
 MOLECULAR DIFFUSION IN SOLIDS
. Permeability Equations for Diffusion in Solids
In many cases the experimental data for diffusion of gases in solids are not given as diffusivities
and solubilities but as permeabilities, PM, in m3 of solute gas A at STP diffusing per second per
m2 cross-sectional area through a solid 1 m thick under a pressure difference of 1 atm
pressure. This can be related to Fick’s equation (6.2) as follows
D AB
N  (C  C ) (6.2)
z 2z
A A1 A2
1

From Eq. (6.2),

𝑆𝑃 𝐴 1 𝑆𝑃 𝐴 2
𝐶𝐴1 = , 𝐶𝐴2 = (6.6)
22 .414 22 .414
Substituting eq. (6.6) into (6.2)
𝐷 𝐴𝐵 𝑆(𝑃 𝐴 1 −𝑃 𝐴 2 ) 𝑃 𝑀 (𝑃 𝐴 1 −𝑃 𝐴 2 )
𝑁𝐴 = = kg mol/s.m2 (6.7)
22.414 (𝑧 2 −𝑧 1 ) 22 .414 (𝑧 2 −𝑧 1 )

where the permeability PM is

𝑚 3 (𝑆𝑇𝑃 )
𝑃𝑀 = 𝐷𝐴𝐵 𝑆 (6.8)
𝑠.𝑚 2 𝐶.𝑆 𝑎𝑡𝑚 /𝑚

For the cgs system, the permeability is given as PM’, cc(STB)/(s.cm2 C.S. atm/cm). This is
related to PM by
𝑃𝑀 = 10−4 PM’ (6.9)
When there are several solids 1,2,3,……., in series and L1,L2,…….,represent the thickness of
each, then Eq.(6.7) becomes
𝑃 𝐴 1 −𝑃 𝐴 2 1
𝑁𝐴 = 𝐿1 𝐿 (6.10)
22 .414
𝑃𝑀1
+𝑃 2 +⋯+
𝑀2
 EXAMPLE

A polyethylene film 0.000 15 m (0.15 mm) thick is being considered for use in packaging a pharmaceutical product at
30°C. If the partial pressure of O 2 outside is 0.21 atm and inside the package it is 0.01 atm, calculate the diffusion flux
of O 2 at steady state. Use permeability data from Table. Assume that the resistances to diffusion outside the film and
inside are negligible compared to the resistance of the film.
6.2. Diffusion on Porous Solids That Depends on Structure
6.2.1 Diffusion of Liquids on Porous Solids

Fig 6. Sketch of a typical porous solid

𝜀 𝐷 𝐴𝐵 (𝐶𝐴 1 −𝐶𝐴 2 )
𝑁𝐴 = (6.11)
𝜏 (𝑧 2 −𝑧1 )
𝜀
𝐷𝐴 𝑒𝑓𝑓 = 𝐷𝐴𝐵 , m2/s (6.12)
𝜏
where 𝜀 is the open void fraction, DAB is the diffusivity of salt in water, and 𝜏 is a factor which
corrects for the path longer than (𝑧2 − 𝑧1 ). Often the terms are combined into an effective
diffusivity.
 EXAMPLE

A sintered solid of silica 2.0 mm thick is porous with a void fraction ε of 0.30 and a tortuosity r of 4.0. The pores are filled
with water at 298 K. At one face the concentration of KCI is held at 0.1 0 g mol/liter, and fresh water flows rapidly by the
other face. Neglecting any other resistances but that in the porous solid, calculate the diffusion of KCl at steady state.
6.2.2 Diffusion of Gases on Porous Solids
If the voids shown in Figure 6 are filled with gases, then a somewhat similar situation exists. If
the pores are very large so that diffusion occurs only by Fickian-type diffusion, then Eq.(6.11)
becomes, for gases,
𝜀𝐷𝐴𝐵 (𝐶𝐴 1 −𝐶𝐴 2 ) 𝜀𝐷𝐴𝐵 (𝑃𝐴 1 −𝑃𝐴 2 )
𝑁𝐴 = = (6.13)
𝜏 (𝑧2 −𝑧1 ) 𝜏𝑅𝑇 (𝑧2 −𝑧1 )
7.DIFFUSION OF GASES IN POROUS
SOLIDS AND CAPILLARIES
In previous section diffusion in porous solids that depends on structure was discussed
for liquids and for gases. For gases it was assumed that the pores were very large and
Fickian-type diffusion occured. However, often the pores are small in diameter and
the mechanism of diffusion is basically changed.
Diffusion of gases in small pores occurs often in heterogeneous catalysis where gases
diffuse through very small pores to react on the surface of the catalyst. In freeze
drying of foods such as turkey meat, gaseous H 2 0 diffuses through very fine pores of
the porous structure.
7.DIFFUSION OF GASES IN POROUS
SOLIDS AND CAPILLARIES
Since the pores or capillaries of porous solids are often small, the diffusion of gases may depend upon the
diameter of the pores. We first define a mean free path λ; which is the average distance a gas molecule travels
before it collides with another gas molecule.
7.DIFFUSION OF GASES IN POROUS
SOLIDS AND CAPILLARIES

Figure 7.1
 Knudsen Diffusion of Gases

Knudsen Diffusion of Gases d t In Fig. 7.1 (a) a gas molecule A at partial pressure PA1 at
the entrance to a capillary is diffusing through the capillary having a diameter of d m.
The total pressure P is constant throughout. The mean free path ), is large compared
to the diameter d. As a result, the molecule collides with the wall and molecule-wall
collisions are important. This type is called Knudsen diffusion.
EXAMPLE

A H2(A)-C2H6(B) gas mixture is diffusing in a pore of a nickel catalyst used for

hydrogenation at 1.01325 x 105 Pa pressure and 373 K. The pore radius is 60 A
(angstrom). Calculate the Knudsen diffusivity DKA of H2 .
 Molecular Diffusion of Gases
As shown in Fig. 7.1 (b )when the mean free path λ is small compared to the pore diameter d or where
NKn ≤1/100, molecule-molecule collisions predominate and molecule-wall collisions are few. Ordinary
molecular or Fickian diffusion holds and Fick's law predicts the diffusion to within about 10%. The error
diminishes as NKn gets smaller since the diffusion approaches more closely the Fickian type.
 Transition-Region Diffusion of Gases
As shown in Fig. 7.1c, when the mean free path λ and pore diameter are intermediate in size between the two limits
given for Knudsen and molecular diffusion, transition-type diffusion occurs where molecule-molecule and molecule-
wall collisions are important in diffusion.
 Flux Ratios for Diffusion of Gases in Capillaries
1. Diffusion in open system. If diffusion in porous solids or channels with no chemical reaction is occurring where the
total pressure P remains constant, then for an open binary counterdiffusing system, the ratio of NA/NB is constant in
all of the three diffusion regimes

2. Diffllsion ill closed system. When molecular diffusion is occurring in a closed system shown in Fig. 7.2 at
constant total pressure P, equimolar counterdiffusion occurs.
EXAMPLE

A gas mixture at a total pressure of 0.10 atm abs and 298 K is composed of N2 (A) and He (B). The mixture is diffusing
through an open capillary 0.010 m long having a diameter of 5 x 10 -6 m. The mole fraction of N2 at one end is
XA1 =0.8 and at the other end is XA2 = 0.2. The molecular diffusivity DAB is 6.98 x 10 -5 m2/s at 1 atm, which is an
average value by several investigators. Calculate the flux NA at steady state