Above Mach 1, You Must Widen to Speed Up — The Laval Nozzle Paradox

On October 14, 1947, Chuck Yeager's X-1 hit Mach 1.06. What surprised the wind-tunnel crews wasn't the sound barrier itself — it was that diffusers behaved like accelerators on the supersonic side, and nozzles behaved like brakes. This post follows the one-line area–Mach identity that explains the swap, shows why a Laval nozzle has to pinch in the middle, and asks where the shock waves come from. A thirty-line NumPy script traces the isentropic curve for one nozzle. Drop the Mach slider below 1, then push it past 2, and the wavefronts fold from circles into a cone that lies behind the source.

Below 1 the waves still run ahead of the source. The moment the slider crosses 1, the envelope becomes a slanted line, and near 2 the half-angle tightens toward 30 degrees. Pause on that line and the region outside the cone stays quiet.

Mach number — one ratio that splits the regimes#

The Mach number (flow speed $U$ over sound speed $a$) is more than a velocity scale. Per-unit-volume inertia goes like $\rho U^2 / l$ and elastic stress like $K l^2 / l^3$, so the square root of their ratio is exactly the Mach number.

Here $K$ is the bulk modulus (stiffness against volume change) and $\rho$ is density. So Mach measures how hard the flow is squeezing on the medium's elasticity. As a rule of thumb, $\mathrm{Ma} < 0.3$ keeps density variation under five percent and lets you pretend the fluid is incompressible. Above that, the flow itself starts compressing what it moves through.

Past the sound barrier — from Doppler to the Mach cone#

A still source emits concentric wavefronts. Move the source and the front shrinks ahead, stretches behind — classical Doppler. Push the source to sonic speed and every front piles onto one plane. Push it past sonic and the source overtakes the fronts it has emitted. Their envelope wraps into a cone, and outside that cone no sound ever arrives.

The half-angle $\mu$ comes from pure geometry.

A Mach 2 jet trails a cone with $\mu \approx 30°$; a Mach 5 vehicle squeezes it to $11.5°$. The slanted line that showed up just after the slider crossed 1 is this half-angle.

The one-line identity — (1−1/Ma²)dρ/ρ + dA/A = 0#

Now to the punchline. Take a steady one-dimensional flow through a duct of varying area $A(x)$. Mass conservation in differential form reads

and the isentropic Euler equation, combined with the sound-speed identity $a^2 = dp/d\rho$, gives

Eliminating $du$ between the two leaves the single line that swallows every paradox.

Sign analysis sweeps everything clean. $\rho$ is density, $A$ is area, $\mathrm{Ma}$ is the Mach number.

Subsonic ($\mathrm{Ma} < 1$, parenthesis negative): shrinking the area ($dA<0$) drops density, so velocity rises. That's the garden hose at your thumb.

Supersonic ($\mathrm{Ma} > 1$, parenthesis positive): every sign flips. A narrower duct now raises density and slows the flow. Trying to squeeze gas faster through a tighter neck actually chokes it. To accelerate, you must let the area grow.

At Mach 0.3, 1/Ma² is 1/0.09, about 11.1, so the parenthesis (1 − 1/Ma²) is about −10.1. Density change is the area change divided by that number. A 10 percent area change moves density by only about 1 percent. Far below sound speed, area and density barely talk to each other.

At Mach 2, 1/Ma² is 0.25 and the parenthesis is +0.75. After the sign flip, density change is about 1.33 times the area change. Open the area by 10 percent and density drops about 13 percent while speed rises. At Mach 1 the parenthesis is zero, so the area change must be zero too or the identity cannot hold. That is the number behind a throat that is a single station.

The Laval nozzle — pinching does not always mean speeding up#

This sign flip dictates the Laval shape. To take subsonic gas all the way to supersonic, you first contract to reach sonic at one slim cross-section, then expand. The throat sits exactly at $\mathrm{Ma}=1$ because the only way the identity stays consistent when $1 - 1/\mathrm{Ma}^2 = 0$ is if $dA = 0$.

Once you fix the throat area $A^*$, every other section follows from Mach alone.

The same area ratio admits two Mach solutions, one subsonic and one supersonic. Which branch the flow picks depends on the back pressure $P_b$. Drop the back pressure low enough and the diverging part rides the supersonic branch with Mach climbing past one. Push back pressure above the design value and a normal shock appears somewhere inside the nozzle, slamming the flow back to subsonic.

Following the curve in NumPy#

The same area–Mach formula plus an isentropic pressure ratio fits in thirty lines. We invert the relation by bisection.

import numpy as np
 
GAMMA = 1.4
 
def area_mach_relation(M, gamma=GAMMA):
    """A/A* — isentropic area-Mach identity."""
    exp = (gamma + 1) / (2 * (gamma - 1))
    inner = (2 / (gamma + 1)) * (1 + 0.5 * (gamma - 1) * M ** 2)
    return (1 / M) * inner ** exp
 
def solve_supersonic_branch(area_ratio, gamma=GAMMA, lo=1.0001, hi=25.0):
    """Area ratio -> supersonic Mach (bisection)."""
    for _ in range(80):
        mid = 0.5 * (lo + hi)
        if area_mach_relation(mid, gamma) < area_ratio:
            lo = mid
        else:
            hi = mid
    return 0.5 * (lo + hi)
 
def pressure_ratio_isentropic(M, gamma=GAMMA):
    """p/p0 — static over stagnation pressure."""
    return (1 + 0.5 * (gamma - 1) * M ** 2) ** (-gamma / (gamma - 1))
 
def normal_shock_p2_p1(M1, gamma=GAMMA):
    """Static pressure ratio across a normal shock."""
    return 1 + (2 * gamma / (gamma + 1)) * (M1 ** 2 - 1)
 
A_ratio = 2.5
M_exit = solve_supersonic_branch(A_ratio)
p_exit = pressure_ratio_isentropic(M_exit)
 
print(f"design exit  Ma = {M_exit:.3f}")
print(f"design exit  p/p0 = {p_exit:.4f}")
print(f"if Ma1=1.5 normal shock, p2/p1 = {normal_shock_p2_p1(1.5):.3f}")

The output:

design exit  Ma = 2.443
design exit  p/p0 = 0.0633
if Ma1=1.5 normal shock, p2/p1 = 2.458

Set the exit-area slider to 2.5 and drop the back-pressure slider under the design static ratio. The design exit Mach on the figure sits near the printed 2.44.

A nozzle with area ratio 2.5 run on design hits Mach 2.44 at the exit and drops static pressure to 6.3 percent of stagnation. Raise the back-pressure ratio above 0.063 and the diverging curve breaks behind the throat, with Mach falling under 1. Embed a Mach 1.5 normal shock inside it and the pressure jumps by a factor 2.46 while the flow snaps back to subsonic.

Put the printed exit Mach 2.443 back into the cone relation. sin μ = 1/2.443 is about 0.409, so the half-angle is about 24 degrees, tighter than the 30 degrees of a Mach 2 cone. The printed 2.458 is not a separate calculation. With γ = 1.4, 2γ/(γ+1) is 2.8/2.4, and Ma² − 1 at Mach 1.5 is 1.25. Add (2.8/2.4)×1.25 to 1 and the result is 2.458.

Shock waves — the limit of compression#

A compression wave that catches up to itself becomes a shock. The trailing edge heats from the compression, sound speed climbs, the rear overtakes the front. The whole pressure jump collapses onto a sheet a few molecular mean free paths thick — micrometers in air.

The normal-shock pressure jump $p_2/p_1 = 1 + \frac{2\gamma}{\gamma+1}(\mathrm{Ma}_1^2 - 1)$ is zero at Mach 1 and grows quadratically. When back pressure exceeds the design value of a supersonic Laval nozzle, the flow plants a shock inside to bridge the pressure mismatch and decelerates itself. That's the over-expanded regime.

The same equation governs the micro-shocks released by collapsing cavitation bubbles, the crack of a whip whose tip just exceeded sound speed, and the heating of meteors entering the atmosphere. Different faces of one identity. The same jump standing on a free surface is the hydraulic jump in a kitchen sink. Only the speed scale changes, from sound speed to a Froude number.

Back pressure above 0.063 kills the exit Mach#

If back pressure sits above the design static ratio 0.0633, the diverging section cannot stay on the supersonic branch. Mach 2.443 and a static pressure of 6.33 percent of stagnation come as a pair. Raise only the back pressure and the exit Mach falls under 1 even though the area ratio is still 2.5. A Mach 1.5 shock standing inside multiplies pressure by 2.458 and closes that gap.

What locks the mass flow once the throat is sonic sits outside this note. Where a converging-diverging nozzle chokes the mass flow takes the flow rate after Mach 1 at the throat. If area ratio and back pressure already match and the mass flow still does not, that is the place to check whether the throat is really sonic.