View source for Choked flow

'''Choked flow''' is a compressible flow effect.  The parameter that becomes "choked" or "limited" is the fluid velocity.

Choked flow is a [[Fluid dynamics|fluid dynamic]] condition associated with the [[Venturi effect]]. When a flowing fluid at a given [[pressure]] and [[temperature]] passes through a constriction (such as the throat of a [[convergent-divergent nozzle]] or a [[valve]] in a [[pipe (material)|pipe]]) into a lower pressure environment the fluid velocity increases. At initially subsonic upstream conditions, the [[conservation of mass]] principle requires the fluid [[velocity]] to increase as it flows through the smaller cross-sectional area of the constriction.  At the same time, the Venturi effect causes the static pressure, and therefore the density, to decrease at the constriction. '''Choked flow''' is a limiting condition where the mass flow will not increase with a further decrease in the downstream pressure environment for a fixed upstream pressure and temperature.

For homogeneous fluids, the physical point at which the choking occurs for [[adiabatic process|adiabatic]] conditions, is when the exit plane velocity is at [[speed of sound|sonic]] conditions; i.e., at a [[Mach number]] of 1.<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. [http://nepis.epa.gov/Exe/ZyNET.exe/10003MK5.TXT?ZyActionD=ZyDocument&Client=EPA&Index=1986+Thru+1990&Docs=&Query=&Time=&EndTime=&SearchMethod=1&TocRestrict=n&Toc=&TocEntry=&QField=pubnumber%5E%22OSWERHCHAP%22&QFieldYear=&QFieldMonth=&QFieldDay=&UseQField=pubnumber&IntQFieldOp=1&ExtQFieldOp=1&XmlQuery=&File=D%3A%5Czyfiles%5CIndex%20Data%5C86thru90%5CTXT%5C00000003%5C10003MK5.TXT&User=ANONYMOUS&Password=anonymous&SortMethod=h%7C-&MaximumDocuments=10&FuzzyDegree=0&ImageQuality=r75g8/r75g8/x150y150g16/i425&Display=p%7Cf&DefSeekPage=x&SearchBack=ZyActionL&Back=ZyActionS&BackDesc=Results%20page&MaximumPages=1&ZyEntry=1&SeekPage=x Handbook of Chemical Hazard Analysis, Appendix B] Click on PDF icon, wait and then scroll down to page 391 of 520 PDF pages.</ref><ref name=PGS2>''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] {{webarchive|url=https://web.archive.org/web/20070809204358/http://www.vrom.nl/pagina.html?id=20725 |date=2007-08-09 }}</ref> At choked flow, the mass flow rate can be increased only by increasing density upstream and at the choke point.

The choked flow of gases is useful in many engineering applications because the mass flow rate is independent of the downstream pressure, and depends only on the temperature and pressure and hence the density of the gas on the upstream side of the restriction. Under choked conditions, valves and calibrated [[orifice plate]]s can be used to produce a desired mass flow rate.

==Choked flow in liquids==
If the fluid is a liquid, a different type of limiting condition (also known as choked flow) occurs when the [[Venturi effect]] acting on the liquid flow through the restriction causes a decrease of the liquid pressure beyond the restriction to below that of the liquid's [[vapor pressure]] at the prevailing liquid temperature. At that point, the liquid will partially [[Flash evaporation|flash]] into bubbles of vapor and the subsequent collapse of the bubbles causes [[cavitation]]. Cavitation is quite noisy and can be sufficiently violent to physically damage valves, pipes and associated equipment. In effect, the vapor bubble formation in the restriction prevents the flow from increasing any further.<ref>[http://www.documentation.emersonprocess.com/groups/public/documents/brochures/d351912x012.pdf Read page 2 of this brochure.]</ref><ref>[http://www.chemicalprocessing.com/Media/MediaManager/control_valves.pdf Control Valve Handbook] Search document for "Choked".</ref>

==Mass flow rate of a gas at choked conditions==

All gases flow from upstream higher pressure sources to downstream lower pressure sources. There are several situations in which choked flow occurs, such as the change of the cross section in a [[de Laval nozzle]] or the flow through an [[orifice plate]]. Here the most important part is where to calculate the choked velocity: at upstream or downstream of a nozzle or orifice. The choked velocity is always observed at upstream of an orifice or nozzle and this velocity is usually less than the speed of sound in air. Another important aspect is that this is the actual velocity of the upstream fluid. Hence, the upstream actual volumetric flow rate, when expanded to downstream pressure, will result in a more actual volumetric flow for the downstream condition. Thus, the overall leakage rate when measured at downstream conditions needs to take care of this fact. When this choked velocity has reached the mass flow rate from upstream to downstream, it can still be increased if the upstream pressure is increased. However, this value of the choked velocity will keep the actual volumetric flow rate (Actual Gas Flow rate, and hence velocity) the same irrespective of the upstream pressure, provided that choked flow conditions prevail. 

===Choking in change of cross section flow===
Assuming ideal gas behaviour, steady-state choked flow occurs when the downstream pressure falls below a critical value <math>p^{*}</math>. That critical value can be calculated from the dimensionless critical pressure ratio equation<ref>Potter & Wiggert, 2010, ''Mechanics of Fluids'', 3rd SI ed., Cengage.</ref>
:<math>\frac{p^*}{p_0} = \left(\frac{2}{\gamma + 1}\right)^\frac{\gamma}{\gamma - 1}</math>,
where <math>\gamma</math> is the [[heat capacity ratio]] <math>c_p/c_v</math> of the gas and where <math>p_0</math> is the total (stagnation) upstream pressure.

For air with a heat capacity ratio <math>\gamma = 1.4</math>, then <math>p^* = 0.528 p_0</math>; other gases have <math>\gamma</math> in the range 1.09 (e.g. butane) to 1.67 (monatomic gases), so the critical pressure ratio varies in the range <math>0.487 < p^*/p_0 < 0.587</math>, which means that, depending on the gas, choked flow usually occurs when the downstream static pressure drops to below 0.487 to 0.587 times the absolute pressure in stagnant upstream source vessel.

When the gas velocity is choked, the equation for the [[mass flow rate]] is:<ref name=Perry/><ref name=Hazards/><ref name=PGS2/>
: <math>\dot{m} = C_d A \sqrt{\gamma \rho_0 P_0 \left(\frac{2}{\gamma + 1}\right)^\frac{\gamma + 1}{\gamma - 1}}</math>

{| border="0" cellpadding="2"
|-
| align=right | Where:
| &nbsp;
|-
! align=right | <math>\dot{m}</math>,
| align=left  | [[mass flow rate]], in kg/s
|-
! align=right | <math>C_d</math>,
| align=left  | [[discharge coefficient]], dimensionless
|-
! align=right | <math>A</math>,
| align=left  | discharge hole cross-sectional area, in m²
|-
! align=right | <math>\gamma</math>,
| align=left  | <math>\frac{c_p}{c_v}</math> of the gas
|-
! align=right | <math>c_p</math>, 
| align=left  | [[specific heat capacity|specific heat]] of the gas at constant pressure
|-
! align=right | <math>c_v</math>,
| align=left  | specific heat of the gas at constant volume
|-
! align=right | <math>\rho_0</math>,
| align=left  | real gas (total) density at total pressure <math>P_0</math> and total temperature <math>T_0</math>, in kg/m³
|-
! align=right | <math>P_0</math>,
| align=left  | absolute upstream total pressure of the gas, in Pa, or kg/m·s²
|-
! align=right | <math>T_0</math>,
| align=left  | absolute upstream total temperature of the gas, in K
|}

The mass flow rate is primarily dependent on the cross-sectional area <math>A</math> of the nozzle throat and the upstream pressure <math>P</math>, and only weakly dependent on the temperature <math>T</math>.  The rate does not depend on the downstream pressure at all.  All other terms are constants that depend only on the composition of the material in the flow. ''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 upstream pressure is increased as this increases the density of the gas entering the orifice.

{{Confusing|date=June 2014}}
The value of <math>C_d</math> can be calculated using the below expression:
: <math>C_d = \frac{\dot{m}}{A \sqrt{2\rho\Delta P}}</math>

{| border="0" cellpadding="2"
|-
| align=right | Where:
| &nbsp;
|-
! align=right | <math>C_d</math>,
| align=left  | discharge coefficient through the constriction (dimensionless)
|-
! align=right | <math>A</math>,
| align=left  | cross-sectional area of flow constriction (unit length squared)
|-
! align=right | <math>\dot{m}</math>,
| align=left  | mass flow rate of fluid through constriction (unit mass of fluid per unit time)
|-
! align=right | <math>\rho</math>,
| align=left  | density of fluid (unit mass per unit volume)
|-
! align=right | <math>\Delta P</math>,
| align=left  | pressure drop across constriction (unit force per unit area)
|}

The above equations calculate the steady state mass flow rate for the pressure and temperature existing in the upstream pressure source.

If the gas is being released from a closed high-pressure vessel, the  above steady state equations may be used to approximate the ''initial'' mass flow rate. Subsequently, the mass flow rate will decrease during the discharge as the source vessel empties and the pressure in the vessel decreases. Calculating the flow rate versus time since the initiation of the discharge is much more complicated, but more accurate. Two equivalent methods for performing such calculations are explained and compared online.<ref>{{Cite web |url=http://www.air-dispersion.com/feature2.html |title=Calculating Accidental Release Rates From Pressurized Gas Systems |access-date=2015-01-01 |archive-url=https://web.archive.org/web/20150101190305/http://www.air-dispersion.com/feature2.html |archive-date=2015-01-01 |url-status=live }}</ref>

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 where M is the molecular weight of the gas.

==Real gas effects==

If the upstream conditions are such that the gas cannot be treated as ideal, there is no closed form equation for evaluating the choked mass flow.  Instead, the gas expansion should be calculated by reference to real gas property tables, where the expansion takes place at constant enthalpy.{{citation needed|date=January 2017}}

==Minimum pressure ratio required for choked flow to occur==

The minimum pressure ratios required for choked conditions to occur (when some typical industrial gases are flowing) are presented in Table 1. The ratios were obtained using the criterion that choked flow occurs when the ratio of the absolute upstream pressure to the absolute downstream pressure is equal to or greater than <math>(2/[\gamma + 1])^{-\gamma/(\gamma - 1)}</math>, where <math>\gamma</math> is the [[heat capacity ratio|specific heat ratio]] of the gas. The minimum pressure ratio may be understood as the ratio between the upstream pressure and the pressure at the nozzle throat when the gas is traveling at Mach 1; if the upstream pressure is too low compared to the downstream pressure, sonic flow cannot occur at the throat.

{| class="wikitable"
|+ Table 1
|-
! Gas
! <math>\gamma = \frac{c_p}{c_v}</math><ref>
{{cite book
 |author1=Perry, Robert H.
 |author2=Green, Don W.
 | title=[[Perry's Chemical Engineers' Handbook]], Table 2-166,
 | edition=6th | publisher=McGraw-Hill Company|year=1984
 | isbn=0-07-049479-7
}}</ref><ref>
{{cite book
 | author=Phillips Petroleum Company
 | title=Reference Data For Hydrocarbons And Petro-Sulfur Compounds
 | publisher=Phillips Petroleum Company
 | edition=Second Printing
 | year=1962
}}</ref>
! {{abbr|Min.|Minimum}} P<sub>u</sub>/P<sub>d</sub> <br/>for choked flow
|-
| Dry air
| 1.400 at 20&nbsp;°C || 1.893
|-
| Nitrogen
| 1.404 at 15&nbsp;°C || 1.895
|-
| Oxygen
| 1.400 at 20&nbsp;°C || 1.893
|-
| Helium
| 1.660 at 20&nbsp;°C || 2.049     
|-
| Hydrogen
| 1.410 at 20&nbsp;°C || 1.899        
|-
| Methane
| 1.307 || 1.837
|-
| Propane
| 1.131 || 1.729
|-
| Butane
| 1.096 || 1.708
|-
| Ammonia
| 1.310 at 15&nbsp;°C || 1.838
|-
| Chlorine
| 1.355 || 1.866
|-
| Sulfur dioxide
| 1.290 at 15&nbsp;°C || 1.826 
|-
| Carbon monoxide
| 1.404 || 1.895 
|-
| Carbon dioxide
| 1.30  || 1.83
|}

Notes:
* P<sub>u</sub>, absolute upstream gas pressure 
* P<sub>d</sub>, absolute downstream gas pressure

==Venturi nozzles with pressure recovery==

The flow through a venturi nozzle achieves a much lower nozzle pressure than downstream pressure.  Therefore, the pressure ratio is the comparison between the upstream and nozzle pressure.  Therefore, flow through a venturi can reach Mach 1 with a much lower upstream to downstream ratio.<ref>https://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/19680011915.pdf</ref>

==Thin-plate orifices==

The flow of real gases through thin-plate orifices never becomes fully choked.  The mass flow rate through the orifice continues to increase as the downstream pressure is lowered to a perfect vacuum, though the mass flow rate increases slowly as the downstream pressure is reduced below the critical pressure.<ref name=thin_plate_1>[http://www.engsoft.co.kr/download_e/steam_flow_e.htm Section 3 -- Choked Flow]</ref> Cunningham (1951) first drew attention to the fact that choked flow will not occur across a standard, thin, square-edged orifice.<ref>Cunningham, R.G., "Orifice Meters with Supercritical Compressible Flow" Transactions of the ASME, Vol. 73, pp. 625-638, 1951.</ref><ref>{{cite book|author=Richard W. Miller|title=Flow Measurement Engineering Handbook|edition=Third|publisher=McGraw Hill|year=1996|isbn=0-07-042366-0}}</ref>

===Vacuum conditions===
In the case of upstream air pressure at atmospheric pressure and vacuum conditions downstream of an orifice, both the air velocity and the mass flow rate becomes choked or limited when sonic velocity is reached through the orifice.

==The flow pattern ==

[[File:Flowpatterns.gif|right|thumb|upright|Figure 1. Flow patterns]]
Figure 1a shows the flow through the nozzle when it is completely subsonic (i.e. the nozzle is not choked). The flow in the chamber accelerates as it converges toward the throat, where it reaches its maximum (subsonic) speed at the throat. The flow then decelerates through the diverging section and exhausts into the ambient as a subsonic jet. Lowering the back pressure, in this state, will increase the flow speed everywhere in the nozzle.<ref name=flow_pattern>[http://www.engapplets.vt.edu/fluids/CDnozzle/cdinfo.html  The flow through the nozzle]</ref>

When the back pressure, p<sub>b</sub>, is lowered enough, the flow speed is Mach 1 at the throat, as in figure 1b. The flow pattern is exactly the same as in subsonic flow, except that the flow speed at the throat has just reached Mach 1. Flow through the nozzle is now choked since further reductions in the back pressure can't move the point of M=1 away from the throat. However, the flow pattern in the diverging section does change as you lower the back pressure further.<ref name="flow_pattern"/>

As p<sub>b</sub> is lowered below that needed to just choke the flow, a region of supersonic flow forms just downstream of the throat. Unlike in subsonic flow, the supersonic flow accelerates as it moves away from the throat. This region of supersonic acceleration is terminated by a normal shock wave. The shock wave produces a near-instantaneous deceleration of the flow to subsonic speed. This subsonic flow then decelerates through the remainder of the diverging section and exhausts as a subsonic jet. In this regime if you lower or raise the back pressure you move the shock wave away from (increase the length of supersonic flow in the diverging section before the shock wave) the throat.<ref name="flow_pattern"/>

If the p<sub>b</sub> is lowered enough, the shock wave will sit at the nozzle exit (figure 1d). Due to the very long region of acceleration (the entire nozzle length) the flow speed will reach its maximum just before the shock front. However, after the shock the flow in the jet will be subsonic.<ref name="flow_pattern"/>

Lowering the back pressure further causes the shock to bend out into the jet (figure 1e), and a complex pattern of shocks and reflections is set up in the jet which will involve a mixture of subsonic and supersonic flow, or (if the back pressure is low enough) just supersonic flow. Because the shock is no longer perpendicular to the flow near the nozzle walls, it deflects the flow inward as it leaves the exit producing an initially contracting jet. This is referred as overexpanded flow because in this case the pressure at the nozzle exit is lower than that in the ambient (the back pressure)- i.e. the flow has been expanded by the nozzle too much.<ref name="flow_pattern"/>

A further lowering of the back pressure changes and weakens the wave pattern in the jet. Eventually the back pressure will be low enough so that it is now equal to the pressure at the nozzle exit. In this case, the waves in the jet disappear altogether (figure 1f), and the jet will be uniformly supersonic. This situation, since it is often desirable, is referred to as the 'design condition'.<ref name="flow_pattern"/>

Finally, if the back pressure is lowered even further we will create a new imbalance between the exit and back pressures (exit pressure greater than back pressure), figure 1g. In this situation (called 'underexpanded') what we call expansion waves (that produce gradual turning perpendicular to the axial flow and acceleration in the jet) form at the nozzle exit, initially turning the flow at the jet edges outward in a plume and setting up a different type of complex wave pattern.<ref name="flow_pattern"/>

==See also==
*[[Accidental release source terms]] includes mass flow rate equations for non-choked gas flows as well.
*[[Orifice plate]] includes derivation of non-choked gas flow equation.
*[[de Laval nozzle]]s are Venturi tubes that produce supersonic gas velocities as the tube and the gas are first constricted and then the tube and gas are expanded beyond the choke plane.
*[[Rocket engine nozzles]] discusses how to calculate the exit velocity from nozzles used in rocket engines.
*[[Hydraulic jump]].

==References==
{{reflist}}

==External links==
*[https://web.archive.org/web/20060501124417/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]
*[http://www.enggcyclopedia.com/calculators/restriction-orifice-sizing Restriction orifice sizing control] Perform orifice plate, restriction orifice sizing calculation for a single phase flow.

[[Category:Flow regimes]]
[[Category:Aerodynamics]]
[[Category:Gas technologies]]