derivation of the advection-dispersion equation (ade) assumptions 1.equivalent porous medium (epm)...
Post on 19-Dec-2015
219 views
TRANSCRIPT
![Page 1: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/1.jpg)
Derivation of theAdvection-Dispersion Equation (ADE)
Assumptions
1. Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely fractured medium with a single network of connected fractures)
2. Miscible flow (i.e., solutes dissolve in water; DNAPL’s and LNAPL’s require a different governing equation. See p. 472, note 15.5, in Zheng and Bennett.)
![Page 2: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/2.jpg)
3. No density effects
Density-dependent flow requiresa different governing equation. SeeZheng and Bennett, Chapter 15.
Figures from Freeze & Cherry (1979)
![Page 3: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/3.jpg)
Derivation of theAdvection-Dispersion Equation (ADE)
s
hhKAQ
12
Darcy’s law:
s
h1
h2
q = Q/A
advective flux fA = q c
s
h1
h2
f = F/A
![Page 4: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/4.jpg)
s
h1
h2fA = advective flux = qc
f = fA + fD
How do we quantify thedispersive flux?
s
ccADF dDiff
12
How about Fick’s law of diffusion? where Dd is the effective
diffusion coefficient.
![Page 5: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/5.jpg)
Fick’s law describes diffusion of ions on a molecular scale as ions diffuse from areas of higher to lower concentrations.
Dual PorosityDomain
Figure from Freeze & Cherry (1979)
![Page 6: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/6.jpg)
We need to introduce a “law” to describedispersion, to account for the deviation ofvelocities from the average linear velocitycalculated by Darcy’s law.
Average linear velocityTrue velocities
![Page 7: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/7.jpg)
We will assume that dispersion followsFick’s law, or in other words, that dispersionis “Fickian”. This is an important assumption;it turns out that the Fickian assumption is notstrictly valid near the source of the contaminant.
s
ccDfD
12
where D is the dispersion coefficient.
porosity
![Page 8: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/8.jpg)
Mathematically, porosity functions as a kind of units conversion factor.
Porosity ()
for example:
q c = v c
Later we will define the dispersion coefficientin terms of v and therefore we insert now:
s
ccDfD
12
![Page 9: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/9.jpg)
Assume 1D flow
qx
s
and a line source
Case 1
![Page 10: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/10.jpg)
cvcx
hhKcqf xxA
][12
Advective flux
x
ccDf xD
12
Dispersive flux
Assume 1D flow
qx
s
D is the dispersion coefficient. It includesthe effects of dispersion and diffusion. Dx is sometimeswritten DL and called the longitudinal dispersion coefficient.
porosityCase 1
![Page 11: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/11.jpg)
Assume 1D flow
qx
s
and a point source
Case 2
![Page 12: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/12.jpg)
fA = qxcAdvective flux
Dx represents longitudinal dispersion (& diffusion);Dy represents horizontal transverse dispersion (& diffusion);Dz represents vertical transverse dispersion (& diffusion).
)(12
x
ccDf xDx
)(12
z
ccDf zDz
Dispersive fluxes )(12
y
ccDf yDy
![Page 13: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/13.jpg)
Figure from Freeze & Cherry (1979)
Continuous point source
Instantaneous point source
Averagelinearvelocity
center of mass
![Page 14: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/14.jpg)
Figure from Wang and Anderson (1982)
InstantaneousPoint Source
transversedispersion
longitudinal dispersion
Gaussian
![Page 15: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/15.jpg)
Derivation of the ADE for1D uniform flow and 3D dispersion(e.g., a point source in a uniform flow field)
f = fA + fD
Mass Balance:Flux out – Flux in = change in mass
vx = a constant vy = vz = 0
![Page 16: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/16.jpg)
Porosity ()
There are two types of porosity in transport problems:total porosity and effective porosity.
Total porosity includes immobile pore water, which contains solute and therefore it should be accounted for when determining the total mass in the system.
Effective porosity accounts for water in interconnected pore space, which is flowing/mobile.
In practice, we assume that total porosity equals effectiveporosity for purposes of deriving the advection-dispersion eqn.See Zheng and Bennett, pp. 56-57.
![Page 17: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/17.jpg)
Definition of the Dispersion Coefficientin a 1D uniform flow field
vx = a constantvy = vz = 0
Dx = xvx + Dd
Dy = yvx + Dd
Dz = zvx + Dd
where x y z are known as dispersivities. Dispersivity is essentially a “fudge factor” to account for the deviations of the true velocities from the average linear velocities calculated from Darcy’s law.
Rule of thumb: y = 0.1x ; z = 0.1y
![Page 18: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/18.jpg)
t
c
x
cv
z
cD
y
cD
x
cD zyx
2
2
2
2
2
2
ADE for 1D uniform flow and 3D dispersion
No sink/source term; no chemical reactions
Question: If there is no source term, how does the contaminant enter the system?
![Page 19: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/19.jpg)
t
c
x
cv
x
cD
2
2Simpler form of the ADE
Uniform 1D flow; longitudinal dispersion;No sink/source term; no chemical reactions
There is a famous analytical solution to this form of the ADE with a continuous line source boundary condition. The solution is called the Ogata & Banks solution.
Question: Is this equation valid for both point and line source boundaries?
![Page 20: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/20.jpg)
Effects of dispersion on the concentration profile
(Zheng & Bennett, Fig. 3.11)
no dispersion dispersion
(Freeze & Cherry, 1979, Fig. 9.1)
t1 t2 t3 t4
![Page 21: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/21.jpg)
Effects of dispersion on the breakthrough curve
![Page 22: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/22.jpg)
Figure from Wang and Anderson (1982)
InstantaneousPoint Source
Gaussian
![Page 23: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/23.jpg)
Breakthroughcurve
Concentrationprofile
long tail
![Page 24: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/24.jpg)
Figure from Freeze & Cherry (1979)
Microscopic or local scale dispersion
![Page 25: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/25.jpg)
Macroscopic Dispersion (caused by the presence of heterogeneities)
Homogeneous aquifer
Heterogeneousaquifers
Figure from Freeze & Cherry (1979)
![Page 26: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/26.jpg)
Dispersivity () is a measure of the heterogeneity present in the aquifer.
A very heterogeneous porous mediumhas a higher dispersivity than a slightlyheterogeneous porous medium.
![Page 27: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/27.jpg)
Dispersion in a 3D flow field
![Page 28: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/28.jpg)
x
z
x’
z’
global local
Kxx Kxy Kxz
Kyx Kyy Kyz
Kzx Kzy Kzz
K’x 0 0
0 K’y 0
0 0 K’z
[K] = [R]-1 [K’] [R]
K =
![Page 29: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/29.jpg)
z
hK
y
hK
x
hKq
z
hK
y
hK
x
hKq
z
hK
y
hK
x
hKq
zzzyzxz
yzyyyxy
xzxyxxx
![Page 30: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/30.jpg)
Dispersion Coefficient (D)
D = D + Dd
Dxx Dxy Dxz
Dyx Dyy Dyz
Dzx Dzy Dzz
D =
In general: D >> Dd
D represents dispersion Dd represents molecular diffusion
![Page 31: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/31.jpg)
z
cD
y
cD
x
cDf
z
cD
y
cD
x
cDf
z
cD
y
cD
x
cDf
zzzyzxDz
yzyyyxDy
xzxyxxDx
![Page 32: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/32.jpg)
In a 3D flow field it is not possible to simplify the dispersiontensor to three principal components. In a 3D flow field, we must consider all 9 components of the dispersion tensor.
The definition of the dispersion coefficient is more complicated for 2D or 3D flow. See Zheng and Bennett, eqns. 3.37-3.42.
Dx = xvx + Dd
Dy = yvx + Dd
Dz = zvx + Dd
Recall, that for1D uniform flow:
![Page 33: Derivation of the Advection-Dispersion Equation (ADE) Assumptions 1.Equivalent porous medium (epm) (i.e., a medium with connected pore space or a densely](https://reader035.vdocument.in/reader035/viewer/2022062714/56649d365503460f94a0eeea/html5/thumbnails/33.jpg)
General form of the ADE:
Expands to 9 terms
Expands to 3 terms
(See eqn. 3.48 in Z&B)