3.1.12Compressible Flow & Aerodynamics

Normal shock properties — M₂, P₂ - P₁, T₂ - T₁, ρ₂ - ρ₁, P₀₂ - P₀₁

1,939 words9 min readdifficulty · medium

1. The three conservation laws (first principles)

Take a control volume straddling the shock. State 1 = upstream (supersonic), state 2 = downstream (subsonic). Steady, 1-D, adiabatic, no area change, no friction work on the CV walls.

Why these three? Mass can't be created, momentum changes only via net force (here the pressure difference), and energy is conserved because no heat is added and the shock does no shaft work. Five unknowns (P2,T2,ρ2,u2,M2P_2,T_2,\rho_2,u_2,M_2), and these laws + equation of state close the system.


2. Deriving the downstream Mach number M2M_2 (the master result)

Everything is cleanest in terms of Mach number. Use u=Mau = Ma, a2=γRTa^2=\gamma R T, so ρu2=ργRTM2=γPM2\rho u^2 = \rho \gamma R T M^2 = \gamma P M^2.

Step A — rewrite momentum (2): Divide nothing, just substitute ρu2=γPM2\rho u^2=\gamma P M^2: P_1(1+\gamma M_1^2) = P_2(1+\gamma M_2^2)\;\Rightarrow\; \frac{P_2}{P_1}=\frac{1+\gamma M_1^2}{1+\gamma M_2^2}\tag{4} Why this step? Converts the momentum balance into a pure pressure ratio in Mach numbers.

Step B — rewrite energy (3): Stagnation temperature T0=T(1+γ12M2)T_0=T\left(1+\frac{\gamma-1}{2}M^2\right) is conserved, so: \frac{T_2}{T_1}=\frac{1+\frac{\gamma-1}{2}M_1^2}{1+\frac{\gamma-1}{2}M_2^2}\tag{5} Why this step? Energy conservation directly gives the temperature ratio once we know both Mach numbers.

Step C — rewrite mass (1): ρu=PRTMa=PMγRTRT=PMγRT\rho u = \frac{P}{RT}Ma = \frac{P M\sqrt{\gamma R T}}{RT}=\frac{PM\sqrt{\gamma}}{\sqrt{RT}}. Setting state 1 = state 2: \frac{P_1 M_1}{\sqrt{T_1}}=\frac{P_2 M_2}{\sqrt{T_2}}\tag{6}

Step D — combine. Put (4) and (5) into (6). After algebra (square it to kill the roots) the PP and TT ratios collapse into a single equation relating M1M_1 and M2M_2, whose non-trivial root is:

Why this is beautiful: M2M_2 depends on M1M_1 only. The trivial root M2=M1M_2=M_1 is the "no shock" solution; equation (7) is the genuine shock branch. For any M1>1M_1>1, equation (7) gives M2<1M_2<1supersonic always becomes subsonic.


3. The property ratios (all in terms of M1M_1)

Substitute (7) back into (4), (5), use mass for density. With γ=1.4\gamma=1.4 unless noted.


4. Stagnation pressure loss P02/P01P_{02}/P_{01} — the entropy fingerprint

T0T_0 is conserved (adiabatic), but P0P_0 is NOT — it drops because entropy rises. From the entropy change of an ideal gas Δs=cplnT2T1RlnP2P1\Delta s = c_p\ln\frac{T_2}{T_1}-R\ln\frac{P_2}{P_1}, and P0=P(1+γ12M2)γ/(γ1)P_0=P\left(1+\frac{\gamma-1}{2}M^2\right)^{\gamma/(\gamma-1)}:

Since Δs>0\Delta s>0 for all M1>1M_1>1, P02/P01<1P_{02}/P_{01}<1 always. This loss is the engineer's enemy: in a supersonic inlet, a strong normal shock destroys stagnation pressure, hurting engine thrust. (Hence multiple weak oblique shocks are preferred.)

Figure — Normal shock properties — M₂, P₂ - P₁, T₂ - T₁, ρ₂ - ρ₁, P₀₂ - P₀₁

5. Worked examples


6. Common mistakes (Steel-man + fix)


Recall Feynman: explain to a 12-year-old

Imagine running so fast that the air can't get out of your way in time. It bunches up in front of you into a super-thin, super-squished wall of air. Crossing that wall, in an instant the air gets hotter, heavier, higher-pressure, and slower than sound. Energy is saved (nothing is burned), but the squishing is messy and wasteful — you can never un-squish it perfectly. That "wasted neatness" is why the air's push-power (stagnation pressure) goes down even though its total energy (stagnation temperature) stays the same.


Flashcards

What is conserved across a normal shock?
Stagnation enthalpy / temperature T0T_0 (adiabatic), plus mass, momentum, energy. NOT entropy, NOT stagnation pressure.
Why is a normal shock not isentropic?
It is irreversible; entropy increases (Δs>0\Delta s>0), so isentropic relations cannot be used across it.
Formula for M22M_2^2 in terms of M1M_1
M22=1+γ12M12γM12γ12M_2^2=\dfrac{1+\frac{\gamma-1}{2}M_1^2}{\gamma M_1^2-\frac{\gamma-1}{2}}
Pressure ratio P2/P1P_2/P_1 across shock
1+2γγ+1(M121)1+\dfrac{2\gamma}{\gamma+1}(M_1^2-1)
Density ratio ρ2/ρ1\rho_2/\rho_1
(γ+1)M12(γ1)M12+2\dfrac{(\gamma+1)M_1^2}{(\gamma-1)M_1^2+2}, equals u1/u2u_1/u_2.
Limiting density ratio as M1M_1\to\infty (air)
γ+1γ1=6\dfrac{\gamma+1}{\gamma-1}=6
How is T2/T1T_2/T_1 obtained from P and ρ ratios?
T2T1=P2/P1ρ2/ρ1\dfrac{T_2}{T_1}=\dfrac{P_2/P_1}{\rho_2/\rho_1} (ideal gas law).
Why does P02/P01<1P_{02}/P_{01}<1?
Because Δs>0\Delta s>0, and P02/P01=eΔs/RP_{02}/P_{01}=e^{-\Delta s/R}.
M2M_2 and P2/P1P_2/P_1 for M1=2M_1=2, air
M2=0.577M_2=0.577, P2/P1=4.5P_2/P_1=4.5.
Can a normal shock occur in subsonic flow?
No — it would require Δs<0\Delta s<0, violating the 2nd law. Needs M1>1M_1>1.

Connections

  • Isentropic Flow Relations — valid before/after, not across the shock.
  • Speed of Sound and Mach Number — why M>1M>1 is the trigger.
  • Oblique Shock Waves — turning + weaker P0P_0 loss; built from normal-shock physics.
  • Rankine–Hugoniot Relations — the PPρ\rho form of these jumps.
  • Entropy and Second Law — source of stagnation-pressure loss.
  • Supersonic Inlets and Diffusers — engineering reason we minimise normal shocks.
  • Stagnation Properties T0 and P0 — what survives vs what's lost.

Concept Map

gives

gives

gives

closes system

substitute rho u sq equals gamma P M sq

T0 conserved

rewrite in Mach

combine

combine

combine

non-trivial root

feeds back

feeds back

density via EOS

entropy rises

Control volume across shock

Mass continuity eq 1

Momentum eq 2

Energy eq 3

Ideal gas law P equals rho R T

Combine 4 5 6

Pressure ratio P2 over P1

Temperature ratio T2 over T1

Mass in Mach form eq 6

Downstream Mach M2 subsonic

Density ratio rho2 over rho1

Stagnation pressure drops P02 over P01

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho, normal shock ek bahut patli si "deewar" hai jahan supersonic flow (M₁ > 1) achanak subsonic (M₂ < 1) ban jaata hai. Kyun? Kyunki gas mein information sound ki speed se travel karti hai. Jab flow sound se tez chal raha ho, to aage wali gas ko "warning" nahi milti ki jagah banao — isliye pressure pile-up ho jaata hai aur ek sudden jump (shock) ban jaata hai. Is shock ke aar-paar P, T aur density sab badh jaate hain, velocity ghat jaati hai.

Saari formulas hum sirf teen basic laws se nikaalte hain: mass, momentum, energy conservation, plus ideal gas law. In tino ko Mach number ki bhasha mein likho to ek master equation milti hai — M2M_2 sirf M1M_1 par depend karta hai. Yaad rakho: P2/P1P_2/P_1, T2/T1T_2/T_1, ρ2/ρ1\rho_2/\rho_1 sab M1M_1 ke functions hain, aur M1=2M_1=2 (air) ke liye yeh approx 4.5, 1.69, 2.67 aate hain, M2=0.577M_2=0.577.

Sabse important concept: T0T_0 (stagnation temperature) constant rehti hai kyunki process adiabatic hai (energy bachi rehti hai), lekin P0P_0 (stagnation pressure) GIR jaati hai kyunki shock irreversible hai, entropy badhti hai. Isi liye P02/P01<1P_{02}/P_{01} < 1 hamesha. Yeh loss engineering mein bada dushman hai — supersonic engine inlets mein strong normal shock se thrust kam ho jaata hai, isliye engineers chhote-chhote oblique shocks use karte hain.

Ek galti se bacho: shock ke across kabhi isentropic relations mat lagao — woh sirf shock se pehle aur baad mein valid hain. Aur ek mast fact: chahe shock kitna bhi strong ho, air ki density max 6 guna hi compress hoti hai (γ+1γ1\frac{\gamma+1}{\gamma-1}), baaki energy temperature aur pressure mein chali jaati hai.

Go deeper — visual, from zero

Test yourself — Compressible Flow & Aerodynamics

Connections