User:Milton Beychok/Sandbox: Difference between revisions

From Citizendium
Jump to navigation Jump to search
imported>Milton Beychok
No edit summary
imported>Milton Beychok
No edit summary
Line 1: Line 1:
The '''AERMOD''' [[atmospheric dispersion modeling]] system is an integrated system that includes three modules:<ref name=EPA-454>[http://www.epa.gov/scram001/7thconf/aermod/aermod_mfd.pdf AERMOD:Description of Model Formulation]</ref><ref name=EP>Prater, E.T. and Midgley, C., ''A new air dispersion modeling system is helping create more accurate industrial source models'', Environmental Protection, Vol. 17, No. 3, Stevens Publishing</ref><ref>[http://www.cleanairinfo.com/regionalstatelocalmodelingworkshop/documents/Brode_AERMOD_TechnicalForum.pdf#search=%22Brode%20AERMOD%20technical%20forum%20%22www.cleanairinfo.com%22%22 Brode, R.W., ''AERMOD Technical Forum'', EPA R/S/L Modelers Workshop], San Diego, California, April 16, 2006</ref>
[[Image:Vap-Liq Separator.png|thumb|right|246px|A typical flash drum]]
'''Flash (or partial) [[evaporation]]''' is the partial [[vaporization]] that occurs when a [[Boiling point|saturated liquid]] stream undergoes a reduction in pressure by passing through a [[thermal expansion valve|throttling valve]] or other throttling device. This process is one of the simplest [[unit operation]]s. If the throttling valve or device is located at the entry into a [[pressure vessel]] so that the flash evaporation occurs within the vessel, then the vessel is often referred to as a [[Vapor-liquid separator|flash drum]].


*A steady-state dispersion model designed for short-range (up to 50 kilometers) dispersion of [[air pollution|air pollutant]] [[emissions]] from [[Air pollution dispersion terminology|stationary industrial sources]].
If the saturated liquid is a single-component liquid (for example,  liquid [[propane]] or liquid [[ammonia]]), a part of the liquid immediately "flashes" into vapor. Both the vapor and the residual liquid are cooled to the [[saturation temperature]] of the liquid at the reduced pressure. This is often referred to as "auto-refrigeration" and is the basis of most conventional [[vapor compression refrigeration]] systems.
* A [[meteorology|meteorological]] data [[preprocessor]] ('''AERMET''') that accepts surface meteorological data, upper air [[Rawinsonde|soundings]], and optionally, data from on-site instrument towers. It then calculates [[Earth's atmosphere|atmospheric]] parameters needed by the dispersion model, such as atmospheric [[Air pollution dispersion terminology|turbulence characteristics]], [[Air pollution dispersion terminology|mixing heights]], [[friction velocity]], [[Air pollution dispersion terminology|Monin-Obukov length]] and [[surface heat flux]].
* A [[terrain]] preprocessor ('''AERMAP''') whose main purpose is to provide a physical relationship between terrain features and the behavior of air pollution plumes. It generates location and height data for each receptor location.  It also provides information that allows the dispersion model to simulate the effects of air flowing over hills or splitting to flow around hills.


AERMOD also includes '''PRIME''' (Plume Rise Model Enhancements)<ref>[http://www.epa.gov/scram001/7thconf/iscprime/tekpapr1.pdf Development and Evaluation of the PRIME Plume Rise and Building Downwash Model]</ref> which is an algorithm  for modeling the effects of [[Air pollution dispersion terminology|downwash]] created by the pollution plume flowing over nearby buildings.
If the saturated liquid is a multi-component liquid (for example, a mixture of [[propane]], [[isobutane]] and normal [[butane]]), the flashed vapor is richer in the more [[volatility|volatile]] components than is the remaining liquid.<br><br>


==History of the development of AERMOD==
==Flash evaporation of a single-component liquid==
The flash evaporation of a single-component liquid is an [[isoenthalpic|isenthalpic]] (i.e., constant [[enthalpy]]) process and is often referred to as an '''[[Adiabatic process|adiabatic]] flash'''. The following equation, derived from a simple heat balance around the throttling valve or device, is used to predict how much of a single-component  liquid is vaporized.<br><br>


AERMOD was developed by the AERMIC ([[American Meteorological Society|AMS]]/[[EPA]] Regulatory Model Improvement Committee), a collaborative working group of scientists from the AMS and the EPA.<ref name=EPA-454/> The AERMIC was initially formed in 1991.
:'''''X = 100 ( H<sub>u</sub><sup>L</sup>  – H<sub>d</sub><sup>L</sup> ) &divide; ( H<sub>d</sub><sup>V</sup>  – H<sub>d</sub><sup>L</sup> )'''''


The AERMIC developed AERMOD in seven steps:
:{| border="0" cellpadding="2"
|-
|align=right|where:
|&nbsp;
|-
!align=right|''X''
|align=left|= &nbsp;weight percent vaporized
|-
!align=right|''H<sub>u</sub><sup>L</sup>''
|align=left|= &nbsp;upstream liquid enthalpy  at upstream temperature and pressure, J/kg
|-
!align=right|''H<sub>d</sub><sup>V</sup>''<br>&nbsp;
|align=left|= &nbsp;flashed vapor enthalpy at downstream pressure and corresponding saturation<br> &nbsp; &nbsp; temperature, J/kg
|-
!align=right|''H<sub>d</sub><sup>L</sup>''<br>&nbsp;
|align=left|= &nbsp;residual liquid enthalpy at downstream pressure and corresponding saturation<br> &nbsp; &nbsp; temperature, J/kg
|}


*Initial model formulation
If the enthalpy data required for the above equation is unavailable, then the following equation may be used.<br><br>
*Developmental evaluation
*Internal [[peer review]] and [[software testing|beta testing]]
*Revised model formulation
*Performance evaluation and sensitivity testing
*External peer review
*Submission to the EPA for consideration as a regulatory model.


On April 21st of 2000, the EPA proposed that AERMOD be adopted as the EPA's preferred regulatory model for both simple and complex terrain.<ref>[http://frwebgate.access.gpo.gov/cgi-bin/getdoc.cgi?dbname=2000_register&docid=00-4235-filed Federal Register: April 21, 2000 (Volume 65, Number 78)] Proposed Rule</ref> On November 9th of 2005, AERMOD was adopted by the EPA and promulgated as their preferred regulatory model, effective as of December 9th of 2005.
:'''''X = 100 &middot; c<sub>p</sub> ( T<sub>u</sub> – T<sub>d</sub> ) &divide; H<sub>v</sub>'''''
<ref>[http://a257.g.akamaitech.net/7/257/2422/01jan20051800/edocket.access.gpo.gov/2005/05-21627.htm Federal Register: November 9, 2006 (Volume 70, Number 216)] Final Rule</ref> The entire developmental and adoption process took 14 years (from 1991 to 2005).


==Features and capabilities of AERMOD==
:{| border="0" cellpadding="2"
|-
|align=right|where:
|&nbsp;
|-
!align=right|''X''
|align=left|= &nbsp;weight percent vaporized
|-
!align=right|''c<sub>p</sub>''
|align=left|= &nbsp;liquid [[specific heat capacity|specific heat]] at upstream temperature and pressure, J/(kg&nbsp;°C)
|-
!align=right|''T<sub>u</sub>''
|align=left|= &nbsp;upstream liquid temperature, °C
|-
!align=right|''T<sub>d</sub>''
|align=left|= &nbsp;liquid [[saturation temperature]] corresponding to the downstream pressure, °C
|-
!align=right|''H<sub>v</sub>''<br>&nbsp;
|align=left|= &nbsp;liquid [[heat of vaporization]] at downstream pressure and corresponding saturation<br> &nbsp; &nbsp; temperature, J/kg
|}


Some of the primary features and capabilities of AERMOD are:<ref name=EP/><ref>[http://www.epa.gov/scram001/7thconf/aermod/aermod_mep.pdf AERMOD:Latest Features and Evaluation Results] (EPA-454/R-03-003)</ref><ref>[http://mcvehil-monnett.com/nwsltr4q2000.pdf McVehil-Monnet Associates ''Air Quality Newsletter'']</ref>
( Note: The words "upstream" and "downstream" refer to before and after the liquid passes through the throttling valve or device.)


*Source types: Multiple [[Air pollution dispersion terminology|point, area and volume sources]]
This type of flash evaporation is used in the [[desalination]] of brackish water or ocean water by "[[Multi-Stage Flash]] Distillation." The water is heated and then routed into a reduced-pressure flash evaporation "stage" where some of the water flashes into steam. This steam is subsequently condensed into salt-free water. The residual salty liquid from that first stage is introduced into a second flash evaporation stage at a pressure lower than the first stage pressure. More water is flashed into steam which is also subsequently condensed into more salt-free water. This sequential use of multiple flash evaporation stages is continued until the design objectives of the system are met. A large part of the world's installed desalination capacity uses multi-stage flash distillation. Typically such plants have 24 or more sequential stages of flash evaporation.
*Source releases: Surface, near surface and elevated [[Air pollution dispersion terminology|emission sources]]
*Source locations: Urban or rural locations. Urban effects are scaled by population.
*Plume types: Continuous, [[Air pollution dispersion terminology|buoyant plumes]]
*Plume deposition: [[Air pollution dispersion terminology|Dry or wet deposition]] of [[Particulate matter|particulates]] and/or [[gas]]es
*[[Atmospheric dispersion modeling|Plume dispersion]] treatment: [[Air pollution dispersion terminology|Gaussian model]] treatment in horizontal and in vertical for stable atmospheres. Non-Gaussian treatment in vertical for unstable atmospheres
*Terrain types: Simple or complex terrain
*Building effects: Handled by PRIME downwash algorithms 
*Meteorology data height levels: Accepts meteorology data from multiple heights
*Meteorological data profiles: Vertical profiles of wind, turbulence and temperature are created


==References==
==Equilibrium flash of a multi-component liquid==
<references/>
 
The '''equilibrium flash''' of a multi-component liquid may be visualized as a simple [[distillation]] process using a single [[equilibrium stage]]. It is very different and more complex than the flash evaporation of single-component liquid. For a multi-component liquid, calculating the amounts of flashed vapor and residual liquid in equilibrium with each other at a given temperature and pressure requires a trial-and-error [[Iterative method|iterative]] solution. Such a calculation is commonly referred to as an equilibrium flash calculation. It involves solving the following '''Rachford Rice equation''':<ref>
[http://faculty.mcneese.edu/pdauten/lec5.ppt McNeese University lecture] (scroll down to Rachford Rice Equation)</ref><ref>[http://pubs.acs.org/cgi-bin/abstract.cgi/iecred/1993/32/i07/f-pdf/f_ie00019a029.pdf?sessid=6006l3 ''Polynomial Objective Functions for Flash Calculations''], John H. Warrent and Michael A. Adewumi, Ind. Eng. Chem. Res. 1993, Vol.32, pages 1528-1530</ref><ref>[http://ascend.cheme.cmu.edu/ftp/pdfPapersRptsSlides/tech_complementarity.pdf "Complementarity Formulation for the Representation of Algebraic Systems Containing Conditional Equations", V. Rico-Ramírez and A. W. Westerberg] Institute for Complex Engineered Systems,
Carnegie Mellon University, Pittsburgh, Pennsylvania, Technical Report ICES 06-243-98</ref><ref>[http://www.infochemuk.com/publicat/poster.pdf Automatic Plotting of Multiple Phase Boundaries and Flash Calculations], Infochem Computer Services, [[United Kingdom]]</ref><ref>[http://www.mathworks.com/matlabcentral/fileexchange/loadFile.do?objectId=9341&objectType=file Flash Calculations using the Soave-Redlich-Kwong equation of state] (view full-size image)</ref>
 
:<math>\sum_i\frac{z_i\, (K_i - 1)}{1 + a\, (K_i - 1)}=0</math>
 
:Overall material balance equation :
 
:<math>F = V + L \,</math>
 
:Material balance equation for any component <math>i</math> :
 
:<math>F\!\cdot\! z_i = V\!\cdot\! y_i + L\!\cdot\! x_i</math>
 
:Equation defining the [[vapor-liquid equilibrium]] constant <math>K_i</math> :


== Further reading==
:<math>K_i = y_i / x_i \,</math>


For those who would like to learn more about this topic, it is suggested that either one of the following books be read:
:{| border="0" cellpadding="2"
|-
|align=right|where:
|&nbsp;
|-
!align=right|''<math>F \,</math>''
|align=left|<math>= moles\; of\; total\; feed\; liquid</math>
|-
!align=right|''<math>V \,</math>''
|align=left|<math>= moles\; of\; flashed\; vapor</math>
|-
!align=right|''<font style="vertical-align:-25%;"><math>a \,</math></font>''
|align=left|<math>= fraction\; of\; feed\; that\; is\; vaporized =  \left(\frac{V}{F}\right)</math>
|-
!align=right|''<math>L \,</math>''
|align=left|<math>= moles\; of\; residual\; liquid</math>
|-
!align=right|''<math>K_i \,</math>''
|align=left|<math>= equilibrium\;constant</math>
|-
!align=right|''<font style="vertical-align:-30%;"><math>z_i \,</math></font>''
|align=left|<math>= mole\; fraction\; of\; component\; i\; in\; the\; feed\; liquid</math>
|-
!align=right|''<font style="vertical-align:-30%;"><math>y_i \,</math></font>''
|align=left|<math>= mole\; fraction\; of\; component\; i\; in\; the\; flashed\; vapor</math>
|-
!align=right|''<font style="vertical-align:-30%;"><math>x_i \,</math></font>''
|align=left|<math>= mole\; fraction\; of\; component\; i\; in\; the\; residual\; liquid</math>
|}


*{{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]
[[Newton's method]] (also known as the Newton-Raphson method) is an efficient iterative algorithm for solving the above Rachford Rice equation. Alternatively, an [[Microsoft Excel|Excel]] spread sheet and the Excel Solver function can be used.


*{{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}} [http://www.air-dispersion.com www.air-dispersion.com]
The equilibrium flash of multi-component liquids is very widely utilized in [[Oil refinery|petroleum refineries]], [[petrochemical]] and [[chemical plant]]s and [[natural gas processing]] plants.


*{{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 }}
(See [[mole fraction]] for a definition of that terminology).


==External links==
==Spray drying==
[[Spray drying]] is the rapid drying of a [[slurry]] of very small solids suspended in a liquid. The slurry is first [[atomization|atomized]] into very small liquid droplets which are then sprayed into a stream of hot dry air. The liquid rapidly evaporates leaving behind dry powder or dry solid granules. The dry powder or solid granules are recovered from the exhaust air by using [[cyclone (industry)|cyclones]], [[Filtration|bag filters]] or [[electrostatic precipitator]]s.


*[http://www.epa.gov/scram001/dispersion_prefrec.htm#aermod The EPA's download site for AERMOD] (Model code, User's Guide and other material)
A brief explanation of spray drying has been included here because some readers may consider spray drying to be a form of flash evaporation. However, although it is a form of liquid evaporation, it is quite different from flash evaporation.
*[http://ams.confex.com/ams/pdfpapers/38497.pdf Brode, R.W., ''Implementation and Evaluation of the AERMOD-PRIME Model''] (AMS Conference, May 21, 2002)


==See also==
==See also==
*[[Bibliography of atmospheric dispersion modeling]]
*[[Evaporation]]
*[[Atmospheric dispersion modeling]]
*[[Evaporator]]
*[[Compilation of atmospheric dispersion models]]
*[[Vapor-liquid separator]]
 
==References==
{{reflist}}
 
==External links==
*[http://www.tlv.com/global/US/steam-theory/flash-steam-and-vapor.html Vapor and Flash Steam] Animation, photos and technical explanation of the difference between Flash Steam and Vapor.
*[http://www.spiraxsarco.com/resources/steam-engineering-tutorials/condensate-recovery/flash-steam.asp Flash Steam Tutorial] The benefits of recovering flash steam, how it is done and typical applications.
*[http://www.escwa.org.lb/information/publications/edit/upload/tech-01-3-e.pdf Water Desalination Technologies] in the Middle East and Western Asia
*[http://www.niro.com/ndk_website/niro/cmsresources.nsf/filenames/spray-drying.pdf/$file/spray-drying.pdf Discussion of spray drying]
*[http://petrochemical.gronerth.com/flash_distillation_program_explorer.htm Flash evaporation program online] Flash distillation of the hydrocarbon compounds.

Revision as of 01:16, 4 July 2008

A typical flash drum

Flash (or partial) evaporation is the partial vaporization that occurs when a saturated liquid stream undergoes a reduction in pressure by passing through a throttling valve or other throttling device. This process is one of the simplest unit operations. If the throttling valve or device is located at the entry into a pressure vessel so that the flash evaporation occurs within the vessel, then the vessel is often referred to as a flash drum.

If the saturated liquid is a single-component liquid (for example, liquid propane or liquid ammonia), a part of the liquid immediately "flashes" into vapor. Both the vapor and the residual liquid are cooled to the saturation temperature of the liquid at the reduced pressure. This is often referred to as "auto-refrigeration" and is the basis of most conventional vapor compression refrigeration systems.

If the saturated liquid is a multi-component liquid (for example, a mixture of propane, isobutane and normal butane), the flashed vapor is richer in the more volatile components than is the remaining liquid.

Flash evaporation of a single-component liquid

The flash evaporation of a single-component liquid is an isenthalpic (i.e., constant enthalpy) process and is often referred to as an adiabatic flash. The following equation, derived from a simple heat balance around the throttling valve or device, is used to predict how much of a single-component liquid is vaporized.

X = 100 ( HuL – HdL ) ÷ ( HdV – HdL )
where:  
X =  weight percent vaporized
HuL =  upstream liquid enthalpy at upstream temperature and pressure, J/kg
HdV
 
=  flashed vapor enthalpy at downstream pressure and corresponding saturation
    temperature, J/kg
HdL
 
=  residual liquid enthalpy at downstream pressure and corresponding saturation
    temperature, J/kg

If the enthalpy data required for the above equation is unavailable, then the following equation may be used.

X = 100 · cp ( Tu – Td ) ÷ Hv
where:  
X =  weight percent vaporized
cp =  liquid specific heat at upstream temperature and pressure, J/(kg °C)
Tu =  upstream liquid temperature, °C
Td =  liquid saturation temperature corresponding to the downstream pressure, °C
Hv
 
=  liquid heat of vaporization at downstream pressure and corresponding saturation
    temperature, J/kg

( Note: The words "upstream" and "downstream" refer to before and after the liquid passes through the throttling valve or device.)

This type of flash evaporation is used in the desalination of brackish water or ocean water by "Multi-Stage Flash Distillation." The water is heated and then routed into a reduced-pressure flash evaporation "stage" where some of the water flashes into steam. This steam is subsequently condensed into salt-free water. The residual salty liquid from that first stage is introduced into a second flash evaporation stage at a pressure lower than the first stage pressure. More water is flashed into steam which is also subsequently condensed into more salt-free water. This sequential use of multiple flash evaporation stages is continued until the design objectives of the system are met. A large part of the world's installed desalination capacity uses multi-stage flash distillation. Typically such plants have 24 or more sequential stages of flash evaporation.

Equilibrium flash of a multi-component liquid

The equilibrium flash of a multi-component liquid may be visualized as a simple distillation process using a single equilibrium stage. It is very different and more complex than the flash evaporation of single-component liquid. For a multi-component liquid, calculating the amounts of flashed vapor and residual liquid in equilibrium with each other at a given temperature and pressure requires a trial-and-error iterative solution. Such a calculation is commonly referred to as an equilibrium flash calculation. It involves solving the following Rachford Rice equation:[1][2][3][4][5]

Overall material balance equation :
Material balance equation for any component  :
Equation defining the vapor-liquid equilibrium constant  :
where:  

Newton's method (also known as the Newton-Raphson method) is an efficient iterative algorithm for solving the above Rachford Rice equation. Alternatively, an Excel spread sheet and the Excel Solver function can be used.

The equilibrium flash of multi-component liquids is very widely utilized in petroleum refineries, petrochemical and chemical plants and natural gas processing plants.

(See mole fraction for a definition of that terminology).

Spray drying

Spray drying is the rapid drying of a slurry of very small solids suspended in a liquid. The slurry is first atomized into very small liquid droplets which are then sprayed into a stream of hot dry air. The liquid rapidly evaporates leaving behind dry powder or dry solid granules. The dry powder or solid granules are recovered from the exhaust air by using cyclones, bag filters or electrostatic precipitators.

A brief explanation of spray drying has been included here because some readers may consider spray drying to be a form of flash evaporation. However, although it is a form of liquid evaporation, it is quite different from flash evaporation.

See also

References

  1. McNeese University lecture (scroll down to Rachford Rice Equation)
  2. Polynomial Objective Functions for Flash Calculations, John H. Warrent and Michael A. Adewumi, Ind. Eng. Chem. Res. 1993, Vol.32, pages 1528-1530
  3. "Complementarity Formulation for the Representation of Algebraic Systems Containing Conditional Equations", V. Rico-Ramírez and A. W. Westerberg Institute for Complex Engineered Systems, Carnegie Mellon University, Pittsburgh, Pennsylvania, Technical Report ICES 06-243-98
  4. Automatic Plotting of Multiple Phase Boundaries and Flash Calculations, Infochem Computer Services, United Kingdom
  5. Flash Calculations using the Soave-Redlich-Kwong equation of state (view full-size image)

External links