Source Terms For Accidental Discharge Flow Rates
Click in table below on desired source term (metric version):
To return to opening page and select United States version, click here
>
Gas Discharge To The Atmosphere From A Pressure Source:^{ 1, 2, 6, 7}
When gas stored under pressure in a closed vessel is discharged to the atmosphere through a hole or other opening, the gas velocity through that opening may be choked (i.e., has attained a maximum) or nonchoked. Choked velocity, which is also referred to as sonic velocity, occurs when the ratio of the absolute source pressure to the absolute downstream ambient pressure is equal to or greater than [ ( k + 1 ) / 2 ]^{ k / ( k  1 )} , where k is the specific heat ratio of the discharged gas. For many gases, k ranges from about 1.09 to about 1.41, and thus [ ( k + 1 ) / 2 ]^{ k / ( k  1 )} 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 ambient atmospheric pressure.
When the gas velocity is choked, the equation for the mass flow rate is:
or this equivalent form:
[ 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 ]^{ k / ( k  1 )}, then the gas velocity is nonchoked (i.e., subsonic)
and the equation for mass flow rate is:
or this equivalent form:
where:
^{ }
_{ }
^{ }
_{ }
^{ }

Q
C
A^{ }
k_{ }
ρ^{ }
P
P_{A}
M
R^{ }
T
Z

= mass flow rate, kg / s
= discharge coefficient (dimensionless, usually about 0.72)
= discharge hole area, m^{ 2}
= c_{p} / c_{v} of the gas = the isentropic expansion coefficient
= (specific heat at constant pressure) / (specific heat at constant volume)
= real gas density, kg / m^{ 3} at P and T
= absolute source or upstream pressure, Pa
= absolute ambient or downstream pressure, Pa _{ }
= gas molecular weight
= the Universal Gas Law Constant = 8314.5 ( Pa · m^{ 3}) / ( kgmol · °K )
= gas temperature, °K
= the gas compressibility factor at P and T (dimensionless)

The above equations calculate the initial instantaneous 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. Click here to learn how such calculations are performed.
When expressed in the customary USA units, the equations above also contain the gravitational conversion factor g_{c} ( 32.17 ft / s^{2} in USA units ). Since g_{c} is 1 ( kg · m ) / ( N · s^{2} ) in the SI metric system of units, the above equations do not include it.
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_{s} which only applies to a specific individual gas. The relationship between the two constants is R_{s} = R / (MW).
Other forms of the above equations have been published by Ramskill. When used properly, Ramskill's equations yield results which are identical to the results yielded by the above equations. Appendix 1 at the end of this article discusses Ramskill's equation for gas flow from a pressurized source.
Notes:
(1) The above equations are for a real gas.
(2) For an ideal gas, Z = 1 and d is the ideal gas density.
(3) kgmol = kilogram mole
^{Return to Top}
Liquid Discharge From A Pressurized Source Vessel:^{ 1, 2}
Initial instantaneous flow through the discharge opening:
(1) Q_{i} = C A [ ( 2 g d^{ 2} H ) + ( 2 d ) ( P  P_{A} ) ]^{ 1/2}
 
Final flow when the liquid level reaches the bottom of the discharge opening:
(2) Q_{f} = C A [ ( 2 d ) ( P  P_{A} ) ]^{ 1/2}
 
Average flow:
(3) Q_{avg} = ( Q_{i} + Q_{f} ) ÷ 2
 
where:
^{ }
^{ }
^{ }
^{ }
_{ }^{ }

Q
C
A^{ }
g^{ }
d^{ }
P^{ }
P_{A}^{ }
H

= mass flow rate, kg / s
= discharge coefficient (dimensionless, usually about 0.62)
= discharge hole area, m^{ 2}
= local gravitational acceleration of 9.807 m / s^{ 2}
= source liquid density, kg / m^{ 3}
= absolute source pressure, Pa ^{ }
= absolute ambient pressure, Pa ^{ }_{ }
= height of liquid above bottom of discharge opening, m

When expressed in the customary USA units, the equations (1) and (2) above also contain the gravitational conversion factor g_{c} ( 32.17 ft / s^{2} in USA units ) associated with the ( P  P_{A} ) term. Since g_{c} is 1 ( kgm ) / ( Ns^{2} ) in the SI metric system of units, equations (1) and (2) above do not include it.
^{Return to Top}
Liquid Discharge From A NonPressurized Source Vessel:^{ 1, 2}
Initial instantaneous flow through the discharge opening:
(1) Q_{i} = C A ( 2 g d^{ 2} H )^{1/2
}  
Final flow when the liquid level reaches the bottom of the discharge opening:
Average flow:
where: ^{ } ^{ } ^{ }

Q
C
A^{ }
g^{ }
d^{ }
H

= mass flow rate, kg / s
= discharge coefficient (dimensionless, usually about 0.62)
= discharge hole area, m^{ 2}
= local gravitational acceleration of 9.807 m / s^{ 2}
= source liquid density, kg / m^{ 3}
= height of liquid above bottom of discharge opening, m

^{Return to Top}
Evaporation From A NonBoiling Liquid Pool:
Three different methods of calculating the rate of evaporation from a nonboiling liquid pool are presented in this
section.
Method developed by the U.S. Air Force:^{ 2}
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.
(1) E = ( 4.161 x 10^{ 5} ) u^{0.75} T_{F} M ( P_{S }/ P_{H} )
 
where:^{ }
_{ }
_{ }
_{ }
_{ }
_{ }

E^{ }
u
T_{A}
T_{F}
T_{P}
M
P_{S}
P_{H}

= evaporation flux, ( kg / minute) / m^{ 2} of pool surface
= windspeed just above the liquid surface, m / s
= ambient temperature, °K_{ }
= pool liquid temperature correction factor_{ }
= pool liquid temperature, °C_{ }
= pool liquid molecular weight
= pool liquid vapor pressure at ambient temperature, mm Hg_{ }
= hydrazine vapor pressure at ambient temperature, mm Hg_{ }

(2) If T_{P} = 0 °C or less, then T_{F} = 1.0
If T_{P} > 0 °C, then T_{F} = 1.0 + 0.0043 T_{P}^{2}
 
(3) P_{H} = 760 exp[ 65.3319  ( 7245.2 / T_{A} )  ( 8.22 ln T_{A} )
+ ( 6.1557 x 10^{  3 }) T_{A} ]
 
Note: The function "ln x" is the natural logarithm (base e) of x, and the function "exp x" is the value of the
constant e (approximately 2.7183) raised to the power x.
Method developed by U.S. EPA:^{ 5, 6}
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 Environmental Protection Agency ( U.S. EPA ) using units which were a mixture of metric usage and United States usage. The nonmetric units
have been converted to metric units for this presentation.
(10.40 ) u^{ 0.78} M^{ 0.667}
A P
(1) E = —————————————
R T
 
where:
^{ }
^{ }

E
u
M
A^{ }
P
T
R^{ }

= evaporation rate, kg / minute
= windspeed just above the pool liquid surface, m / s
= molecular weight of the pool liquid
= surface area of the pool liquid, m^{ 2}
= vapor pressure of the pool liquid at the pool temperature, kPa
= pool liquid temperature, °K
= the Universal Gas Law constant = 82.05 ( atm · cm^{ 3} ) / ( gmol · °K )

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_{s} which only applies to a specific individual gas. The relationship between the two constants is R_{s} = R / (MW).
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:
(2) A = ( cubic meters of pool liquid ) / ( 0.01 m )
 
Notes:
1 kPa = 0.0102 kg / cm^{ 2} = 0.01 bar
gmol = gram mole.
atm = atmosphere
Method developed by Stiver and Mackay:^{ 3}
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.
(1) E = k P M / ( R T_{A })
(2) k = 0.002 u ^{ }
 
where:^{ }
_{ }
^{ }

E^{ }
k
T_{A}
M
P
R^{ }
u

= evaporation flux, ( kg / s ) / m^{ 2} of pool surface
= mass transfer coefficient, m / s
= ambient temperature, °K _{ }
= pool liquid molecular weight
= pool liquid vapor pressure at ambient temperature, Pa
= the Universal Gas Law constant = 8314.5 ( Pa · m^{ 3} ) / ( kgmol · °K )
= windspeed just above the liquid surface, m / s

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_{s} which only applies to a specific individual gas. The relationship between the two constants is R_{s} = R / (MW).
Note: kgmol = kilogram mole
^{Return to Top}
Evaporation From A Boiling Pool Of Cold Liquid:^{ 2}
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 zero degrees Centigrade or less ).
(1) E = ( 0.0001 ) ( 7.7026  0.0288 B ) ( M ) e^{  ( 0.0077 B )  0.1376}
 
where: ^{ }

E^{ }
B
M
e

= evaporation flux, ( kg / minute) / m^{ 2} of pool surface
= atmospheric boiling point of pool liquid, °C
= molecular weight of pool liquid
= 2.7183 ( the number that is the base of the natural logarithm system )

^{Return to Top}
Discharge Of Flashing Saturated Liquid:^{ 2, 4}
(1) Q = 0.7584 D^{ 2} P [ ln ( P / 101,325 ) ] ( T_{B} / T )
( T / c_{p} )^{1/2} ( T  T_{B} )^{  1}
 
where:
^{ } _{ }
_{ }

Q
D
P^{ }
T
T_{B}
c_{p}

= initial instantaneous mass flow, kg / s
= discharge hole diameter, m
= absolute source pressure, Pa
= source liquid temperature, °K
= atmospheric boiling point of source liquid, °K_{ }
= source liquid specific heat, J / kg / °C_{ }

Note: ln = natural logarithm (base e)
^{Return to Top}
Discharge of Flashing SubCooled Liquid:^{ 4}
(1) _{ }
(2) _{ }

Calculate the singlephase flow component ( Q_{S }) for the source liquid by using
the same equation as for a liquid discharge from a pressurized source, except
substitute the source pressure minus the source liquid vapor pressure for the
source pressure.
Calculate the flashing flow component ( Q_{F }) by using the same equation as for
a flashing saturated liquid.


(3) Q_{ Total} = ( Q_{S}^{2} + Q_{F}^{2} )
^{ 1/2}
 
where:

Q_{ Total}

= initial instantaneous mass flow, kg / s

^{Return to Top}
Adiabatic Flash of a Liquified Gas Release Into Atmosphere:
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 "adiabatic flashing" and
the following equation, derived from a simple heat balance, is used to predict how much of the liquified gas is vaporized.
(1) X = 100 ( H_{s}^{L}  H_{a}^{L} ) ÷
( H_{a}^{V}  H_{a}^{L} )
 
where:
_{ }^{ }
_{ }^{ }
_{ }^{ }

X
H_{s}^{L}
H_{a}^{V}
H_{a}^{L}

= weight percent vaporized
= source liquid enthalpy at source temperature and pressure, J / kg_{ }^{ }
= flashed vapor enthalpy at atmospheric boiling point and pressure, J / kg_{ }^{ }
= residual liquid enthalpy at atmospheric boiling point and pressure, J / kg_{ }^{ }

If the enthalpy data required for the above equation is unavailable, then the following equation may be used.
(2) X = 100 [ c_{p} ( T_{s}  T_{b} ) ] ÷ H
 
where:
_{ }
_{ }
_{ }

X
c_{p}
T_{s}
T_{b}
H

= weight percent vaporized
= source liquid specific heat, J / kg / °C _{ }
= source liquid temperature, °K _{ }
= source liquid atmospheric boiling point, °K _{ }
= source liquid heat of vaporization at atmospheric boiling point, J / kg

^{Return to Top}
References:
(1) "Perry's Chemical Engineers' Handbook, Sixth Edition, McGrawHill Co., 1984
(2) "Handbook of Chemical Hazard Analysis Procedures", Federal Emergency Management Agency, U.S. Dept. of
Transportation, and U.S. Environmental Protection Agency, 1989 (available at http://nepis.epa.gov/Adobe/PDF/10003MK5.PDF) provides references to (2a), (2b) and (2c) below
(2a) Clewell, H.J., "A Simple Method For Estimating the Source Strength Of Spills Of Toxic Liquids", Energy
Systems Laboratory, ESLTR8303, 1983 (Available at Air Force Weather Technical Library, Asheville, North Carolina)
(2b) Ille, G. and Springer, C., "The Evaporation And Dispersion Of Hydrazine Propellants From Ground Spills", Civil
and Environmental Engineering Development Office, CEEDO 7127830, 1978 (Available at Air Force Weather Technical
Library, Asheville, North Carolina)
(2c) Kahler, J.P., Curry, R.C. and Kandler, R.A., "Calculating Toxic Corridors", Air Force Weather Service, AWS
TR80/003, 1980 (Available at Air Force Weather Technical Library, Asheville, North Carolina)
(3) Stiver, W. and Mackay, D., "A Spill Hazard Ranking System For Chemicals", Environment Canada First
Technical Spills Seminar, Toronto, Canada, 1983
(4) Fauske, Hans K., "Flashing Flows: Some Guidelines For Emergency Releases", Plant/Operations Progress,
July 1985
(5) "Technical Guidance For Hazards Analysis", U.S, EPA and U.S. FEMA, December 1987 [ Equation (7),
Section G2, Appendix G. Available at http://www.epa.gov/OEM/docs/chem/tech.pdf
(6) "Risk Management Program Guidance For Offsite Consequence Analysis", U.S. EPA publication
EPA550B99009, April 1999. [ Equation (D1), Section D.2.3, and Equation (D7), Section D.6, Appendix D. Available at http://www.epa.gov/emergencies/docs/chem/ocaall.pdf ]
(7) "Methods For The Calculation Of Physical Effects Due To Releases Of Hazardous Substances (Liquids and Gases)", CPR 14E, Third Edition Second Revised Print, The Netherlands Organization Of Applied Scientific Research, The Hague, 2005. [ Equations (2.22) and (2.25) on page 2.68. ]
APPENDIX 1: THE RAMSKILL EQUATION FOR THE NONCHOKED FLOW OF GAS:
The first section in this article, entitled "Gas Discharge To The Atmosphere From A Pressure Source", includes two equivalent forms of an equation for the nonchoked flow of a gas from a pressurized source. Another form of the equation has been published by Ramskill. When used properly, Ramskill's equation yields results which are identical to the two forms in the above first section. Ramskill's equation is shown below as Equation (1):
The ideal gas density at the downstream conditions of temperature and pressure in Ramskill's equation is ρ_{A} and it is defined in Equation (2) using the ideal Gas Law:
Since the downstream temperature T_{A} is not known, the adiabatic expansion Equation (3) below is used to determine T_{A} in terms of the known upstream temperature T:
Combining Equations (2) and (3) results in Equation (4) which defines ρ_{A} in terms of the known upstream temperature T:
Using Equations (1) and (4) to determine nonchoked mass flow rates for ideal gases gives identical results to the results obtained using equation (5) below (which is one of the two equivalent equations in the "Gas Discharge To The Atmosphere From A Pressure Source" section of this article):
where:
^{ }
_{ }
^{ }
_{ }
^{ }

Q
C
A^{ }
k_{ }
ρ^{ }
P
P_{A}
M
R^{ }
T

= mass flow rate, kg / s
= discharge coefficient (dimensionless, usually about 0.72)
= discharge hole area, m^{ 2}
= c_{p} / c_{v} of the gas = the isentropic expansion coefficient
= (specific heat at constant pressure) / (specific heat at constant volume)
= gas density, kg / m^{ 3} at upstream source conditions P and T
= absolute source or upstream pressure, Pa
= absolute ambient or downstream pressure, Pa _{ }
= gas molecular weight
= the Universal Gas Law Constant = 8314.5 ( Pa · m^{ 3}) / ( kgmol · °K )
= upstream source gas temperature, °K

References for Equation (1):
Gierer, Conrad and Hyatt, Nigel,"Using Source Term Analysis Software for Calculating Fluid Flow Release Rates", Dyadem International Ltd., CACHE Newsletter No.48, Spring 1999,Austin, Texas (www.che.utexas.edu/cache/newsletters/Spr_99.pdf)
Ramskill, P.K., "Discharge Rate Calculation Methods for Use In Plant Safety Assessments", Safety and Reliability Directory, 1986, United Kingdom Atomic Energy Authority
References for Equation (5):
Handbook of Chemical Hazard Analysis Procedures" Appendix B, Federal Emergency Management Agency, U.S. Dept. of Transportation, and U.S. Environmental Protection Agency, 1989.
"Risk Management Program Guidance For Offsite Consequence Analysis", U.S. EPA publication EPA550B99009, April 1999.
"Methods For The Calculation Of Physical Effects Due To Releases Of Hazardous Substances (Liquids and Gases)", PGS2 CPR 14E, Chapter 2, Section 2.5.2.3, The Netherlands Organization Of Applied Scientific Research, The Hague, 2005.
^{Return to Top}
>> Was this article useful? My book "Fundamentals of Stack Gas Dispersion" is better yet. <<
>> Click on "How To Buy" (see below) for price and method of payment. <<
Click 
HERE 
to see the front cover of the book or 

click any of these links to visit other parts of this site:
Or click on one of these Feature Technical Articles: