CME 301 - Mass Transfer Differential MT Equations For Unsteady-State Molecular Diffusion (USS-MD)
CME 301 - Mass Transfer Differential MT Equations For Unsteady-State Molecular Diffusion (USS-MD)
CME 301 - Mass Transfer Differential MT Equations For Unsteady-State Molecular Diffusion (USS-MD)
2
Assumptions:
3
Model:
∂c A
−∇N A + RA =
∂t ∂ 2 c A ∂c A
DAB =
∂c A ∂z 2
∂t
NA =− DAB + c AV
∂z
IC:
cA(z,0) = cA0 @ t = 0 & all z
BCs:
cA(0,t) = cAs @ z = 0 & all t
cA(∞,t) = cA0 @ z = ∞ & all t 4
Solution:
c As − c A z
=
erf erf (φ )
c As − c A 0 2 D t
AB
NOTE:
Er r or Funct ion
erf(0) = 0
erf(∞) = 1
cA − cA 0
1 − erf (φ ) =
= erfc (φ )
c As − c A 0
Compl ement ar y
Er r or Funct ion
5
Flux at z = 0:
∂c A
NA = −D AB
z =0
∂z z =0
− [c As − c A 0 ]
= −D AB
πD t
AB
D AB
= [c As − c A 0 ]
πt
6
2. USS-MD in a Medium with Finite
Dimensions and Simple Geometry
Assumptions:
1. In solids & liquids, bulk flow = 0 cAs NA NA cAs
cA.V = 0
2. Dilute A UMD approaches EMCD
L/2 L/2
(same as point 1)
3. No Reaction (RA = 0)
z=0 z=L
4. Constant T and PT cT = constant
5. Binary Mixture
7
Solution
∂c A
−∇N A + RA =
∂t ∂ 2 c A ∂c A
DAB =
∂c A ∂z 2
∂t
NA =− DAB + c AV
∂z
IC: cA(z,0) = cA0
BCs:
cA(0,t) = cAs
cA(L,t) = cAs
∂cA(L/2,t)/∂z = 0
8
Concentration profile of A, i.e., cA(z, t):
c A − c As 4 ∞ 1 z 2 2 DAB t
Y == ∑ sin nπ exp −n π
c A0 − c As π n =1 n L L2
n=
1,3,
∂c A 4 DAB ∞
nπ 2 2 DAB t
NA L =
− DAB = ( cAs − cA0 ) ∑ cos
exp −n π 2
=0
z=
2 ∂z L n =1 2 L
∂c A 4 DAB ∞
2 2 DABt
N A ( 0, t ) =
− DAB =− ( cAs cA0 ) ∑ exp −n π 2 ; n=
1,3,
∂z z = 0 L n =1 L
9
Typical Concentration Profiles for
slab with Finite Boundaries
NA NA
cAs cA0
t
t
cA0 cAs
10
Time-Concentration Charts
(Appendix F) for USS-MD
11
Using Time-Concentration Charts:
Relative
4. Define n = x/x1 position
12
How to measure x1 and x?
Slab Cylinder (radial Sphere (radial
symmetry) symmetry)
cAs
x = x1 x1 = R
x
x1 = L/2 x1 = R
L
x=0 cAs
x1 = L/2 x=r cAs
x = x1 x=r
cAs
cAs
x = x1
x
x1 = L L
x=0 13
What is the relative resistance (m)?
DAB 1 kc x1 convective mass transfer resistance
=
m = =
kc x1 1 DAB molecular mass transfer resistance
m=6
Log scale! Y
m=2
m=1
m=0, m=0
n=1
X
Variable n
15
16
17
18
Figures F.4 to F.9 can also be used by changing α by
DAB and (T-T∞)/ (T0-T∞) by Y
m = 100
(T-T∞)/ (T0-T∞)
m=0
n=0
α ⋅t x
2
1 19
20
21
22
23
What if you have an object with more than
one symmetry plane?
Solution:
Solve for each symmetry alone, then:
Yoverall = Ysym1 Ysym2 Ysym3
24