Вы находитесь на странице: 1из 9

# MIT OpenCourseWare

http://ocw.mit.edu
1.061 / 1.61 Transport Processes in the Environment
Fall 2008
3.DiffusionofanInstantaneous,
PointRelease
This section follows on from Section 1, exploring the Gaussian
concentration profile that develops after an instantaneous, point release.
Several important properties of the Gaussian distribution (and their
utility) are emphasized. The spread of a point release in all 3
dimensions is discussed, noting that the profile of concentration along
any straight line cut through the patch will have a Gaussian distribution.
The animation shows the spread of an instantaneous, point release with
anisotropic diffusivity, highlighting the dependence of length scales of
the Gaussian distribution on the square root of D.
1
3. Diffusion of an Instantaneous Point Source
The equation of conservation of mass is also known as the transport equation, because it
describes the transport of scalar species in a fluid systems. In this and subsequent
sections we consider analytical solutions to the transport equation that describe the fate of
mass under different boundary and initial conditions. For incompressible flow with
isotropic and homogeneous diffusion coefficients, we can begin with the following form
of the transport equation. At this point we will also ignore the source/sink term.
oC oC oC oC | o
2
C o
2
C o
2
C |
(1) + u + v + w = D + +
ot ox oy oz |
|
ox
2
oy
2
oz
2
|
|
The first solutions describes the evolution of a slug of mass, M, released instantaneously
at time t = 0 into a stationary fluid (u,v,w = 0). We consider a one-dimensional system in
the x-direction, i.e. the system is uniform in y and z, such that /y = 0 and /z = 0. A
model for such a system is shown below. A slab of mass (gray) is released into a long,
narrow tube. The mass is distributed uniformly in the y-z plane and with negligible
dimension in x, such that the initial concentration is C(x) = M (x), where () is the Dirac
delta function. The transport equation for this system is given in (2).
oC o
2
C
(2) = D
2
, M released at t = 0, x = 0.
ot ox
The solution, C(x,t), can be found as follows. From (2) we anticipate that the
concentration will be a function of M, D [L
2
T
-1
], x[L], and t[T], with dimensions
expressed in brackets in terms of mass [M], length [L], and time [T]. In a one-
dimensional system, concentration has units of mass per length, C[M L
-1
]. From
dimensional analysis (e.g. Buckingham Pi Theorem) the above parameter set forms the
following dimensionless groups.
C x

(3) =f
M 4Dt

4Dt

2
The factors 4 and 4 are added now for convenience, because they appear in the final
solution. Originally, these factors were found through trial and error. Defining
=x 4Dt , and then evaluating (2) with (3),
d

df
(4) +2f =0
d d
Setting the terms within the bracket to zero provides one solution to (4),
df
(5) +2f =0 ,
d
which in turn has the solution
(6) f =A
0
exp(
2
) ,
with A
0
an arbitrary constant. (6) gives us a specific form for (3) that is a solution for (2).
A
o
x
2

## (7) C(x, t) = exp -

4Dt 4Dt
Since the total mass must be conserved over time, for all time
(8) M = Cdx ,
x=
so that A
0
= M. In summary we have the following solution.
One-dimensional, Instantaneous Point Source
M released at t = 0 and x = 0.
Initial condition, C(x) = M (x).
(9) C(x, t) =
M
exp

-
x
2

=[M/L]
4Dt 4Dt

3
Note that the above solution is correct for a one-dimensional mathematical system, in
which C[M/L]. In real (three-dimensional) systems concentration must have units of
M/L
3
. To convert the mathematical solution given in (9) to real space, divide by the
neglected dimensions, here the cross-sectional area of the system in the y-z plane, A
yz
.
(10) C(x, t) =
M
exp

-
x
2

=[M/L
3
]
A 4Dt 4Dt
yz
Recall that the same solution was derived from a statistical model of diffusion, as
described in Conceptual Model of Diffusion.
Moments of the Concentration Distribution.
It is useful to characterize the geometry of the evolving distribution of mass based on its
moments. The zeroth moment gives the total mass, M, within the cloud,
zeroth moment = M
0
= M = C(x, t)dx .
Higher moments indicate the centroid and variance of the mass distribution.
first moment = M
1
= xC(x, t)dx
second moment = M
2
= x
2
C(x, t)dx
The mean indicates the centroid of mass for the cloud,
M
1
Mean = = .
M
0
The variance indicates the width of the distribution.

(x - )
2
C(x, t)dx

2
=

=
M
2

2
.

C(x, t)dx
M
0
Evaluating the mean and variance for the distribution given in (10), we find that the
mean, = 0 , i.e. the centroid of the patch remains at x = 0 for all time. The variance is
(11)
2
= 2Dt .

|
|

|
|

|
\
j
J
4
If the coefficient of diffusion, D, is constant, the variance of the mass distribution grows
linearly in time. Taking the derivative of (11),
1 d
2
1
2
2
-
1
2
(12) D = = .
2 dt 2 t
2
- t
1
This provides a useful tool for estimating the diffusion coefficient by tracking the
variance of a diffusing tracer over time. It is useful to note that for a Gaussian
distribution, the concentration at one standard deviation from the center of mass is,
M 2Dt
(13) C(x = = 2Dt ) = C
max
(t)exp(-1/2) = 0.61C = exp - ,
A 4Dt
max
4Dt
yz
where C
max
is the maximum concentration in the patch at time t, and is given by the
bracketed coefficient to the exponential term in (13), C
max
= M / (A
yz
(4Dt)
1/2
)).
Finally, as described in Chapter 1, the standard deviation, , provides a useful length-
scale for the diffusing patch. Noting that 95% of the total mass is contained within -2
of the centroid, a conventional length-scale for a patch is
(14) L =4.
With this definition for patch size, the edges of the patch are defined by the concentration
C(x = 2) = 0.14 C
max
.
Instantaneous Point Source in Two-Dimensions
A slug of mass, M, is released at x = y = z = 0 in a stationary fluid (u = v = w = 0). The
flow conditions are uniform in z (/z = 0). The mass slug is initially distributed
uniformly in z and with negligible dimension in y and x, such that the initial
concentration is C(x, y, t=0) = M (x)(y), where () is the Dirac delta function.
- -

5
For generality we will assume anisotropic diffusion, D
x
D
y
. The transport equation for
this system is then,
oC o
2
C o
2
C
(15) = D
2
+ D
y
.
ot
x
ox oy
2
From Ficks Law and by inspection of (15), the diffusion in x (1
st
term on right hand side)
depends only on the distribution in x and the diffusion in y (2
nd
term on right hand side)
depends only on the distribution in y. This leads us to propose that the solution to (15)
may be the product of two uncoupled solutions describing the distribution in x and y.
(16) C(x, y, t) = MC
1
(x, t)C
2
(y, t) ,
where, explicitly, C
1
is not a function of y and C
2
is not a function of x. Inserting the
proposed solution form into (15),
(17)
o
(C
1
C
2
) =C
1
oC
2
+ C
2
oC
1
=D
x
C
2
o
2
C
2
1
+ D
y
C
1
o
2
C
2
2
ot ot ot ox oy
The above equation can be rewritten to isolate the functions C
1
and C
2
.
| oC
1
o
2
C
1
| | oC
2
o
2
C
2
|
(18) C
2
|
- D
x 2
|
+ C
1
|
- D
y 2
|
= 0 .
| ot ox | | ot ox |
The trivial, but useless, solution is C
1
= C
2
= 0. Alternatively each bracketed set of terms
could be zero. For example, setting the first bracketed set to zero yields,
oC o
2
C
(19)
1
= D
x 2
1
,
ot ox
which is the 1-D, diffusion equation (2). From (7), we can write the solution to (19).
2
4D
A
1
(20) C
1
(x, t) = exp
t
x
-
4D t
x x
A similar form is found for C
2
by setting the second bracketed term in (18) to be zero.
The full solution will then have the form,

2
A
1
A
2
y
2
exp
4 t D
x
D
y
4D
x
t 4D
y
t
,
x
(21) C(x, y, t) = C
1
(x, t)C
2
(y, t) =
with A
1
and A
2
constants. The product A
1
A
2
= M is found by applying the condition

6
(22) M = Cdxdy .
y= x =
Finally, the two-dimensional solution (21) has units [M L
-2
], because the third dimension
has been neglected. To recover the correct concentration in real, three-dimensional
space, divide the two-dimensional solution by the length of the missing dimension, e.g.
L
z
in the figure given above. In summary we have the following solution.
Instantaneous Point Source in Two-Dimensions
M released at t = 0 and x = y = 0.
Initial condition, C = M (x) (y)
(23) C(x, y, t) =
M
exp

-
x
2
-
y
2

= [ML
-3
].
L 4 t D
x
D
y
4D
x
t 4D
y
t

z
From (23) the profile of concentration, C(x,y) along any straight line cut through the
patch will have a Gaussian distribution. Obvious examples are the profiles C(x, y = 0)
and C(x = 0, y). This is the direct result of Ficks Law, which requires that the flux in
any direction be proportional only to the gradients in that direction. If the diffusion
coefficients are anisotropic, the cloud will grow anisotripically, increasing in length more
quickly along the axis of maximum diffusion rate. The length scale of the cloud along
any axes will be proportional to the diffusion coefficient along that axes.
(24)
L = 4 = 4 2D t
x x x
L = 4 = 4 2D t
y y y
With this length scale definition, the edge of the patch is located 2 from the centroid
along a concentration contour C = 0.14 C
max
, similar to the one-dimensional case
described above. The basic characteristics of a two-dimensional slug release are depicted
in the following animation.
Animation - Anisotropic Diffusion in Two-Dimensions - see chapter homepage
This animation depicts the diffusion of a discrete mass released at (x=0, y=0, t=0). The
diffusion is anisotropic, D
x
= 4 D
y
. As indicated in (24), the length scales grow in
proportion to the square root of the diffusion, such that the dimensions of the cloud are
anisotropic, with L
x
/L
y
= D
x
/D
y
. In this system, L
x
= 2L
y
. Note that the profiles of
concentration along the x- and y-axes are Gaussian in shape. Finally, noting that the 2
contour reaches x = L at t = 600s provides an estimate for the diffusion coefficient, D
x
.
Specifically, use the constraint L
x
= 2L at t = 600s in (24) to estimate D
x
.
The above analyses are easily extended to three dimensions.
- - -

7
Instantaneous, Point Release in Three Dimensions
M released at t = 0 and x = y = z = 0.
Initial Condition, C = M (x)(y)(z).
M

2 2 2
y x z
= [ML
-3
] (25) C(x, y, z, t) = exp
( 4 t )
3/2
D
x
D
y
D
z
4D
x
t 4D
y
t 4D
z
t
Dilution Rate:
The dilution rate of a patch may be described by the rate of change in the maximum
concentration within the patch. For a slug release the maximum concentration, C
max
(t),
occurs at the centroid and corresponds to the coefficient leading the exponential in (10),
(23), and (25).
M
1- D: C
max
=
(4D t)
1/ 2
A
yz
M
(26) 2 - D: C
max
=
L 4 t D D
z x y
M
3 - D: C =
max
( 4 t )
3/2
D
x
D
y
D
z
The three-dimensional system provides the greatest dilution rate, C
max
~ t
-3/2
, because
diffusion can occur in all directions. In the 1-D and 2-D systems, diffusive flux can only
occur in 1 or 2 directions, respectively, such that the rate of dilution is diminished to t
-1/2
and t
-1
, respectively. The difference in dilution rate is depicted below. A slug of mass is
released into a 1-D, 2-D, and 3-D system, such that at some small time the maximum
concentration in each cloud is the same. Subsequently, the maximum concentration
decays most rapidly in the least constrained system, i.e. the 3-D system.
Cmax [M L-3]
1.2
1
0.8
0.6
0.4
0.2
0
3-D
1-D
2-D
0 2 4 6 8
Time [T]
The text passage is a detailed explanation of the mathematical modeling of the diffusion of an instantaneous point release in one, two, and three dimensions. It includes the derivation of the concentration profile over time and space, and discusses the properties of the Gaussian distribution in these contexts.
no_question

the passage does not contain any explicit mathematical questions or word problems that require solving. it is primarily a descriptive and analytical text explaining the mathematical concepts and derivations related to diffusion processes.
Given a text passage, determine:
    1. If it contains a mathematical question (including word problems or calculations)
    2. If an answer to that mathematical question is provided in the text

First, describe the content structure you found in 1-2 sentences.

Then provide your analysis in the following format:
TYPE: <one of the following values>
- NO_QUESTION
- SINGLE_QUESTION
- SINGLE_QUESTION_WITH_ANSWER
- MULTIPLE_QUESTIONS
- MULTIPLE_QUESTIONS_WITH_ANSWERS
system
Вы находитесь на странице: 1из 9

# MIT OpenCourseWare

http://ocw.mit.edu
1.061 / 1.61 Transport Processes in the Environment
Fall 2008
3.DiffusionofanInstantaneous,
PointRelease
This section follows on from Section 1, exploring the Gaussian
concentration profile that develops after an instantaneous, point release.
Several important properties of the Gaussian distribution (and their
utility) are emphasized. The spread of a point release in all 3
dimensions is discussed, noting that the profile of concentration along
any straight line cut through the patch will have a Gaussian distribution.
The animation shows the spread of an instantaneous, point release with
anisotropic diffusivity, highlighting the dependence of length scales of
the Gaussian distribution on the square root of D.
1
3. Diffusion of an Instantaneous Point Source
The equation of conservation of mass is also known as the transport equation, because it
describes the transport of scalar species in a fluid systems. In this and subsequent
sections we consider analytical solutions to the transport equation that describe the fate of
mass under different boundary and initial conditions. For incompressible flow with
isotropic and homogeneous diffusion coefficients, we can begin with the following form
of the transport equation. At this point we will also ignore the source/sink term.
oC oC oC oC | o
2
C o
2
C o
2
C |
(1) + u + v + w = D + +
ot ox oy oz |
|
ox
2
oy
2
oz
2
|
|
The first solutions describes the evolution of a slug of mass, M, released instantaneously
at time t = 0 into a stationary fluid (u,v,w = 0). We consider a one-dimensional system in
the x-direction, i.e. the system is uniform in y and z, such that /y = 0 and /z = 0. A
model for such a system is shown below. A slab of mass (gray) is released into a long,
narrow tube. The mass is distributed uniformly in the y-z plane and with negligible
dimension in x, such that the initial concentration is C(x) = M (x), where () is the Dirac
delta function. The transport equation for this system is given in (2).
oC o
2
C
(2) = D
2
, M released at t = 0, x = 0.
ot ox
The solution, C(x,t), can be found as follows. From (2) we anticipate that the
concentration will be a function of M, D [L
2
T
-1
], x[L], and t[T], with dimensions
expressed in brackets in terms of mass [M], length [L], and time [T]. In a one-
dimensional system, concentration has units of mass per length, C[M L
-1
]. From
dimensional analysis (e.g. Buckingham Pi Theorem) the above parameter set forms the
following dimensionless groups.
C x

(3) =f
M 4Dt

4Dt

2
The factors 4 and 4 are added now for convenience, because they appear in the final
solution. Originally, these factors were found through trial and error. Defining
=x 4Dt , and then evaluating (2) with (3),
d

df
(4) +2f =0
d d
Setting the terms within the bracket to zero provides one solution to (4),
df
(5) +2f =0 ,
d
which in turn has the solution
(6) f =A
0
exp(
2
) ,
with A
0
an arbitrary constant. (6) gives us a specific form for (3) that is a solution for (2).
A
o
x
2

## (7) C(x, t) = exp -

4Dt 4Dt
Since the total mass must be conserved over time, for all time
(8) M = Cdx ,
x=
so that A
0
= M. In summary we have the following solution.
One-dimensional, Instantaneous Point Source
M released at t = 0 and x = 0.
Initial condition, C(x) = M (x).
(9) C(x, t) =
M
exp

-
x
2

=[M/L]
4Dt 4Dt

3
Note that the above solution is correct for a one-dimensional mathematical system, in
which C[M/L]. In real (three-dimensional) systems concentration must have units of
M/L
3
. To convert the mathematical solution given in (9) to real space, divide by the
neglected dimensions, here the cross-sectional area of the system in the y-z plane, A
yz
.
(10) C(x, t) =
M
exp

-
x
2

=[M/L
3
]
A 4Dt 4Dt
yz
Recall that the same solution was derived from a statistical model of diffusion, as
described in Conceptual Model of Diffusion.
Moments of the Concentration Distribution.
It is useful to characterize the geometry of the evolving distribution of mass based on its
moments. The zeroth moment gives the total mass, M, within the cloud,
zeroth moment = M
0
= M = C(x, t)dx .
Higher moments indicate the centroid and variance of the mass distribution.
first moment = M
1
= xC(x, t)dx
second moment = M
2
= x
2
C(x, t)dx
The mean indicates the centroid of mass for the cloud,
M
1
Mean = = .
M
0
The variance indicates the width of the distribution.

(x - )
2
C(x, t)dx

2
=

=
M
2

2
.

C(x, t)dx
M
0
Evaluating the mean and variance for the distribution given in (10), we find that the
mean, = 0 , i.e. the centroid of the patch remains at x = 0 for all time. The variance is
(11)
2
= 2Dt .

|
|

|
|

|
\
j
J
4
If the coefficient of diffusion, D, is constant, the variance of the mass distribution grows
linearly in time. Taking the derivative of (11),
1 d
2
1
2
2
-
1
2
(12) D = = .
2 dt 2 t
2
- t
1
This provides a useful tool for estimating the diffusion coefficient by tracking the
variance of a diffusing tracer over time. It is useful to note that for a Gaussian
distribution, the concentration at one standard deviation from the center of mass is,
M 2Dt
(13) C(x = = 2Dt ) = C
max
(t)exp(-1/2) = 0.61C = exp - ,
A 4Dt
max
4Dt
yz
where C
max
is the maximum concentration in the patch at time t, and is given by the
bracketed coefficient to the exponential term in (13), C
max
= M / (A
yz
(4Dt)
1/2
)).
Finally, as described in Chapter 1, the standard deviation, , provides a useful length-
scale for the diffusing patch. Noting that 95% of the total mass is contained within -2
of the centroid, a conventional length-scale for a patch is
(14) L =4.
With this definition for patch size, the edges of the patch are defined by the concentration
C(x = 2) = 0.14 C
max
.
Instantaneous Point Source in Two-Dimensions
A slug of mass, M, is released at x = y = z = 0 in a stationary fluid (u = v = w = 0). The
flow conditions are uniform in z (/z = 0). The mass slug is initially distributed
uniformly in z and with negligible dimension in y and x, such that the initial
concentration is C(x, y, t=0) = M (x)(y), where () is the Dirac delta function.
- -

5
For generality we will assume anisotropic diffusion, D
x
D
y
. The transport equation for
this system is then,
oC o
2
C o
2
C
(15) = D
2
+ D
y
.
ot
x
ox oy
2
From Ficks Law and by inspection of (15), the diffusion in x (1
st
term on right hand side)
depends only on the distribution in x and the diffusion in y (2
nd
term on right hand side)
depends only on the distribution in y. This leads us to propose that the solution to (15)
may be the product of two uncoupled solutions describing the distribution in x and y.
(16) C(x, y, t) = MC
1
(x, t)C
2
(y, t) ,
where, explicitly, C
1
is not a function of y and C
2
is not a function of x. Inserting the
proposed solution form into (15),
(17)
o
(C
1
C
2
) =C
1
oC
2
+ C
2
oC
1
=D
x
C
2
o
2
C
2
1
+ D
y
C
1
o
2
C
2
2
ot ot ot ox oy
The above equation can be rewritten to isolate the functions C
1
and C
2
.
| oC
1
o
2
C
1
| | oC
2
o
2
C
2
|
(18) C
2
|
- D
x 2
|
+ C
1
|
- D
y 2
|
= 0 .
| ot ox | | ot ox |
The trivial, but useless, solution is C
1
= C
2
= 0. Alternatively each bracketed set of terms
could be zero. For example, setting the first bracketed set to zero yields,
oC o
2
C
(19)
1
= D
x 2
1
,
ot ox
which is the 1-D, diffusion equation (2). From (7), we can write the solution to (19).
2
4D
A
1
(20) C
1
(x, t) = exp
t
x
-
4D t
x x
A similar form is found for C
2
by setting the second bracketed term in (18) to be zero.
The full solution will then have the form,

2
A
1
A
2
y
2
exp
4 t D
x
D
y
4D
x
t 4D
y
t
,
x
(21) C(x, y, t) = C
1
(x, t)C
2
(y, t) =
with A
1
and A
2
constants. The product A
1
A
2
= M is found by applying the condition

6
(22) M = Cdxdy .
y= x =
Finally, the two-dimensional solution (21) has units [M L
-2
], because the third dimension
has been neglected. To recover the correct concentration in real, three-dimensional
space, divide the two-dimensional solution by the length of the missing dimension, e.g.
L
z
in the figure given above. In summary we have the following solution.
Instantaneous Point Source in Two-Dimensions
M released at t = 0 and x = y = 0.
Initial condition, C = M (x) (y)
(23) C(x, y, t) =
M
exp

-
x
2
-
y
2

= [ML
-3
].
L 4 t D
x
D
y
4D
x
t 4D
y
t

z
From (23) the profile of concentration, C(x,y) along any straight line cut through the
patch will have a Gaussian distribution. Obvious examples are the profiles C(x, y = 0)
and C(x = 0, y). This is the direct result of Ficks Law, which requires that the flux in
any direction be proportional only to the gradients in that direction. If the diffusion
coefficients are anisotropic, the cloud will grow anisotripically, increasing in length more
quickly along the axis of maximum diffusion rate. The length scale of the cloud along
any axes will be proportional to the diffusion coefficient along that axes.
(24)
L = 4 = 4 2D t
x x x
L = 4 = 4 2D t
y y y
With this length scale definition, the edge of the patch is located 2 from the centroid
along a concentration contour C = 0.14 C
max
, similar to the one-dimensional case
described above. The basic characteristics of a two-dimensional slug release are depicted
in the following animation.
Animation - Anisotropic Diffusion in Two-Dimensions - see chapter homepage
This animation depicts the diffusion of a discrete mass released at (x=0, y=0, t=0). The
diffusion is anisotropic, D
x
= 4 D
y
. As indicated in (24), the length scales grow in
proportion to the square root of the diffusion, such that the dimensions of the cloud are
anisotropic, with L
x
/L
y
= D
x
/D
y
. In this system, L
x
= 2L
y
. Note that the profiles of
concentration along the x- and y-axes are Gaussian in shape. Finally, noting that the 2
contour reaches x = L at t = 600s provides an estimate for the diffusion coefficient, D
x
.
Specifically, use the constraint L
x
= 2L at t = 600s in (24) to estimate D
x
.
The above analyses are easily extended to three dimensions.
- - -

7
Instantaneous, Point Release in Three Dimensions
M released at t = 0 and x = y = z = 0.
Initial Condition, C = M (x)(y)(z).
M

2 2 2
y x z
= [ML
-3
] (25) C(x, y, z, t) = exp
( 4 t )
3/2
D
x
D
y
D
z
4D
x
t 4D
y
t 4D
z
t
Dilution Rate:
The dilution rate of a patch may be described by the rate of change in the maximum
concentration within the patch. For a slug release the maximum concentration, C
max
(t),
occurs at the centroid and corresponds to the coefficient leading the exponential in (10),
(23), and (25).
M
1- D: C
max
=
(4D t)
1/ 2
A
yz
M
(26) 2 - D: C
max
=
L 4 t D D
z x y
M
3 - D: C =
max
( 4 t )
3/2
D
x
D
y
D
z
The three-dimensional system provides the greatest dilution rate, C
max
~ t
-3/2
, because
diffusion can occur in all directions. In the 1-D and 2-D systems, diffusive flux can only
occur in 1 or 2 directions, respectively, such that the rate of dilution is diminished to t
-1/2
and t
-1
, respectively. The difference in dilution rate is depicted below. A slug of mass is
released into a 1-D, 2-D, and 3-D system, such that at some small time the maximum
concentration in each cloud is the same. Subsequently, the maximum concentration
decays most rapidly in the least constrained system, i.e. the 3-D system.
Cmax [M L-3]
1.2
1
0.8
0.6
0.4
0.2
0
3-D
1-D
2-D
0 2 4 6 8
Time [T]
user
The text passage is a detailed explanation of the mathematical modeling of the diffusion of an instantaneous point release in one, two, and three dimensions. It includes the derivation of the concentration profile over time and space, and discusses the properties of the Gaussian distribution in these contexts.

TYPE: NO_QUESTION

The passage does not contain any explicit mathematical questions or word problems that require solving. It is primarily a descriptive and analytical text explaining the mathematical concepts and derivations related to diffusion processes.