Fact-checked by Grok 2 weeks ago

Choked flow

Choked flow, also known as critical flow, is a phenomenon in compressible fluid dynamics where the mass flow rate through a restriction, such as a nozzle throat or valve orifice, reaches a maximum value and becomes independent of downstream pressure conditions.^[1] This limiting condition arises when the fluid velocity at the constriction accelerates to the local speed of sound, resulting in a Mach number of unity (M = 1), beyond which further reductions in downstream pressure do not increase the flow rate.^[2] The effect is prominent in gases and vapors but can also occur in two-phase flows, where compressibility plays a key role in limiting throughput.^[3] In a converging-diverging nozzle, choked flow typically establishes at the throat when the upstream-to-downstream pressure ratio falls below a critical value, approximately 0.528 for air under isentropic conditions (γ = 1.4).^[2] The choked mass flow rate can be calculated using the formula \dot{m} = A^* \cdot p_t \cdot \sqrt{\frac{\gamma}{R T_t}} \cdot \left( \frac{2}{\gamma + 1} \right)^{\frac{\gamma + 1}{2(\gamma - 1)}}, where A^* is the throat area, p_t and T_t are the stagnation pressure and temperature, \gamma is the specific heat ratio, and R is the gas constant; this expression highlights how the flow is dictated solely by upstream stagnation properties once choking occurs.^[1] Downstream of the throat, the flow may expand supersonically or form shock waves if the back pressure is not matched, but these adjustments do not affect the throat conditions or mass flow.^[2] Choked flow is fundamental in engineering applications, including rocket propulsion systems where it ensures predictable thrust in nozzles, as the exhaust achieves sonic velocity at the throat before accelerating to supersonic speeds.^[4] It also governs performance in jet engines, steam turbines, and pressure relief valves, where precise control of maximum flow rates is essential for safety and efficiency, such as in nuclear reactor cooling or industrial gas handling to prevent overpressurization.^[5] In two-phase systems, like cryogenic propellants in hybrid rockets, choking influences injector design and mass flow stability, impacting overall system reliability.

Introduction and Fundamentals

Definition and Basic Principles

Choked flow represents a critical limiting condition in fluid dynamics, occurring when the velocity of a fluid at a constriction, such as the throat of a nozzle, reaches the local speed of sound, resulting in a Mach number of 1. At this point, the mass flow rate through the restriction becomes independent of any further reductions in downstream pressure, as the flow cannot accelerate beyond sonic conditions due to the propagation limits of pressure disturbances. This phenomenon is fundamental to compressible flows, where significant density variations arise from high velocities and pressure drops.^[1]^[6] The concept of choked flow was first explored in the 19th century amid advancements in steam engines and nozzles, where engineers observed flow limitations in high-speed steam expansion. A pivotal contribution came from Swedish inventor Gustaf de Laval, who in the 1880s developed the convergent-divergent nozzle design for steam turbines, enabling efficient supersonic flow and highlighting choking effects at the throat. This innovation laid the groundwork for understanding sonic limitations in nozzles, influencing later applications in propulsion systems.^[7]^[8] Choked flow presupposes knowledge of compressible flow principles, where fluids exhibit density changes under acceleration, contrasting with incompressible approximations valid at low speeds. Key physics include isentropic expansion, an idealized reversible adiabatic process that conserves entropy during flow through a constriction. The behavior is governed by three conservation laws: the continuity equation, ensuring mass balance across varying cross-sections; the momentum equation, relating velocity changes to pressure gradients; and the energy equation, accounting for total enthalpy preservation in adiabatic conditions. These equations collectively explain how sonic conditions emerge at the minimum area without requiring detailed derivations here.^[9]^[10] For ideal gases, the local speed of sound, which defines the choking threshold, is expressed as

a = \sqrt{\gamma R T}

where \gamma is the specific heat ratio (ratio of specific heats at constant pressure and volume), R is the specific gas constant for the gas, and T is the absolute temperature. This formula underscores the temperature dependence of sonic velocity, critical for predicting when flow chokes in gaseous systems.^[11]

Occurrence in Gases and Liquids

In compressible gases, choked flow arises from the inherent compressibility effects that govern the fluid's behavior under acceleration. As the gas flows through a restriction, such as a nozzle throat, its velocity increases while density decreases due to expansion, ultimately reaching the local speed of sound (Mach number of 1) at the minimum cross-sectional area.^[1] This sonic condition limits the mass flow rate, preventing further increases even if downstream pressure is reduced, as upstream conditions cannot be influenced beyond the sonic barrier.^[12] The phenomenon requires assumptions of compressible, isentropic flow, where thermodynamic properties like pressure and temperature vary significantly along the flow path.^[1] In contrast, choked flow in liquids occurs through cavitation, a process driven by phase change rather than compressibility. When the local static pressure in the liquid drops below its vapor pressure at the prevailing temperature—typically in a constriction like a valve or pipe bend—vapor bubbles form and expand, creating a two-phase mixture that effectively blocks the flow path and restricts throughput.^[13] This rate-limiting effect is independent of achieving sonic speeds in the liquid phase alone, as liquids are nearly incompressible; instead, the vapor cavities disrupt the continuity of the flow, leading to a maximum discharge similar to gaseous choking.^[14] The prerequisite is a sufficient pressure differential to induce vaporization, often exacerbated by high velocities in hydraulic components.^[13] The key differences between these occurrences lie in their underlying physics: in gases, choking stems from thermodynamic expansion and the propagation limit imposed by the speed of sound, whereas in liquids, it involves a phase transition to vapor and the resulting two-phase flow dynamics that impede liquid motion.^[15] For instance, air flowing through a converging nozzle experiences sonic choking due to compressibility, while water in a pipe restriction undergoes cavitation choking when pressures fall below the vapor threshold, potentially causing erosion from collapsing bubbles.^[16] Experimental observations confirm cavitation's role in liquid choking within hydraulic systems, such as control valves and pumps, where flow rates plateau once the downstream-to-upstream pressure ratio drops below a critical value tied to vapor pressure (e.g., around 0.45 for certain orifice configurations under high upstream pressures).^[13] In pump inducers and centrifugal impellers, choked conditions manifest as performance breakdowns at low cavitation numbers, with vapor cavities extending along blades and calculable via vapor pressure thresholds to predict maximum flow limits.^[14]

Theoretical Foundations for Gases

Mass Flow Rate at Choked Conditions

In choked flow conditions for an ideal gas, the mass flow rate reaches its maximum value when the flow velocity at the throat of a converging nozzle equals the local speed of sound, corresponding to a Mach number of 1.^[1] This maximum rate is independent of downstream pressure as long as the pressure ratio remains below the critical value, and it serves as the quantitative limit for gas flow through restrictions.^[17] The derivation assumes isentropic flow, a perfect gas law (p = \rho R T), and constant specific heat ratio \gamma = c_p / c_v.^[12] The choked mass flow rate \dot{m} is derived from the continuity equation combined with isentropic relations and the area-velocity relation in compressible flow. Start with the mass flow rate expression: \dot{m} = \rho^* V^* A^*, where \rho^*, V^*, and A^* are the density, velocity, and cross-sectional area at the sonic throat, respectively.^[1] At sonic conditions, V^* = a^* = \sqrt{\gamma R T^*}, the speed of sound.^[17] Using the isentropic temperature relation from stagnation conditions, T^* / T_0 = 2 / (\gamma + 1), and the corresponding pressure and density ratios, p^* / p_0 = \left[2 / (\gamma + 1)\right]^{\gamma / (\gamma - 1)} and \rho^* / \rho_0 = \left[2 / (\gamma + 1)\right]^{1 / (\gamma - 1)}, substitute \rho^* = p^* / (R T^*).^[12] This yields the choked mass flow rate:

\dot{m} = A^* P_0 \sqrt{\frac{\gamma}{R T_0}} \left( \frac{2}{\gamma + 1} \right)^{\frac{\gamma + 1}{2(\gamma - 1)}}

where P_0 and T_0 are the upstream stagnation pressure and temperature.^[1] The area-velocity relation, dA / A = (M^2 - 1) dV / V, confirms that the minimum area (throat) occurs at M = 1, maximizing \dot{m} for fixed stagnation conditions.^[17] The mass flow rate depends on upstream stagnation conditions (P_0, T_0), throat geometry (A^*), the gas constant R, and \gamma. For diatomic gases like air, \gamma = 1.4 at standard conditions, yielding a numerical factor of approximately 0.6847 in the formula when expressed as \dot{m} / A^* = 0.6847 P_0 / \sqrt{R T_0}.^[18] Higher \gamma (e.g., for monatomic gases near 1.67) increases the choked rate slightly due to the functional dependence.^[18] In SI units, \dot{m} is expressed in kg/s, with P_0 in Pa, T_0 in K, R in J/(kg·K), and A^* in m². For air (R = 287 J/(kg·K), \gamma = 1.4) at standard sea-level stagnation conditions (P_0 = 101325 Pa, T_0 = 288.15 K) and a throat area A^* = 10^{-4} m² (1 cm²), the choked mass flow rate is approximately 0.0241 kg/s.^[19] This example illustrates the scale: the flow rate scales linearly with A^* and P_0, and inversely with \sqrt{T_0}, enabling practical computations for nozzle design.^[1]

Critical Pressure Ratio for Choking

The critical pressure ratio for choking in ideal gas flow is the ratio of the throat pressure to the stagnation pressure at sonic conditions, given by \frac{p^*}{p_0} = \left( \frac{2}{\gamma + 1} \right)^{\frac{\gamma}{\gamma - 1}}.^[2] For diatomic gases like air with \gamma = 1.4, this evaluates to approximately 0.528. Choking occurs when the downstream-to-upstream pressure ratio falls below this value, such that further reductions in back pressure do not increase the mass flow rate.^[1] The derivation follows from isentropic relations: at M = 1, T^* = T_0 \frac{2}{\gamma + 1}, and p^* = p_0 \left( \frac{T^*}{T_0} \right)^{\frac{\gamma}{\gamma - 1}}. This ratio depends only on \gamma; for monatomic gases (\gamma = 1.67), it is about 0.487, while for polyatomic gases (\gamma \approx 1.3), it increases to around 0.546.^[2] In a converging nozzle, subsonic flow throughout requires the back pressure ratio p_b / p_0 > p^* / p_0. When p_b / p_0 < p^* / p_0, the throat becomes sonic, establishing choked flow. This threshold is crucial for predicting flow regimes in nozzles and orifices under varying pressure conditions.^[17]

Real Gas Effects

The ideal gas model for choked flow assumes a compressibility factor Z = 1, but at high pressures or low temperatures, real gases exhibit Z \neq 1, leading to deviations in density and thermodynamic properties that alter the speed of sound and consequently the mass flow rates at choked conditions.^[20] This deviation arises from intermolecular forces, causing the actual gas density to differ from the ideal \rho = P / (RT), which impacts the sonic velocity a = \sqrt{\gamma Z R T} and the point at which flow chokes.^[20] For instance, Z < 1 at moderate temperatures near the critical point increases density relative to ideal predictions, while Z > 1 at high temperatures reduces it.^[20] To adjust for these real gas effects, equations of state such as the van der Waals or Redlich-Kwong models are applied, which account for finite molecular volume and attractive forces through parameters a and b.^[20] These require numerical solutions, often via integration of the conservation equations along the flow path, as no closed-form expression for choked mass flow exists unlike the ideal case.^[21] Alternatively, tabulated thermophysical properties from authoritative databases enable iterative calculations of enthalpy, entropy, and density during isentropic expansion to the throat.^[22] In practical examples, such as supercritical CO₂ flow in turbine or compressor nozzles, real gas effects can reduce the choke margin by up to 9% compared to ideal gas assumptions, primarily due to variations in the isentropic exponent and compressibility near the critical point.^[23] Similarly, for high-pressure steam in turbine stages, real gas modeling reveals discrepancies in predicted flow rates under extreme conditions like those in rocket propulsion analogs, highlighting the need for property corrections.^[20] These adjustments are critical in systems where ideal predictions overestimate capacity, potentially leading to operational inefficiencies.^[20] Recent advancements (post-2020) integrate real gas equations of state into through-flow and CFD simulations for high-pressure choked flows, such as in turbines handling CO₂ or natural gas mixtures, where look-up tables or databases provide accurate inputs for compressibility and transport properties to resolve numerical instabilities.^[24] These methods enable precise prediction of flow profiles without simplifying assumptions, improving design reliability in applications like carbon capture systems.^[24]

Flow Configurations and Devices

Converging-Diverging Nozzles

Converging-diverging nozzles, also known as de Laval nozzles, feature a geometry consisting of a converging section that narrows to a throat of minimum cross-sectional area, followed by a diverging section that expands the flow path. In the converging section, subsonic flow accelerates toward the throat, reaching sonic conditions (Mach number = 1) at this minimum area under choked flow. The diverging section then allows the flow to expand further, accelerating to supersonic velocities (Mach > 1) if the throat is choked and the back pressure is sufficiently low. This configuration requires the flow to be choked at the throat to achieve supersonic exit conditions, enabling efficient conversion of thermal energy to kinetic energy in high-speed applications.^[25] The choking mechanism in these nozzles occurs at the throat, where the minimum area limits the mass flow rate to a maximum value independent of downstream conditions once sonic velocity is attained. Upstream of the throat, pressure and density decrease as velocity increases, but perturbations in the diverging section do not propagate back through the sonic throat due to the speed of sound barrier, isolating the upstream flow. This allows isentropic expansion in the diverging section without altering the choked mass flow at the throat, provided the nozzle is designed for the specific pressure ratio.^[25]^[26] Key design parameters include the ratio of the exit area to the throat area (A_e / A_t), which determines the achievable exit Mach number based on the area-Mach number relation for isentropic flow:

\frac{A}{A^*} = \frac{1}{M} \left[ \frac{2}{\gamma + 1} \left(1 + \frac{\gamma - 1}{2} M^2 \right) \right]^{\frac{\gamma + 1}{2(\gamma - 1)}}

Here, A^* is the throat area, M is the local Mach number, and \gamma is the specific heat ratio of the gas. For rocket nozzles, area ratios are often selected to achieve exit Mach numbers exceeding 3, such as A_e / A_t \approx 10 for \gamma = 1.2 to produce high exhaust velocities on the order of 3-4 km/s.^[27]^[25] This nozzle design was invented by Swedish engineer Gustaf de Laval in 1888 for use in impulse steam turbines, where it enabled efficient supersonic expansion of steam to drive turbine blades. The principle later became foundational for rocketry, adapting the geometry for combustion gases to generate thrust in vacuum or atmospheric conditions.^[28]

Venturi Nozzles with Pressure Recovery

Venturi nozzles, also known as Venturi tubes, feature a converging inlet section that accelerates the flow to a minimum cross-sectional area at the throat, followed by a diverging outlet section designed to decelerate the flow and facilitate pressure recovery. In subsonic choked flow conditions, the flow reaches sonic velocity (Mach 1) at the throat, where the mass flow rate becomes independent of downstream pressure, while the diverging diffuser converts a portion of the kinetic energy back into static pressure, achieving partial recovery typically up to 80-90% of the inlet pressure. This configuration is particularly suited for precise flow metering in gases, as the smooth geometry minimizes turbulence and boundary layer separation compared to abrupt restrictions.^[29]^[30] The flow characteristics in Venturi nozzles under choked conditions differ from those in simple converging nozzles due to the pressure recovery in the diffuser, which allows choking to occur at a lower critical pressure ratio—defined as the inlet-to-outlet pressure ratio required for sonic conditions at the throat—typically ranging from 1.3 to 1.5 for air, compared to approximately 1.89 for an ideal isentropic converging nozzle without recovery. This lower ratio enables stable choked operation over a broader range of downstream pressures, with the throat pressure maintained near the critical value of about 0.528 times the inlet stagnation pressure for diatomic gases like air. In metering applications, the differential pressure across the inlet and throat is measured to infer flow rate, leveraging the compressible flow expansion factor to account for density variations.^[31] The mass flow rate through a choked Venturi nozzle is given by the adapted discharge equation:

\dot{m} = C_d A_t \sqrt{2 \rho_1 (P_1 - P_2)}

where C_d is the discharge coefficient, approximately 0.98 for well-designed Venturi nozzles with smooth contours and optimal diffuser angles (typically 7-15 degrees), A_t is the throat area, \rho_1 is the inlet density, and \Delta P = P_1 - P_2 is the pressure differential between inlet and throat. For fully choked compressible flow, this is often refined using the isentropic choked flow relation multiplied by C_d, ensuring accuracy within 0.5-1% for Reynolds numbers above 10^5. The high C_d value reflects the nozzle's efficiency in achieving near-ideal flow attachment.^[32]^[29] A key advantage of Venturi nozzles over thin-plate orifices in choked flow metering is the minimal permanent pressure loss, often limited to 10-20% of the total differential pressure due to effective recovery in the diffuser, which reduces energy dissipation and operational costs in continuous flow systems such as gas pipelines or engine testing. This recovery also enhances measurement repeatability, with discharge coefficients remaining stable across Mach numbers from 0.2 to 1.0 and varying minimally with Reynolds number in turbulent regimes.^[29]

Thin-Plate Orifices

Thin-plate orifices are simple flow restriction devices consisting of a flat plate, typically with a thickness-to-diameter ratio of 0 to 0.015, featuring a central sharp-edged hole inserted perpendicularly into a conduit. Upon passing through the orifice, the fluid streamlines converge and contract downstream, forming a vena contracta—the narrowest cross-section of the jet where velocity is maximum and pressure is minimum—due to the abrupt geometry that prevents immediate expansion. For a circular hole in such a thin plate, the coefficient of contraction C_c, defined as the ratio of the vena contracta area to the orifice area, is approximately 0.61 at high Reynolds numbers, reflecting the inertial dominance over viscous effects in the flow.^[33]^[34] Unlike converging-diverging nozzles, choked flow in thin-plate orifices is rarely fully realized because the free jet downstream allows expansion waves to propagate, preventing the sonic condition from being strictly confined to the orifice plane. As downstream pressure decreases below the critical ratio, the vena contracta shifts upstream toward the orifice but does not coincide with it, enabling a slight increase in mass flow rate even into near-vacuum conditions, rather than plateauing at a maximum. This partial choking behavior, first experimentally demonstrated for compressible flows, results in flow rates that approach but do not attain the theoretical choked limits of isentropic nozzle flow.^[33] The mass flow rate \dot{m} through a thin-plate orifice under non-choked or near-choked conditions is commonly estimated using the formula

\dot{m} = C_d A \sqrt{2 \rho (P_1 - P_2)},

where C_d \approx 0.6 is the discharge coefficient (product of C_c and the velocity coefficient C_v \approx 0.98, accounting for frictional losses), A is the geometric orifice area, \rho is the upstream fluid density, and P_1 and P_2 are the upstream and downstream pressures, respectively. This expression, derived from Bernoulli's principle with empirical corrections, provides a practical approximation for both liquids and gases when the pressure differential \Delta P = P_1 - P_2 is significant but subcritical, though compressible effects require adjustments for gases near choking.^[33] For accurate flow measurement using thin-plate orifice meters, calibration follows international standards such as ISO 5167-2, which specifies geometry, installation (e.g., flange or corner tappings), and Reynolds number ranges (> 5,000 for reliability) for differential pressure devices in subsonic, turbulent flows with P_2 / P_1 \geq 0.75. In high \Delta P scenarios approaching choked limits, these standards guide uncertainty estimates (typically ±1-2% for discharge coefficients), but extrapolations beyond subsonic applicability demand experimental validation or advanced modeling to account for incomplete choking.^[33]

Applications and Special Cases

Vacuum and Low-Pressure Conditions

In vacuum systems, choked flow plays a key role during the initial stages of evacuation, where gas accelerates to sonic velocity through restrictions such as valves or orifices when the upstream-to-downstream pressure ratio exceeds the critical value of approximately 0.528 for air at 20°C. This limits the mass flow rate, influencing pump-down times in rough vacuum conditions (1000–1 mbar). For instance, in leak detection, choked flow through calibrated orifices enables precise gas flow measurements using evacuated chambers as reference standards. In diffusion pumps, maintaining the backing pressure below the critical forepressure (typically around 0.1 mbar) is essential to avoid disrupting the vapor jet, ensuring reliable compression from ultrahigh vacuum levels (10^{-9} mbar) to higher pressures. These considerations are vital for applications in semiconductor processing and space simulation chambers, where precise control of low-pressure gas dynamics prevents inefficiencies or contamination.^[35]

Industrial and Aerospace Applications

In industrial settings, choked flow is critical for flow control in chemical plants, where control valves manage high-pressure gases and liquids to prevent overpressurization and ensure process stability. For instance, in processes involving compressible fluids, choked conditions limit the maximum mass flow rate through valves, requiring precise sizing to avoid cavitation-induced erosion and excessive noise exceeding 100 dB.^[15] In steam turbines, choked flow occurs at the exit of the last blade stage when the back pressure results in sonic velocity, optimizing efficiency but necessitating careful design to handle the transition to supersonic flow in the exhaust.^[36] Cryogenic valves in chemical processing plants, such as those handling liquefied natural gas or oxygen, often operate under choked conditions to regulate low-temperature fluids, where rapid pressure drops can induce flashing and two-phase flow.^[37] Additionally, safety relief valves in these facilities are designed per API Standard 520 to accommodate choked flow during overpressure events, ensuring the relieving pressure does not exceed the critical pressure ratio for gases, typically calculated using the omega method for real fluids to determine the effective discharge area.^[38] In aerospace applications, choked flow is fundamental to rocket engines, where the throat of converging-diverging nozzles achieves sonic conditions to maximize mass flow and thrust, as seen in liquid rocket engines like those using RP-1 and LOX propellants.^[39] The SpaceX Merlin engine exemplifies this, with its choked throat enabling high-pressure combustion chamber operation up to 9.7 MPa, facilitating reliable ignition and sustained performance during ascent.^[40] Jet engine exhausts rely on choked flow in the nozzle throat to accelerate exhaust gases to supersonic speeds, contributing to thrust generation in afterburning turbofans.^[41] Wind tunnels simulating aerospace conditions often induce choked flow through variable throat nozzles to achieve Mach numbers up to 5, allowing precise replication of flight aerodynamics.^[42] Furthermore, choked flow enables thrust vectoring in fluidic systems, where secondary air injection shifts the effective throat or creates oblique shocks in the divergent section, achieving vector angles up to 18° with minimal efficiency loss (thrust coefficient 0.86–0.98).^[43] Recent developments in the 2020s have integrated additive manufacturing for nozzles in hypersonic vehicles, producing complex geometries that maintain choked flow at throats while enhancing thermal management for speeds exceeding Mach 5.^[44] In CO2 sequestration pipelines, choked flow models predict temperature drops to -75°C during depressurization through valves, informing designs to prevent solid CO2 formation and ensure safe transport in carbon capture systems.^[45] Challenges in applying choked flow include scaling models for real gases, where deviations from ideal behavior require compressibility factors (Z) and specific heat ratios (k) in sizing equations, as outlined in API 520 for non-ideal fluids.^[38] Erosion at throats from high-velocity particles or cavitation poses significant risks, with computational studies showing maximum rates increasing with inlet pressure and particle mass flow.^[46] Mitigation strategies employ materials like Inconel 718 or 625 for nozzle throats, offering superior oxidation and creep resistance to withstand thermal and erosive stresses in rocket applications.^[47]

Flow Characteristics

Velocity and Mach Number Profiles

In choked flow through a converging section, such as in a de Laval nozzle, the subsonic flow accelerates progressively toward the throat, with velocity increasing as the cross-sectional area decreases, consistent with the area-velocity relation for compressible fluids.^[25] This acceleration raises the Mach number from near-zero at the inlet to exactly 1 at the throat, where sonic conditions are achieved and the flow chokes, limiting the mass flow rate.^[48] Downstream of the throat in a diverging section under choked conditions, the flow transitions to supersonic speeds and continues to accelerate if the expansion is isentropic, with the Mach number increasing as the area expands.^[25] However, in overexpanded nozzles where the exit pressure is lower than the ambient pressure, oblique or normal shocks may form within the divergent section, abruptly reducing the Mach number and velocity to match external conditions.^[49] The Mach-area relation governs these profiles in isentropic choked flow, plotting the area ratio A/A^* (where A^* is the throat area) against Mach number M. At the minimum area (A/A^* = 1), M = 1; for A/A^* > 1, the curve exhibits two branches—subsonic (M < 1, decreasing toward the throat) and supersonic (M > 1, increasing away from the throat)—symmetric in shape but representing distinct flow regimes.^[48] Schlieren imaging visualizes these profiles by capturing density gradients, revealing the sonic line as a sharp boundary at the throat and expansion fans or shock waves in supersonic regions downstream, as seen in experimental studies of nozzle plumes.^[50]

Pressure and Density Distributions

In isentropic choked flow through a converging-diverging nozzle, the static pressure decreases from the stagnation value P_0 upstream of the throat to the critical pressure P^* at the sonic throat, where P^*/P_0 = \left[2/(\gamma + 1)\right]^{\gamma/(\gamma - 1)}. For air with \gamma = 1.4, this ratio is approximately 0.528, marking the point where the flow reaches Mach 1 and becomes choked, limiting further mass flow increases despite lower downstream pressures.^[9] Downstream of the throat in the diverging section, the pressure continues to decrease along the isentropic relation P/P_0 = \left[1 + (\gamma - 1)/2 \cdot M^2 \right]^{-\gamma/(\gamma - 1)} for supersonic expansion, reaching the exit pressure P_e if matched to the back pressure; however, if the nozzle is underexpanded (P_e > P_b), expansion fans form outside, while in overexpanded conditions (P_e < P_b), oblique or normal shocks adjust the flow.^[9] In cases with pressure recovery, such as subsonic flow or post-shock deceleration, the pressure rises toward the back pressure P_b, but this is absent in fully supersonic isentropic choked flow.^[51] The density distribution follows the isentropic relation \rho / \rho_0 = (P / P_0)^{1/\gamma}, resulting in a sharp drop from stagnation density \rho_0 to critical density \rho^* at the throat. For air (\gamma = 1.4), \rho^*/\rho_0 \approx 0.634, reflecting the sonic condition where density decreases as velocity increases to sonic speed.^[9] More generally, \rho / \rho_0 = \left[1 + (\gamma - 1)/2 \cdot M^2 \right]^{-1/(\gamma - 1)}, so density continues to fall in the supersonic diverging section until the exit, with the drop most pronounced near the throat due to the accelerating flow.^[9] This density reduction supports the choked mass flow rate \dot{m} = A^* P_0 \sqrt{\gamma / (R T_0)} \cdot \left[2/(\gamma + 1)\right]^{(\gamma + 1)/(2(\gamma - 1))}, but spatial variations emphasize the thermodynamic compression upstream and expansion downstream.^[52] In overexpanded choked flows, where the back pressure exceeds the design exit pressure, a normal shock may form in the diverging section, causing an abrupt pressure rise across the shock wave while ending the supersonic region. The pressure ratio across the normal shock is P_2 / P_1 = \left[2 \gamma M_1^2 - (\gamma - 1)\right] / (\gamma + 1), where M_1 is the upstream Mach number; for example, at M_1 = 2 and \gamma = 1.4, this yields P_2 / P_1 \approx 4.5, rapidly increasing pressure to match the higher back pressure and transitioning the flow to subsonic downstream.^[53] This shock-induced pressure jump also causes a corresponding density increase post-shock, as \rho_2 / \rho_1 = \left[(\gamma + 1) M_1^2\right] / \left[(\gamma - 1) M_1^2 + 2\right], approximately 2.67 for the M_1 = 2 case, though total pressure losses occur due to the non-isentropic nature of the shock.^[53] Such shocks position based on back pressure, moving downstream as P_b decreases until exiting the nozzle for perfectly matched conditions.^[51] For air in a choked converging-diverging nozzle, pressure-temperature-entropy (P-T-s) diagrams illustrate these distributions: starting from stagnation conditions (P_0, T_0, s_0), the isentropic path traces a vertical line on the T-s plane from T_0 to T^* = T_0 \cdot 2/(\gamma + 1) \approx 0.833 T_0 at the throat (constant entropy), with pressure dropping to $0.528 P_0.^[9] Further expansion to the exit follows another isentropic segment, lowering T and P while density falls to \rho_e / \rho_0 = (P_e / P_0)^{1/\gamma}; if a normal shock occurs, the path deviates rightward on the T-s diagram (entropy increase), with post-shock pressure recovery shown as a subsonic compression line back toward P_b.^[54] These diagrams highlight the thermodynamic efficiency loss from shocks, with numerical examples for air confirming the throat values as benchmarks for design.^[51]

References

[1]
Mass Flow Choking
The maximum mass flow rate through the system occurs when the flow is choked at the smallest area. This location is called the throat of the nozzle.
[2]
[PDF] Choked Flow
Thus the flow is said to be choked in the sense that it does not matter what the downstream conditions are as long as they are such that the Mach number in the ...Missing: definition | Show results with:definition
[3]
[PDF] THE RELAP5 CHOKED FLOW MODEL AND APPLICATION TO A ...
Choked flow is defined as the conditions wherein the mass flow rate becomes independent of the downstream conditions, that is, that point at which further ...
[4]
Rocket Thrust Equations
The smallest cross-sectional area of the nozzle is called the throat of the nozzle. The hot exhaust flow is choked at the throat, which means that the Mach ...Missing: valves | Show results with:valves
[5]
[PDF] Module 13: Critical Flow Phenomenon Joseph S. Miller, PE
• Critical flow (“choked flow”) involves situation where fluid moves from higher pressure volume at speed limited only by speed of sound for fluid – such as ...
[6]
Sonic Flow - an overview | ScienceDirect Topics
When the fluid flow velocity in a choke reaches the traveling velocity of sound in the fluid under the in situ condition, the flow is called sonic flow. Under ...
[7]
de Laval nozzle – Knowledge and References - Taylor & Francis
The convergent–divergent nozzle design is also known as the de Laval nozzle after the Swedish inventor Gustaf de Laval who developed it in the 19th century. As ...
[8]
Converging Diverging Nozzle
This device was invented by Carl de Laval toward the end of the l9th century and is thus often referred to as the 'de Laval' nozzle. This applet is intended to ...
[9]
Isentropic Flow Equations
On this slide we have collected many of the important equations which describe an isentropic flow. We begin with the definition of the Mach number.
[10]
[PDF] Compressible Flow
BASIC EQUATIONS OF COMPRESSIBLE FLUID FLOW. The basic equations of compressible fluid flow are : (i) Continuity equation, (ii) Momentum equation, (iii) ...
[11]
Speed of Sound - HyperPhysics
R = the universal gas constant = 8.314 J/mol K, · T = the absolute temperature · M = the molecular weight of the gas in kg/mol · γ = the adiabatic constant, ...Missing: γRT | Show results with:γRT
[12]
[PDF] Notes on Thermodynamics, Fluid Mechanics, and Gas Dynamics
Dec 15, 2021 · Because the flow is choked, the throat (exit) pressure will be the sonic pressure: (3). (4). The mass flow rate will be the choked flow mass ...
[13]
[PDF] Characterization of high-pressure cavitating flow through a thick ...
If the local static pressure becomes lower than the vapor pressure corresponding to the fluid temperature, the liquid phase will be vaporized, thus inducing ...
[14]
Chapter 7 - Hydrodynamics of Pumps - Christopher E. Brennen
Changes in the temperature of the liquid being pumped will clearly affect the vapor pressure, pV, and therefore the NPSH or cavitation number. This effect ...
[15]
Choked Flow in Control Valves
In liquid applications, choking is a result of the reduction in pressure through the control elements. Figure 2 shows the instantaneous pressure as liquid moves ...
[16]
Critical Flow vs. Choked Flow: Simple Explanation for Valve Sizing
Since gases are compressible, they can reach sonic speeds at a restriction. Liquids don't normally do this, so true choked flow doesn't occur. Instead, liquids ...
[17]
Mass Flow Rate Equations | Glenn Research Center - NASA
Nov 21, 2023 · This equation is shown in the red box on this slide and relates the mass flow rate to the flow area A, total pressure pt and temperature Tt of ...Missing: ideal | Show results with:ideal
[18]
Specific Heat - Cp & Cv | Glenn Research Center - NASA
Jun 30, 2025 · For air, gamma = 1.4 for standard day conditions. “Gamma” appears in many fluids equations including the equation relating pressure, temperature ...
[19]
[PDF] The Ideal Gas Constant
Jan 6, 2014 · • The specific gas constant (R, sometimes Rgas) is not universal ... • For air in S.I. units, air. kJ. 8.3143. kJ. J kmol K. 0.2870. 287.0 kg.
[20]
[PDF] Compressible flow at high pressure with linear equation of state
This paper analyzes compressible flow at high pressures, where real-gas behavior differs from ideal-gas, using a linearization of the equation of state.
[21]
[PDF] Flow Characteristics of High-Pressure Hydrogen Gas in the Critical ...
In the present study, Redlich-Kwong's equation of state is selected to simulate the high-pressure hydrogen gas through a critical nozzle, since it bears a ...
[22]
[PDF] COMPARISON OF FIVE NATURAL GAS EQUATIONS OF STATE ...
addition to the compressibility factor. However, only the results for the ... [15] Johnson, R. C., Calculations of Real-Gas. Effects in Flow Through ...
[23]
[PDF] An Investigation of Real Gas Effects in Supercritical CO2 Centrifugal ...
Feb 25, 2015 · A reduction of up to 9% in the choke margin and an increase of up to 5% in the swing between inlet and outlet corrected mass flow per unit area ...
[24]
Method for considering real gas equations of state in through-flow ...
Jul 1, 2024 · The jump in the compressibility factor from 0.60 to 1.06 further verifies the robustness of the calculation for a strong change in the ...
[25]
Converging/Diverging (CD) Nozzle
The flow in the throat is sonic which means the Mach number is equal to one in the throat. Downstream of the throat, the geometry diverges and the flow is ...
[26]
[PDF] Converging-Diverging Nozzle
shows that M = f(A/A*), i.e., the Mach number at any location in the duct is a function of the ratio of the local duct area to the sonic throat area. • Note ...Missing: de parameters
[27]
[PDF] NOZZLES
Notice that, to get the same high Mach number, e.g. M=3, the area ratio needed is A*/A=0.33 for γ=1.67 and A*/A=0.15 for γ=1.20, i.e. more than double exit area ...
[28]
[PDF] Analysis of the Ping-Pong Ball Gun-Theoretical & Computational ...
Jul 3, 2014 · This converging- diverging nozzle was invented by 'Gustaf de Laval' in 1888 for use in steam turbines. The De Laval nozzle relies on the ...
[29]
[PDF] PERFORMANCE OF A VENTURI METER WITH SEPA~RABLE ...
In situations requiring low pressure loss, the venturi meter is one of the most com- monly used devices for the measurement of flow rate. Considerable work has ...
[30]
Sonic Nozzles - FlowMaxx Engineering
The Sonic Nozzle is operated by either pressurizing the inlet (P1) or evacuating the exit (P3), to achieve a pressure ratio of 1.2 to 1 or greater, inlet to ...<|separator|>
[31]
CEESI Critical Flow Venturi/Sonic Nozzle Publications
The choking pressure ratio (CPR) for a Critical Flow Venturi (CFV) is defined as the ratio of the maximum permissible exit pressure to the inlet pressure ...Missing: tube | Show results with:tube
[32]
Discharge Coefficients of Critical Venturi Nozzles for CO2 and SF6
The discharge coefficients of critical Venturi nozzles were measured for CO2 and SF6 on the Reynolds number range from 3×103 to 2×105. The results showed.
[33]
[PDF] Modeling Choked Flow Through an Orifice
Jun 23, 2010 · AFT Standard air is used as the fluid with the orifice upstream pressure constant at 50 psia @ 70F and downstream pressure varying from 45.
[34]
Flow Through an Orifice - Richard Fitzpatrick
For instance, if the orifice is a circular hole punched in a thin plate then the contraction coefficient is observed to take the value $ 0.62$ (Batchelor ...
[35]
[PDF] Calculation of the Choked Back Pressure for Steam Turbines with ...
The turbine choked back pressure refers to the back pressure under which the steam flow speed at the exit of the turbine last blade is close to the local sound ...
[36]
Computational study of compressible flow through choke valve
This study investigates compressible flow through a choke valve, quantifying vorticity, Mach number, pressure, temperature, and density, and observing vortex ...Missing: plants | Show results with:plants
[37]
[PDF] PRESSURE RELIEF VALVE ENGINEERING HANDBOOK | Emerson
... pressure relief valve orifice area and maximum available flow. This sizing program is a powerful tool, yet easy to use. its many features include quick and ...
[38]
Rocket Engines – Introduction to Aerospace Flight Vehicles
... critical in achieving optimal thrust performance from the rocket engine. ... A throat where the flow becomes sonic or choked; 3. A diverging or “bell ...
[39]
Review of Thrust Regulation and System Control Methods of ... - MDPI
Liquid rocket engines, as critical power units in aerospace propulsion ... In cavitating-choked flow, mechanical energy is primarily consumed to enlarge ...
[40]
Wind Tunnel Installation Effects on a High-Speed Exhaust Flow ...
May 19, 2025 · A jet vectoring effect is identified due to the pylon, which is intensified under choked tunnel operation. Additionally, a methodology is ...
[41]
[PDF] An Evaluation of Wind Tunnel Test Techniques for Aircraft Nozzle ...
thrust vectoring will require a more extensive jet-effects testing effort to determine system performance. For several years AEDC research and technology ...
[42]
[PDF] SUMMARY OF FLUIDIC THRUST VECTORING RESEARCH ...
The desire for such integrated exhaust nozzles was the catalyst for new fluidic control techniques; including throat area control, expansion control, and thrust ...
[43]
3D printing technologies for creating models and nozzles in an ...
Aug 6, 2025 · Modern additive technologies allow the production of models. The basic properties of the materials and the advantages of 3D printing methods are ...
[44]
An experimental study on the choked flow characteristics of CO2 ...
Zhao et al. [14] developed a choked flow model based on the adiabatic process assumption to calculate the temperature drop of pressurized CO2 pipeline ...
[45]
(PDF) CFD Simulations of Flow Jetting Impact and High Erosion ...
In this study, computational fluid dynamics (CFD) is used to investigate erosion on a production choke valve and its downstream spool. There are three main ...
[46]
Comparative Analysis of Inconel 718 and Inconel 625 for Rocket ...
Aug 10, 2025 · Inconel 718 and Inconel 625 are widely used nickel-based superalloys due to their high temperature strength, oxidations resistance, and ...
[47]
[PDF] Isentropic Flow with Area Change Mach Number Relation
Area Ratio. • For isentropic flow, look at effect of area change on M by comparing A at any point to area at sonic point (A*).
[48]
[PDF] Dimensional, Convergent= Divergent Exhaust Nozzle
When a nozzle is overexpanded (Po> Ps,,), the flow often separates or shocks form that cause nonuniformities in exit flow velocity and pressure profiles.
[49]
[PDF] Exhaust Nozzle Plume Effects on Sonic Boom Test Results for ...
Schlieren images are shown in Figure 16, and also demonstrate the same flow features as presented in CFD contour plots of Mach number. A description of the ...
[50]
[PDF] Notes on Thermodynamics, Fluid Mechanics, and Gas Dynamics
Dec 15, 2021 · To determine the exit plane pressure and initial mass flow rate through the nozzle, first determine whether or not the flow is choked. Determine ...
[51]
Critical Flow - Thermopedia
Sep 7, 2010 · Also, sharp-edged inlet geometries generate vena contracta, reducing the area of the choked plane and hence mass flow rates, whereas well- ...
[52]
Normal Shock Wave Equations
A normal shock occurs in front of a supersonic object if the flow is turned by a large amount and the shock cannot remain attached to the body. The detached ...
[53]
[PDF] AME 436 - Paul D. Ronney
... Flow. Isentropic nozzle flow. ➢ When can choking occur? If M ≥ 1 or so need ... T-s diagram for the various diabatic flows. (reference state Mach = 1).<|control11|><|separator|>