Atmospheric Dispersion Modeling Article Index for
Atmospheric
Website Links For
Atmospheric
 

Information About

Atmospheric Dispersion Modeling




The dispersion models require the input of data which includes:

  • Meteorological conditions such as wind speed and direction, the amount of atmospheric Turbulence (as characterized by what is called the "stability class"), the ambient air temperature and the height to the bottom of any Inversion aloft that may be present.

  • Emissions parameters such as source location and height, source vent stack diameter and exit Velocity , exit temperature and Mass Flow Rate .

  • Terrain elevations at the source location and at the receptor location.

  • The location, height and width of any obstructions (such as buildings or other structures) in the path of the emitted gaseous plume.


Many of the modern, advanced dispersion modeling programs include a pre-processor module for the input of meteorological and other data, and many also include a post-processor module for graphing the output data and/or plotting the area impacted by the air pollutants on maps.

The atmospheric dispersion models are also known as atmospheric diffusion models, air dispersion models, air quality models, and air pollution dispersion models.




GAUSSIAN AIR POLLUTANT DISPERSION EQUATION


The technical literature on air pollution dispersion is quite extensive and dates back to the 1930's and earlier. One of the early air pollutant plume dispersion equations was derived by Bosanquet and Pearson in 1936. Their equation did not assume Gaussian Distribution nor did it include the effect of ground reflection of the pollutant plume.

Sir Graham Sutton derived an air pollutant plume dispersion equation in 1947 which did include the assumption of Gaussian distribution for the vertical and crosswind dispersion of the plume and also included the effect of ground reflection of the plume.

Under the stimulus provided by the advent of stringent Environmental Control Regulations , there was an immense growth in the use of air pollutant plume dispersion calculations between the late 1960s and today. A great many computer programs for calculating the dispersion of air pollutant emissions were developed during that period of time and they were called "air dispersion models". The basis for most of those models was the Complete Equation For Gaussian Dispersion Modeling Of Continuous, Buoyant Air Pollution Plumes shown below:



C = rac{\;Q}{u}\cdot rac{\;f}{\sigma_y\sqrt{2\pi}}\;\cdot rac{\;g_1 + g_2 + g_3}{\sigma_z\sqrt{2\pi}}



where:

f =\; crosswind\;\, dispersion\;\, parameter

f =\;\exp\; {Link without Title}

g =\;vertical\;\,dispersion\;\,parameter =\,g_1 + g_2 + g_3

g_1 =\;vertical\;\,dispersion\;\, with\;\, no \;\, reflections

g_1 =\; \exp\; - H)^2/\,(2\;\sigma_z^2\;)\;

g_2 =\;vertical\;\,dispersion\;\,for\;\,reflection\;\, from\;\,the\;\, ground

g_2=\;\exp\; + H)^2/\,(2\;\sigma_z^2\;)\;

g_3 =\;vertical\;\,dispersion\;\,for\;\,reflection\;\, from\;\,inversion\;\,aloft

g_3 =\;\sum_{m=1}^\infty\;\big\{\exp\; - H - 2mL)^2/\,(2\;\sigma_z^2\;)\;

:::+\, \exp\; + H + 2mL)^2/\,(2\;\sigma_z^2\;)\;
:::+\, \exp\; + H - 2mL)^2/\,(2\;\sigma_z^2\;)\;
:::+\, \exp\; - H + 2mL)^2/\,(2\;\sigma_z^2\;)\; \big\}



C =\;concentration\;\,of\;\,emissions,\;\, \mbox{g}/\mbox{m}^3, \;\,at\;\,any\;\,receptor\;\,located:

:::x\;\,meters\;\,downwind\;\,from\;\,the\;\,emission\;\,source\;\,point

:::y\;\,meters\;\,crosswind\;\,from\;\,the\;\,emission\;\,plume\;\,centerline

:::z\;\,meters\;\,above\;\,ground\;\,level

Q\; =\;source\;\,pollutant\;\,emission\;\,rate, \;\,\mbox{g}/\mbox{s}

 u\; =\;horizontal\;\,wind\;\,velocity\;\,along \;\,the\;\,plume\;\,centerline,\;\, \mbox{m}/\mbox{s}

H\, =\;\,height\;\,of\;emission\;\, plume\;\,centerline\;\, above\;\,ground\;\,level,\,\; \mbox{m}

\sigma_z =\;\,vertical\;\,standard\;\, deviation\;\,of\;\,emission\;\,distribution,\,\; \mbox{m}

\sigma_y =\;\,horizontal\;\,standard\;\, deviation\;\,of\;\,emission\;\,distribution,\,\; \mbox{m}

L\; =\;\,height\;\,from\;\,ground\;\,level\;\, to\;\,bottom\;\,of\;\,inversion\;\,aloft,\,\; \mbox{m}



The above equation not only includes upward reflection from the ground, it also includes downward reflection from the bottom of any inversion lid present in the atmosphere.



The sum of the four exponential terms in ''g''3 converges to a final value quite rapidly. For most cases, the summation of the series with ''m'' = 1, ''m'' = 2 and ''m'' = 3 will provide an adequate solution.



It should be noted that σz and σy are functions of the atmospheric stability class (i.e., a measure of the turbulence in the ambient atmosphere) and of the downwind distance to the receptor. The two most important variables affecting the degree of pollutant emission dispersion obtained are the height of the emission source point and the degree of atmospheric turbulence. The more turbulence, the better the degree of dispersion.

The resulting calculations for air pollutant concentrations are often expressed as an Air Pollutant concentration Contour Map in order to show the spatial variation in contaminant levels over a wide area under study. In this way the contour lines can overlay sensitive receptor locations and reveal the spatial relationship of air pollutants to areas of interest.


THE BRIGGS' PLUME RISE EQUATIONS


The Gaussian air pollutant dispersion equation (discussed above) requires the input of ''H'' which is the pollutant plume's centerline height above ground level—and H
is the sum of ''H''s (the actual physical height of the pollutant plume's emission source point) plus Δ''H'' (the plume rise due the plume's buoyancy).

To determine Δ''H'', many if not most of the air dispersion models developed and used between the late 1960s and the early 2000s used what are known as "the Brigg's equations". Briggs published his first plume rise model observations and comparisons in 1965. In 1968, at a symposium sponsored by CONCAWE (a Dutch organization), he compared many of the plume rise models then available in the literature. In that same year, Briggs also wrote the section of the publication edited by Slade dealing with the comparative analyses of plume rise models. That was followed in 1969 by his classical critical review of the entire plume rise literature, in which he proposed a set of plume rise equations which have became widely known as "the Briggs equations". Subsequently, Briggs modified his 1969 plume rise equations in 1971 and in 1972.

Briggs divided air pollution plumes into these four general categories:

  • Cold jet plumes in calm ambient air conditions

  • Cold jet plumes in windy ambient air conditions

  • Hot, buoyant plumes in calm ambient air conditions

  • Hot, buoyant plumes in windy ambient air conditions


Briggs considered the trajectory of cold jet plumes to be dominated by their initial velocity momentum, and the trajectory of hot, buoyant plumes to be dominated by their buoyant momentum to the extent that their initial velocity momentum was relatively unimportant. Although Briggs proposed plume rise equations for each of the above plume categories, it is important to emphasize that "the Briggs equations" which become widely used are those that he proposed for bent-over, hot buoyant plumes.

In general, Briggs' equations for bent-over, hot buoyant plumes are based on observations and data involving plumes from typical combustion sources such as the Flue Gas Stacks from steam-generating boilers burning Fossil Fuel s in large power plants. Therefore the stack exit velocities were probably in the range of 20 to 100 ft/s (6 to 30 m/s) with exit temperatures ranging from 250 to 500 °F (120 to 260 °C).

A logic diagram for using the Briggs equations to obtain the plume rise trajectory of bent-over buoyant plumes is presented below:


FURTHER READING


For those who would like to learn more about this topic, it is suggested that either one of the following books be read:




REFERENCES