|
|
| Line 1: |
Line 1: |
| '''Accidental release source terms''' are the mathematical equations that quantify the flow rate at which accidental releases of air [[pollutant]]s into the ambient [[Natural environment|environment]] can occur at industrial facilities such as [[Oil refinery|petroleum refineries]], [[petrochemical]] plants, [[natural gas]] processing plants, oil and gas transportation [[Pipeline transport|pipelines]], chemical plants, and many other industrial activities. Governmental regulations in many countries require that the probability of such accidental releases be analyzed and their quantitative impact upon the environment and human health be determined so that mitigating steps can be planned and implemented. | | The writer's title is Christy Brookins. Her family members lives in Alaska but her spouse wants them to transfer. My spouse doesn't like it the way I do but what I truly accurate psychic predictions [[http://cspl.postech.ac.kr/zboard/Membersonly/144571 cspl.postech.ac.kr]] like doing is caving but I don't have the time recently. My working day occupation is an invoicing officer but I've currently utilized for an additional one.<br><br>Check out my web-site; free psychic; [http://www.weddingwall.com.au/groups/easy-advice-for-successful-personal-development-today/ www.weddingwall.com.au], online psychics - [http://srncomm.com/blog/2014/08/25/relieve-that-stress-find-a-new-hobby/ please click the next webpage] - |
| | |
| There are a number of mathematical calculation methods for determining the flow rate at which gaseous and liquid pollutants might be released from various types of accidents. Such calculational methods are referred to as ''source terms'', and this article on accidental release source terms explains some of the calculation methods used for determining the [[mass flow rate]] at which gaseous pollutants may be accidentally released.
| |
| | |
| ==Accidental release of pressurized gas==
| |
| When gas stored under [[pressure]] in a closed vessel is discharged to the [[Earth's atmosphere|atmosphere]] through a hole or other opening, the gas [[velocity]] through that opening may be choked (i.e., it has attained a maximum) or it may be non-choked.
| |
| | |
| Choked velocity, also referred to as sonic velocity, occurs when the ratio of the absolute source pressure to the absolute downstream pressure is equal to or greater than '''[(''k'' + 1) ÷ 2 ]<sup> ''k''÷(''k'' - 1 )</sup>''', where ''k'' is the [[heat capacity ratio|specific heat ratio]] of the discharged gas (sometimes called the [[adiabatic index|isentropic expansion factor]] and sometimes denoted as <math>\gamma</math>).
| |
| | |
| For many gases, ''k'' ranges from about 1.09 to about 1.41, and therefore '''[(''k'' + 1) ÷ 2 ]<sup> ''k''÷(''k'' - 1 )</sup>''' ranges from 1.7 to about 1.9, which means that choked velocity usually occurs when the absolute source vessel pressure is at least 1.7 to 1.9 times as high as the absolute downstream ambient atmospheric pressure.
| |
| | |
| When the gas velocity is choked, the equation for the [[mass flow rate]] in SI metric units is:<ref name=Perry>''[[Perry's Chemical Engineers' Handbook]]'', Sixth Edition, McGraw-Hill Co., 1984.</ref><ref name=Hazards>''Handbook of Chemical Hazard Analysis Procedures'', Appendix B, Federal Emergency Management Agency, U.S. Dept. of Transportation, and U.S. Environmental Protection Agency, 1989. Also provides the references below:<br>– Clewell, H.J., ''A Simple Method For Estimating the Source Strength Of Spills Of Toxic Liquids'', Energy Systems Laboratory, ESL-TR-83-03, 1983.<br>– Ille, G. and Springer, C., ''The Evaporation And Dispersion Of Hydrazine Propellants From Ground Spill'', Environmental Engineering Development Office, CEEDO 712-78-30, 1978.<br>– Kahler, J.P., Curry, R.C. and Kandler, R.A.,''Calculating Toxic Corridors'' Air Force Weather Service, AWS TR-80/003, 1980.<br>[http://nepis.epa.gov/Adobe/PDF/10003MK5.PDF Handbook of Chemical Hazard Analysis, Appendix B] Scroll down to page 391 of 520 PDF pages.</ref><ref name=Risk>[http://nepis.epa.gov/Exe/ZyPDF.cgi?Dockey=1000383V.PDF "Risk Management Program Guidance For Offsite Consequence Analysis"] U.S. EPA publication EPA-550-B-99-009, April 1999. (See derivations of equations D-1 and D-7 in Appendix D)</ref><ref name=Netherlands>"Methods For The Calculation Of Physical Effects Due To Releases Of Hazardous Substances (Liquids and Gases)", PGS2 CPR 14E, Chapter 2, The Netherlands Organization Of Applied Scientific Research, The Hague, 2005. [http://vrom.nl/pagina.html?id=20725 PGS2 CPR 14E]</ref><br><br>
| |
| | |
| <math>Q\;=\;C\;A\;\sqrt{\;k\;\rho\;P\;\bigg(\frac{2}{k+1}\bigg)^{(k+1)/(k-1)}}</math>
| |
| | |
| or this equivalent form:
| |
| | |
| <math>Q\;=\;C\;A\;P\;\sqrt{\bigg(\frac{\;\,k\;M}{Z\;R\;T}\bigg)\bigg(\frac{2}{k+1}\bigg)^{(k+1)/(k-1)}}</math><br><br>
| |
| | |
| For the above equations, '''it is important to note that although the gas velocity reaches a maximum and becomes choked, the mass flow rate is not choked'''. The mass flow rate can still be increased if the source pressure is increased.
| |
| | |
| Whenever the ratio of the absolute source pressure to the absolute downstream ambient pressure is less than
| |
| '''[ ( ''k'' + 1 ) ÷ 2 ]<sup> ''k'' ÷ ( ''k'' - 1 )</sup>''', then the gas velocity is non-choked (i.e., sub-sonic) and the equation for mass flow rate is:<br><br>
| |
| | |
| <math>Q\;=\;C\;A\;\sqrt{\;2\;\rho\;P\;\bigg(\frac{k}{k-1}\bigg)\Bigg[\,\bigg(\frac{\;P_A}{P}\bigg)^{2/k}-\;\,\bigg(\frac{\;P_A}{P}\bigg)^{(k+1)/k}\;\Bigg]}</math>
| |
| | |
| or this equivalent form:
| |
| | |
| <math>Q\;=\;C\;A\;P\;\sqrt{\bigg(\frac{2\;M}{Z\;R\;T}\bigg)\bigg(\frac{k}{k-1}\bigg)\Bigg[\,\bigg(\frac{\;P_A}{P}\bigg)^{2/k}-\;\,\bigg(\frac{\;P_A}{P}\bigg)^{(k+1)/k}\;\Bigg]}</math><br><br>
| |
| | |
| {| border="0" cellpadding="2"
| |
| |-
| |
| |align=right|where:
| |
| |
| |
| |-
| |
| !align=right| ''Q''
| |
| |align=left|= [[mass flow rate]], kg/s
| |
| |-
| |
| !align=right| ''C''
| |
| |align=left|= discharge coefficient, dimensionless (usually about 0.72)
| |
| |-
| |
| !align=right| ''A''
| |
| |align=left|= discharge hole area, m²
| |
| |-
| |
| !align=right| ''k''
| |
| |align=left|= c<sub>p</sub>/c<sub>v</sub> of the gas
| |
| |-
| |
| !align=right| ''c<sub>p</sub>''
| |
| |align=left|= [[specific heat capacity|specific heat]] of the gas at constant pressure
| |
| |-
| |
| !align=right| ''c<sub>v</sub>''
| |
| |align=left|= specific heat of the gas at constant volume
| |
| |-
| |
| !align=right| ''<math>\rho</math>''
| |
| |align=left|= [[ideal gas|real gas]] [[density]] at '''''P''''' and '''''T''''', kg/m³
| |
| |-
| |
| !align=right| ''P''
| |
| |align=left|= absolute upstream pressure, Pa
| |
| |-
| |
| !align=right| ''P<sub>A</sub>''
| |
| |align=left|= absolute ambient or downstream pressure, Pa
| |
| |-
| |
| !align=right| ''M''
| |
| |align=left|= the gas [[molecular mass]], kg/kmol (also known as the molecular weight)
| |
| |-
| |
| !align=right| ''R''
| |
| |align=left|= the [[Gas constant|Universal Gas Law Constant]] = 8314.5 Pa·m³/(kmol·K)
| |
| |-
| |
| !align=right| ''T''
| |
| |align=left|= absolute upstream gas temperature, K
| |
| |-
| |
| !align=right| ''Z''
| |
| |align=left|= the gas [[Compressibility factor#Thermodynamics|compressibility factor]] at ''P'' and ''T'', dimensionless
| |
| |}
| |
| | |
| The above equations calculate the '''<u>initial instantaneous</u>''' mass flow rate for the pressure and temperature existing in the source vessel when a release first occurs. The initial instantaneous flow rate from a leak in a pressurized gas system or vessel is much higher than the average flow rate during the overall release period because the pressure and flow rate decrease with time as the system or vessel empties. Calculating the flow rate versus time since the initiation of the leak is much more complicated, but more accurate. Two equivalent methods for performing such calculations are presented and compared at [http://www.air-dispersion.com/feature2.html www.air-dispersion.com/feature2.html].
| |
| | |
| The technical literature can be very confusing because many authors fail to explain whether they are using the universal gas law constant ''R'' which applies to any [[ideal gas]] or whether they are using the gas law constant ''R''<sub>s</sub> which only applies to a specific individual gas. The relationship between the two constants is ''R''<sub>s</sub> = ''R''/''M''.
| |
| | |
| Notes:<br>
| |
| * The above equations are for a real gas.
| |
| * For an ideal gas, '''''Z''''' = 1 and '''''ρ''''' is the ideal gas density.
| |
| * 1 kilomole (kmol) = 1000 [[mole (unit)|moles]] = 1000 gram-moles = kilogram-mole.
| |
| | |
| ===Ramskill's equation for non-choked mass flow===
| |
| | |
| P.K. Ramskill's equation <ref>[http://www.che.utexas.edu/cache/newsletters/Spr_99.pdf CACHE Newsletter No.48, Spring 1999] Gierer, C. and Hyatt, N.,''Using Source Term Analysis Software for Calculating Fluid Flow Release Rates'' Dyadem International Ltd.</ref><ref>Ramskill, P.K. (1986), ''Discharge Rate Calculation Methods for Use In Plant Safety Assessments'', Safety and Reliability Directory, United Kingdom Atomic Energy Authority</ref> for the non-choked flow of an ideal gas is shown below as equation (1):<br><br>
| |
| | |
| (1) <math>Q = C \;\rho_A\;A\;\sqrt{\frac{\;\,2\;P}{\rho}\cdot\frac{k}{k-1}\cdot{\Bigg[\; 1 - {\bigg(\frac{P_A}{P}\bigg)^{(k-1)/k)}}\Bigg]}}</math><br><br>
| |
| | |
| The gas density, '''''<math>\rho</math><sub>A</sub>''''', in Ramskill's equation is the ideal gas density at the downstream conditions of temperature and pressure and it is defined in equation (2) using the [[ideal gas law]]: <br><br>
| |
| | |
| (2) <math>\rho_A = \frac{M\;P_A}{R \;T_A}</math><br><br>
| |
| | |
| Since the downstream temperature '''''T<sub>A</sub>''''' is not known, the isentropic expansion equation below <ref>[http://www.grc.nasa.gov/WWW/K-12/airplane/compexp.html Isentropic Compression or Expansion]</ref> is used to determine '''''T<sub>A</sub>''''' in terms of the known upstream temperature '''''T''''':
| |
| | |
| (3) <math>\frac{T_A}{T} = \bigg(\frac{P_A}{P}\bigg)^{(k-1)/k}</math><br><br>
| |
| | |
| Combining equations (2) and (3) results in equation (4) which defines '''''<math>\rho</math><sub>A</sub>''''' in terms of the known upstream temperature '''''T''''':<br><br>
| |
| | |
| (4) <math>\rho_A = \frac{M \;P^{\;(k-1)/k}}{R \;T \;P_A^{\ -1/k}}</math>
| |
| | |
| Using equation (4) with Ramskill's equation (1) to determine non-choked mass flow rates for ideal gases gives identical results to the results obtained using the non-choked flow equation presented in the previous section above.
| |
| | |
| ==Evaporation of non-boiling liquid pool==
| |
| Three different methods of calculating the rate of evaporation from a non-boiling liquid pool are presented in this section. The results obtained by the three methods are somewhat different.
| |
| | |
| ===The U.S. Air Force method===
| |
| The following equations are for predicting the rate at which liquid evaporates from the surface of a pool of liquid which is at or near the ambient temperature. The equations were derived from field tests performed by the U.S. Air Force with pools of liquid hydrazine. <ref name=Hazards/><br><br>
| |
| | |
| <math>E=(4.161\times{10^{-5}})\cdot u^{0.75}\cdot T_F\cdot M\cdot (P_S / P_H)</math>
| |
| | |
| {| border="0" cellpadding="2"
| |
| |-
| |
| |align=right|where:
| |
| |
| |
| |-
| |
| !align=right| ''E''
| |
| |align=left|= evaporation flux, (kg/min)/m² of pool surface
| |
| |-
| |
| !align=right| ''u''
| |
| |align=left|= windspeed just above the liquid surface, m/s
| |
| |-
| |
| !align=right| ''T<sub>A</sub>''
| |
| |align=left|= absolute ambient temperature, K
| |
| |-
| |
| !align=right| ''T<sub>F</sub>''
| |
| |align=left|= pool liquid temperature correction factor, dimensionless
| |
| |-
| |
| !align=right| ''T<sub>P</sub>''
| |
| |align=left|= pool liquid temperature, °C
| |
| |-
| |
| !align=right| ''M''
| |
| |align=left|= pool liquid molecular weight, dimensionless
| |
| |-
| |
| !align=right| ''P<sub>S</sub>''
| |
| |align=left|= pool liquid vapor pressure at ambient temperature, [[mmHg]]
| |
| |-
| |
| !align=right| ''P<sub>H</sub>''
| |
| |align=left|= hydrazine vapor pressure at ambient temperature, mmHg (see equation below)
| |
| |}<br>
| |
| | |
| '''If ''T''<sub>P</sub> = 0 °C or less, then ''T''<sub>F</sub> = 1.0'''<br>
| |
| '''If ''T''<sub>P</sub> > 0 °C, then ''T''<sub>F</sub> = 1.0 + 0.0043 ''T''<sub>P</sub><sup>2</sup>'''<br><br>
| |
| | |
| <math>P_H = 760\cdot e^{65.3319-(7245.2/T_A)-(8.22\;\ln T_a) + (0.0061557\;T_A)}</math>
| |
| | |
| {| border="0" cellpadding="2"
| |
| |-
| |
| |align=right|where:
| |
| |
| |
| |-
| |
| |align=right| <math>e</math>
| |
| |align=left|= 2.7183, the base of the natural logarithm system
| |
| |-
| |
| |align=right| <math>\ln</math>
| |
| |align=left|= natural logarithm
| |
| |}
| |
| | |
| ===The U.S. EPA method===
| |
| The following equations are for predicting the rate at which liquid evaporates from the surface of a pool of liquid which is at or near the ambient temperature. The equations were developed by the United States [[United States Environmental Protection Agency|Environmental Protection Agency]] using units which were a mixture of metric usage and United States usage.<ref name=Risk/> The non-metric units have been converted to metric units for this presentation.<br><br>
| |
| | |
| <math>E = \frac{0.1288\cdot A\cdot P\cdot M^{0.667}\cdot u^{0.78}}{T}</math><br><br> NB The constant used here is 0.284 from the mixed unit formula/2.205 lb/kg. The 82.05 become 1.0 = (ft/m)^2 x mmHg/kPa.
| |
| | |
| {| border="0" cellpadding="2"
| |
| |-
| |
| |align=right|where:
| |
| |
| |
| |-
| |
| !align=right| ''E''
| |
| |align=left|= evaporation rate, kg/min
| |
| |-
| |
| !align=right| ''u''
| |
| |align=left|= windspeed just above the pool liquid surface, m/s
| |
| |-
| |
| !align=right| ''M''
| |
| |align=left|= pool liquid molecular weight, dimensionless
| |
| |-
| |
| !align=right| ''A''
| |
| |align=left|= surface area of the pool liquid, m²
| |
| |-
| |
| !align=right| ''P''
| |
| |align=left|= vapor pressure of the pool liquid at the pool temperature, kPa
| |
| |-
| |
| !align=right| ''T''
| |
| |align=left|= pool liquid absolute temperature, K
| |
| |}<br>
| |
| | |
| The U.S. EPA also defined the pool depth as 0.01 m (i.e., 1 cm) so that the surface area of the pool liquid could be calculated as:<br><br>
| |
| | |
| :'''''A''''' = (pool volume, in m³)/(0.01)<br><br>
| |
| | |
| Notes:<br>
| |
| * 1 kPa = 0.0102 [[kgf]]/cm² = 0.01 bar
| |
| * mol = [[mole (unit)|mole]]
| |
| * atm = atmosphere
| |
| | |
| ===Stiver and Mackay's method===
| |
| The following equations are for predicting the rate at which liquid evaporates from the surface of a pool of liquid which is at or near the ambient temperature. The equations were developed by Warren Stiver and Dennis Mackay of the Chemical Engineering Department at the University of Toronto. <ref name=Stiver>Stiver, W. and Mackay, D., ''A Spill Hazard Ranking System For Chemicals'', Environment Canada First Technical Spills Seminar, Toronto, Canada, 1993.</ref><br><br>
| |
| | |
| <math>E = \frac{k\cdot P\cdot M}{R\cdot T_A}</math>
| |
| | |
| {| border="0" cellpadding="2"
| |
| |-
| |
| |align=right|where:
| |
| |
| |
| |-
| |
| !align=right| ''E''
| |
| |align=left|= evaporation flux, (kg/s)/m² of pool surface
| |
| |-
| |
| !align=right| ''k''
| |
| |align=left|= mass transfer coefficient, m/s = 0.002 '''''u'''''
| |
| |-
| |
| !align=right| ''T<sub>A</sub>''
| |
| |align=left|= absolute ambient temperature, K
| |
| |-
| |
| !align=right| ''M''
| |
| |align=left|= pool liquid molecular weight, dimensionless
| |
| |-
| |
| !align=right| ''P''
| |
| |align=left|= pool liquid vapor pressure at ambient temperature, Pa
| |
| |-
| |
| !align=right| ''R''
| |
| |align=left|= the universal gas law constant = 8314.5 Pa·m³/(kmol·K)
| |
| |-
| |
| !align=right| ''u''
| |
| |align=left|= windspeed just above the liquid surface, m/s
| |
| |}<br>
| |
| | |
| ==Evaporation of boiling cold liquid pool==
| |
| The following equation is for predicting the rate at which liquid evaporates from the surface of a pool of cold liquid (i.e., at a liquid temperature of about 0 °C or less). <ref name=Hazards/><br><br>
| |
| | |
| <math>E = ( 0.0001\;M) (7.7026-0.0288\;B)\;e^{-(0.0077\;B)-0.1376}</math>
| |
| | |
| {| border="0" cellpadding="2"
| |
| |-
| |
| |align=right|where:
| |
| |
| |
| |-
| |
| !align=right| ''E''
| |
| |align=left|= evaporation flux, (kg/min)/m² of pool surface
| |
| |-
| |
| !align=right| ''B''
| |
| |align=left|= pool liquid [[atmospheric boiling point]], °C
| |
| |-
| |
| !align=right| ''M''
| |
| |align=left|= pool liquid molecular weight, dimensionless
| |
| |-
| |
| !align=right| ''e''
| |
| |align=left|= the base of the natural logarithm system = 2.7183
| |
| |}
| |
| | |
| ==Adiabatic flash of liquified gas release==
| |
| Liquified gases such as ammonia or chlorine are often stored in cylinders or vessels at ambient temperatures and pressures well above atmospheric pressure. When such a liquified gas is released into the ambient atmosphere, the resultant reduction of pressure causes some of the liquified gas to vaporize immediately. This is known as [[flash evaporation|"adiabatic flashing"]] and the following equation, derived from a simple heat balance, is used to predict how much of the liquified gas is vaporized.<br><br>
| |
| | |
| <math>X = 100\;\frac{H_s^L - H_a^L}{H_a^V - H_a^L}</math>
| |
| | |
| {| border="0" cellpadding="2"
| |
| |-
| |
| |align=right|where:
| |
| |
| |
| |-
| |
| !align=right| ''X''
| |
| |align=left|= weight percent vaporized
| |
| |-
| |
| !align=right| ''H<sub>s</sub><sup>L</sup>''
| |
| |align=left|= source liquid [[enthalpy]] at source temperature and pressure, J/kg<br>
| |
| |-
| |
| !align=right| ''H<sub>a</sub><sup>V</sup>''
| |
| |align=left|= flashed vapor enthalpy at atmospheric boiling point and pressure, J/kg
| |
| |-
| |
| !align=right| ''H<sub>a</sub><sup>L</sup>''
| |
| |align=left|= residual liquid enthalpy at atmospheric boiling point and pressure, J/kg
| |
| |}<br>
| |
| | |
| If the enthalpy data required for the above equation is unavailable, then the following equation may be used.<br><br>
| |
| | |
| <math>X = 100\cdot c_p\cdot (T_s - T_b)/H</math>
| |
| | |
| {| border="0" cellpadding="2"
| |
| |-
| |
| |align=right|where:
| |
| |
| |
| |-
| |
| !align=right| ''X''
| |
| |align=left|= weight percent vaporized
| |
| |-
| |
| !align=right| ''c<sub>p</sub>''
| |
| |align=left|= source liquid [[specific heat]], J/(kg °C)
| |
| |-
| |
| !align=right| ''T<sub>s</sub>''
| |
| |align=left|= source liquid absolute temperature, K
| |
| |-
| |
| !align=right| ''T<sub>b</sub>''
| |
| |align=left|= source liquid absolute atmospheric boiling point, K
| |
| |-
| |
| !align=right| ''H''
| |
| |align=left|= source liquid [[heat of vaporization]] at atmospheric boiling point, J/kg
| |
| |}
| |
| | |
| ==See also==
| |
| *[[Choked flow]]
| |
| *[[Orifice plate]]
| |
| *[[Flash evaporation]]
| |
| | |
| ==References==
| |
| {{reflist}}
| |
| | |
| ==External links==
| |
| *[http://www.che.utexas.edu/cache/newsletters/Spr_99.pdf Ramskill's equations] are presented and cited in this pdf file (use search function to find "Ramskill").
| |
| *More release source terms are available in the feature articles at [http://www.air-dispersion.com www.air-dispersion.com]
| |
| *[http://www.okcc.com/PDF/Choked%20Flow%20of%20Gases%20pg.48.pdf Choked flow of gases]
| |
| *[http://www.qub.ac.uk/qc/webpages/whatwedo/researchgroups/environmentalmodelling/ia/documents/chapter5.pdf Development of source emission models]
| |
| | |
| [[Category:Atmospheric dispersion modeling]]
| |
| [[Category:Air pollution]]
| |
| [[Category:Chemical engineering]]
| |