1 Introduction To Mass Transfer and Diffusion PDF
1 Introduction To Mass Transfer and Diffusion PDF
MASS TRANSFER
1
TOPIC 1: DIFFUSION
COURSE OUTCOME
2
Course Activities:
Lectures; in-class examples; tutorials and assignments
EKC 217: Introduction to Mass Transfer and Diffusion
Gas-Gas
Gas-Liquid
Gas-Solid
Liquid-Liquid
Liquid-Solid
Solid-Solid
1) Gas-Gas
The phase contact between gas-gas is not practically
realized as most gasses are completely soluble in each
other.
2) Gas-Liquid
This is on of the most common phase contact in the
industries.
Among some of the industrial process that has this two
phase contact are distillation, absorption,
desorption/stripping.
a) Distillation
The main function is to separate a
liquid mixture, of miscible and
volatile substances into individual
components by vaporization;
Utilizes the differences in volatility
for separation while the vapor phase
is created from the liquid by
application of heat;
Example of industrial process:
Separation of a mixture of alcohol
and water into its component and
the separation of a mixture of
benzene and toluene into its
component (Figure 1).
EKC 217: Introduction to Mass Transfer and Diffusion
10
Benzene
Mixture of
benzene &
toluene
Toluene
Figure 1
b) Gas Absorption
In gas absorption, a soluble vapor
is absorbed by means of a liquid in
which the solute gas is more or less
soluble using selective absorbent.
For instance, if a mixture of air
and ammonia is in contact with
liquid water, a large portion of
ammonia but essentially no air will
dissolve in the liquid
(Figure 2).
Through this, the air/ammonia
mixture can be separated.
The process need more than one
column in order to achieve 100%
separation.
EKC 217: Introduction to Mass Transfer and Diffusion
11
Air (outlet)
Water (inlet)
Absorption
column
Figure 2
c) Gas Desorption/Stripping
Similar to gas absorption, but the only difference is purely
in the direction of solute transfer;
For example, if air is brought into contact with an
ammonia-water solution, some of the ammonia will leave
the liquid phase and enters the gas phase.
3) Gas-Solid
If a solid is moistened with a volatile liquid and is exposed
to a relatively dry gas, the liquid leaves the solid and
diffuses into the gas, an operation generally known as
DRYING.
12
4) Liquid-Liquid
Separations involving the contact of two insoluble
liquid phase are known as liquid-liquid extraction.
A simple and familiar example is the separation of acetone
from acetone-water solution using carbon tetrachloride. If the
acetone-water solution is shaken in a separatory funnel with
carbon tetrachloride and the liquids are allowed to settle, a
large portion of the acetone will be found in the carbon
tetrachloride-rich phase and will thus have been separated
from the water.
Water
Acetone-water
solution
Carbon
tetrachloride-acetone
solution
Carbon tetrachloride
solution
13
5) Liquid-Solid
Crystallization is one of the process that involves liquid and solid
phase. It is used to obtain materials in attractive and uniform
crystals of good purity, separating a solute from the melt or a
solution and leaving impurities behind .
Leaching is a process whereby liquid solvent is used to dilute
selective component in a solid mixture. For example, the leaching of
gold from its ore by cyanide solution and cotton seed oil from the
seeds by hexane.
6) Solid-Solid
E.g.: diffusion of carbon into iron during
case-hardening, doping of semiconductors for transistors, migration
of doped molecules in semiconductors
at high temperature.
EKC 217: Introduction to Mass Transfer and Diffusion
14
DIFFUSION
15
GAS
LIQUID
is assumed.
EKC 217: Introduction to Mass Transfer and Diffusion
16
MOLECULAR DIFFUSION
A container is filled
with dye solution
Water molecules
Dye molecules
EKC 217: Introduction to Mass Transfer and Diffusion
17
18
the lower region than the upper, there will be a net rate
of mass transfer of dye molecules will take place from
the lower region to the upper region.
After a longer period, the concentration of dye will be
uniform throughout the solution.
19
20
FICKS LAW
21
J A DAB
dc A
dz
(for component A)
A and B
------ (1.1)
where:
JA
DAB
cA
z
J A cDAB
------ (1.3)
dx A
dz
where:
c = Total molar concentration of component A and B
xA = Molar fraction of component A
NOTE: cxA = cA
EKC 217: Introduction to Mass Transfer and Diffusion
22
(since yA = pA/pT)
p A p A pT
p
yA T
RT pT RT
RT
where:
pA = partial pressure of A
pT = total pressure of the system
Therefore:
J A DAB
dc A
p dy
DAB T A
dz
RT dz
23
------ (1.4)
REVIEW:
Different expressions of the concentrations in a binary mixtures of component
A and B.
24
Exercise:
The composition of air is often given in terms of only the two
principal species in the gas mixture:
Oxygen, O2
Nitrogen, N2
yO2 = 0.21
yN2 = 0.79
25
J He DHeN 2
dcHe
dz
26
cHe 2
z1
cHe1
J He dz DHeN 2
J He
dcHe
pHe1 nHe
cHe1
RT
V
27
28
P
c A cB c
RT
------ (1.5)
dcA dcB dc 0
EKC 217: Introduction to Mass Transfer and Diffusion
------ (1.6)
J Az J Bz 0
------ (1.7)
dc A
J A DAB
dz
and
dcB
J B DBA
dz
30
------ (1.8)
dc A
dcB
DBA
0
dz
dz
------ (1.9)
DAB DBA
------ (1.10)
31
m kg mole
J A (kg mole A/m 2 s) Ad c A
3
m
s
33
------ (1.11)
Now lets consider what happens when the whole fluid is moving in bulk
or convective flow to the right.
The molar average velocity of the whole fluid relative to a stationary
point is M (m/s). Component A is still diffusing to the right, but now
its diffusion velocity, Ad , is measured relative to the moving fluid.
To a stationary observer, A is moving faster than the bulk of the phase,
since its diffusion velocity, Ad , is added to that of the bulk phase, M .
Expressed mathematically, the velocity of A relative to the stationary
point is the sum of the diffusion velocity and the average or convective
velocity:
------ (1.12)
A Ad M
34
Expressed pictorially:
A
Ad
A = actual velocity of A
Ad = diffusional velocity of A
M = velocity of bulk
c A A c A Ad c AM
------ (1.13)
Each of the three terms represents a flux. The first term, cAA, can
be represented by the flux NA kg mole A/m2s. This is the total flux
of A relative to the stationary point. The second term is JA, the
diffusion flux relative to the moving fluid. The third term is the
convective flux of A relative to the stationary point.
EKC 217: Introduction to Mass Transfer and Diffusion
35
N A J A cAM
------ (1.14)
Let N be the total convective flux of the whole stream relative to the
stationary point. Then,
------ (1.15)
N cM N A N B
Solving for M (bulk velocity) gives:
N A NB
c
------ (1.16)
cA
N A J A (N A NB )
c
EKC 217: Introduction to Mass Transfer and Diffusion
------ (1.17)
36
N A cDAB
dx A c A
(N A NB )
dz
c
------ (1.18)
N B cDBA
dxB cB
(N A NB )
dz
c
------ (1.19)
To solve Eq. (1.18) or (1.19), the relation between the flux NA and NB
must be known. Eq. (1.18) and (1.19) hold for diffusion in gas, liquid or
solid.
37
Eq. (1.18) and (1.19) can also be used in different forms. For
example, since N = cM and cxA = cA, thus:
N A c A M DAB
dc A
dz
------ (1.20)
38
39
Equimolar counterdiffusion
40
In equimolar counterdiffusion,
the molar fluxes of A and B is
equal, but in opposite direction
or the net volumetric and molar
flows are zero.
A typical example of this case is
the diffusion of A and B in the
vapor phase for distillation that
have constant molar overflow.
Since the net volumetric and molar flows are zero, thus Eq.
(1.18) can be used with the convective term is set to zero, as
shown below:
------ (1.21)
dx A
J A cDAB
dz
zT
DAB c dx A J A dz
x Ai
------ (1.22)
41
------ (1.23)
42
Example 2:
Ammonia gas (A) is diffusing through a uniform tube 0.10 m
long containing N2 gas (B) at 1.0132 x 105 Pa pressure and 298 K.
At a point 1, pA1 = 1.013 x 104 Pa and at a point 2, pA2 = 0.507 x 104
Pa. The diffusivity DAB = 0.230 x 10-4 m2/s. Calculate the flux JA
at steady state and repeat for JB.
Solution:
Given: P = 1.0132 x 105 Pa
z2 z1 = 0.10 m
T = 298 K
43
0.23 x 10 1.013 x 10
JA
0.507 x 10 4
8314(298)(0.10 0)
4
0.23 x 10 9.119 x 10
JB
9.625 x 10 4
8314(298)(0.10 0)
4
44
Unimolecular Diffusion
45
In unimolecular diffusion,
mass transfer of
component A occurs
through stagnant
component B, NB = 0.
Therefore, the total flux to
or away from the interface,
N is the same as NA.
N A cDAB
dx A
xA N A
dz
------ (1.24)
46
Rearranging:
dx A
N A (1 x A ) cDAB
dz
NA
1
or
dz
dx
DAB c
1 xA
------ (1.25)
Integrating:
x A dx
N A zT
1 xA
dz
ln
x Ai 1 x
DAB c 0
1 x Ai
A
------ (1.26)
Or:
DAB c 1 x A
NA
ln
zT
1 x Ai
------ (1.27)
47
dx A c A
(N A NB )
dz
c
------ (1.28)
Since NB = 0, therefore:
N A cDAB
dx A c A
(N A )
dz
c
------ (1.29)
49
PT
c
,
RT
------ (1.30)
p A x A PT
Hence:
DAB dp A p A
NA
NA
RT dz PT
------ (1.31)
Rearranging:
pA
DAB dp A
N A 1
RT dz
PT
------ (1.32)
Integrating:
NA
z2
z1
DAB p A 2 dp A
dz
RT p A1 1 p A / PT
------ (1.33)
50
DAB PT
PT p A2
NA
ln
RT ( z2 z1 ) PT p A1
------ (1.34)
Since:
PT pA2 = pB2 ,
PT pA1 = pB1 ,
pB2 pB1 = pA1 pA2 , then:
DAB PT
p A1 p A2
pB 2
NA
ln
RT ( z2 z1 ) pB 2 pB1
pB1
------ (1.35)
pBM
pB 2 pB1
ln( pB 2 / pB1 )
------ (1.36)
51
DAB PT
NA
( p A1 p A2 )
RT ( z2 z1 ) pBM
------ (1.37)
52
cA M y A
and
N
M
M
where:
M = molar density (kgmole/m3)
= 1/22.41 kgmole/m3 (at standard conditions, 0C & 1 atm)
yA = mole fraction of component A in vapor phase
N = total convective flux of the whole stream relative to
the stationary point (kgmole/m2s)
M = molar average velocity (ms-1)
cA = molar concentration of component A (kgmole/m3)
EKC 217: Introduction to Mass Transfer and Diffusion
53
dy A
N A y A N DAB M
dz
------ (1.38)
dy A
N A y A N A DAB M
dz
54
------ (1.39)
dy A
N A (1 y A ) DAB M
dz
z2
yA2
z1
y A1
N A dz DAB M
------ (1.40)
dy A
(1 y A )
DAB M 1 y A2
NA
ln
z2 z1 1 y A1
------ (1.41)
------ (1.42)
55
Similarly,
yB1 = 1 yA1
yB2 = 1 yA2
yB2 yB1 = yA1 yA2
Then,
DAB M y A1 y A2 yB 2
NA
ln
z2 z1 yB 2 yB1 yB1
The logarithmic mean of
yB1 and yB2 is given by:
yBM
yB 2 yB1
ln( yB 2 / yB1 )
------ (1.43)
------ (1.44)
DAB M
NA
( y A1 y A2 )
( z2 z1 ) yBM
EKC 217: Introduction to Mass Transfer and Diffusion
56
------ (1.45)
Example 3:
Water in the bottom of a narrow metal
tube is held at a constant temperature of
293 K. The total pressure of air (assume
dry) is 1.01325 x 105 Pa (1.0 atm) and the
temperature is 293 K (20C). Water
evaporates and diffuses through the air in
the tube and the diffusion path z2 z1 is
0.1524 m (0.5 ft) long. The diagram is
similar to the shown figure. Calculate
the rate of evaporation at steady state in
lb mol/ft2 h and kgmole/m2 s. The
diffusivity of water vapor at 293 K and 1
atm pressure is 0.250 x 10-4 m2/s. Assume
that the system is isothermal. Use SI and
English units.
EKC 217: Introduction to Mass Transfer and Diffusion
57
pA2 = 0
Air (B)
0.5 ft
NA
pA1 = 0.0231 atm
Water (A)
T = 293 K
Solution:
The diffusivity is converted to ft2/h by using the conversion factor (refer
Appendix 1, McCabe, Smith and Harriott).
2345.6
p A1
0.0231 atm
5
1.01325 x 10
and
p A2 0 (pure air)
58
pB 2 pB1
1.00 0.9769
Since pB1 is close to pB2 , the linear mean (pB1 +pB2)/2 could be
used and would be very close to pBM .
59
DAB PT
( p A1 p A 2 )
RT ( z 2 z1 ) pBM
0.969(1.0)(0.0231 0)
0.730(528)(0.5)(0.988)
1.175 x 10 4 lbmole/ft 2 h
NA
DAB PT
( p A1 p A2 )
RT ( z2 z1 ) pBM
8314(293)(0.1524)(1.001 x 105 )
1.595 x 10 7 kgmole/m 2 s
60
Exercise:
a) For the diffusion of solute A through a layer of gas to an
absorbing liquid with yAi = 0.20 and yA = 0.10, calculate the
rate transfer for unimolecular diffusion compared to that for
equimolar counter diffusion.
b) What is the value of yA halfway through the layer for
unimolecular diffusion?
[Ans: yA = 0.1515]
61
63
------ (1.46)
cav
M av
1
2
M1 M 2
------ (1.47)
where:
cav
M1
1
64
NA = constant, NB = 0
DAB
NA
( z2 z1 ) xBM
where:
( x A1 x A2 )
M av
xBM
------ (1.48)
xB 2 xB1
ln( xB 2 / xB1 )
DAB
NA
(c A1 c A2 )
( z2 z1 )
EKC 217: Introduction to Mass Transfer and Diffusion
------ (1.49)
65
Example 4:
Calculate the rate of diffusion of acetic acid (A) across a film of
nondiffusing water (B) solution 1 mm thick at 17C when the
concentrations on opposite sides of the film are 9 and 3 wt %,
respectively. The diffusivity of acetic acid in the solution is
0.95 x 10-9 m2/s.
Solution:
Given:
(z2 z1) = 0.001 m
MA = 60.03 kg/kmole
MB = 18.02 kg/kmole
At 17 C: Density of the 9% solution = 1012 kg/m3
Density of the 3% solution = 1003.2 kg/m3
EKC 217: Introduction to Mass Transfer and Diffusion
66
0.09 / 60.03
0.0015
1012
52.7 kmole/m 3
M 1 19.21
67
1
19.21 kg/kmole
0.0520
Similarly, at point 2:
x A2
0.03 / 60.03
0.0005
1
18.42 kg/kmole
0.0543
1003.2
54.5 kmole/m 3
M2
18.42
Then,
M av
52.7 54.5
53.6 kmole/m 3
2
68
xBM
xB 2 xB1
0.9908 0.9712
0.980
ln( xB 2 / xB1 ) ln(0.9908 / 0.9712)
0.95 x 10 9
53.6(0.0288 0.0092)
NA
(0.001)(0.980)
1.018 x 10-6 kmole/m 2 s
69