HW1 Solution
HW1 Solution
HW1 Solution
(0.0.1)
a) Write all boundary and initial conditions needed to solve this equation for C(z, t).
b) Find the steady state solution (infinite time) for C(z).
c) Plot the initial and steady state solutions on a C versus z plot. Draw a plausible guess for an
intermediate solution (at some time between t=0 and infinity).
Solution:
a) The boundary condition is defined by time (t) and space (z ). From the information given,at t=0,
before the diffusion event, the concentration within the barrier is 0mol/L. So we have:
C(z, 0) = 0, z > 0.
(0.0.2)
At any given time, the concentration in the pool and ground water is constant, we have:
C(0, t) = 2mol/L
(0.0.3)
C(50cm, t) = 0mol/L
(0.0.4)
(0.0.5)
dC(z)
d2 C(z)
d dC(z)
=D
=D (
)=0
2
dt
dz
dz dz
(0.0.6)
d dC(z)
(
)dz =
dz dz
1
Z
0dz
(0.0.7)
Direction of diffusion
-
Pit
Barrier
c=2mol/L
Groundwater
c=0mol/L
z = 0cm
z = 50cm
Z
dC(z)
dz = C1 dz
dz
C(z) = C1 z + C2
(0.0.8)
(0.0.9)
Therefore,
C1 = 0.04molL1 cm1
and C2 = 2molL1 .
(0.0.10)
(0.0.11)
C(z)
C(z)
t=0
2mol/L
steady
6
s
2mol/L
state
6
s
@
@
@
@
@
@
@
@
@
s -
0mol/L
0cm
@
@
@s -
0mol/L
0cm
50cm
50cm
Figure 2: left) initial concentration profile, right) steady state concentration profile.
When second order derivative is positive, the curve concaves up ^ instead of concaves down _ .
An reasonable intermediate state solution should has a similar shape with the curves shown in Figure
3.
C(t, z)
6
2mol/L s
intermediate states
t increases
s -
0mol/L
0cm
50cm
Problem 2
A balloon, initially at a volume of 0.3 cm3 and filled with an ideal gas at a pressure of 1 atm, is to be
filled by flowing in helium at a rate of 1 cm3 /s and a pressure of 25 atm. The balloon will expand as
helium enters to maintain a pressure inside the balloon of 1 atm. Helium may be treated as an ideal
gas, and both the flow in and the temperature in the balloon are maintained at 25 C).
a) If the inlet of the balloon has an area of 0.7 mm2 , determine the velocity (in cm/s) and molar flux
(in moles/cm2 s) of helium entering the balloon.
b) Once filled to a volume of 1200 cm3 , the balloon leaks due to the diffusion of helium through the
balloon material. The flux of helium across the balloon surface can be expressed as kC, where C is
the concentration of gas inside the balloon and k is a mass transfer coefficient with a value of 1x10-7
cm/s. Determine how long it takes for the balloon to shrink to a volume of 600 cm3 . You may treat
the balloon as spherical.
Useful expressions (where r is the sphere radius):
Volume of a sphere = 4/3 r3
Area of a sphere = 4 r2
dr3 /dr = 3r2
c) Write the mole balance that would be used to re-examine the filling process if the diffusion of
helium out of the balloon was considered during filling. Do you anticipate the filling time changing
substantially due to helium diffusion? You do not need to resolve for the filling time.
Solution:
a) From the problem statement, we have:
Fin = 1.0cm3 /s,
Pin = 25atm,
A = 0.7mm2 ,
T = 298K(25 C).
(0.0.13)
The velocity:
v=
1.0cm3 /s 100mm2
Fin
=
= 142.9cm/s.
A
0.7mm2 cm2
(0.0.14)
Pin
142.9cm/s 25atm
=
= 0.146mol cm2 s.
RT
82.06cm3 mol1 K 1 atm 298K
(0.0.15)
b) Since the balloon is filled, the inlet term in molar balance is zero. The outlet flux of helium is
expressed by kC. The molar balance is:
dCV
= Jdif f usion A = kCA.
(0.0.16)
dt
We know the pressure and temperature remain constant inside the balloon, from the ideal gas law,
C =P /RT, the concentration C stays constant. Thus, from (0.0.16), we have:
dV
= kA.
dt
(0.0.17)
Both V and A are dependent on the radius of the sphere r. With V = 4/3 r3 and A = 4 r2 . Plug
these expressions into equation (0.0.17):
d 4 3
( r ) = k4r2 .
dt 3
(0.0.18)
d 4 3
4 dr3
4 dr3 dr
4
dr
dr
( r ) =
=
= 3r2 = 4r2 = k4r2
dt 3
3 dt
3 dr dt
3
dt
dt
(0.0.19)
dr
= k.
(0.0.20)
dt
The initial volume V0 =1200cm3 , the final volume V =600cm3 . So we can calculate the initial and
final diameter of the sphere (r 0 and r ):
r0 = (
3V0 1/3
) = 6.59cm,
4
r = (
3V 1/3
) = 5.23cm.
4
kdt
r|rr0 = kt|0
r0 r
.
k
(0.0.21)
(0.0.22)
The time it takes for the balloon to shrink to a volume of 600cm3 is:
=
r0 r
= 1.36 107 = 157.4 days.
k
(0.0.23)
c) If we consider the helium gas flowing in, the mole balance becomes:
dCV
= Fin Cin Jdif f usion A = Fin Cin kCA.
dt
(0.0.24)
(0.0.25)
k
=
(
|{z}
2P
2
dt |4r
4r
{z } outlet
(0.0.26)
inlet
The change of radius comes from both inlet helium gas and outlet diffusion. However, the inlet
contribution is much larger than the outlet. So the filling time wont change substantially due to
helium diffusion.
Another way to understand this is from class example, it takes less than 1 second to fill the balloon.
From problem 2b, the time it takes to allow half of the helium gas diffuse out of the balloon is 157.4
days. So diffusion is a much slower process.