MIT OpenCourseWare

http://ocw.mit.edu

1.061 / 1.61 Transport Processes in the Environment

Fall 2008

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
Answer 3.1

Estimate the coefficient of diffusivity within this region of atmosphere.

Hint 1 - Can you assume uniform concentration in any direction?

The gas is effectively released into an unbounded domain, so that one cannot expect it to
mix rapidly to uniform conditions in any direction, i.e. ∂C/∂y ≠ 0, ∂C/∂z ≠ 0, ∂C/∂y ≠ 0.
All three dimensions must be retained in the governing equation.

Hint 2 - What assumption can you make about the air currents?
The atmosphere is stagnant, which implies that there are no ambient air currents. You
may assume that u=v=w=0. The passage of each plane will create air movement. This
movement may locally and temporarily enhance the dilution of the cloud. However, this
effect dies out within a few minutes, so that over an hour time frame, the effect is
negligible.

Hint 3 - What assumption must you make about the coefficient of diffusion?
Since you are only given one value of concentration, you can solve uniquely for only one
value of diffusion. Thus, to make the problem tractable, you must assume that the
diffusion is isotropic.

With the above assumptions, the transport equation reduces to,

∂C  ∂2 C ∂2 C ∂2 C 
=D  + + .
∂t  ∂x 2 ∂y 2 ∂z 2 

Letting the initial release point be at (x, y, z)=0, the initial condition can be written,

C(x,y,z,t=0) = M δ(x) δ(y) δ(z).

The solution to the above equation and initial condition is Equation 25 in Chapter 3.

M  x 2 + y 2 + z2 
C(x,y,z, t) = 
exp -  .
( 4 π D t )3/2  4Dt 

The maximum concentration within the cloud is given in Equation 26 in Chapter 3.

M
C MAX = .
( 4 π D t )3/2

Using the information given in the problem statement, we can find D,

D = (M C MAX ) 4πt = (10kg 3 × 10 kgm )

2/3 -5 −3 2/3 2 −1
(4 π 3600 s) = 0.1m s
Answer 3.2

What governing equation describes the evolution of the gas concentration in the hall?
For isotropic diffusion, the governing equation is:

∂C ∂C ∂C ∂C  ∂2 C ∂2 C ∂2 C 
+u +v +w =D 2 + 2 + 2
∂t ∂x ∂y ∂z  ∂x ∂y ∂z 

Align the coordinate x to the length of the hallway and define the y-z plane as the cross-section.
Since the gas "mixes rapidly" in the vertical and horizontal, we assume ∂C/∂y = ∂C/∂z = 0. The
problem statement gives no information about air currents in the hallway, so we assume they are
negligible, u = 0. With these assumptions, the governing equation is reduced to,

∂C ∂2 C
=D 2 .
∂t ∂x

If the spill location is taken as x = 0, the initial condition is then, M = δ(x)

At what time after the spill do you smell the gas?

The governing equation and initial condition above describe an instantaneous, point release
diffusing in one-dimension. The concentration field is described by Equation 10 in chapter 3.

M  x2 
exp  -
3
3.10 C(x, t) = = [M/L ]
A yz 4πDt  4Dt 

Use 3.10 to find the time at which the concentration at your door, C(x=20m, t) is 10 µgl-1. As
time appears in both the exponential and leading terms, it is simpler to use a graphical solution.

µg/l Concentration at X = 20m

16
14
12
10
8
6
4
2
0
0 200 400 600 800 1000 1200 1400
TIME [SECONDS]

You will smell the gas at your office approximately 1100 seconds after the spill.
When does the smell, as perceived by humans, disappear from the hallway?
To answer this question we need information on the end conditions of the hallway. Lets first
consider that the hallway is open at both ends, so that odor can diffuse beyond the end of the
hall. Then, we need only consider the evolution of the maximum concentration, located at the
site of the spill, x = 0. We seek the time for which

M −1
C max (t) = < 10µgl ,
Ayz 4πDt

or simply

 10g 2
t >  
 2m × 3m × 0.01gm −3 
(4 π × 0.05m s ) = 44,232 seconds = 12.3 hrs
2 −1

At this time, the length of the cloud will be 4 σ = 4 2Dt = 266m , indicating that the cloud has
diffused beyond the length of the hallway. The above time scale is correct, only if the hallway is
open at both ends. If the hallway is shut off by fire doors at both ends, then in reality the cloud
cannot diffuse beyond the length of the hall. Under these conditions, the final concentration in
the hallway set the maximum possible dilution, which is determined by distributing the total
mass released over the total volume of the hallway (2m x 3m x 100m).

C final = 10g (2m × 3m × 100m ) = 0.016gm = 16µgl .

−3 −1

Since the final concentration is above the detection limit, the smell will not disappear until the
fire doors are opened.
Answer 3.3
Hint 1: Make sketch that defines the diameter of each cloud

a. While both clouds are fully visible (C > 10-g l-1), which cloud will appear larger, and by
how much?
Since the size of the cloud increases in proportion to the diffusion coefficient, the blue cloud will
grow more rapidly, and thus appear bigger, than the red cloud. Specifically, the length scale of
each cloud, as defined in chapter 3, equation 26, is

L B = 4σ B = 4 2D B t and L R = 4σ R = 4 2DR t .
The ratio of dye cloud diameters is then, L B L R = DB D R = 2 .

b. At what time and at what location will the two dye clouds first appear to touch?
Hint 2: Simplify governing equation with assumptions.
Begin with the full transport equation that governs the evolution of both dye drops.

∂C ∂C ∂C ∂C ∂ ∂C ∂ ∂C ∂ ∂C
+u +v +w = Dx + Dy + Dz ±S
∂t ∂x ∂y ∂z ∂x ∂x ∂y ∂y ∂z ∂z
(1) (2) (3) (4) (5) (6) (7) (8)

Simplifying assumptions:
• The dyes do not interact, so we can drop the source/sink term 8.
• The fluid is stagnant, thus u=v=w=0, and we can drop terms 2,3,4
• Molecular diffusion is homogeneous and isotropic, Dx=Dy=Dz=D, and D ≠ f(x,y,z).
This reduces the diffusion terms to: D (∂2C/dx2 + ∂2C/dy2 + ∂2C/dz2).
• The dyes mix rapidly between the plates, such that ∂C/∂z = 0, and we can drop term 7.

∂C  ∂2 C ∂2 C 
Simplified governing equation: = D  2 + 2 .
∂t  ∂x ∂y 

Initial conditions: Blue Dye: CB(x,y,t=0) = M δ(x) δ(y)

Red Dye: CR(x,y,t=0) = M δ(x-5) δ(y)

The diffusion of each dye is thus described by the two-dimensional, instantaneous, point release
solution. See equation 23 in Chapter 3, with Dx=Dy=D. Note, here LZ is the plate gap, 5mm.
 M  x2 + y2 
C B(x,y, t) = exp - 
L z 4 π DB t  4DB t 

M  (x - 5)2 + y2 
C R (x,y,t) = 
exp - 
L z 4π DR t  4DRt 

Note that the exponential is written to be one at the center of the red cloud, i.e. at (x=10,y=0),
with the position given in units of centimeters.

Hint 3: Write mathematical criteria for condition when clouds first touch
From geometry, the two clouds will first touch along the line y = 0. We wish to find the time, t,
and location, x, for which CB(x, y=0, t) = CR(x, y=0, t) = 10 g/l. This could be tackled
analyticaly by first finding the position xBR(t) at which CB=CR. Then solve for t using this
position to constrain x, e.g. solve CB(x = xBR, y=0, t). However, a simpler approach is to graph
CB and CR in an interactive graphing package, such as Excel, and then vary time until the
intersection of the two concentration curves lies at C = 10 g/l

Solution: - At t = 20500s, the intersection of the blue and red concentration curves corresponds
to C = 10 g/l and is located at x = 3.1 cm. The clouds will first appear to touch and x = 3.1 cm.
Two-dimensional Dye Diffusion
Concentration in g/l
Time, t = 20500s.
20

15

10

0
-5 0 X [ cm ] 5 10

Make a rough estimate of the location using your result from part a?
Based on the definition used in a. and the definition sketch shown above, the two clouds first
appear to touch when, (LB/2)+(LR/2) = 5. Additionally, LR = (LB/2), so the edge of the blue
cloud will be at x = (LB/2) = 5/1.5=3.3cm, when it first touches the red cloud.
c. At what time will the line connecting the release points be completely purple?
Hint 4: - Define a mathematical criteria for this to occur?
This condition requires that CB and CR > 10 g/l in the region x = 0 to 10 cm. Use the
spreadsheet created in part b to interrogate the concentration field over a range of time.

Solution - Graph CB and CR in an interactive graphing package, such as Excel and vary time
until the above criteria is met. You will find that the criteria is never met. Between 0<x<5 cm,
when CR ≥ 10g/l then CB < 10 g/l, and when CB≥10g/l, CR < 10 g/l.

Two-dimensional Dye Diffusion

Concentration in g/l
Time, t = 660000s.
30
25
20
15
10
5
0
-5 0 X [ cm ] 5 10

Two-dimensional Dye Diffusion

Concentration in g/l
Time, t = 170000s.
30
25
20
15
10
5
0
-5 0 X [ cm ] 5 10