3.1.17 · Physics › Compressible Flow & Aerodynamics
Jab ek supersonic flow apne aap se door mud-ti hai (ek convex corner par), toh woh ek single shock se nahi mud sakti — yeh physics ke rules tod-ti hai. Balki woh dheere-dheere ek infinite fan of weak Mach waves se mud-ti hai, har ek flow ko thoda sa mod-ti hai aur thoda sa speed up karti hai. Prandtl–Meyer function ν ( M ) ek bookkeeping device hai: yeh batati hai ki ek flow M = 1 se kisi given Mach number M tak accelerate hone mein kitna total angle "kharch" kar chuka hai . Flow ko angle Δ θ se mod-ne par simply ν utni hi amount se badal jaati hai.
Definition Prandtl–Meyer function
ν ( M ) woh angle (radians ya degrees mein) hai jis se ek initially sonic (M = 1 ) flow ko isentropically expand karke Mach number M tak pahunchaya jaata hai . Isse aise define kiya gaya hai ki ν ( 1 ) = 0 .
Ek isentropic expansion turn ke liye jiska magnitude θ ho:
θ = ν ( M 2 ) − ν ( M 1 )
jahan M 2 > M 1 (flow accelerate hoti hai) aur turn flow se door hota hai.
Yeh key practical baat hai: turning angle, ν ke change ke barabar hota hai. Koi shock tables nahi, koi iteration on shock relations nahi — bas do lookups.
Intuition Expansion ≠ compression kyun
Ek compression corner (flow apni taraf mud-ti hai) Mach waves ko ek saath stack kar deta hai → woh ek shock mein mil jaate hain (entropy badhti hai). Ek expansion corner (flow door mud-ti hai) Mach waves ko alag-alag phailaata hai → woh kabhi cross nahi karti → koi discontinuity nahi → isentropic . Kyunki yeh isentropic hai, hum smooth differential relations use kar sakte hain aur unhe integrate kar sakte hain. Woh integral hi ν ( M ) hai.
Har infinitesimal Mach wave local flow ke saath Mach angle μ banaati hai,
μ = arcsin ( M 1 ) .
Hum ek tiny flow deflection d θ aur ek Mach wave ke across tiny speed change ke beech relation banate hain.
Step 1 — Mach wave ke across velocity triangle.
Ek Mach wave itni weak hoti hai ki sirf wave ke normal velocity component hi change hota hai; tangential component preserved rehta hai. Velocity V ko wave ke along aur across resolve karke aur ek small turn d θ apply karke, velocity triangle ki geometry deti hai:
d θ = V M 2 − 1 d V .
Yeh step kyun? Normal Mach number M sin μ = 1 hai, aur tangential balance iss exact factor M 2 − 1 = M 2 sin 2 μ ⋅ s i n 2 μ 1 − 1 ko force karta hai… concretely cot μ = M 2 − 1 triangle se aata hai.
Step 2 — d V / V ko d M / M mein convert karo.
V = M a aur isentropic temperature relation a 2 a 0 2 = 1 + 2 γ − 1 M 2 se, log lo aur differentiate karo:
V d V = 1 + 2 γ − 1 M 2 1 M d M .
Yeh step kyun? V depend karta hai M par directly bhi aur local speed of sound a ke through bhi; accelerating flow mein temperature drop V ki growth ko slow kar deta hai, isliye denominator aata hai.
Step 3 — Combine karo.
d θ = 1 + 2 γ − 1 M 2 M 2 − 1 M d M .
Step 4 — M = 1 (jahan θ = 0 ) se M tak integrate karo:
ν ( M ) = ∫ 1 M 1 + 2 γ − 1 M 2 M 2 − 1 M d M .
Step 5 — Closed form (standard integral, M 2 − 1 trig substitute karo):
ν m a x ko samajhna
ν m a x ≈ 130.4 5 ∘ woh maximum turn hai jo ek flow expansion se le sakti hai M = ∞ (pressure → 0) tak pahunchne se pehle. Agar kisi corner ko ν m a x − ν ( M 1 ) se zyada turn chahiye, toh flow aur wall ke beech ek vacuum / void ban jaata hai.
Worked example Example 1 — Corner ke upar Expansion
Air (γ = 1.4 ) M 1 = 2.0 par ek θ = 1 5 ∘ convex corner se mud-ti hai. M 2 nikalo.
Step 1. ν ( M 1 ) compute karo. M = 2 formula mein plug karne par ν ( 2 ) = 26.3 8 ∘ milta hai.
Kyun? Yeh "angle budget jo already kharch ho chuka hai" M = 2 tak sonic se pahunchne mein hai.
Step 2. Turn add karo: ν ( M 2 ) = ν ( M 1 ) + θ = 26.3 8 ∘ + 1 5 ∘ = 41.3 8 ∘ .
Kyun? Expansion badhati hai ν ko exactly turn angle se (boxed relation).
Step 3. ν ( M 2 ) = 41.3 8 ∘ ⇒ M 2 ≈ 2.60 invert karo.
Kyun? ν , M mein monotonic hai, isliye ek ν ↔ ek M .
Worked example Example 2 — Fan ke across Pressure ratio
Example 1 ke liye, p 2 / p 1 nikalo.
Step 1. Flow isentropic hai, isliye p 0 constant hai: p 1 p 2 = p 1 / p 0 p 2 / p 0 .
Step 2. Use karo p 0 p = ( 1 + 2 γ − 1 M 2 ) − γ / ( γ − 1 ) .
M 1 = 2 ke liye: p 1 / p 0 = 0.1278 . M 2 = 2.6 ke liye: p 2 / p 0 = 0.0501 .
Step 3. p 2 / p 1 = 0.0501/0.1278 ≈ 0.392 .
Kyun? Expansion pressure girata hai — flow accelerate hone mein kaam karta hai, thanda hota hai aur rarefied hota hai.
Worked example Example 3 — Maximum turn check
Flow M 1 = 3 par enter karti hai (ν = 49.7 6 ∘ ). Sabse bada possible turn kya hai?
θ m a x = ν m a x − ν ( M 1 ) = 130.4 5 ∘ − 49.7 6 ∘ = 80.6 9 ∘ .
Kyun? Isse aage, M → ∞ , p → 0 : physically vacuum ban jaata hai.
Common mistake "Compression corners ke liye bhi
ν use karo."
Kyun sahi lagta hai: ν itna convenient hai ki tum ise har jagah use karna chahte ho, aur ek compression "negative expansion" jaisa lagta hai.
The fix: Compression corners shocks create karte hain, jo non-isentropic hote hain. ν ( M ) isentropic flow assume karke derive kiya gaya tha. Compression ke liye oblique-shock relations use karo. (Ek smooth isentropic compression , jaise concave wall, ν ko decrease karte hue use kar sakta hai — lekin ek sharp corner nahi kar sakta.)
Common mistake "Turn angle seedha
M mein add karo."
Kyun sahi lagta hai: Dono angles/numbers hain, M 2 = M 1 + θ likhna tempting lagta hai.
The fix: Tum θ , ν mein add karte ho, M mein nahi. ν aur M nonlinearly related hain; tumhe M 1 → ν 1 → ν 2 = ν 1 + θ → M 2 jaana hoga.
θ aur ν ke units bhool jaao."
Degrees aur radians mix karna sab kuch kharaab kar deta hai. Tables usually ν degrees mein deti hain; θ bhi degrees mein rakho. Closed-form integral naturally radians mein deta hai.
Recall Feynman: ek 12-saal ke bachche ko explain karo
Socho tum ek curving wall ke andar ke edge ke paas bahut tez daud rahe ho, sound se bhi tez. Jab wall suddenly tumse door mud-ti hai, tumhe thodi extra jagah milti hai, toh tum speed up karte ho aur phail jaate ho — jaise paani ek ramp ke end se rush karta hai. Yeh ek achanak "bump" se nahi hota; yeh bahut saari choti fan-shaped pushes se hota hai. Prandtl–Meyer number bas ek score hai jo count karta hai ki tumhari flow ne kitna "spreading-out turning" kiya hai jab se woh barely supersonic thi. Agar tum wall ko 15 degrees door mod-do, tera score 15 se badh jaata hai. Score se tum padh sakte ho ki tum ab kitni tez ja rahe ho.
"E.A.V." — E xpansion → ν A dd hoti hai (ν 2 = ν 1 + θ ) → V elocity (Mach) badhti hai. Compression mein S ubtract hota, lekin S hock use karo uske liye.
Kaun sa reference state ν = 0 define karta hai?
Expansion fan isentropic kyun hai lekin compression corner nahi?
Air ke liye ν m a x kya hai aur wahan physically kya hota hai?
Prandtl–Meyer function ν(M) physically kya represent karta hai? Woh angle jis se ek sonic (M=1) flow ko isentropically expand karke Mach number M tak pahunchaya jaata hai; ν(1)=0.
Expansion ke liye turning angle aur ν ka relation? θ = ν(M₂) − ν(M₁), jisme M₂ > M₁.
Expansion corner smooth fan kyun form karta hai shock ki jagah? Mach waves alag-alag phailti hain (kabhi coalesce nahi hoti), isliye flow isentropic rehti hai — koi discontinuity nahi.
M ke terms mein local Mach angle μ? μ = arcsin(1/M).
Derivation ko drive karne waala differential relation? dθ = [√(M²−1) / (1 + (γ−1)/2·M²)] · dM/M.
ν(M) ka closed form? ν = √((γ+1)/(γ−1))·arctan√((γ−1)/(γ+1)·(M²−1)) − arctan√(M²−1).
γ=1.4 ke liye ν_max aur uska matlab? ≈130.45°; M→∞ (p→0) ke liye max turn; zyada turn ⇒ vacuum ban jaata hai.
Sharp compression corner ke liye ν(M) use kyun nahi kar sakte? Woh ek shock produce karta hai, jo non-isentropic hai; ν isentropic flow assume karke derive kiya gaya tha.
Prandtl–Meyer fan ke across static pressure badhta hai ya girta hai? Girta hai (expansion: flow accelerate, cool, aur rarefy hoti hai); p₀ constant rehta hai.
M₁ aur turn θ diye ho toh M₂ nikalne ka procedure? M₁→ν₁; ν₂=ν₁+θ; ν₂→M₂ invert karo.
Oblique Shock Waves — compression counterpart (non-isentropic).
Mach Waves and Mach Angle — fan ke building blocks.
Isentropic Flow Relations — p / p 0 , T / T 0 deta hai M 2 nikalne ke baad use karne ke liye.
Method of Characteristics — ν ( M ) ± θ ko Riemann invariants ke roop mein use karta hai.
Expansion Fan / Centered Rarefaction — woh physical structure jise ν describe karta hai.
Nozzle Design (Supersonic) — contour design Prandtl–Meyer expansion par rely karta hai.
gives cot mu = sqrt of M2-1
from V=Ma and isentropic T
Expansion corner turns away
Infinite fan of weak Mach waves
Prandtl-Meyer function v of M
Mach angle mu = arcsin 1/M
d-theta from velocity triangle
Combined differential relation
Turn angle = v of M2 minus v of M1