Problems 0
Problems 0
Cedric J. Gommes
February 6, 2014
Contents
Ficks laws: a back to basics
3 Lifetime of an aerosol
11
15
(1)
at a surface
(4)
You may also know the flux of molecules J0 through a surface. In this case,
the condition is
Dn C = J0 at a surface
(5)
where n is the normal vector to the surface. You may of course have (slightly)
more complicated conditions, such as
Dn C = kC
at a surface
(6)
Think about what the latter equation means. When you can expect to encounter this type of boundary condition?
3
The situation where the concentration does not depend on the time, but
only on the position, is referred to as being stationary. It is important that
you understand that you can have transport phenomena going on although
the concentrations are stationary. It is the gradients that matter for
the transport, not the time derivatives. In a stationary situation, the
concentration obeys simply
C = 0
(7)
In Cartesian coordinates, the Laplacian operator is
2
2
2
C
C + zz
C + yy
C = xx
(8)
1
r (r2 r C)
r2
(9)
(10)
You are asked to estimate how long it takes for sparkling water to go flat
when it is kept in a PET bottle. The initial pressure in the bottle is 2 bars,
and one asks when the pressure will reach the limit 1.5 bars.
The volume of the bottle is 1.5 L. It is made of PET with thickness 250
microns, and the area of the bottle is 800 cm2 . The solubility of CO2 in PET
obeys Henrys law P = kc with k P ET = 29.4 L bar/mol. The solubility of
W
CO2 in water also obeys Henrys law with 1/kH
= 1.7m3 (ST P )/(m3 bar).
The diffusion coefficient of CO2 in PET is 2 1013 m2 /s.
Consider that the system reaches rapidly a stationary state. You can also
neglect the amount of CO2 inside the bottle that is in the gas phase.
CCO2
CW
0.25 mm
CPET
water
PET
air
Figure 3: Sketch of the CO2 concentration profile along a direction perpendicular to the bottle wall.
Solution
We shall first express the two Henry constants in the same units. They are
1/k W = 1.7/22.4 = 0.0759 mol/L/bar
1/k P ET = 0.034 mol/L/bar
where we have used the known volume of a mole of gas in STP conditions.
The CO2 concentration in the water is related to the pressure through
CW =
P
kW
(11)
P
k P ET
(12)
The value 1/k P ET being different from 1/k W , there is a concentration discontinuity at the water/PET interface.
The PET wall being very thin compared to the radius of the bottle,
the PET wall can be assumed to be flat for all practical purposes. When
the stationary state is reached, the concentration profile in the PET wall is
therefore linear, ranging from the value given in Eq. (12) to 0 at the outer
surface of the PET wall. The general situation is sketched in Fig. 1.
6
The flux per unit surface area is calculated through Ficks first law as
JD = D
P
CP ET
= D P ET
l
lk
(13)
where l is the thickness of the wall (250 microns). To calculate the number
of moles of CO2 that leave the bottle per unit time, one has to multiply that
expression by the total area A (800 cm2 ).
To calculate the evolution of the pressure in teh bottle, one has to relate
it to the total number of CO2 molecules in it. The relation is
NCO2 = cW VW = P
VW
kW
(14)
(15)
VW l k P ET
AD k W
(16)
Using the numerical values given here, one has ' 542 days.
The solution of Eq. (15) is an exponential of the form
P (t) = P (0)et/
(17)
Imagine you have an infinite flat surface of water from which molecules evaporate. You are asked to calculate the time-dependent water vapour concentration profile over the surface, assuming that there are initially no water
molecules in the air.
1. Write the equation that has to be solved to calculate C(z, t), where z is
the height above the surface and t is the time, as well as the boundary
conditions. Do not attempt to solve the equation yet.
2. Use dimensional analysis to find the general form of the solution.
3. Using the general form of the solution obtained from dimensional analysis, solve the equation you wrote at point 1. Express the solution in
terms of the error function erf(x) defined as
Z x
2
2
et dt
erf (x) =
0
which takes the value erf() = 1.
Solution
1. The equation is Ficks second law
2
t C = Dzz
C
(18)
(19)
C
C0
(20)
(21)
This means that the solution of the problem has to be of the type
z
C
=F
(22)
C0
Dt
which is the answer to the question.
A few comments. From the boundary conditions, we know that F (0) =
1 and F () = 0. The function F () characterises the shape of the
water-vapour concentration front. Of course at this stage, we do not
know the shape of the front. For instance, we do not known whether
it is sharp of blurry, etc. However, Eq. 22 tells already a lot. For
example, let 1/2 be the value of where F = 1/2. We dont know
the value of 1/2 but this is not important. The diffusion front arrives
9
at the height z when z/ Dt = 1/2 . In other words, the front moves
according to
z = 1/2 Dt
(23)
This tells you that if you have to wait four times longer to see the front
go two times further. We havent solved any equation yet, and look
what we already know! Isnt that beautiful?
3. Using the general form obtained from dimensional analysis
t C =
Dz 1 dF
2 Dt d
2
zz
C=
1 d2 F
Dt d 2
(24)
Introducing this is Ficks second law, one finds the following ODE
dF
d2 F
=
2
d
2 d
(25)
t2 /4
0
F () = A
e
dt + B = A erf
+B
2
0
(26)
(27)
10
Lifetime of an aerosol
11
Solution
1. Thats Ficks second law. The stationary version of it, i.e. c = 0. In
spherical coordinates, it takes the form
1 d
2 dc
r
=0
(29)
r2 dr
dr
With boundary conditions c = c0 at r = R (i.e. saturating vapour
pressure at the surface of the liquid) and c = c1 far from the particle
r where c1 is the ambient humidity.
2. The evaporation rate is by definition the number of water molecules
that diffuse away from the droplet. The rate is therefore calculated by
Ficks first law. All we have to do is chose a surface that surrounds the
droplet and integrate the flux on it
Z
Dn c dA
(30)
N=
surf ace
(31)
N = 4R D
dr r=R
3. In dimensional form, the relation we are looking for is
N = f (R, D, c0 , c1 )
(32)
There are 5 variables and 3 dimensions. Thats two dimensionless numbers. You have, however, some freedom in the way to write them. A
convenient choice is
N
c1
0 =
1 =
(33)
DRc0
c0
The general function is therefore
N
c1
=F
(34)
DRc0
c0
12
If you are really on the ball, you may have noticed at this stage that
c0 and c1 are not really two independent numbers, because you expect
their difference c0 c1 to be the important variable. You are therefore
left with only one dimensionless number, namely
0 =
N
DR(c0 c1 )
(35)
(37)
N = 4DR(c0 c1 )
(38)
c = (c0 c1 )
Introducing this in Eq. 31 leads to
(40)
13
(43)
R(0)2
2vm D(c0 c1 ))
(44)
R(t) = R(0) 1
with
=
(42)
The quantity is the time at which the particle reaches the size R = 0.
In other words, is its lifetime.
6. Let us first put all relevant data in consistent units. The partial pressure 5 kPa is converted to mol/m3 using the perfect gas law c =
P/(RG T ) where RG = 8.314 J/mol/K is the gas constant. This leads
to
c0 c1 =
5 103
(1 0.3) ' 1.39 103 mol/m3
8.314 (273 + 30)
(45)
(46)
14
15