Atmospheric Dispersion (AD)
Atmospheric Dispersion (AD)
Atmospheric Dispersion (AD)
Matus Martini
1
Outline
• Eulerian approach
• Lagrangian approach
eqns for mean concentration, solutions for
instantaneous and continuous source
• Gaussian plume eqn
AD parameterizations (P-G curves), plume
rise
2
Air pollution dispersion models
• Box model air pollutants inside the box are
homogeneously distributed
• Gaussian model is perhaps the oldest (circa 1936)
• Lagrangian model - statistics of the trajectories of a
large number of the pollution plume parcels.
• Eulerian model - fixed three-dimensional Cartesian grid
• Dense gas model
• Hybrids (Plume in Grid model)
3
Leonhard Euler (1707-1783) v. Joseph Louis Lagrange (1736-1813)
4
Lagrange v. Euler
•PROS over Eulerian models:
– no Courant number restrictions
– no numerical diffusion/dispersion
– easily track air parcel histories
– invertible with respect to time
•CONS:
– need very large # points for statistics
– inhomogeneous representation of domain
– convection is poorly represented
– nonlinear chemistry is problematic
c i ci
2
( u j ci ) Di
t x j x j x j
Ri ( c1 , ... , c N , T )
S i ( x1 , x 2 , x 3 , t )
8
Eulerian approach
10
Lagrangian approach
• behavior of representative fluid particles
• consider a single particle located at location x’ at time t’ in a
turbulent fluid
• trajectory of the subsequent motion: X[x’,t’;t] at any later time t
• probability density function
( x, t ) Q( x , t | x' , t' ) ( x' , t' ) dx'
11
Lagrangian approach
• Ensemble mean concentration:
c( x , t )
• General formula for the mean concentration:
i ( x , t )
t
c( x , t ) Q( x , t | x0 , t 0 ) c( x 0 , t 0 ) dx0 Q( x , t | x' , t' )S ( x' , t' )dt' dx'
particles present at t0 added from sources between t0 tand
0
t
12
• Eulerian statistics are readily measurable (fixed network of
anemometers)
• can include detailed chemical mechanisms
• serious mathematical obstacle: closure problem
13
Mean concentration – K theory
• Eulerian approach: approx AD eqn
• Molecular diffusion is negligible compared with turbulent
diffusion
2 c
Di u j' c'
x j x j x j
• linearization: incompressible atmosphere
• Ri almost always nonlin, the most obvious approx:
Ri ( c1 , ... , c N ,T ) Ri ( c1 , ... , c N , T )
c i c ci
uj i K jj Ri ( c1 , ... , c N , T ) S i ( x, t )
t x j x j x
j
• reaction processes are slow compared w/turbulent transport
14
• distribution of sources is “smooth” (violated near strong isolated sources)
Mean concentration – statistical theory
• Lagrangian approach: stationary, homogeneous Gausian flow
• In highly idealized example: u(t) is a random variable
depending only on time, and is stationary, Gaussian
random proces, u(t) pdf 1 ( u u )2
pu ( u ) exp 2
2 u 2 u
• stationarity implies that the statistical properties of u at
two different times depend only on t–and not on t and
individually, transition probabilty density
Q( x , t | x' , t' ) Q( x x' , t t' )
• Then mean concentration is itself Gaussian (!):
1 ( x u t )2
c( x , t ) exp
2 x ( t ) 2 2 ( t )
x
15
Instantaneous point source
• Eulerian approach
• eddy diffusivities Kxx, Kyy, Kzz = const
c c 2 c 2 c 2 c
u K xx K yy K zz
t x x 2
y 2
z 2
c( x , y , z , 0 ) S ( x ) ( y ) ( z )
c( x , y , z , t ) 0 x , y , z
Solution:
S ( x u t )2 y 2
z 2
c( x , y , z , t ) exp
8 t
3/2
K xx K yy K zz 4 K t 4 K t 4 K t
xx yy zz
16
Instantaneous point source
• If we define:
x2 2 K xx t y2 2 K yv t z2 2 K zz t
17
Continuous source, steady state
• Lagrangian approach
• Source began emitting at t = 0, mean concentration
achieves a steady state (independent of time), and source
strength q in [g s-1]: S ( x , y , z ,0 ) q ( x ) ( y ) ( z )
t
c( x , t ) lim c( x , t ) lim
t t Q( x , t | 0 , t' ) q dt'
0
Q( x , t | 0 , t' ) Q( x' , t t' | 0 , 0 )
• Eulerian approach
c 2 c 2 c 2 c
u K q ( x ) ( y ) ( z )
x x 2
y 2
z 2
q u ( r x )
where r x y z
2 2 2 2
Solution: c( x , y , z , t ) exp
4 K r 2K
2 2 K yy x 2 2 K zz x
Lagrangian and Eulerian expressions are identical if y , z
u u
19
Recapitulation
x2 2 K xx t y2 2 K yv t z2 2 K zz t
23
Plume rise h
27
Dispersion parameters in Gaussian models
• Derived from concentrations measured in actual atmospheric diffusion
experiments
z w t Fz
y v t Fy
Pasquill
stability
classes
30
Pasquill-Gifford (P-G) curves
Horizontal y and vertical z dispersion coeff
q H2
c( x , 0 , 0 ) exp
u y z 2
2
z
33
Critical downwind distance
as a function of source height and a stability class
(plume that has reached its final height)
34
Summary 3
Gaussian expressions
fail near the surface, since no vertical shear is present
no chemical reactions, either
Eulerian approach
AD eqn provides more general approach (special cases: uniform wind
speed and constant eddy diffusivities), key problem is to choose the
functional forms of the wind speeds and the eddy diffusivities