Lecture 4: Diffusion: Fick's Second Law: DX X DC) A
Lecture 4: Diffusion: Fick's Second Law: DX X DC) A
Lecture 4: Diffusion: Fick's Second Law: DX X DC) A
Todays topics
Learn how to deduce the Ficks second law, and understand the basic meaning, in
comparison to the first law.
Learn how to apply the second law in several practical cases, including homogenization,
interdiffusion in carburization of steel, where diffusion plays dominant role.
Continued from last lecture, we will learn how to deduce the Ficks second law, and understand
the meanings when applied to some practical cases.
C
dx A
j(x) j(x+dx)
A
x x+dx
dx
We can define the local concentration and diffusion flux (through a unit area, A) at position x
as:
c (x,t), J(x)
we have
[ J ( x) J ( x + dx)]dtA
dc(x) = , J(x+dx) = J(x) + dJ
Adx
1
dc ( x, t ) dJ
then, we have =-
dt dx
or, rewrite it in this format with replacing d
c ( x, t ) J
=-
t x
dc( x) c( x)
from the first law: J = -D = -D
dx x
then, we have
c( x, t ) J 2c
=- =D 2 (1)
t x x
This is the Ficks second law.
So, Ficks first law can be considered as a specific (simplified) format of the second law when
applied to a steady state.
Now, lets consider two real practical cases, and see how to solve the Ficks second law in these
specific cases.
Consider a composition profile as superimposed sinusoidal variation as shown below, where the
solid line represents the initial concentration profile (at t=0), and the dashed line represents the
profile after time .
2
C
t=0
_ l
C
x
t=
x
At t=0, c = c + 0 sin
l
c 2c
with the Ficks second law, =D 2, where D is the diffusion coefficient, a constant.
t x
x t
At time t, c(x,t) = c + 0 sin e
l
Where = l2/2D, is defined as the relaxation time.
l x
The largest Deviation of concentration at x = , where sin = 1, the maximum.
2 l
It is an exponential decay, the longer the wavelength (l), the longer the relaxation time (), then
the slower decay. Short wavelength dies fast. Thats why shaking always helps speed up the
dispersion, because it enables wide spreading (smaller l) of the stuff (like particles) you try to
disperse.
Consider a thin layer of B deposited onto A, through annealing at high temperature, we will be
able to get the concentration profile at different times, from there then we can determine the
diffusion coefficient, D
3
B A
0 x
c 2c
Infinite integration of the Ficks second law, =D 2
t x
We have
2
c(x,t) = ex / 4 Dt
(2)
t
As the B diffuses into A, the total amount of B is fixed
c( x, t )dx = N = constant
0
Then,
x
( )2 x
x 2 / 4 Dt
e dx = 2 D e 2 Dt
d( )=N
0 t 0 Dt
x
To solve the above equation, lets define y = then we have,
2 Dt
2
2 D e y dy = N
0
2
since e y dy =
0 2
then,
N
=
D
so Eq. (2) can now be written as
4
N 2
c(x,t) = e x / 4 Dt (3)
Dt
as determined by this diffusion kinetics equation, the concentration profile of carbon at various
times will be like this
C t1 t3 > t2 > t1
t2
t3
t2 > t1 = N
t1 2 D
t2
N 2
c (x,t) = e x / 4 Dt , x (-, )
2 Dt
Situation b): Doping with a fixed surface concentration (e.g. carburization of steel)
Fe
c
CH4/CO
x=0
5
The solution of the Ficks second law can be obtained as follows, the surface is in contact with an
infinite long reservoir of fixed concentration of Cs. For x < 0, choose a coordinate system u.
diffuse
Cs u
du x=0 x
X
The fixed amount of dopant per area is Cs du=N, which diffuse toward right.
Then using Eq. (3) above, the slab du contributes to the concentration at x is
Cs du u 2 / 4 Dt
dc(x, t) = e
Dt
u
Defining y = , then,
2 Dt
2c s 2
c (x, t) = e y dy ,
x / 2 Dt
2c s 2 x / 2 Dt 2
= [ e y dy - e y dy ]
0 0
2c s x / 2 Dt 2
= [ - e y dy ]
2 0
x
= cs [ 1- erf ( )]
2 Dt
2 z
y2
Where error function erf (z) = e dy
0
6
Considering boundary conditions:
c (x = 0) = cs , constant, fixed.
c (x = ) = c0, corresponding to the original concentration of carbon existing in the phase, c0
remains constant in the far bulk phase at x = .
x
c (x, t) = cs (cs c0)erf ( )
2 Dt
the concentration profile shown above follows this diffusion equation.
Now lets consider Interdiffusion as shown below, which represents more general cases.
c1 > c2
C
c1 t=0
C1 C2 c1 + c2
2 t
c2
x=0 0