MT Lecture 3 - Diffusion in Gases
MT Lecture 3 - Diffusion in Gases
Overview
Diffusion of gases A & B + convection
Diffusion of A through stagnant, non- diffusing B
Diffusion through a varying cross-sectional area
Diffusivities of gases
J A DAB
dC A
dz
.1
The molar rate of transfer of A, per unit area due to molecular motion is
given by:
dx
J A* C.DAB
dz
.2
We have considered Ficks law for diffusion in a stationary fluid; i.e. there has been no net
movement or convective flow of the binary mixture of A and B.
The flux NA can be converted to velocity of diffusion of A by:
m kgmolA
J * A (kgmolA / s.m 2 ) v Ad C A
3
s
m
.3
The average molecular velocity of the whole fluid relative to a stationary point is
vM (m/s).
Expressed mathematically, the overall velocity of A relative to a stationary point
is the sum of the diffusion velocity and the convective velocity:
v A v Ad vM
.4
VA
VAD
VM
v A v Ad vM
v Ac A v Ad c A vM c A
.5
N AT J
vM c A
.6
Let NT be the total convective flux of the whole stream relative to a stationary point.
NT cvM N AT N BT
Then solve it for VM
N AT N BT
vM
C
N AT
.8
N AT J * A vM c A
CA
N AT N BT
J A
C
*
.7
We get
.9
J C.DAB
*
A
N AT C.DAB
dx A
dz
dx A C A
N AT N BT .10a
dz
C
or
N AT DAB
dC A C A
N AT N BT
dz
C
.10b
Diffusion plus convection when the flux NA is used relative to a stationary point.
for Component B
dxB C B
N AT N BT
N BT C.DBA
dz
C
dC B C B
N AT N BT
N BT DBA
dz
C
.11a
.11b
nT
V
CT =
P
=
RT
(12b)
Total pressure
DAB
RT
dpA
dz
pA (N + N )
A
B
+
P
(..13)
10
In terms of concentration of A:
NA = -DAB dCA
dz
CA (N + N )
A
B
+
CT
(..14)
(.16)
11
Air (B)
12
z2 z1
1
Liquid
Benzene
(A)
13
DAB
RT
dpA
dz
pA (N + 0)
A
+
P
(.17)
NA (1 - pA/P) = -
DAB
RT
z2
N dz =
A
z1
NA =
dpA
dz
pA2
DAB
RT
pA1
dpA
(1 - pA/P)
DAB P
P - pA2
ln
RT(z2 z1)
P pA1
(.18)
14
pB,LM =
(pB2 pB1 )
ln(pB2 /pB1 )
(P pA2 ) (P pA1 )
ln[(P - pA2 )/ (P - pA1 )]
(pA1 pA2 )
ln[(P - pA2 )/ (P - pA1 )]
DAB P
(pA1 - pA2 )
RT(z2 z1) pB,LM
(.19)
15
(xB2 xB1 )
ln(xB2 /xB1 )
(1 xA2 ) (1 xA1 )
(xA1 xA2 )
ln[(1 - xA2 )/ (1 - xA1 )]
(.21)
16
Example:
Diffusion of water through stagnant, non-diffusing air:
Water in the bottom of a narrow metal tube is held at a constant
temperature of 293 K. The total pressure of air (assumed to be dry)
is 1 atm and the temperature is 293 K. Water evaporates and
diffuses through the air in the tube, and the diffusion path is 0.1524
m long. Calculate the rate of evaporation at steady state. The
diffusivity of water vapour at 1 atm and 293 K is 0.250 x 10-4 m2/s.
Assume that the vapour pressure of water at 293 K is 0.0231 atm.
Answer: 1.595 x 10-7 kmol/m2.s
17
Air (B)
where
pB,LM =
z2 z1
(pA1 pA2 )
m2/s;
DAB = 0.250 x
P = 1 atm; T = 293 K;
Water (A)
z2 z1 = 0.1524 m;
pA1 = 0.0231 atm (saturated vapour pressure);
pA2 = 0 atm (water vapour is carried away by air at point 2)
R = 8314.34 m3.Pa/kg mol.K
18
NA =
(0.0231 0 )
ln[(1 - 0 )/ (1 0.0231 )]
= 0.988 atm
19
20
21
In the below Fig we see a sphere of fixed radius r1 (m) in an infinite gas
medium. Component (A) at partial pressure pA1 at the surface is diffusing
into the surrounding stagnant medium (B), where pA2 = 0 at some large
distance away. Steady-state diffusion will be assumed.
N AT N AT
A
4r 2
N AT
NAT = kg.mol/s
.22
Since this is a case of A diffusing through stagnant B, NAT will be equated to Eq. (22),
giving:
N AT
N AT
DAB
dPA
2
4r
RT 1 PA P dr
.23
Note that dr was substituted for dz. Rearranging and integrating between r1
and r2 (at a point a large distance away),
N AT
4
r2
r1
N AT
4
dr
DAB
2
r
RT
PA1
dPA
P 1 PA P
A2
1 1 DAB P P PA2
ln
RT
P PA1
r1 r2
.24
.25
23
N AT
4
1 1 DAB P P PA2
ln
RT
P PA1
r1 r2
.26
. 27
This equation can be simplified further. If PA1 is small compared to P (a dilute gas
phase), PBM P.
Also, setting 2r1 = D1 diameter and CA1 = PA1/RT, we obtain
24
2 DAB
N AT 1
PA1 PA2
RTD1
2 DAB
N AT 1
C A1 C A2
r1
. 28
. 29
This equation can also be used for liquids, where DAB is the diffusivity of A in the
liquid.
10/3/2016
Solution:
DAB = 6.92 x 10-6 m2/s
P = 1 atm = 1.01325 x 105 Pa
PA1 = 0.555 mm Hg = 74.0 Pa
PA2 = 0
R = 8314.34 m3.Pa/kg mol.K
PB1 = P PA1=1.01251 x 105 Pa
PB2 = P PA2=1.01325 x 105 Pa
Since the PB1 and PB2 are close to each other,
PBM
PB 2 PB1
1.0129 105 Pa
2
DAB P PA1 PA 2
N AT 1
9.68 10 8 kgmolA / s.m 2
RTr1
PBM
27
Explanation
DAB a T3/2
DAB a (1/P)
DAB a (1/MA)1/2
DAB a (1/MB)1/2
1.8583 10 7 T
2
PW D , ABs AB
3
2
1
1
MA MB
Units:
DAB = Diffusion coefficient (m2/s)
T = Temperature (K)
P = Absolute Pressure (atm)
MA, MB = Molecular weights of A and B, respectively
sAB = Collision diameter () = (sA + sB)/2
WD,AB = Diffusive collision integral (dimensionless)
sAB and WD = Lennard-Jones (L-J) parameters
2. Fuller-Schettler-Giddings
Correlation
When L-J parameters are not available then use
this method:
1/ 2
1.00 x 10
D AB
P v A
1.75
1/ 3
1
1
MA MB
vB
1/ 3 2
Diffusivity of Gases
Example
Calculate the diffusivity of Butanol (A)
gas in air (B) at 1 atm abs using the
Fuller et al. method.
a) 0oC
b) 26oC
Ans 7.73x10-6 and 9.05x10-6 m2/s
Ans:
MA (butanol) = 74.1 (C4H10O)
MB (air) = 29
= 4 16.5 + 10 1.98 + 1 5.48 = 91.28
= 20.1
1/ 2
1.00 x 10
DAB
P v A
1.75
1/ 3
1
1
MA MB
vB
1/ 3 2
Diffusivities of Gases
37
THANK YOU~