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 4.1 - Describe the concentration field inside the box from t = 0 to t = L2/D.

Hint 1 - When will the mass be mixed uniformly in the vertical?

From equation 4.24 a mass released mid-way between two parallel boundaries a distance
LZ apart will be mixed to a uniform concentration between those boundaries in time t =
LZ2/4D, where LZ is the distance between the boundaries. The time for the concentration
to become well-mixed in the vertical is then, t = (0.01L)2/D = 0.0001 L2/D.

Hint 2 - Estimate when the mass will reach each vertical wall in the box.
We estimate the time that the cloud will first touch a boundary based on the standard
deviation of the mass distribution. The cloud will touch the top and bottom walls when,
3σ = 3 2Dt = L, or at the time t = L2 18D . The cloud will touch the far wall when
3σ = 3 2Dt = 99L, or at the time t ≈ 545L2 D .

Hint 3 - How will each boundary impact the solution in the time t = 0 to L2/D?
Use the time scales determined in hint 1 and 2. Because the concentration field becomes
uniform in the vertical very rapidly, within 1/10,000th of the time of interest, we will
assume that the concentration is instantly uniform in z, i.e. ∂C/∂z = 0 for all time. In the
time of interest, the cloud will never reach the far walls (545L2/D >> L2/D), so these
boundaries do not impact the solution. The cloud will reach the near walls (L2/18D <
L2/D), and a correction must be made to satisfy the no flux condition at these walls.

Hint 4 - Place images sources to satisfy the no-flux boundary condition.

Image sources are needed at the following locations. Image 2 balances the loss of mass
from Image 3 across the x-axis and the loss of mass from Image 1 across the y-axis.

Solution - The concentration in the box is described by a superposition of two-

dimensional, instantaneous, slug releases (equation 3.23) at each of the above sources.
M
C(x, y, t) = •
L z 4π Dt
(x − L)2 + (y − L)2 (x − L) 2 + (y + L) 2 (x + L)2 + (y + L)2 (x + L) 2 + (y − L) 2
(exp(- ) + exp(- ) + exp(- ) + exp(- ))
4Dt 4Dt 4Dt 4Dt

real source image 1 image 2 image3

Solution 4.2

(a) An image source is needed at x = -150 m to account for the no-flux boundary at the
end of the canal. The concentration of fuel in the canal is described by the superposition
of these two instantaneous, 1-D sources.

M   − x2   − (x + 150) 2 
C ( x, t ) = exp  + exp 
A 4πDt   4Dt   4Dt 
REAL IMAGE

Therefore, at x = -50 m, t = 10 hrs (= 36000 s),

1000 g   − (−50 m) 2   − (100 m) 2 

C=  exp 2 −1

 + exp
 4(0.01 m 2 s −1 )(36000 s) 
5 m2 4π (0.01 m 2 s −1 )(36000 s)   4(0.01 m s )(36000 s)   
= 2.97 gm -3 (0.176 + 0.001) = 0.53 gm - 3

(b) This is best solved graphically. The plot of:

1000 g   − (−50 m) 2   − (100 m) 2 
C(-50 m, t) = exp
   2 −1

 + exp 2 −1

5 m 2 4π (0.01 m 2 s −1 )t   4(0.01 m s )t   4(0.01 m s )t 
is shown below.

1.4
1.2
1
C (g/m3)

0.8
0.6
0.4
0.2
0
0.E+00 2.E+05 4.E+05 6.E+05 8.E+05 1.E+06
t (s)

From this plot, we can see that the maximum concentration at the house is approximately
1.3 gm-3 (when t ≈ 2.4 x 105 s = 67 hrs)
(c) The easiest ways to solve the equation in (b) for C = 0.2 gm-3 are graphically, or by
trial and error. As the above plot lacks detail in the early stages, we will use trial and
error in our spreadsheet to solve
1000 g   − (−50 m) 2   − (100 m) 2 
0.2 gm =-3
exp 2 −1

 + exp 4(0.01 m 2 s −1 )t 
5 m 4π (0.01 m s )t  
2 2 −1 4( 0 . 01 m s )t   

for t. This yields a solution of t = 21060 s = 5.8 hrs.

Compare this solution to that of Problem 1.5 to see the effect that the no-flux boundary
has on the concentration of fuel observed at your house.

(d) An image sink is needed at x = -150 m to account for the perfectly absorbing
boundary at the end of the canal (see Chapter 4 notes, p 4)

M   − x2   − ( x + 150) 2 
C ( x, t ) = exp  − exp 
A 4πDt   4 Dt   4 Dt 
REAL IMAGE

Repeating (a) œ (c) is relatively simple, the results being:

• C(x = -50 m, t = 36000 s) = 0.52 gm-3 (i.e. at this short time, the boundary has little/no
effect).

• Cmax (x = -50 m) = 0.82 gm-3 (i.e. at t ≈ 83000 s = 23 hrs).

0.9
0.8
0.7
0.6
C (g/m3)

0.5
0.4
0.3
0.2
0.1
0
0.E+00 2.E+05 4.E+05 6.E+05 8.E+05 1.E+06
t (s)
 • C = 0.2 gm-3 when t = 21060 s. This is the same answer as in (c) œ before the cloud
encounters the boundary (i.e. for 2 2Dt < 75 m ⇒ t < 70300 s ), the boundary has
little/no effect.

Answer 4.3
a) Estimate the time scale, T, at which the dye becomes uniformly distributed in x.
The diffusion of dye released at the end wall (top figure) will be similar to the diffusion of dye
released mid-way between end walls placed twice as far apart (bottom figure). This is because
the diffusion in both systems proceeds as a function of exp(-x2/4Dt), which is symmetric about x
= 0. From this similarity we expect that the dye will be well mixed in both systems in the time-
scale already established for the bottom system. Specifically, from eq. 4.24 applied to the
bottom system, T =(2L)2/4D= L2/D= 107sec.

b) Confirm your estimate by plotting C(x) at the times t = T/10, T/4, T/2, T.
Theoretically an infinite number of images is needed to satisfy a no-flux boundary condition at
two parallel boundaries. For simplicity we initially consider only one image position for each
boundary. The boundary at x=0 requires image I1 at x=0, i.e. co-located with the source. The
boundary at x=L requires image I2 at x=2L for the real source, as well as image I3 for the image
source I1. The concentration field is,

M
C(x, t) =
A yz 4π Dt
(2exp(-x 2
)
4Dt ) + 2exp(-(x - 2L)2 4Dt )
 Real + I1 I2+I3

0.200
T = L^2/10D
0.150
C[mg/cm3] T = L^2/4D
0.100
T = L^2/2D
0.050
0.000 T = L^2/D

0 5 10
X [cm]

The solution indicates that three images are not sufficient, because mass is not conserved in the
time of interest. The final concentration should be 1mg/(10cm3) = 0.1 mg cm-3, but the above
solution yields 0.66 mg cm-3 at t =L2/D. In addition, the gradient of concentration should be zero
at the boundaries, ∂C/∂x = 0, to satisfy the no-flux condition. This condition is not met at x =
10cm. For comparison, we now consider a solution with seven images, each denoted by I.

€ I1 balances the real source at boundary x = 0.

I2 and I3 balance the real source and I1 at boundary x = L.
I4 and I5 balance I2 and I3 across the boundary at x = 0.
I6 and I7 balance I4 and I5 across the boundary at x = L.

C(x, t) =
M
A yz 4π Dt
(2exp(-x 2 4Dt ) + 2exp(-(x - 2L)2 4Dt ) + 2exp(-(x + 2L)2 4Dt ) + 2exp(-(x - 4L)2 4Dt ) )
 Real+I1 I2 + I3 I4 + I5 I6+I7

0.200

T = L^2/10D
C[mg/cm3] 0.150

T = L^2/4D
0.100
T=L^2/2D

0.050 T=L^2/D

0.000
0 5 X [cm] 10

With seven images the concentration is correct at T = L2/D. At longer times mass will still be
lost due to unbalanced images and the concentration will decline. However, longer times hold
€
no interest, because once the system is well mixed in x (∂C/∂x = 0 for all x), the solution may be
ended as no further evolution of the true profile will occur. Note, the shape of the concentration
profiles are similar with 2 and 6 images, and both solutions indicate a uniform distribution in x is
achieved by T = L2/D.
 Answer 4.4:
a) Estimate the time at which the concentration at A and A’ begin to diverge?
The concentrations at A and A’ diverge when the boundary impacts the solution in System 1.
This occurs when the diffusing cloud reaches the boundary. Estimate this time by equating the
edge of the cloud with the length scale 3 σ . That is, the cloud will touch the boundary when
2
3σ = 3 2Dt = 50cm, such that t = (50cm ) 18 •2cm 2 s−1 = 70 seconds .

b) What is the final concentration at A (A’), and when is this concentration achieved?
The final concentration in System 1 will be C = (100g)/ (1m x 1m x 0.1m) = 1000 ppm.
It is achieved when the mass is fully mixed across the domain. Using the largest dimension to
estimate this time-scale, t = (L)2/(4 D) = 1250 s. As noted in the text, this is a conservative
estimate, and t = (L)2/(8 D) = 625 s, is also a reasonable estimate. Because System 2 is
unbounded, infinite dilution is possible and the final concentration is C = 0 ppm, but
theoretically this will take infinite time. Because the probe has a detection limit of 10 ppm, zero
concentration will be recorded for any concentration C < 10 ppm, which occurs in a finite time.

c) Describe the evolution of the concentration field in each system, i.e. C(x,y,z,t).
In both systems the concentration is uniform in z by t = (10cm)2/(4x2cm2s-1)=12.5 s after release.
Because this is short relative to other time scales of interest (1250 sec), we neglect the 3-D phase
and use a two-dimensional solution. For System 1, an infinite number of image sources is
needed to satisfy the no-flux boundaries.

M   x2 + y2  ∞  x 2 + (y ± 2nL)2  ∞  (x ± 2nL)2 + y 2 
C1(x, y, t) = exp−  + ∑ exp-  + ∑ exp- 
L z 4π Dt   4Dt  n=1  4Dt  n=1  4Dt 
real source images along y axis images along x axis

M ∞  (x ± 2nL)2 + (y ± 2nL)2 
+ ∑ exp - 
L z 4π Dt n=1  4Dt 
corner images

The images in the second line are located at the corners

(see sketch), to offset the losses of y-axis images across
€ the boundaries at x ± L/2; and the loss of x-axis images
across the boundaries at y = ± L/2. In practice an infinite
number of images is not needed. For time less than
required to reach a well-mixed condition between the
boundaries (t < L2/4D), three images per boundary is
sufficient to approximate the full solution with infinite
images. Beyond this time, the concentration is steady
and uniform, and the detailed solution above is no longer
needed. System 2 is described by a simple two-
dimensional slug-release,

M  x2 + y2 
C2(x, y, t) = exp - .
L z 4πD t  4Dt 