4.1 Compressible Gas Flow, Sonic Velocity, and Choked Flow

4.1 Compressible Gas Flow, Sonic Velocity, and Choked Flow

In chemical manufacturing plants, gases frequently flow under conditions involving large pressure gradients, extreme velocities, or rapid depressurization. Examples include emergency relief valves (PRVs), flare headers, restriction orifices, gas pipeline ruptures, high-pressure steam turbines, and vacuum ejectors. When density changes across a piping element become significant, the constant-density assumptions of incompressible fluid mechanics fail completely. Mastering compressible gas dynamics is vital for the NCEES PE Chemical exam, where questions routinely evaluate sonic speed, choking thresholds, relief valve capacities, and stagnation conditions.

---

1. Fundamentals of Compressible Flow and Compressibility Criteria

A fluid flow is classified as compressible when variations in fluid density $\rho$ significantly influence the momentum and energy balances throughout the flow field. For gases, density is coupled to local static pressure $P$ and absolute temperature $T$ via the thermodynamic equation of state:

$$\rho = \frac{P M}{Z R T}$$

where $M$ is molecular weight, $R = 8,314.5\text{ J/(kmol}\cdot\text{K)} = 10.731\text{ psia}\cdot\text{ft}^3/(\text{lbmol}\cdot^\circ\text{R})$ is the universal gas constant, and $Z$ is the compressibility factor.

The Engineering Threshold for Incompressible vs. Compressible Modeling

Two practical diagnostic criteria dictate when compressible flow equations must be deployed on the PE exam:

1. The Mach Number Criterion: If the fluid velocity $v$ exceeds $30\%$ of the local speed of sound ($Ma = v/c > 0.30$), inertial compression effects cause density changes exceeding $5\%$. Below $Ma = 0.30$, gases can be modeled as incompressible fluids with less than $1\%$ error.
2. The Fractional Pressure Drop Criterion: For flow through a pipeline, valve, or restriction, if the overall pressure drop exceeds $10\%$ of the absolute upstream pressure: $$\frac{\Delta P}{P_1} = \frac{P_1 - P_2}{P_1} > 0.10$$ the fluid density variation is too large for standard incompressible Darcy-Weisbach formulations. If $\Delta P / P_1 \le 0.10$, an engineer may use incompressible formulas evaluated at the arithmetic mean density $\rho_{avg} = (\rho_1 + \rho_2)/2$.

---

2. Speed of Sound (Sonic Velocity) and the Mach Number

The speed of sound ($c$) is the velocity at which an infinitesimal acoustic pressure wave propagates through a fluid medium. Because acoustic wave propagation occurs with negligible heat transfer and minimal viscous dissipation, the disturbance travels as an isentropic process (reversible and adiabatic):

$$c = \sqrt{\left(\frac{\partial P}{\partial \rho}\right)_s}$$

For an ideal gas obeying the isentropic relation $P \rho^{-k} = \text{constant}$ (where $k = C_p / C_v$ is the heat capacity ratio / isentropic expansion exponent):

$$\left(\frac{\partial P}{\partial \rho}\right)_s = \frac{k P}{\rho} = \frac{k R T}{M}$$

Yielding the fundamental sonic velocity formula:

$$c = \sqrt{\frac{k R T}{M}} = \sqrt{k R_{sp} T}$$

where $R_{sp} = R / M$ is the specific gas constant (in $\text{J/(kg}\cdot\text{K)}$ or $\text{ft}\cdot\text{lb}_f/(\text{lb}_m\cdot^\circ\text{R})$).

In US Customary units, apply the gravitational conversion constant $g_c = 32.174\text{ lb}_m\cdot\text{ft}/(\text{lb}_f\cdot\text{s}^2)$:

$$c = \sqrt{\frac{k g_c R T}{M}} = \sqrt{k g_c R_{sp} T}$$

Flow Classification by Mach Number

The Mach number ($Ma$) is the dimensionless ratio of bulk flow velocity $v$ to the local acoustic velocity $c$:

$$Ma = \frac{v}{c}$$

• Subsonic Flow ($Ma < 1.0$): Pressure disturbances travel faster than the fluid stream. Upstream fluid elements sense downstream flow restrictions (such as valves or bends) and adjust smoothly prior to arrival.
• Sonic Flow ($Ma = 1.0$): Fluid velocity exactly matches the speed of sound. Small pressure disturbances cannot propagate upstream against the flow.
• Supersonic Flow ($Ma > 1.0$): Fluid velocity exceeds the speed of sound. Upstream fluid cannot receive acoustic warning of downstream obstacles, leading to shock waves.

---

3. Stagnation State and Isentropic Flow Relations

The stagnation state (denoted by subscript $0$) is the thermodynamic state a fluid would achieve if brought to rest reversibly and adiabatically ($v \to 0$). For a flowing ideal gas with constant specific heats, the energy balance equates enthalpy and kinetic energy:

$$h_0 = h + \frac{v^2}{2} \implies C_p T_0 = C_p T + \frac{v^2}{2}$$

Dividing by $C_p T$ and substituting $C_p = \frac{k R}{M (k - 1)}$ and $c^2 = \frac{k R T}{M}$ yields the stagnation temperature ratio:

$$\frac{T_0}{T} = 1 + \left(\frac{k - 1}{2}\right) Ma^2$$

For an isentropic expansion, combining the temperature ratio with the polytropic isentropic path $P / T^{\frac{k}{k-1}} = \text{constant}$ yields the stagnation pressure ratio and stagnation density ratio:

$$\frac{P_0}{P} = \left[1 + \left(\frac{k - 1}{2}\right) Ma^2\right]^{\frac{k}{k-1}}$$ $$\frac{\rho_0}{\rho} = \left[1 + \left(\frac{k - 1}{2}\right) Ma^2\right]^{\frac{1}{k-1}}$$

---

4. One-Dimensional Variable-Area Duct Dynamics

Combining the differential continuity equation ($d\rho/\rho + dA/A + dv/v = 0$), differential Euler momentum equation ($dP = -\rho v dv$), and sonic definition ($c^2 = dP/d\rho$) yields the area-velocity relationship (Hugoniot equation):

$$\frac{dA}{A} = (Ma^2 - 1) \frac{dv}{v}$$

This fundamental equation reveals why compressible gas flow behaves counter-intuitively compared to liquid flow:

1. Subsonic Regime ($Ma < 1.0$): Since $Ma^2 - 1 < 0$, a decrease in area ($dA < 0$) causes an increase in velocity ($dv > 0$). Fluid accelerates through a converging nozzle, exactly like an incompressible liquid.
2. Supersonic Regime ($Ma > 1.0$): Since $Ma^2 - 1 > 0$, an increase in area ($dA > 0$) causes an increase in velocity ($dv > 0$). Fluid accelerates through a diverging duct.
3. Sonic Condition ($Ma = 1.0$): At $Ma = 1$, $dA/A = 0$, meaning sonic velocity can only occur at a local minimum cross-sectional area—the throat ($dA = 0$).

---

5. Choked Flow and the Critical Pressure Ratio

When a gas expands from a high-pressure reservoir ($P_0, T_0$) through an orifice, nozzle, control valve, or relief valve into a downstream region of lower pressure ($P_b$), the flow rate increases as backpressure $P_b$ is reduced. However, once the throat velocity reaches the local speed of sound ($Ma_t = 1.0$), the system reaches choked flow (or critical flow).

Evaluating the stagnation pressure formula at the sonic throat condition ($Ma = 1.0$) yields the critical pressure ratio ($r_c = P^* / P_0$):

$$\frac{P^*}{P_0} = \left(\frac{2}{k + 1}\right)^{\frac{k}{k - 1}}$$

where $P^*$ is the critical throat pressure (the static pressure at the throat when $Ma = 1$).

Typical Values of the Critical Pressure Ratio

• Diatomic gases and air ($k = 1.40$): $$\frac{P^*}{P_0} = \left(\frac{2}{1.4 + 1}\right)^{\frac{1.4}{0.4}} = (0.8333)^{3.5} = 0.52828 \approx 0.528$$
• Triatomic gases, methane, and steam ($k = 1.30$ to $1.31$): $$\frac{P^*}{P_0} = \left(\frac{2}{1.31 + 1}\right)^{\frac{1.31}{0.31}} = (0.8658)^{4.226} = 0.5439 \approx 0.544$$
• Monatomic gases ($k = 1.667$): $$\frac{P^*}{P_0} = \left(\frac{2}{1.667 + 1}\right)^{\frac{1.667}{0.667}} = (0.750)^{2.5} = 0.4871 \approx 0.487$$

The Choking Criterion

• Unchoked (Subcritical) Flow ($P_b / P_0 > P^* / P_0$): Throat pressure equals the receiver backpressure ($P_t = P_b$). The throat Mach number is subsonic ($Ma_t < 1$), and mass flow rate depends on both $P_0$ and $P_b$.
• Choked (Critical) Flow ($P_b / P_0 \le P^* / P_0$): Throat pressure is locked at the critical pressure ($P_t = P^*$). The throat Mach number is identically $1.0$. Acoustic pressure disturbances in the receiver propagate at speed $c$, which cannot travel upstream against the sonic throat velocity. Consequently, lowering the downstream backpressure further has zero effect on the throat pressure, velocity, or mass flow rate.

---

6. Choked Mass Flow Rate Formulations and Relief Sizing

For an ideal gas undergoing isentropic choked flow through a restriction of throat area $A_t$ with discharge coefficient $C_d$, the maximum mass flow rate is governed by:

$$\dot{m}_{choked} = C_d A_t \rho^* c^*$$

Substituting the critical density $\rho^* = \rho_0 \left(\frac{2}{k+1}\right)^{\frac{1}{k-1}}$ and critical sonic velocity $c^* = \sqrt{\frac{k R T_0}{M} \left(\frac{2}{k+1}\right)}$ leads to the classic choked mass flow equation:

$$\dot{m}_{choked} = C_d A_t P_0 \sqrt{\frac{k M}{R T_0}} \left(\frac{2}{k + 1}\right)^{\frac{k + 1}{2(k - 1)}}$$

In US Customary units ($P_0$ in $\text{lb}_f/\text{ft}^2$, $A_t$ in $\text{ft}^2$, $T_0$ in $^\circ\text{R}$, $g_c = 32.174$, $R = 1,545.3\text{ ft}\cdot\text{lb}_f/(\text{lbmol}\cdot^\circ\text{R})$):

$$\dot{m}_{choked} = C_d A_t P_0 \sqrt{\frac{k M g_c}{R T_0}} \left(\frac{2}{k + 1}\right)^{\frac{k + 1}{2(k - 1)}}$$

API Standard 520 Formulation for Safety Relief Valves (PRVs)

In the chemical process industries and on the NCEES PE Chemical exam, vapor relief valve sizing adheres to API 520 Part I. For critical vapor flow, the required effective discharge area $A$ is expressed as:

$$W = C K_d P_1 A K_w K_c \sqrt{\frac{M}{T Z}}$$

where:

• $W$ is mass flow rate ($\text{lb/h}$),
• $C$ is the gas expansion coefficient: $C = 520 \sqrt{k \left(\frac{2}{k+1}\right)^{\frac{k+1}{k-1}}}$,
• $K_d$ is the effective discharge coefficient (typically $0.975$ for certified gas nozzles),
• $P_1$ is upstream relieving pressure in absolute units ($\text{psia}$), equal to set pressure plus allowable overpressure plus atmospheric pressure,
• $K_w$ is the backpressure correction factor ($K_w = 1.0$ for conventional PRVs discharging to atmospheric or balanced bellows valves with backpressure $< 30\%$),
• $K_c$ is the rupture disk combination factor ($1.0$ without disk, $0.9$ if disk installed upstream),
• $T$ is relieving temperature ($^\circ\text{R} = ^\circ\text{F} + 459.67$).

---

7. Fanno Flow vs. Isothermal Flow in Process Piping

When compressible gas flows through a long conduit of constant cross-sectional area, friction drives changes in fluid properties:

Fanno Flow (Adiabatic Pipe Flow with Wall Friction)

In insulated process piping or short relief header lines, heat transfer is negligible ($q = 0$), so stagnation temperature remains constant ($T_0 = \text{constant}$):

• Friction accelerates a subsonic gas flow ($dMa/dx > 0$). As velocity increases, static pressure $P$ and static temperature $T$ drop.
• The maximum possible length of straight pipe before choking occurs at the exit ($Ma = 1.0$) is the limiting length ($L^*$): $$\frac{f_D L^*}{D} = \frac{1 - Ma^2}{k Ma^2} + \left(\frac{k + 1}{2k}\right) \ln\left[ \frac{(k + 1) Ma^2}{2 + (k - 1) Ma^2} \right]$$
• If the actual pipe length exceeds $L^*$, the upstream pressure $P_0$ must increase or the inlet mass flow rate will decrease (choking at the pipe exit).

Isothermal Flow (Gas Flow with Heat Transfer Through Uninsulated Walls)

In long, uninsulated buried transmission pipelines, rapid heat conduction through the pipe wall maintains the gas at constant temperature ($T = T_{ambient} = \text{constant}$):

• As pressure drops along the pipeline due to wall friction, gas expands, density decreases, and velocity rises.
• Unlike adiabatic Fanno flow, the limiting maximum velocity in isothermal flow is not sonic ($Ma = 1$), but rather: $$v_{limit} = \sqrt{\frac{R T}{M}} = \frac{c}{\sqrt{k}} \implies Ma_{limit} = \frac{1}{\sqrt{k}}$$ For air ($k = 1.4$), isothermal flow chokes when $Ma = 1/\sqrt{1.4} = 0.845$.

---

8. Summary Table: Compressible Flow Regimes and Governing Equations

---

9. Step-by-Step Worked Numerical Example: High-Pressure Blowdown Orifice

Problem Statement

An emergency depressuring restriction orifice is installed in an emergency blowdown line connected to a high-pressure natural gas treator. The treator contains dry methane gas ($M = 16.04\text{ kg/kmol}$, $k = 1.31$, compressibility factor $Z \approx 1.0$) at an upstream stagnation pressure of $P_0 = 4.00\text{ MPa(a)} = 4,000\text{ kPa(a)}$ and a stagnation temperature of $T_0 = 310.0\text{ K}$ ($36.85^\circ\text{C}$). The orifice discharges into an atmospheric flare collection header maintained at a backpressure of $P_b = 300.0\text{ kPa(a)}$ ($3.00\text{ bar(a)}$).

The restriction orifice bore diameter is $d_t = 25.0\text{ mm} = 0.0250\text{ m}$, and the calibrated discharge coefficient is $C_d = 0.950$. The universal gas constant is $R = 8,314.5\text{ J/(kmol}\cdot\text{K)}$.

Determine:

1. The critical pressure ratio $P^* / P_0$ and the critical pressure $P^*$.
2. Whether the restriction orifice is choked.
3. The static temperature $T^*$ and sonic speed $c^*$ at the orifice throat.
4. The gas density at the throat $\rho^*$.
5. The choked mass flow rate $\dot{m}_{choked}$ in kilograms per second ($\text{kg/s}$) and kilograms per hour ($\text{kg/h}$).

Solution

Step 1: Calculate the critical pressure ratio ($P^* / P_0$) and critical pressure ($P^*$). Using the isentropic sonic relation for $k = 1.31$:

$$\frac{P^*}{P_0} = \left(\frac{2}{k + 1}\right)^{\frac{k}{k - 1}} = \left(\frac{2}{1.31 + 1}\right)^{\frac{1.31}{1.31 - 1}} = \left(\frac{2}{2.31}\right)^{\frac{1.31}{0.31}} = (0.865801)^{4.2258}$$

Evaluating the power:

$$\ln(0.865801) = -0.144101$$ $$4.2258 \times (-0.144101) = -0.60894$$ $$\frac{P^*}{P_0} = e^{-0.60894} = 0.54393 \approx 0.544$$

Calculate the absolute critical pressure $P^*$:

$$P^* = 0.54393 \times P_0 = 0.54393 \times 4,000\text{ kPa} = 2,175.7\text{ kPa(a)} \approx 2.176\text{ MPa(a)}$$

Step 2: Determine if the flow is choked. Compare the downstream flare header pressure $P_b$ to the critical throat pressure $P^*$:

$$P_b = 300.0\text{ kPa(a)} < P^* = 2,175.7\text{ kPa(a)}$$ $$\frac{P_b}{P_0} = \frac{300.0}{4,000.0} = 0.0750 \ll 0.544$$

Because the receiver backpressure is substantially below the critical choking pressure, the orifice is completely choked ($Ma_t = 1.0$). The static pressure at the orifice throat is fixed at $P_t = P^* = 2,175.7\text{ kPa(a)}$, regardless of the $300\text{ kPa}$ header pressure.

Step 3: Calculate static temperature ($T^*$) and sonic velocity ($c^*$) at the throat. The throat temperature under choked isentropic expansion is:

$$T^* = T_0 \left(\frac{2}{k + 1}\right) = 310.0\text{ K} \times \left(\frac{2}{2.31}\right) = 310.0 \times 0.865801 = 268.40\text{ K} \quad (-4.75^\circ\text{C})$$

The specific gas constant for methane is:

$$R_{sp} = \frac{R}{M} = \frac{8,314.5\text{ J/(kmol}\cdot\text{K)}}{16.04\text{ kg/kmol}} = 518.36\text{ J/(kg}\cdot\text{K)}$$

The sonic velocity at the throat is:

$$c^* = \sqrt{k R_{sp} T^*} = \sqrt{1.31 \times 518.36\text{ J/(kg}\cdot\text{K)} \times 268.40\text{ K}} = \sqrt{182,258\text{ m}^2/\text{s}^2} = 426.92\text{ m/s}$$

Step 4: Calculate gas density at the throat ($\rho^*$). Using the ideal gas equation of state at throat conditions ($P^*, T^*$):

$$\rho^* = \frac{P^* M}{R T^*} = \frac{2,175.7 \times 10^3\text{ Pa} \times 16.04\text{ kg/kmol}}{8,314.5\text{ J/(kmol}\cdot\text{K)} \times 268.40\text{ K}} = \frac{34,898,228}{2,231,612} = 15.638\text{ kg/m}^3$$

Step 5: Calculate choked mass flow rate ($\dot{m}_{choked}$). The cross-sectional throat area is:

$$A_t = \frac{\pi}{4} d_t^2 = \frac{\pi}{4} (0.0250\text{ m})^2 = 4.9087 \times 10^{-4}\text{ m}^2$$

The mass flow rate is:

$$\dot{m}_{choked} = C_d A_t \rho^* c^* = 0.950 \times (4.9087 \times 10^{-4}\text{ m}^2) \times 15.638\text{ kg/m}^3 \times 426.92\text{ m/s}$$ $$\dot{m}_{choked} = 0.950 \times 3.2771\text{ kg/s} = 3.113\text{ kg/s}$$

Converting to kilograms per hour:

$$\dot{m}_{choked} = 3.1132\text{ kg/s} \times 3,600\text{ s/h} = 11,208\text{ kg/h} \approx 11,210\text{ kg/h}$$

---

10. Common PE Exam Traps in Compressible Flow

1. Using Gauge Pressure in Choking Equations: The critical pressure ratio $P^*/P_0$ and the choked mass flow equation require absolute pressure ($P_{abs} = P_{gauge} + P_{atm}$). Using gauge pressure results in severe calculation errors and incorrect choking determinations.
2. Using Static Temperature Instead of Stagnation Temperature: The total energy of the gas includes its bulk velocity head. In equations for relief valves or blowdown orifices, $T_0$ must be the stagnation temperature of the upstream vessel or reservoir, not the expanded throat temperature.
3. Assuming Downstream Pressure Drops Increase Choked Flow: Once a nozzle or orifice is choked ($P_b \le P^*$), further decreases in receiver backpressure do NOT increase mass flow rate. The throat velocity remains fixed at $Ma = 1.0$.
4. Confusing Specific Gas Constant ($R_{sp}$) with Universal Constant ($R$): In metric units, $R = 8,314.5\text{ J/(kmol}\cdot\text{K)}$. If using specific gas constant $R_{sp} = R/M$, ensure units are in $\text{J/(kg}\cdot\text{K)}$. In US Customary, forgetting $g_c = 32.174$ in $c = \sqrt{k g_c R T / M}$ causes an answer that is off by a factor of $\sqrt{32.174} = 5.67$!