3.1.6 · D2Compressible Flow & Aerodynamics

Visual walkthrough — Area-Mach number relation A - A - = f(M) — isentropic flow

1,880 words9 min readBack to topic

Step 0 — The picture we start from

WHAT. A tube whose width changes. Air enters from a big calm reservoir on the left and threads through the tube. Somewhere the tube pinches to its narrowest point — the throat.

WHY. Everything we prove is a statement about this streamtube, so we must draw it first and name its parts before any symbol is allowed to appear.

PICTURE.

Figure — Area-Mach number relation A - A -  = f(M) — isentropic flow

Everything below is a conversation between these four quantities.


Step 1 — Mass cannot pile up (continuity)

WHAT. Whatever mass of air crosses one doorway per second must cross every doorway per second. We write that fixed rate as (kilograms per second):

WHY this tool. We need one thing that is the same everywhere so we can compare two different stations. Conservation of mass (continuity) gives exactly that anchor. Nothing fancier is required yet.

PICTURE. Two doorways, wildly different sizes, but the number of air-molecules crossing per second (the red flux) is identical.

Figure — Area-Mach number relation A - A -  = f(M) — isentropic flow

Step 2 — Invent the yardstick

WHAT. Pick the one imaginary doorway where the air would be moving at exactly the speed of sound, . Call its area , its density , its speed (the sonic speed). Because is the same everywhere, the sonic-station mass flow equals any station's:

WHY. A raw area in cm² tells you nothing on its own. We need a reference size to compare against — like measuring your height in "heads." The sonic doorway is the perfect ruler because it is constant along an isentropic streamtube (same mass flow, same reservoir). See Choked Flow & Maximum Mass Flow for why it is a true constant.

PICTURE. The real station (left) and the imaginary sonic station (right, in red), both fed the same reservoir, both passing the same .

Figure — Area-Mach number relation A - A -  = f(M) — isentropic flow

Now divide both sides by ... actually, rearrange to isolate the ratio we want:


Step 3 — The density ratio

WHAT. We route both densities through the stagnation density (the density in the calm reservoir, where ). The isentropic stagnation relation says, for local Mach :

WHY this tool. We can't compare to directly — they live at different points. But both are tied to the same reservoir . So we insert as a common reference and cancel it:

=\frac{\left(1+\frac{\gamma-1}{2}\cdot 1^2\right)^{-\frac{1}{\gamma-1}}} {\left(1+\frac{\gamma-1}{2}M^{2}\right)^{-\frac{1}{\gamma-1}}}$$ The top just sets $M=1$ (that's what the star means). Term by term: - $\gamma$ — the ratio of specific heats (for air $\approx 1.4$); it controls how compressible the gas is. - The exponent $\frac{1}{\gamma-1}$ comes straight from the isentropic law $p\propto\rho^\gamma$. - Setting $M=1$ in the numerator freezes it to the sonic value. **PICTURE.** A "bridge" diagram: both $\rho$ and $\rho^*$ climb up to the shared reservoir $\rho_0$ and cross over — the reservoir is the common landing they meet on. ![[deepdives/dd-physics-3.1.06-d2-s04.png]] --- ## Step 4 — The speed ratio $a^*/V$ **WHAT.** Our actual speed is $V=M\,a$, where $a=\sqrt{\gamma R T}$ is the local sound speed. So $a^*/V$ splits into an explicit $1/M$ and a temperature piece: $$\frac{a^{*}}{V}=\frac{a^{*}}{M\,a}=\frac{1}{M}\,\frac{a^{*}}{a} =\frac{1}{M}\sqrt{\frac{T^{*}}{T}}$$ **WHY this tool — why a square root?** Because sound speed obeys $a=\sqrt{\gamma R T}$: it grows with the *square root* of temperature. So the ratio of two sound speeds is the square root of the temperature ratio. That's the only reason the $\sqrt{\;}$ appears — it is inherited from the [[Speed of Sound a = sqrt(gamma R T)|definition of $a$]]. Route the temperatures through the reservoir $T_0$ the same way we did densities: $$\frac{T^{*}}{T}=\frac{T^{*}/T_0}{T/T_0} =\frac{\left(1+\frac{\gamma-1}{2}\cdot 1^2\right)^{-1}}{\left(1+\frac{\gamma-1}{2}M^{2}\right)^{-1}}$$ Term by term: - $\frac{1}{M}$ — the raw "how many Mach" factor; big $M$ makes $V$ big, shrinking $a^*/V$. - $\sqrt{T^*/T}$ — corrects for the fact that sound itself travels at different speeds at different temperatures. **PICTURE.** $V$ shown as $M$ copies of the local sound-speed arrow $a$; alongside, the sonic arrow $a^*$ tied to the throat temperature $T^*$. ![[deepdives/dd-physics-3.1.06-d2-s05.png]] --- ## Step 5 — Multiply the two ratios and collect exponents **WHAT.** Put Step 3 and Step 4 into the boxed line of Step 2: $$\frac{A}{A^{*}}=\underbrace{\frac{\rho^{*}}{\rho}}_{\text{Step 3}}\cdot\underbrace{\frac{a^{*}}{V}}_{\text{Step 4}} =\frac{1}{M}\left(1+\frac{\gamma-1}{2}\right)^{-\frac{1}{\gamma-1}+\left(-\frac12\right)} \left(1+\frac{\gamma-1}{2}M^{2}\right)^{\frac{1}{\gamma-1}+\frac12}$$ **WHY.** The two constants $\left(1+\frac{\gamma-1}{2}\right)$ sit side by side, so we simply add their exponents. Same for the two $M$-dependent brackets. Adding the exponents: $$-\frac{1}{\gamma-1}-\frac{1}{2}=-\frac{\gamma+1}{2(\gamma-1)},\qquad \frac{1}{\gamma-1}+\frac{1}{2}=+\frac{\gamma+1}{2(\gamma-1)}$$ The constant term with the *negative* exponent flips up as $\frac{2}{\gamma+1}$ raised to the *positive* exponent, and combines with the $M$-bracket: > [!formula] The Area–Mach relation, assembled > $$\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)}}$$ > - $\frac{1}{M}$ — the speed factor; it dominates at small $M$, blowing the ratio up. > - The bracket — the temperature/density factor; it dominates at large $M$, also blowing it up. > - The tug-of-war between them creates a **minimum**. **PICTURE.** The exponent arithmetic laid out visually: two stacked exponents merging into one. ![[deepdives/dd-physics-3.1.06-d2-s06.png]] --- ## Step 6 — Every case: read the finished curve **WHAT.** Plot $A/A^*$ against $M$ for $\gamma=1.4$ and check all regimes. **WHY.** A formula you can't sanity-check is a liability. We verify the throat, both limits, and the two-root feature — so the reader never meets a case we skipped. **PICTURE.** ![[deepdives/dd-physics-3.1.06-d2-s07.png]] - **$M\to 0$ (dead-slow subsonic):** $1/M\to\infty$ ⇒ $A/A^*\to\infty$. Slow air needs a *huge* pipe to pass the same mass — it's thin and lazy. - **$M=1$ (throat):** plug in — the bracket becomes $\left[\frac{2}{\gamma+1}\cdot\frac{\gamma+1}{2}\right]=1$ and $1/M=1$, so $A/A^*=1$ exactly. The **minimum**. - **$M\to\infty$ (wildly supersonic):** the bracket grows without bound ⇒ $A/A^*\to\infty$ again. To keep accelerating supersonically the pipe must *keep opening up* ([[Area-Velocity Relation dA-A = (M^2-1) dV-V|because $dA/A=(M^2-1)\,dV/V$]] flips sign). - **The two-root fact:** every horizontal line at height $>1$ cuts the U in **two** places — one subsonic (left), one supersonic (right). Geometry alone can't say which; the back-pressure does (see [[Converging-Diverging (de Laval) Nozzle]]). > [!mistake] "Any area gives one Mach number." > **Why it feels right:** the formula looks like a plain function of $M$. > **The fix:** it is a function *from* $M$, but not one-to-one. Reading it *backwards* (area → > Mach) hits the U-curve twice. Only a normal shock or a chosen back-pressure resolves the > ambiguity — see [[Normal Shock Waves]]. --- ## The one-picture summary ![[deepdives/dd-physics-3.1.06-d2-s08.png]] This single frame carries the whole story: mass conservation ties two doorways together, the sonic doorway $A^*$ (red) is the yardstick, stagnation relations convert densities and temperatures into $M$, and the finished U-curve turns *any* area ratio into a Mach number (subsonic or supersonic branch). > [!recall]- Feynman retelling of the whole walkthrough > Picture air marching through a pipe that pinches and swells. Rule one: nobody vanishes — > the same amount of air passes every doorway each second ($\rho A V$ constant). Now imagine > one special doorway, exactly the size where the air hits the speed of sound; call it $A^*$ > and use it as a ruler. Comparing any real doorway to that ruler turns "constant flow" into > "a ratio of packing $\times$ a ratio of speeds." We rewrite the packing ratio using the calm > reservoir as a middleman, and the speed ratio splits into a plain $1/M$ times a square root > of temperatures (square root because sound speed grows like $\sqrt{T}$). Multiply, add the > exponents, and out pops one clean formula. Plot it and it's a valley: bottom at Mach 1 where > area equals the ruler, walls shooting to infinity both for crawling-slow and screaming-fast > flow. And the sneaky part — the same width fits a slow crowd *or* a supersonic sprint, so you > must know which side of the throat you're on. --- ## Active recall > [!recall]- Why is $\dot m=\rho A V$ constant along the tube? > Because mass is conserved — no air is created or destroyed, so the kilograms-per-second > crossing every doorway must be identical. Which single quantity lets us compare two different stations? ::: The mass flow $\dot m=\rho A V$, which is the same everywhere. Why does a square root appear in the $a^*/V$ ratio? ::: Because sound speed $a=\sqrt{\gamma R T}$, so a ratio of sound speeds is a square root of the temperature ratio. What common reference do we route both densities (and temperatures) through? ::: The stagnation reservoir values $\rho_0$, $T_0$, which cancel when taking the ratio. At $M=1$, what does the bracket $\left[\frac{2}{\gamma+1}(1+\frac{\gamma-1}{2}M^2)\right]$ equal? ::: Exactly 1, giving $A/A^*=1$. Why does $A/A^*\to\infty$ at both $M\to0$ and $M\to\infty$? ::: At small $M$ the $1/M$ factor blows up; at large $M$ the bracket factor blows up.