15.07.2013 Views

Atmospheric Air Pollutant Dispersion

Atmospheric Air Pollutant Dispersion

Atmospheric Air Pollutant Dispersion

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

<strong>Atmospheric</strong> <strong>Air</strong> <strong>Pollutant</strong><br />

• <strong>Atmospheric</strong> Stability<br />

<strong>Dispersion</strong><br />

Estimate air pollutant concentrations<br />

downwind of emission point sources.<br />

Downwind <strong>Air</strong> Pollution Concentrations<br />

Are a function of:<br />

• <strong>Air</strong> Temperature Lapse Rates<br />

• <strong>Atmospheric</strong> <strong>Air</strong> Inversions<br />

• <strong>Atmospheric</strong> Mixing Height<br />

• <strong>Dispersion</strong> from Point Emission Sources<br />

• <strong>Dispersion</strong> Coefficients<br />

In CEE490/ENVH 461 class we will<br />

Primarily use EPA SCREEN software


<strong>Atmospheric</strong> <strong>Air</strong> Vertical Stability<br />

Dry Adiabatic<br />

=<br />

Lapse Rate<br />

ΔTemp<br />

ΔAltitude<br />

ΔT<br />

ΔZ<br />

• Adiabatic vertical air movement causes a<br />

change in pressure and temperature:<br />

– dT/dz = g/C p = γ d (dry adiabatic lapse rate)<br />

γ d = 9.8 K/km<br />

--Stable lapse rate: γ < γ d<br />

– Unstable lapse rate: γ > γ d<br />

=<br />

o<br />

9.76K<br />

=<br />

=<br />

1000meters<br />

o<br />

5.4F<br />

1000ft


<strong>Atmospheric</strong> Stability<br />

Characterized by vertical temperature<br />

gradients (Lapse Rates)<br />

– Dry adiabatic lapse rate (Γ) = 0.976 o C/100 m ~<br />

1 o C/100 m<br />

– International standard lapse rate = 0.0066 o C/m<br />

Does the air temperature lapse rate have anything to<br />

do with air quality?<br />

Yes, because it is related to amount of vertical mixing<br />

of emitted air pollutants.


• First Law of Thermodynamics<br />

= 0 for adiabatic expansion<br />

dq = dh −υdP<br />

= C pdT<br />

• Barometric Equation<br />

dP<br />

= −ρg<br />

dZ<br />

1<br />

⇒ C pdT<br />

= dP = −gdZ<br />

ρ<br />

⇒<br />

dT<br />

dZ<br />

=<br />

−<br />

g<br />

C<br />

p<br />

−<br />

Lapse Rate<br />

1<br />

dP<br />

ρ


100 m<br />

Elevation<br />

(m)<br />

<strong>Air</strong> Temperature Lapse Rates<br />

Superadiabatic<br />

Superadiabatic<br />

Dry Adiabatic Lapse Rate<br />

Subadiabatic<br />

Subadiabatic<br />

T<br />

Temperature ( oC) Inversion Inversion


Stability Conditions<br />

Adiabatic lapse rate<br />

Actual <strong>Air</strong> Temperature lapse rate


z<br />

1 km<br />

Mixing<br />

depth<br />

0<br />

NIGHT<br />

Diurnal Cycle of Surface<br />

MIDDAY<br />

Heating /Cooling<br />

Subsidence<br />

inversion<br />

MORNING<br />

T<br />

NIGHT MORNING AFTERNOON


Actual Temperature Sounding


Superadiabatic Lapse Rates (Unstable air)<br />

• Temperature decreases are greater than -10 o C/1000 meters<br />

• Occur on sunny days<br />

• Characterized by intense vertical mixing<br />

• Excellent dispersion conditions


<strong>Atmospheric</strong> Stability<br />

Superadiabatic – Strong Lapse Rate<br />

Unstable Conditions


Neutral <strong>Air</strong> Temp Lapse Rates<br />

• Temperature decrease with altitude is similar to the adiabatic lapse rate<br />

• Results from:<br />

– Cloudy conditions<br />

– Elevated wind speeds<br />

– Day/night transitions<br />

• Describes OK dispersion conditions<br />

Isothermal Lapse Rates (Weakly Stable)<br />

• Characterized by no temperature change with height<br />

• Atmosphere is somewhat stable<br />

• <strong>Dispersion</strong> conditions are moderate


<strong>Atmospheric</strong> Stability<br />

Subadiabatic – Weak Lapse Rate<br />

Stable Conditions


• 2 major types of inversion:<br />

Subsidence: descent of a layer of air within a high pressure<br />

air mass (descending air increases pressure & temp.)<br />

Radiation: thermal radiation at night from the earth’s<br />

surface into the clear night sky<br />

Radiation<br />

Inversion


Inverted <strong>Air</strong> Temp Lapse Rates (Strongly Stable)<br />

Characterized by increasing air temperature with height<br />

Does it occur during the day or at night? (Both)<br />

Is it associated with high or low air pressure systems? (High)<br />

Does it improve or deteriorate air quality? (deteriorate)<br />

www.ew.govt.nz/enviroinfo/air/weather.htm<br />

Inversion<br />

www.co.mendocino.ca.us/aqmd/Inversions.htm


• Inversion: <strong>Air</strong> Temperature increases with altitude


http://www.co.mendocino.ca.us/aqmd/pages/Inversion-Art-(web).jpg


Radiation Inversions<br />

• Result from radiational cooling of the ground<br />

• Occur on cloudless nights and clear sky – nocturnal<br />

• Are intensified in valleys (heavier cooled air descends to valley floor)<br />

• Cause air pollutants to be “trapped” (poor vertical transport)<br />

www.co.mendocino.ca.us/aqmd<br />

/Inversions.htm<br />

What happens to inversion when sun rises?


Radiation Inversions<br />

• Inversion Breaks up after sunrise<br />

• Breakup results in elevated ground level concentrations<br />

• Breakup described as a fumigation<br />

Red Line is <strong>Air</strong> Temperature<br />

Blue Line is adiabatic air temp lapse


Radiation Inversions<br />

• Elevated inversions are formed over urban areas<br />

– Due to heat island effect


Subsidence Inversion<br />

• Associated with atmospheric high-pressure systems<br />

• Inversion layer is formed aloft due to subsiding air<br />

• Persists for days


Subsidence Inversion<br />

• Migrating high-pressure systems: contribute to the hazy<br />

summer conditions<br />

• Semi-permanent marine high-pressure systems<br />

– Results in a large<br />

number of sunny calm<br />

days<br />

– Inversion layer closest to<br />

the ground on<br />

continental side<br />

– Responsible for air<br />

stagnation over<br />

Southern California<br />

www.oceansatlas.org/.../datard.htm


• Advective - warm air flows over a cold surface


• Mixing Height = Height of air that is mixed and where dispersion occurs<br />

What is the Mixing Height in a radiational inversion?<br />

When does the max MH occur during a day? Min MH?<br />

Which season has the max MH? Min MH?<br />

Why does Phoenix have a larger MH than New Orleans?<br />

Why is agricultural burning allowed only during daytime?


<strong>Air</strong> <strong>Pollutant</strong> <strong>Dispersion</strong> from Point Sources<br />

• Plume rise affects dispersion and transport<br />

– Affects maximum ground level concentrations<br />

– Affects distance to maximum ground level conc.


Lapse Rates and <strong>Atmospheric</strong> Stability<br />

Weak Lapse Condition (Coning)<br />

Wind<br />

Z = altitude<br />

T = Temp<br />

Γ = adiabatic lapse rate


Stack Plume: Coning<br />

Strong wind, no turbulence<br />

What is the stability class? (dashed line is adiabatic<br />

lapse rate) C<br />

Is there good vertical mixing? OK<br />

On sunny or cloudy days? Partly cloudy<br />

Good for dispersing pollutants? OK


Lapse Rates and <strong>Atmospheric</strong> Stability<br />

Inversion Condition (Fanning)<br />

Wind


Stack Plume: Fanning<br />

What is the stability class? (solid line is actual air<br />

temperature with altitude air temp lapse rate)<br />

What is the top view of the plume?<br />

http://www.med.usf.edu/~npoor/4


Stack Plume: Fumigation<br />

Why can’t the pollutants be dispersed upward?<br />

Plume trapped by inversion above stack height.<br />

Does it happen during the day or night? Morning


Lapse Rates and <strong>Atmospheric</strong> Stability<br />

Inversion Below, Lapse Aloft (Lofting)<br />

Wind


Stack Plume: Lofting<br />

Why can’t the pollutants be dispersed downward?<br />

What time of the day or night does this happen?<br />

Evening – night as radiation inversion forms


Lapse Rates and <strong>Atmospheric</strong> Stability<br />

Weak Lapse Below, Inversion Aloft (Trapping)<br />

Wind


Stack Plume: Trapping<br />

What weather conditions cause plume trapping?<br />

Radiation inversion at ground level, subsidence<br />

inversion at higher altitude (evening – night)


Stack Plume: Looping<br />

Strong turbulence<br />

http://www.med.usf.edu/~npoor/3<br />

Is it at stable or unstable condition? Unstable<br />

High or low wind speed? Low wind speed.<br />

Does it happen during the day or night? Day<br />

Is it good for dispersing pollutants? Yes


Dilution of <strong>Pollutant</strong>s in the Atmosphere<br />

• <strong>Air</strong> movement can dilute and remove<br />

pollutants (removal by absorption and<br />

deposition by snow, rain, & to surfaces)<br />

• <strong>Pollutant</strong> dilution is variable, from quite<br />

good to quite poor, according to the wind<br />

velocity and the air stability (lapse rate).


Characteristics of <strong>Dispersion</strong> Models<br />

• The accuracy of air pollutant dispersion models<br />

varies according to the complexity of the<br />

terrain and the availability of historic<br />

meteorological data.<br />

• The acceptability of the results of dispersion<br />

models varies with the experience and<br />

viewpoint of the modeler, the regulator and the<br />

intervener.


<strong>Air</strong> Quality Modeling<br />

Gaussian <strong>Dispersion</strong> Model<br />

EPA <strong>Air</strong> Quality Models (SCREEN TSCREEN,ISC, AERMOD, etc.)<br />

are computer software with equation parameters that include<br />

pollutant emission rate Q (gms/sec), stack height (meters), stack<br />

inside diameter at exit, stack gas temp, stack gas exit velocity<br />

(m/s) ambient air temp, receptor height (m), topography, etc.<br />

and calculate the downwind air pollutant concentrations.<br />

The EPA dispersion software models are used to:<br />

1. Evaluate compliance with NAAQS & prevention of<br />

significant air quality deterioration (PSD required for permit<br />

to construct)<br />

2. Find pollutant emission reductions required.<br />

3. Review permit to construct applications.


Is the <strong>Air</strong> Quality OK?<br />

What is the level of my exposure to these emissions?<br />

Is my family safe?<br />

Where is a safe location with regards to air quality?<br />

How about the adverse impact on the environment<br />

(plants, animals, buildings)?<br />

How to predict the impact of air pollutant emissions<br />

resulting from population growth?<br />

Where is the cleanest air; in city center, in rural area?<br />

Note: Children health (respiratory effects) has been<br />

correlated to the distance from their home to nearest<br />

highway or busy street (diesel engine emissions)


When are model applications required for regulatory<br />

purposes?<br />

• SCREEN3, TSCREEN,<br />

• ISC (Industrial Source Complex),<br />

• AERMOD<br />

AERMOD stands for American merican Meteorological Society<br />

Environmental nvironmental Protection Agency Regulatory egulatory Model Model<br />

Formally Proposed as replacement for ISC in 2000<br />

Adopted as Preferred Model November 9, 2005


Regulatory Application of Models<br />

• PSD: Prevention of Significant Deterioration of <strong>Air</strong> Quality<br />

in relatively clean areas (e.g. National Parks, Wilderness<br />

Areas, Indian Reservations)<br />

• SIP: State Implementation Plan revisions for existing<br />

sources and for New Source Reviews (NSR)


Classifications of <strong>Air</strong> Quality Models<br />

• Developed for a number of air pollutant types and<br />

time periods<br />

– Short-term models – for a few hours to a few<br />

days; worst case episode conditions<br />

– Long-term models – to predict seasonal or<br />

annual average concentrations; health effects<br />

due to exposure<br />

• Classified by<br />

– Non-reactive models – pollutants such as SO 2 and CO<br />

– Reactive models – pollutants such as O 3 , NO 2 , etc.


<strong>Air</strong> Quality Models<br />

• Classified by coordinate system used<br />

– Grid-based<br />

• Region divided into an array of cells<br />

• Used to determine compliance with NAAQS<br />

– Trajectory<br />

• Follow plume as it moves downwind<br />

• Classified by sophistication level<br />

– Screening: simple estimation use preset, worstcase<br />

meteorological conditions to provide<br />

conservative estimates.<br />

– Refined: more detailed treatment of physical and<br />

chemical atmospheric processes; require more<br />

detailed and precise meteorological and<br />

topographical input data.<br />

http://www.epa.gov/scram001/i<br />

mages/smokestacks.jpg


US EPA <strong>Air</strong> Quality Models<br />

• Screening models available at:<br />

www.epa.gov/scram001/tt22.htm#screen<br />

• Preferred models available at:<br />

http://www.epa.gov/scram001/tt22.htm#rec<br />

– A single model found to outperform others<br />

• Selected on the basis of other factors such as past use,<br />

public familiarity, cost or resource requirements and<br />

availability<br />

• No further evaluation of a preferred model is required


Gaussian <strong>Dispersion</strong> Models<br />

• Most widely used<br />

• Based on the assumption<br />

– plume spread results primarily by diffusion<br />

– horizontal & vertical pollutant concentrations in<br />

the plume have double Gaussian distribution)<br />

Q<br />

u<br />

H<br />

z<br />

x<br />

y


Gaussian Model Assumptions<br />

• Gaussian dispersion modeling based on a number of<br />

assumptions including<br />

– Source pollutant emission rate = constant (Steady-state)<br />

– Constant Wind speed, wind direction, and atmospheric<br />

stability class<br />

– <strong>Pollutant</strong> Mass transfer primarily due to bulk air<br />

motion in the x-direction<br />

– No pollutant chemical transformations occur<br />

– Wind speeds are >1 m/sec.<br />

– Limited to predicting concentrations > 50 m downwind


Gaussian <strong>Dispersion</strong> Model


Characteristics of <strong>Pollutant</strong> Plume<br />

Horizontal (y) and vertical (z) dispersion, is<br />

caused by eddies and random shifts of wind<br />

direction.<br />

• Key parameters are:<br />

– Physical stack height (h) – Plume rise (Δh)<br />

– Effective stack height (H) – Wind speed (ux )


Plume <strong>Dispersion</strong> Coordinate System


The Gaussian Model<br />

• C = C(x, y, z, stability)<br />

C<br />

⎛ 2 ⎞<br />

2<br />

Q<br />

⎧ ⎛ ⎞ ⎛<br />

⎜<br />

1 y<br />

⎟<br />

1 ( z − H ) 1 ( z + H )<br />

exp − ⎨exp<br />

⎜<br />

⎜−<br />

⎟ + exp ⎜<br />

⎜−<br />

⎜ 2 ⎟<br />

2<br />

2πσ yσ<br />

z ⎝ 2 σ y ⎠⎩<br />

⎝ 2 σ z ⎠ ⎝ 2 σ z<br />

= 2<br />

• σ y and σ z depend on the atmospheric conditions<br />

• <strong>Atmospheric</strong> stability classifications are defined<br />

in terms of surface wind speed, incoming solar<br />

radiation and cloud cover<br />

2<br />

⎞⎫<br />

⎟<br />

⎟⎬<br />

⎠⎭


C<br />

Gaussian <strong>Dispersion</strong> Equation<br />

( x,<br />

y,<br />

z)<br />

⎡ ⎛ 2<br />

Q 1<br />

⎢ ⎜<br />

y<br />

= exp − +<br />

πσ σ ⎢<br />

⎜ 2<br />

2<br />

2 y zu<br />

⎣<br />

2 ⎝ σ y σ z<br />

( z − H )<br />

σ y & σ z = f(downwind distance x & atmos stability)<br />

• Q = pollutant emission rate (grams/sec)<br />

• H = effective stack height (meters) = stack height + plume rise<br />

• u = wind speed (m/sec)<br />

σ y = horizontal crosswind dispersion coefficient (meters)<br />

σ z= vertical dispersion coefficient (meters)<br />

2<br />

⎞⎤<br />

⎟<br />

⎥<br />

⎠⎥⎦


Plume <strong>Dispersion</strong> Equations<br />

General Equation – Plume with Reflection for Plume Height H<br />

2<br />

( z H ) ⎞ ⎛ ( z + H )<br />

2<br />

Q ⎡ ⎛ y ⎞ ⎧ ⎛<br />

C(<br />

x,<br />

y,<br />

z;<br />

H ) = • ⎢exp⎜<br />

⎟<br />

⎜<br />

− • ⎨ ⎜<br />

⎟<br />

exp −<br />

2<br />

2πuσ σ ⎢⎣<br />

⎝ 2σ<br />

y z<br />

y ⎠ ⎩ ⎝<br />

−<br />

2<br />

2σ<br />

z<br />

⎟ + exp ⎜<br />

⎜−<br />

⎠ ⎝<br />

2<br />

2σ<br />

z<br />

Ground Level Concentration – Plume at Height H<br />

C(<br />

x,<br />

y,<br />

0;<br />

H )<br />

2<br />

2<br />

Q ⎡ ⎛ ⎞ ⎛ ⎞⎤<br />

⎢ ⎜<br />

y H<br />

• exp ⎟<br />

⎜<br />

−<br />

⎟<br />

• exp⎜<br />

− ⎟⎥<br />

πuσ<br />

σ ⎣⎢<br />

⎝ 2σ<br />

⎠ ⎝ 2σ<br />

y z<br />

y<br />

z ⎠⎥⎦<br />

= 2<br />

2<br />

Ground Level Center Line Conc (y = 0) – Plume Height H<br />

C(<br />

x,<br />

0,<br />

0;<br />

H )<br />

2<br />

Q ⎡ ⎛ H ⎞⎤<br />

• ⎢exp⎜<br />

− ⎟⎥<br />

πuσ<br />

σ ⎣ ⎝ 2σ<br />

y z<br />

z ⎠⎦<br />

= 2<br />

Ground Level Center Line – Ground Point Source (y = 0, H = 0)<br />

Q<br />

C(<br />

x,<br />

0,<br />

0;<br />

0)<br />

=<br />

πuσ<br />

σ<br />

y<br />

z<br />

2<br />

⎞⎫⎤<br />

⎟<br />

⎟⎬⎥<br />

⎠⎭⎥⎦


Gaussian <strong>Dispersion</strong> Equation<br />

If the emission source is at ground level with no<br />

effective plume rise then<br />

C<br />

( x,<br />

y,<br />

z)<br />

⎡ ⎛ 2 2<br />

Q 1<br />

⎢ ⎜<br />

y z<br />

= exp − +<br />

πσ σ ⎢<br />

⎜ 2 2<br />

y zu<br />

⎣<br />

2 ⎝ σ y σ z<br />

Plume Rise<br />

• H is the sum of the physical stack height and plume rise.<br />

H<br />

=<br />

Δh<br />

plume rise<br />

+<br />

h<br />

actual<br />

stack<br />

⎞⎤<br />

⎟<br />

⎥<br />

⎠⎥⎦


Key to Stability Categories


<strong>Atmospheric</strong> Stability Classes


Horizontal <strong>Dispersion</strong> Coefficient σ y


Vertical <strong>Dispersion</strong> Coefficient σ z


Maximum Downwind Ground-Level<br />

Concentration (C max )<br />

• <strong>Pollutant</strong>s require time (and distance) to reach the<br />

ground.<br />

• C max decreases as effective plume height H<br />

increases.<br />

• Distance to C max increases as H increases.


Maximum Concentration (C max ) and Distance to<br />

C max (x max )


Maximum Ground Level Concentration<br />

Under moderately stable to near neutral conditions,<br />

σ y = 1<br />

k σ<br />

z<br />

The ground level concentration at the center line is<br />

2<br />

Q ⎡ H ⎤<br />

C x,<br />

0,<br />

0 = exp⎢−<br />

2<br />

2<br />

πk1σ<br />

z u ⎣ 2σ<br />

z<br />

( ) ⎥ ⎦<br />

The maximum occurs at<br />

dC / dσz<br />

= 0 ⇒ σz<br />

=<br />

H<br />

2<br />

Once σ z is determined, x can be known and subsequently C.<br />

C<br />

Q<br />

= = 0.<br />

1171<br />

πσ σ u<br />

σ<br />

( x,<br />

0,<br />

0)<br />

exp[<br />

−1]<br />

y<br />

z<br />

Q<br />

σ<br />

y<br />

z<br />

u


Plume Rise<br />

• H is the sum of the physical stack height and plume rise.<br />

H<br />

=<br />

Δh<br />

plumerise<br />

+<br />

h<br />

actualstack


<strong>Pollutant</strong> Plume Rise<br />

Plume rises from stack exit.


Plume Rise<br />

Buoyant plume: Initial buoyancy >> initial momentum<br />

Forced plume: Initial buoyancy ~ initial momentum<br />

Jet: Initial buoyancy 55 m / s )<br />

where buoyancy flux is<br />

2<br />

F = gV d ( T −T<br />

) / 4T<br />

s<br />

u<br />

a<br />

S<br />

)<br />

V s: Stack exit velocity, m/s<br />

d: top inside stack diameter, m<br />

T s: stack gas temperature, K<br />

T a: ambient temperature, K<br />

g: gravity, 9.8 m/s 2


Carson and Moses: vertical momentum & thermal<br />

buoyancy, based on 615 observations involving 26 stacks.<br />

Δh<br />

Δh<br />

Δh<br />

plume rise<br />

plume rise<br />

plume rise<br />

=<br />

=<br />

3.<br />

47<br />

0.<br />

35<br />

Vsd<br />

u<br />

Vsd<br />

u<br />

Vsd<br />

= −1.<br />

04<br />

u<br />

( T T )<br />

Q = m&<br />

C −<br />

h<br />

m&<br />

=<br />

πd<br />

4<br />

2<br />

p<br />

s<br />

P<br />

Vs<br />

RT<br />

s<br />

a<br />

+<br />

+<br />

5.<br />

15<br />

2.<br />

64<br />

+<br />

2.<br />

24<br />

Q<br />

u<br />

Q<br />

u<br />

h<br />

h<br />

Q<br />

u<br />

h<br />

(unstable)<br />

(neutral)<br />

(stable)<br />

(heat emission rate, kJ/s)<br />

(stack gas mass flow rate. kg/s)


C<br />

( x,<br />

y,<br />

z)<br />

Q<br />

2πσ<br />

σ<br />

⎡ y<br />

exp⎢−<br />

⎢⎣<br />

2σ<br />

2<br />

Wark & Warner, “<strong>Air</strong> Pollution: Its Origin & Control”<br />

⎤⎪⎧<br />

⎡<br />

⎥⎨exp⎢−<br />

⎥⎦<br />

⎪⎩ ⎣<br />

2 ( z − H ) ⎤ ⎡ ( z + H )<br />

= 2<br />

2<br />

2<br />

y zu<br />

y 2σz<br />

2σz<br />

⎥ + exp⎢−<br />

⎦ ⎣<br />

Plume is Reflected when it touches ground surface.<br />

2<br />

⎤⎪⎫<br />

⎥⎬<br />

⎦⎪⎭


Ground level concentration<br />

C<br />

Q<br />

πσ σ<br />

⎡<br />

exp⎢−<br />

⎢⎣<br />

y<br />

2σ<br />

⎤ ⎡<br />

⎥ exp⎢−<br />

⎥⎦<br />

⎣<br />

H<br />

2σ<br />

= 2<br />

2<br />

y zu<br />

y<br />

z<br />

2<br />

2<br />

⎤<br />

⎥<br />


<strong>Dispersion</strong> of SO 2 from Stack<br />

• An industrial boiler is burning at 12 tons of 2.5% sulfur coal/hr<br />

with an emission rate of 151 g/s. The following exist : H = 120<br />

m, u = 2 m/s, y = 0. It is one hour before sunrise, and the sky<br />

is clear. Find the downwind ground level SO2 concentration at<br />

x =2 km, y = 0, and z = 0.<br />

Stability class =<br />

σy =<br />

σz =<br />

C(x=2 km, y=0, z=0) =<br />

C<br />

Q<br />

πσ σ<br />

⎡ 2<br />

y<br />

exp⎢−<br />

⎢⎣<br />

2σ<br />

⎤ ⎡ H<br />

⎥exp⎢−<br />

⎥⎦<br />

⎣ 2σ<br />

= 2<br />

2<br />

y zu<br />

y<br />

z<br />

2<br />

⎤<br />

⎥<br />


<strong>Dispersion</strong> of Ground Level <strong>Air</strong> <strong>Pollutant</strong> Emissions<br />

• <strong>Air</strong> pollutant emissions are from a ground level source<br />

with H = 0, u = 4 m/s, Q = 100 g/s, and the stability class<br />

= B, what is downwind concentration at x = 200 m, y = 0,<br />

and z = 0?<br />

At 200 m:<br />

σ y =<br />

σ z =<br />

C(x=200 m, y=0, z=0) =


• Calculate H using plume rise equations for an 80 m high stack (h) with a stack<br />

diameter = 4 m, stack gas velocity = 14 m/s, stack gas temperature = 90o C (363<br />

K), ambient temperature = 25 oC (298 K), 10 meter high wind speed u at 10 m<br />

= 4m/s, and stability class = B. Find the Max Ground Level air pollutant<br />

concentration & its (location downwind distance x).<br />

F =<br />

Δh plume rise =<br />

H =<br />

σz =<br />

σy =<br />

Cmax =<br />

Plume Rise and Max Conc<br />

Residents around the Rock Cement Plant are complaining that its emission<br />

are in violation and the ground level concentrations exceed the air quality<br />

standards. The plant has its facility within 0.2 km diameter fence. Its<br />

effective stack height is 50 m. You are a government agency environmental<br />

engineer. Where are you going to locate your air quality monitors? Why?


Complex Horizontal, Vertical, and Temporal Wind Structure<br />

Winds aloft have no continuous measurements


Meteorology


In most of cases terrain is not flat terrain – topographic complexity<br />

Complex horizontal, vertical, and temporal dispersion


Plume <strong>Air</strong> <strong>Pollutant</strong> Building Downwash


Building Downwash


Building Downwash


Building Downwash – Short Stack


Building Downwash – Taller Stack (Not GEP)


Building Downwash – Tallest Stack<br />

H GEP = 2.5 H Bldg Prior to 1979 H GEP = H Bldg + 1.5 L


Cavity and Wakes


Cavity and Far Wake<br />

Cavity


Building Downwash for 2 Identical Stack Emissions at Different Locations<br />

Building<br />

Cavity<br />

The stack on the left is located on top of a building and this structure affects the<br />

wind-flow which, in turn, affects the plume trajectory, pulling it down into the<br />

cavity zone (near wake) behind the building. The stack on the right is located far<br />

enough downwind of the building to be unaffected by the cavity (near wake)<br />

effects and is only affected by the air flow in the far wake.


Building Downwash<br />

Cavity


Z R<br />

Cavity<br />

X R<br />

2Y R<br />

Cavity and Wakes<br />

Z R<br />

2Y R


Aerodynamic Wake<br />

• Region where local air velocities are different from<br />

the free stream values<br />

• Streamline Separation – at an object<br />

– Eddy Recirculation (generally lower velocity region)<br />

– Turbulent Shear Region (generally higher velocity)<br />

• Reattachment of Streamlines<br />

• Near Wake (Cavity)<br />

– Usually on the ‘lee’ side of the object.<br />

• Far Wake<br />

– Effect of another object on the separated streamlines<br />

• Estimation of Wake/Cavity Boundary<br />

– Effect of Building geometry


Two Points of Concern<br />

• <strong>Dispersion</strong> of <strong>Air</strong> <strong>Pollutant</strong> Plume in the presence<br />

of Buildings<br />

– What happens to the plume from an existing stack and<br />

existing buildings?<br />

• Design of Stacks in the Presence of Building<br />

– What are the design guidelines if a new emission source<br />

is being proposed in an area full of buildings?


Good Engineering Practice GEP Analysis<br />

• H GEP = good engineering practice height<br />

• H GEP = H + 1.5 L<br />

– H = height of adjacent or nearby structure<br />

– L = lesser dimension height or proj. width<br />

– 5L = region of influence<br />

• The structure with the greatest influence<br />

is then used in the model to evaluate<br />

wake effects and downwash.


Good Engineering Practice Stack Height<br />

1985 Regulations -<br />

GEP Stack height – greater of the following<br />

- a) 65 m from the base of the stack<br />

-b) H GEP = 2.5 H<br />

for stacks in existence before Jan, 1979.<br />

For All Other Stacks:<br />

H GEP = H + 1.5 L<br />

Where H is the height of the nearest building and<br />

L = lesser dimension height or projected building width


GEP stack height<br />

• Building downwash can occur when<br />

H Stack = H S < H b + 1.5L<br />

H Stack = Height of Stack<br />

H b = Height of Building<br />

L = lesser of H b or PBW<br />

PBW = Maximum Projected Building Width<br />

Screen models will do this calculation when the<br />

building downwash option is used. If HS > Hb +<br />

1.5L, then building downwash will not be shown in<br />

SCREEN results


Getting Started – Screen 3<br />

• Convert all lengths and distances to meters<br />

• Convert temperatures to degrees Kelvin<br />

• Identify building contributions to air<br />

dispersion (stack emissions)<br />

• Screen 3 should be run in regulatory<br />

default mode<br />

• Note that SCREEN3 differs from TSCREEN (SCREEN3<br />

can select atmospheric stability classes or run Full<br />

Meteorology, will not calculate concentrations at other<br />

averaging times, etc.)


Information Required to run<br />

Screen Models for Point Source<br />

• Emission rate (g/s)<br />

• Stack Height (m)<br />

• Shortest distance to property line<br />

• Stack velocity (or volumetric airflow SCREEN3)<br />

• Stack gas temperature ( o K)<br />

• Stack Inside Diameter<br />

• Building Height, Length, Width


SCREEN Example – Stack emissions<br />

• Emission rate (g/s) .01 g/s<br />

• Stack Height 15.24m<br />

• Building Height 6.096m<br />

• Shortest Distance to property line 91.44m<br />

• Stack airflow in acfm 20,000 acfm<br />

• Stack gas temperature 294.3° K<br />

• Stack inside diameter 1.143m<br />

• Building dimensions 30.48m L, 30.48m W, 10.67m H


SCREEN Example – Stack emissions<br />

Note that this is the output from SCREEN3 software (not TSCREEN)


SCREEN Example – Stack emissions


SCREEN Example – Stack emissions


SCREEN Example – Stack emissions


SCREEN Example – Stack emissions


SCREEN Example – Stack emissions


SCREEN Example – Stack emissions<br />

The most conservative scenario gives a maximum 1-hr 1 hr<br />

concentration of 2.595 ug/m 3 at a distance of 91 meters


Flare Emissions Elevated point source


Screen Model for Flare Source<br />

• Emission Rate<br />

• Flare Stack Height<br />

• Total Heat Release Rate<br />

• Shortest Distance to property line<br />

• Influential Building Dimensions


SCREEN Example – Flare emissions


SCREEN Example – Flare emissions


SCREEN Example – Flare emissions


SCREEN Example – Flare emissions


Screen Model for Area Source<br />

• Emission Rate<br />

• Source Release Height<br />

• Larger Side Length of Rectangular Area<br />

• Smaller Side Length of Rectangular Area<br />

• Shortest Distance to property line


SCREEN Example – Area Source


SCREEN Example – Flare emissions


Screen Model for Volume Source<br />

• Emission Rate<br />

• Source Release Height<br />

• Initial Lateral Dimension<br />

• Initial Vertical Dimension<br />

• Shortest Distance to Property Line


Volume Source<br />

• Source Release Height is the center of the<br />

Volume Source:<br />

If the Source is from a building, the release<br />

height is set equal to one half of the building<br />

height.<br />

• Volume sources are modeled as a square in<br />

Screen3. If the source is not square, the width<br />

should be set to the minimum length.


Volume Source<br />

Initial Lateral Dimension (σ y0 )<br />

Single Volume Source<br />

σ y0 = length of side divided by 4.3<br />

Line Source composed of several volume sources<br />

σ y0 = length of side divided by 2.15<br />

Line Source composed of separated volume sources<br />

σ y0 = center to center distance divided by 2.15


Volume Source<br />

Initial Vertical Dimension (σ z0 )<br />

Surface-Based Source<br />

σ z0 = vertical dimension of source divided by 2.15<br />

Elevated Source on or adjacent to a building<br />

σ z0 = building height divided by 2.15<br />

Elevated Source not on or adjacent to a building<br />

σ z0 = vertical dimension of source divided by 4.3


Vertical Dimension is<br />

height of building<br />

divided by 2.15<br />

Example – Volume Source<br />

Volume Source from a Building<br />

Lateral Dimension is<br />

minimum length of<br />

building divided by 4.3<br />

Release Height is ½ of<br />

Building Height


SCREEN Example – Volume Source


SCREEN Example – Volume Source<br />

Initial Lateral Dimension obtained by taking<br />

building length of 30.48m divided by 4.3<br />

Using a building as the volume source, so the<br />

initial vertical dimension is the height of the<br />

building divided by 2.15


SCREEN Example – Volume Source


SCREEN Example – Volume Source


Screen model results<br />

• Screen model results are all<br />

maximum 1-hr concentrations (except<br />

for complex terrain and if SCREEN is<br />

run inside of TSCREEN can obtain<br />

concentrations in other averaging<br />

times).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!