Air pollution dispersion modeling: Difference between revisions
imported>Milton Beychok m (Added a CZ link) |
Pat Palmer (talk | contribs) mNo edit summary |
||
| (28 intermediate revisions by 7 users not shown) | |||
| Line 1: | Line 1: | ||
{{subpages}} | {{subpages}} | ||
{{Image|AirPollutionSource.jpg|right|160px|Industrial air pollution sources}} | |||
'''Air pollution dispersion modeling''' is the [[Mathematics|mathematical]] [[Computer simulation|simulation]] of how air pollutants disperse in the ambient atmosphere. It is performed with computer programs, called [[Air pollution dispersion terminology#Dispersion models|dispersion models]], that solve the mathematical equations and [[algorithm]]s which simulate the pollutant dispersion. The dispersion models are used to estimate or to predict the downwind concentration of air pollutants emitted from [[Air pollution dispersion terminology#Air pollutant emission sources|emission sources]] such as industrial plants and vehicular traffic. Such models are important to governmental agencies tasked with protecting and managing the ambient air quality. The models are typically employed to determine whether existing or proposed new industrial facilities are or will be in compliance with the National Ambient Air Quality Standards (NAAQS) in the [[United States]] and similar standards in other nations. The models also serve to assist in the design of effective control strategies to reduce emissions of harmful air pollutants. | '''Air pollution dispersion modeling''' is the [[Mathematics|mathematical]] [[Computer simulation|simulation]] of how air pollutants disperse in the ambient atmosphere. It is performed with computer programs, called [[Air pollution dispersion terminology#Dispersion models|dispersion models]], that solve the mathematical equations and [[algorithm]]s which simulate the pollutant dispersion. The dispersion models are used to estimate or to predict the downwind [[concentration]] of air pollutants emitted from [[Air pollution dispersion terminology#Air pollutant emission sources|emission sources]] such as industrial plants and vehicular traffic. Such models are important to governmental agencies tasked with protecting and managing the ambient [[air quality]]. The models are typically employed to determine whether existing or proposed new industrial facilities are or will be in compliance with the [[National Ambient Air Quality Standards]] (NAAQS) in the [[United States of America]] and similar standards in other nations. The models also serve to assist in the design of effective control strategies to reduce emissions of harmful air pollutants. | ||
The dispersion models require the input of data which includes: | The dispersion models require the input of data which includes: | ||
* Meteorological conditions such as [[wind]] speed and direction, the amount of [[Air pollution dispersion terminology|atmospheric turbulence]] (as characterized by what is referred to the [[Air pollution dispersion terminology|''stability class'']], the ambient air temperature and the height to the bottom of any temperature [[Air pollution dispersion terminology|inversion]] that may be present aloft. | * Meteorological conditions such as [[wind]] speed and direction, the amount of [[Air pollution dispersion terminology|atmospheric turbulence]] (as characterized by what is referred to as the [[Air pollution dispersion terminology|''stability class'']], the ambient air temperature and the height to the bottom of any temperature [[Air pollution dispersion terminology|inversion]] that may be present aloft. | ||
* Emissions parameters such as source location and height, source vent stack diameter and exit velocity, exit temperature and mass flow rate. | * 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. | * 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 [[Gas|gaseous]] [[Air pollution dispersion terminology|emission plume]]. | * The location, height and width of any obstructions (such as buildings or other structures) in the path of the [[Gas|gaseous]] [[Air pollution dispersion terminology|emission 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. | 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. Currently, the [[AERMOD air pollution dispersion model]] is the preferred regulatory model of the [[U.S. Environmental Protection Agency]]. | ||
The air pollution dispersion models are also known as atmospheric dispersion models, atmospheric diffusion models, air dispersion models and air quality models. | The air pollution dispersion models are also known as atmospheric dispersion models, atmospheric diffusion models, air dispersion models and air quality models. | ||
| Line 16: | Line 16: | ||
==Gaussian air pollutant dispersion equation== | ==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.<ref>Bosanquet, C.H. and Pearson, J.L., "The spread of smoke and gases from chimneys", Trans. Faraday Soc., 32:1249, 1936</ref> Their equation did not assume [[Gaussian distribution]] nor did it include the effect of ground reflection of the pollutant plume. | 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.<ref>Bosanquet, C.H. and Pearson, J.L., "The spread of smoke and gases from chimneys", Trans. Faraday Soc., 32:1249, 1936</ref> Their equation did not assume [[Normal distribution|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<ref>Sutton, O.G., "The problem of diffusion in the lower atmosphere", QJRMS, 73:257, 1947 and "The theoretical distribution of airborne pollution from factory chimneys", QJRMS, 73:426, 1947</ref> 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. | Sir Graham Sutton derived an air pollutant plume dispersion equation in 1947<ref>Sutton, O.G., "The problem of diffusion in the lower atmosphere", QJRMS, 73:257, 1947 and "The theoretical distribution of airborne pollution from factory chimneys", QJRMS, 73:426, 1947</ref> 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. Since personal computers also came into existence during that period, a great many computer programs for calculating the dispersion of air pollutant emissions were developed in that same period. The basis for most of those models was the '''Complete Equation For Gaussian Dispersion Modeling Of Continuous, [[Buoyant]] Air Pollution Plumes | 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. Since personal computers also came into existence during that period, a great many computer programs for calculating the dispersion of air pollutant emissions were developed in that same period. The basis for most of those models was the '''Complete Equation For Gaussian Dispersion Modeling Of Continuous, [[Buoyant]] Air Pollution Plumes''' shown below: <ref>{{cite book|author=Turner, D.B.|title=[[Workbook of Atmospheric Dispersion Estimates]]| edition=2nd Edition|publisher=CRC Press|year=1994|id=ISBN 1-56670-023-X}} [http://www.crcpress.com/shopping_cart/products/product_detail.asp?sku=L1023&parent_id=&pc= www.crcpress.com]</ref><ref name=Beychok>{{cite book|author=Beychok, Milton R.|title=[[Fundamentals of Stack Gas Dispersion]]|edition=4th Edition|publisher=author-published|year=2005|id=ISBN 0-9644588-0-2}}</ref><br><br> | ||
<math>C = \frac{\;Q}{u}\cdot\frac{\;f}{\sigma_y\sqrt{2\pi}}\;\cdot\frac{\;g_1 + g_2 + g_3}{\sigma_z\sqrt{2\pi}}</math><br><br> | <math>C = \frac{\;Q}{u}\cdot\frac{\;f}{\sigma_y\sqrt{2\pi}}\;\cdot\frac{\;g_1 + g_2 + g_3}{\sigma_z\sqrt{2\pi}}</math><br><br> | ||
| Line 63: | Line 63: | ||
|- | |- | ||
!align=right| | !align=right| | ||
|align=left| <math>+\, \exp\;[-\,(z - H + 2mL)^2/\,(2\;\sigma_z^2\;)\;]</math> | |align=left| <math>+\, \exp\;[-\,(z - H + 2mL)^2/\,(2\;\sigma_z^2\;)\;]\big\}</math> | ||
|- | |- | ||
!align=right|<math>C</math> | !align=right|<math>C</math> | ||
| Line 99: | Line 99: | ||
|} | |} | ||
The above equation not only includes upward reflection of the pollution plume from the ground, it also includes downward reflection from the bottom of any temperature inversion lid present in the atmosphere. | The above equation not only includes upward reflection of the pollution plume from the ground, it also includes downward reflection from the bottom of any temperature inversion lid present in the atmosphere. | ||
The sum of the four exponential terms in <math>g_3</math> 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. | The sum of the four exponential terms in <math>g_3</math> 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 <math>\sigma_z</math> and <math>\sigma_y</math> 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. | It should be noted that <math>\sigma_z</math> and <math>\sigma_y</math> 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. | ||
| Line 109: | Line 109: | ||
==The Briggs plume rise equations== | ==The Briggs plume rise equations== | ||
The Gaussian air pollutant dispersion equation (discussed above) requires the input of ''H'' (also known as the effective plume height, ''H<sub>e</sub>'') which is the pollutant plume's centerline height above ground level. ''H<sub>e</sub>'' at any distance from the pollutant plume's source is the sum of ''H''<sub>s</sub> (the actual physical height of the pollutant plume's source point) plus Δ''H'' (the plume rise due the plume's buoyancy) at that distance | The Gaussian air pollutant dispersion equation (discussed above) requires the input of ''H'' (also known as the effective plume height, ''H<sub>e</sub>'') which is the pollutant plume's centerline height above ground level. ''H<sub>e</sub>'' at any distance from the pollutant plume's source is the sum of ''H''<sub>s</sub> (the actual physical height of the pollutant plume's source point) plus Δ''H'' (the plume rise due the plume's buoyancy) at that distance. | ||
[[Image:Gaussian Plume.png|thumb|right|333px|{{#ifexist:Template:AirPollutionSource.jpg/credit|{{Gaussian Plume.png/credit}}<br/>|}}Visualization of a buoyant Gaussian air pollutant dispersion plume]] | [[Image:Gaussian Plume.png|thumb|right|333px|{{#ifexist:Template:AirPollutionSource.jpg/credit|{{Gaussian Plume.png/credit}}<br/>|}}Visualization of a buoyant Gaussian air pollutant dispersion plume]] | ||
To determine Δ''H'', many if not most of the air dispersion models developed between the late 1960s and the early 2000s used what are known as "the Briggs equations." G.A. Briggs first published his plume rise observations and comparisons in 1965.<ref>Briggs, G.A., "A plume rise model compared with observations", JAPCA, 15:433-438, 1965</ref> In 1968, at a symposium sponsored by CONCAWE (a Dutch organization), he compared many of the plume rise models then available in the literature.<ref>Briggs, G.A., "CONCAWE meeting: discussion of the comparative consequences of different plume rise formulas", Atmos. Envir., 2:228-232, 1968</ref> In that same year, Briggs also wrote the section of the publication edited by Slade<ref>Slade, D.H. (editor), "Meteorology and atomic energy 1968", Air Resources Laboratory, U.S. Dept. of Commerce, 1968</ref> 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,<ref>Briggs, G.A., "Plume Rise", USAEC Critical Review Series, 1969</ref> 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.<ref>Briggs, G.A., "Some recent analyses of plume rise observation", Proc. Second Internat'l. Clean Air Congress, Academic Press, New York, 1971</ref><ref>Briggs, G.A., "Discussion: chimney plumes in neutral and stable surroundings", Atmos. Envir., 6:507-510, 1972</ref> | To determine Δ''H'', many if not most of the air dispersion models developed between the late 1960s and the early 2000s used what are known as "the Briggs equations." G.A. Briggs first published his plume rise observations and comparisons in 1965.<ref>Briggs, G.A., "A plume rise model compared with observations", JAPCA, 15:433-438, 1965</ref> In 1968, at a symposium sponsored by CONCAWE (a Dutch organization), he compared many of the plume rise models then available in the literature.<ref>Briggs, G.A., "CONCAWE meeting: discussion of the comparative consequences of different plume rise formulas", Atmos. Envir., 2:228-232, 1968</ref> In that same year, Briggs also wrote the section of the publication edited by Slade<ref>Slade, D.H. (editor), "[[Meteorology]] and atomic energy 1968", Air Resources Laboratory, U.S. Dept. of Commerce, 1968</ref> 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,<ref>Briggs, G.A., "Plume Rise", USAEC Critical Review Series, 1969</ref> 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.<ref>Briggs, G.A., "Some recent analyses of plume rise observation", Proc. Second Internat'l. Clean Air Congress, Academic Press, New York, 1971</ref><ref>Briggs, G.A., "Discussion: chimney plumes in neutral and stable surroundings", Atmos. Envir., 6:507-510, 1972</ref> | ||
Briggs divided air pollution plumes into these four general categories: | Briggs divided air pollution plumes into these four general categories: | ||
| Line 122: | Line 122: | ||
* Hot, buoyant plumes in windy 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'''''. | 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 (as depicted in the adjacent diagram of a plume)'''''. | ||
In general, Briggs's equations for bent-over, hot buoyant plumes are based on observations and data involving plumes from typical combustion sources such as the [[flue gas stack]]s from [[steam-generating boiler]]s burning [[fossil fuel]]s in large [[power plant]]s. 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). | In general, Briggs's equations for bent-over, hot buoyant plumes are based on observations and data involving plumes from typical combustion sources such as the [[flue gas stack]]s from [[steam-generating boiler]]s burning [[fossil fuel]]s in large [[power plant]]s. 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). | ||
| Line 138: | Line 138: | ||
|- | |- | ||
!align=right| F<sup> </sup> <!-- The HTML is needed to line up characters. Do not remove.--> | !align=right| F<sup> </sup> <!-- The HTML is needed to line up characters. Do not remove.--> | ||
|align=left|= buoyancy factor, in m<sup>4</sup>/s<sup>3</sup> | |align=left|= [[buoyancy]] factor, in m<sup>4</sup>/s<sup>3</sup> | ||
|- | |- | ||
!align=right| x | !align=right| x | ||
| Line 169: | Line 169: | ||
* RIMPUFF: Developed in [[Denmark]] | * RIMPUFF: Developed in [[Denmark]] | ||
The European Topic Centre on Air and Climate Change, which is part of the European Environment Agency (EEA), maintains an online Model Documentation System (MDS)that includes descriptions and other information for almost all of the dispersion models developed by the countries of Europe. The MDS currently contains | The European Topic Centre on Air and Climate Change, which is part of the [[European Environment Agency]] (EEA), maintains an online Model Documentation System (MDS) that includes descriptions and other information for almost all of the dispersion models developed by the countries of Europe. The MDS currently contains about 140 models developed in [[Europe]] (excluding the United Kingdom).<ref>[http://pandora.meng.auth.gr/mds/strquery.php?wholedb The Model Documentation System (MDS)]</ref> | ||
==References== | ==References== | ||
| Line 178: | Line 178: | ||
For those who would like to learn more about this topic, it is suggested that either one of the following books be read: | For those who would like to learn more about this topic, it is suggested that either one of the following books be read: | ||
*{{cite book | author=Turner, D.B. | title=Workbook of | *{{cite book | author=Turner, D.B. | title=[[Workbook of Atmospheric Dispersion Estimates]]| edition=2nd Edition| publisher=CRC Press | year=1994 | id=ISBN 1-56670-023-X }} [http://www.crcpress.com/shopping_cart/products/product_detail.asp?sku=L1023&parent_id=&pc= www.crcpress.com] | ||
*{{cite book|author=Beychok, Milton R.|title=Fundamentals of Stack Gas Dispersion|edition=4th Edition|publisher=author-published|year=2005|id=ISBN 0-9644588-0-2}} | *{{cite book|author=Beychok, Milton R.|title=[[Fundamentals of Stack Gas Dispersion]]|edition=4th Edition|publisher=author-published|year=2005|id=ISBN 0-9644588-0-2}} | ||
*{{cite book | author=Schnelle, Jr., Karl B. and Dey, Partha R. | title=Atmospheric Dispersion Modeling Compliance Guide| publisher=McGraw-Hill | year=2000 | id=ISBN 0-07-058059-6 }} | *{{cite book | author=Schnelle, Jr., Karl B. and Dey, Partha R. | title=Atmospheric Dispersion Modeling Compliance Guide| publisher=McGraw-Hill | year=2000 | id=ISBN 0-07-058059-6 }}[[Category:Suggestion Bot Tag]] | ||
[[Category: | |||
Latest revision as of 11:23, 1 October 2024
Air pollution dispersion modeling is the mathematical simulation of how air pollutants disperse in the ambient atmosphere. It is performed with computer programs, called dispersion models, that solve the mathematical equations and algorithms which simulate the pollutant dispersion. The dispersion models are used to estimate or to predict the downwind concentration of air pollutants emitted from emission sources such as industrial plants and vehicular traffic. Such models are important to governmental agencies tasked with protecting and managing the ambient air quality. The models are typically employed to determine whether existing or proposed new industrial facilities are or will be in compliance with the National Ambient Air Quality Standards (NAAQS) in the United States of America and similar standards in other nations. The models also serve to assist in the design of effective control strategies to reduce emissions of harmful air pollutants.
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 referred to as the stability class, the ambient air temperature and the height to the bottom of any temperature inversion that may be present aloft.
- 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 gaseous emission 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. Currently, the AERMOD air pollution dispersion model is the preferred regulatory model of the U.S. Environmental Protection Agency.
The air pollution dispersion models are also known as atmospheric dispersion models, atmospheric diffusion models, air dispersion models and air quality 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.[1] 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[2] 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. Since personal computers also came into existence during that period, a great many computer programs for calculating the dispersion of air pollutant emissions were developed in that same period. The basis for most of those models was the Complete Equation For Gaussian Dispersion Modeling Of Continuous, Buoyant Air Pollution Plumes shown below: [3][4]
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle C = \frac{\;Q}{u}\cdot\frac{\;f}{\sigma_y\sqrt{2\pi}}\;\cdot\frac{\;g_1 + g_2 + g_3}{\sigma_z\sqrt{2\pi}}}
| where: | |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f} | = crosswind dispersion parameter |
|---|---|
| = Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \exp\;[-\,y^2/\,(2\;\sigma_y^2\;)\;]} | |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g} | = vertical dispersion parameter = Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \,g_1 + g_2 + g_3} |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g_1} | = vertical dispersion with no reflections |
| = Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \; \exp\;[-\,(z - H)^2/\,(2\;\sigma_z^2\;)\;]} | |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g_2} | = vertical dispersion for reflection from the ground |
| = Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \;\exp\;[-\,(z + H)^2/\,(2\;\sigma_z^2\;)\;]} | |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g_3} | = vertical dispersion for reflection from an inversion aloft |
| = Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sum_{m=1}^\infty\;\big\{\exp\;[-\,(z - H - 2mL)^2/\,(2\;\sigma_z^2\;)\;]} | |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle +\, \exp\;[-\,(z + H + 2mL)^2/\,(2\;\sigma_z^2\;)\;]} | |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle +\, \exp\;[-\,(z + H - 2mL)^2/\,(2\;\sigma_z^2\;)\;]} | |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle +\, \exp\;[-\,(z - H + 2mL)^2/\,(2\;\sigma_z^2\;)\;]\big\}} | |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle C} | = concentration of emissions, in g/m³, 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 | |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Q} | = source pollutant emission rate, in g/s |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle u} | = horizontal wind velocity along the plume centerline, m/s |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle H} | = height of emission plume centerline above ground level, in m |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma_z} | = vertical standard deviation of the emission distribution, in m |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma_y} | = horizontal standard deviation of the emission distribution, in m |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle L} | = height from ground level to bottom of the inversion aloft, in m |
| Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \exp}
|
= the exponential function Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle e}
, which equals approximately 2.71828, and is also known as Euler's Number |
The above equation not only includes upward reflection of the pollution plume from the ground, it also includes downward reflection from the bottom of any temperature inversion lid present in the atmosphere.
The sum of the four exponential terms in Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 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 Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma_z} and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma_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 (also known as the effective plume height, He) which is the pollutant plume's centerline height above ground level. He at any distance from the pollutant plume's source is the sum of Hs (the actual physical height of the pollutant plume's source point) plus ΔH (the plume rise due the plume's buoyancy) at that distance.
To determine ΔH, many if not most of the air dispersion models developed between the late 1960s and the early 2000s used what are known as "the Briggs equations." G.A. Briggs first published his plume rise observations and comparisons in 1965.[5] In 1968, at a symposium sponsored by CONCAWE (a Dutch organization), he compared many of the plume rise models then available in the literature.[6] In that same year, Briggs also wrote the section of the publication edited by Slade[7] 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,[8] 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.[9][10]
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 (as depicted in the adjacent diagram of a plume).
In general, Briggs's 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 fuels 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[4] to obtain the plume rise trajectory of bent-over buoyant plumes is presented below:
where: Δh = plume rise, in m F = buoyancy factor, in m4/s3 x = downwind distance from plume source, in m xf = downwind distance from plume source to point of maximum plume rise, in m u = windspeed at actual stack height, in m/s s = stability parameter, in s-2
The above parameters used in the Briggs' equations are discussed in much more detail by Beychok.[4]
Major air pollution dispersion models in current use
These are some of the major air pollution dispersion models currently being used. There are literally dozens of other models as well.
- ADMS 3: Developed in the United Kingdom
- AERMOD: Developed in the United States
- AUSPLUME: Developed in Australia
- CALPUFF: Developed in the United States
- DISPERSION2:Developed in Sweden
- ISC3: Developed in the United States
- LADM: Developed in Australia
- NAME: Developed in the United Kingdom
- MERCURE: Developed in France
- RIMPUFF: Developed in Denmark
The European Topic Centre on Air and Climate Change, which is part of the European Environment Agency (EEA), maintains an online Model Documentation System (MDS) that includes descriptions and other information for almost all of the dispersion models developed by the countries of Europe. The MDS currently contains about 140 models developed in Europe (excluding the United Kingdom).[11]
References
- ↑ Bosanquet, C.H. and Pearson, J.L., "The spread of smoke and gases from chimneys", Trans. Faraday Soc., 32:1249, 1936
- ↑ Sutton, O.G., "The problem of diffusion in the lower atmosphere", QJRMS, 73:257, 1947 and "The theoretical distribution of airborne pollution from factory chimneys", QJRMS, 73:426, 1947
- ↑ Turner, D.B. (1994). Workbook of Atmospheric Dispersion Estimates, 2nd Edition. CRC Press. ISBN 1-56670-023-X. www.crcpress.com
- ↑ 4.0 4.1 4.2 Beychok, Milton R. (2005). Fundamentals of Stack Gas Dispersion, 4th Edition. author-published. ISBN 0-9644588-0-2.
- ↑ Briggs, G.A., "A plume rise model compared with observations", JAPCA, 15:433-438, 1965
- ↑ Briggs, G.A., "CONCAWE meeting: discussion of the comparative consequences of different plume rise formulas", Atmos. Envir., 2:228-232, 1968
- ↑ Slade, D.H. (editor), "Meteorology and atomic energy 1968", Air Resources Laboratory, U.S. Dept. of Commerce, 1968
- ↑ Briggs, G.A., "Plume Rise", USAEC Critical Review Series, 1969
- ↑ Briggs, G.A., "Some recent analyses of plume rise observation", Proc. Second Internat'l. Clean Air Congress, Academic Press, New York, 1971
- ↑ Briggs, G.A., "Discussion: chimney plumes in neutral and stable surroundings", Atmos. Envir., 6:507-510, 1972
- ↑ The Model Documentation System (MDS)
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:
- Turner, D.B. (1994). Workbook of Atmospheric Dispersion Estimates, 2nd Edition. CRC Press. ISBN 1-56670-023-X. www.crcpress.com
- Beychok, Milton R. (2005). Fundamentals of Stack Gas Dispersion, 4th Edition. author-published. ISBN 0-9644588-0-2.
- Schnelle, Jr., Karl B. and Dey, Partha R. (2000). Atmospheric Dispersion Modeling Compliance Guide. McGraw-Hill. ISBN 0-07-058059-6.
- Pages using ISBN magic links
- Editable Main Articles with Citable Versions
- CZ Live
- Engineering Workgroup
- Earth Sciences Workgroup
- Chemical Engineering Subgroup
- Environmental Engineering Subgroup
- Articles written in American English
- Advanced Articles written in American English
- All Content
- Engineering Content
- Earth Sciences Content
- Chemical Engineering tag
- Environmental Engineering tag
- Pages with math errors
- Pages with math render errors