Topic 2 Non Steady State Diffusion
Topic 2 Non Steady State Diffusion
Topic 2 Non Steady State Diffusion
( Netintothe
rate of molar flow rate − Net rate of molar flow rate = Rate of accummulation of species A
control volume ) ( out of the control volume ) ( ¿ thecontrol volume )
The molar flow rate of species A by diffusion at the plane Z=Z is given by Fick's law:
∂C A
W Az =−D AB S ( )
∂z z
−−−1
The molar flow rate of species A by diffusion at the plane Z=Z+ΔZ, is:
W Az =−D AB S ( ∂C∂ z )
A
z +∆ z
−−−2
∂C A
S ∆ z−−−3
∂z
∂ CA ∂ CA ∂ CA
−D AB S ( )
∂z z
+ D AB S
∂z ( ) z +∆ z
=S
∂z
∆ z−−−4
1
Rearranging and simplifying:
( ∂ C A /∂ z )z +∆ z−( ∂C A /∂ z )z ∂ CA
D AB ( ∆z )
=
∂t
−−−5
In the limit, as ΔZ→0
∂C A ∂2 C A
=D AB −−−6 ¿
∂t ∂ z2
The more general form, for three dimensions rectangular coordinates, is:
∂C A ∂2 C A ∂2 C A ∂2 C A
∂t
=D AB
∂x
2
+
∂y
2
+ (
∂z
2
−−−7 )
For one-dimensional diffusion in the redial direction only, for cylindrical and spherical
coordinates, Fick's second law becomes, respectively:
∂C A D AB ∂ ∂ C A
∂t
=
r ∂r
r
∂r (
−−−8 )
And
∂C A D AB ∂ 2 ∂ C A
∂t
= 2
r ∂r
r( ∂r) −−−9
Analytical solutions to these partial differential equations are available for a variety of boundary
conditions.
Only two of the most common useful solutions will be considered.
2
To solve the above partial differential equation, either the method of combination of variables or
the Laplace method is applicable. The result, in terms of the fractional accomplished
concentration change (θ), is:
C A−C A 0 z
θ=
C A 1−C A 0
=erfc
(2 √ D AB t)−−−10
∂ CA C AS−C A 0 −z 2
n A =−D AB S ( )
∂z z =0 (
=D AB S
√ π D AB t
exp
) (
4 D AB t )z=0
−−11
Thus:
D AB
n A|z =0=
√ πt
S ( C As−C A 0 )−−−12
We can determine the total number of moles of solute, NA, transferred into the semi-infinite
medium by integrating equation 12 with respect to time:
t t
D AB dt D AB t
0 √
N A =∫ n A|z=0 dt =
πt
S ( C As−C A 0 ) ∫ =2 S ( C As−C A 0 )
0 √t π √
−−13
3
Semi-infinite Medium correlations
4
5