Critical Mach number — onset of local supersonic flow

3.1.24 · Physics › Compressible Flow & Aerodynamics

WHY do we care?

WHY it matters: $M_{cr}$ se upar, wing par ek supersonic flow ka pocket bann jaata hai. Yeh pocket aksar ek shock wave mein khatam hota hai → drag mein achanak badhoti (drag divergence) → woh badnaam "sound barrier" buffeting jo $M_{\infty} < 1$ par hoti hai.
WHAT it is NOT: Yeh nahi hai ki $M_{\infty} = 1$ . Plane abhi bhi subsonic hai; sirf ek local patch sonic hai.
The 80/20: Ek hi idea master karo — "low surface pressure ⇒ high local speed ⇒ pehla point jo sonic hoga" — aur poori derivation khud niklegi.

HOW: derive $M_{cr}$ from first principles

Humein do physics facts aur ek geometry fact chahiye.

Fact 1 — Isentropic pressure–Mach relation

Steady, adiabatic, frictionless (isentropic) compressible flow ke liye, stagnation pressure $p_{0}$ ek streamline ke saath constant rehti hai. Energy aur perfect-gas relations se: $\frac{p _{0}}{p} = (1 + \frac{γ - 1}{2} M^{2})^{\frac{γ}{γ - 1}}$

Why this step? Energy conservation kehti hai ki total temperature constant hai: $T_{0} / T = 1 + \frac{γ - 1}{2} M^{2}$ . Isentropy $p$ aur $T$ ko $p / T^{γ / (γ - 1)} = const$ se link karti hai. Dono combine karo → upar wala formula.

Fact 2 — Definition of pressure coefficient

$C_{p} \equiv \frac{p - p _{\infty}}{\frac{1}{2} ρ _{\infty} V _{\infty}^{2}}$ $\frac{1}{2} ρ_{\infty} V_{\infty}^{2} = \frac{γ}{2} p_{\infty} M_{\infty}^{2}$ use karke (kyunki $a^{2} = γ p / ρ$ aur $V = M a$ ): $C_{p} = \frac{2}{γ M _{\infty}^{2}} (\frac{p}{p _{\infty}} - 1)$

Why this step? Yeh dynamic pressure ko purely $M_{\infty}$ ke terms mein rewrite karta hai — exactly wohi jo humein chahiye.

Putting it together — the critical pressure coefficient $C_{p, cr}$

Usi streamline par same stagnation pressure se:

= \left(\frac{1+\frac{\gamma-1}{2}M_\infty^2}{1+\frac{\gamma-1}{2}M_{local}^2}\right)^{\frac{\gamma}{\gamma-1}}$$ **Local** point ko sonic set karo, $M_{local}=1$, aur free stream ko $M_\infty=M_{cr}$: $$\boxed{\,C_{p,cr} = \frac{2}{\gamma M_{cr}^2}\left[\left(\frac{1+\frac{\gamma-1}{2}M_{cr}^2}{1+\frac{\gamma-1}{2}}\right)^{\frac{\gamma}{\gamma-1}} - 1\right]\,}$$ > [!formula] Yeh curve kya hai > $C_{p,cr}$ vs $M_{cr}$ ek **universal** curve hai (sirf $\gamma$ par depend karti hai). Yeh kehti hai: *"free-stream Mach $M_{cr}$ par local flow ko sonic banane ke liye, minimum surface pressure exactly isi $C_p$ tak dip karni chahiye."* Yeh **airfoil se independent** hai. ### The geometry fact — airfoil ke saath loop close karna Airfoil ka apna **minimum-pressure** curve hota hai: surface ka sabse low $C_p$ kaise badhta hai jab $M_\infty$ badhti hai. Ek simple model hai **Prandtl–Glauert** compressibility correction: $$C_{p} = \frac{C_{p,0}}{\sqrt{1-M_\infty^2}}$$ jahan $C_{p,0}$ incompressible minimum pressure coefficient hai (wing ki shape se ek fixed negative number). > [!intuition] Intersection HI answer hai > $M_\infty$ ke against do curves plot karo: > 1. $C_{p,cr}(M_\infty)$ — **gas-dynamics** requirement (universal). > 2. $C_{p,\min}(M_\infty)=C_{p,0}/\sqrt{1-M_\infty^2}$ — **airfoil** ki actual lowest pressure. > > Jahan yeh **cross** karte hain, wing ka minimum pressure point *abhi abhi* sonic reach kiya hai. Woh crossing $M_\infty$ **hi** $M_{cr}$ hai. ![[3.1.24-Critical-Mach-number-—-onset-of-local-supersonic-flow.png]] --- ## Worked examples > [!example] Example 1 — Kisi diye $M_{cr}$ ke liye $C_{p,cr}$ compute karo > $M_{cr}=0.70$, $\gamma=1.4$ ke liye critical pressure coefficient nikalo. > > **Step 1.** Inner ratio: > $\dfrac{1+0.2(0.49)}{1+0.2} = \dfrac{1.098}{1.2}=0.915$. *Why?* Yeh sonic local point par $\frac{p_\infty/p_0}{p^*/p_0}$ hai. > **Step 2.** $\gamma/(\gamma-1)=3.5$ tak raise karo: $0.915^{3.5}=0.728$. *Why?* Isentropic $p$ ratio. > **Step 3.** $C_{p,cr}=\dfrac{2}{1.4(0.49)}(0.728-1)=\dfrac{2}{0.686}(-0.272)=-0.793$. > > **Answer:** $C_{p,cr}\approx -0.79$. Ek wing jiska minimum $C_p$ $M_\infty=0.70$ par $-0.79$ tak pahunche, woh apne critical Mach par hai. > [!example] Example 2 — Ek real airfoil ke liye $M_{cr}$ nikalo > Ek airfoil ka incompressible minimum $C_{p,0}=-0.55$ hai. $M_{cr}$ estimate karo. > > **Step 1.** Airfoil curve: $C_{p,\min}=\dfrac{-0.55}{\sqrt{1-M_\infty^2}}$. *Why?* Suction peak ki Prandtl–Glauert scaling. > **Step 2.** $C_{p,\min}(M_\infty)=C_{p,cr}(M_\infty)$ numerically solve karo. > - $M=0.70$ try karo: airfoil $=-0.55/\sqrt{0.51}=-0.770$; required $C_{p,cr}=-0.793$. Airfoil abhi itna low nahi. > - $M=0.71$ try karo: airfoil $=-0.55/\sqrt{0.496}=-0.781$; $C_{p,cr}(0.71)\approx-0.762$. Ab airfoil curve se *neeche* hai → already sonic ke past. > **Step 3.** Crossing approximately $M_{cr}\approx 0.705$ par hai. > > **Answer:** $M_{cr}\approx 0.70$. *Why this matters:* ek **thinner** airfoil → kam suction → chhota $|C_{p,0}|$ → **zyada** $M_{cr}$. Isliye fast jets ki wings thin hoti hain. > [!example] Example 3 — Qualitative forecast > **Forecast:** Wing ko sweep back karo. $M_{cr}$ badhega ya ghattega? > **Verify:** Sweep leading edge ke *normal* velocity component ko reduce karta hai jo airfoil "dekhti" hai, toh effective $M_\infty$ ban jaata hai $M_\infty\cos\Lambda$. Wing ek *zyada* true $M_\infty$ par sonic reach karti hai. → **$M_{cr}$ badhta hai.** ✅ (Swept wings sound barrier ko kyun beat karti hain.) --- ## Common mistakes > [!mistake] "$M_{cr}=1$ — yahi woh point hai jab plane sound ki speed hit karta hai." > **Why it feels right:** "Sonic" naturally poora craft Mach 1 par suggest karta hai. > **The fix:** Hawa curvature ke upar *accelerate* karti hai, isliye **local** flow free stream se aage nikal jaata hai. Plane apni *skin par kahi* sonic hota hai jabki $M_\infty\approx 0.7$–$0.85$ par fly kar raha hota hai. $M_{cr}<1$ hamesha. > [!mistake] "Zyada negative (lower) $C_{p,cr}$ ⇒ sonic reach karna mushkil." > **Why it feels right:** Bade numbers "zyada extreme" lagte hain. > **The fix:** $C_{p,cr}$ *kam* negative hota jaata hai jab $M_{cr}$ badhta hai. Airfoil ko sirf universal curve *match* karni hoti hai. Relevant comparison yeh hai ki suction peak kaunsi curve se milti hai — **intersection** dekho, sirf magnitude nahi. > [!mistake] Prandtl–Glauert ko $M_{cr}$ tak bilkul exact use karna. > **Why it feels right:** Yeh woh simple closed form hai jo sablog quote karte hain. > **The fix:** P–G ek *linearized* estimate hai; $M_{cr}$ ke paas yeh suction under-predict karta hai. Better corrections (Karman–Tsien, Laitone) thoda *lower* $M_{cr}$ deti hain. P–G ko intuition ke liye use karo, final design ke liye nahi. --- > [!recall]- Feynman: 12-saal ke bachche ko samjhao > Socho tum hose par apna angutha dabate ho — paani narrow gap se **tezi se nikalta** hai. Hawa bhi wing ki **curved top** par yahi karti hai: woh speed up hoti hai. Agar plane itna tez ude, toh wing ke sabse tez spot par hawa sound ki speed tak pahunch jaati hai *bhale hi poora plane abhi bhi awaaz se slow ho*. Jis exact plane-speed par yeh pehli baar hota hai, uska ek naam hai: **critical Mach number**. Iske baad, wing par ek chota "boom zone" banta hai aur plane ko achanak bahut zyada extra drag feel hoti hai — yahi woh "sound barrier" hai jisse early pilots ladte the. > [!mnemonic] Yaad rakho > **"CRITICAL = Curvature Rushes Inner-flow To 1 — Causes A Local boom."** > Aur: *thin & swept ⇒ high $M_{cr}$* → **"Skinny swept wings sneak past sonic."** --- ## #flashcards/physics Critical Mach number $M_{cr}$ kya hai? ::: Woh free-stream Mach number jis par body ke kisi point par flow pehli baar local Mach 1 tak pahunchti hai. $M_{cr}$ 1 se kam, barabar, ya zyada hota hai? ::: 1 se kam — plane abhi bhi subsonic hai jabki sirf ek local patch sonic hai. Wing ki upper surface par flow kyun accelerate karti hai? ::: Curvature streamtube ko constrict karti hai; continuity se flow speed up hoti hai aur local pressure drop hoti hai (suction peak). Isentropic stagnation-to-static pressure ratio batao. ::: $p_0/p = (1+\frac{\gamma-1}{2}M^2)^{\gamma/(\gamma-1)}$. $C_p$ ko $p/p_\infty$ aur $M_\infty$ ke terms mein likho. ::: $C_p = \frac{2}{\gamma M_\infty^2}(p/p_\infty - 1)$. Critical pressure coefficient formula likho. ::: $C_{p,cr}=\frac{2}{\gamma M_{cr}^2}\big[\big(\frac{1+\frac{\gamma-1}{2}M_{cr}^2}{1+\frac{\gamma-1}{2}}\big)^{\gamma/(\gamma-1)}-1\big]$. Kya $C_{p,cr}$ curve airfoil par depend karti hai? ::: Nahi — yeh sirf $\gamma$ par depend karti hai; yeh ek universal gas-dynamics curve hai. $M_{cr}$ graphically kaise nikaalte hain? ::: Universal $C_{p,cr}(M_\infty)$ curve ko airfoil ki minimum-pressure curve $C_{p,0}/\sqrt{1-M_\infty^2}$ se intersect karo. Thinner airfoil ka $M_{cr}$ par kya effect hota hai? ::: Kam suction (chhota $|C_{p,0}|$) ⇒ zyada $M_{cr}$. Wing sweep ka $M_{cr}$ par kya effect hota hai? ::: Badhata hai; airfoil sirf $M_\infty\cos\Lambda$ "dekhti" hai, isliye sonic onset ke liye zyada true $M_\infty$ chahiye. $M_{cr}$ se thoda upar kya hota hai? ::: Ek local supersonic pocket banta hai, jo aksar ek shock se khatam hota hai → drag divergence (sound barrier). $M_{cr}=0.70$, $\gamma=1.4$ ke liye $C_{p,cr}$ kya hai? ::: Lagbhag $-0.79$. --- ## Connections - [[Isentropic Flow Relations]] — $p_0/p$ aur $T_0/T$ provide karta hai. - [[Pressure Coefficient Cp]] — $C_{p,cr}$ banane mein use hone wali definition. - [[Prandtl–Glauert Compressibility Correction]] — airfoil ki minimum-pressure curve. - [[Drag Divergence Mach Number]] — $M_{cr}$ se *thoda upar* kya hota hai. - [[Shock Waves & Normal Shocks]] — local supersonic pocket ko khatam karte hain. - [[Swept Wings & Transonic Design]] — $M_{cr}$ ka engineering response. - [[Area Rule (Transonic)]] — $M_{cr}$ ke baad wave drag minimize karna. ## 🖼️ Concept Map ```mermaid flowchart TD A[Subsonic flight] -->|air speeds over curved wing| B[Local flow accelerates] B -->|first reaches M_local=1| C[Critical Mach M_cr] C -->|above M_cr| D[Supersonic pocket on wing] D -->|ends in| E[Shock wave] E -->|causes| F[Drag divergence] F -->|felt as| G[Sound barrier buffeting] H[Isentropic p0-M relation] -->|combined with| J[Pressure coefficient Cp] J -->|set M_local=1| K[Critical Cp,cr formula] C -->|defines condition for| K K -->|depends only on gamma| L[Universal curve] K -->|intersect airfoil min-pressure| C ```