Visual walkthrough — Mass flow rate ṁ and its relation to throat area

3.3.7 · D2 · Physics › Rocket Propulsion › Mass flow rate ṁ and its relation to throat area

---

Step 1 — "Mass flow rate" asal mein kya count karta hai

KYA HAI. Ek seedha pipe socho. Uspar ek vertical line khiincho — ek cross-section. Ek simple sawaal pucho: har second kitne kilograms gas us line se cross karti hai? Woh number mass flow rate hai, likha jaata hai $\dot{m}$ (dot ka matlab hai "per second"), measure hota hai $\text{kg/s}$ mein.

KYU HAI. Rocket ek aisi machine hai jo mass ko peeche fenkhti hai. Isse pehle ki hum baat karein ki throat gas ko kitni tezi se bahar nikalta hai, humein ek clean count chahiye ki har second ek line se kitni gas pass hoti hai. Woh count hi poora khel hai.

PICTURE. Figure mein shaded slug of gas dekho. Ek chote time $dt$ mein gas right taraf $vdt$ length drift karti hai, jahan $v$ uski speed hai. Woh slug ek cylinder hai: area $A$ cross mein, length $vdt$ along mein. (One-dimensional assumption: wahi $v$ face ke har jagah.)

Slug ka volume hai $A\times (vdt)$. Uski mass hai density times volume, $\rho Avdt$. Time $dt$ se divide karo:

$\dot{m}=\underset{\text{kitna packed}}{\underbrace{\rho }}\underset{\text{kitna choda}}{\underbrace{A}}\underset{\text{kitna fast}}{\underbrace{v}}$

Teen honest factors — pack, width, speed. $\dot{m}$ bas yahi hota hai, hamesha.

---

Step 2 — Har station par wahi $\dot{m}$ (continuity)

KYA HAI. Ab pipe ko ek funnel ki tarah draw karo jo narrow aur wide hoti hai. Do lines chuno: station 1 (wide) aur station 2 (narrow). Hum claim karte hain ki dono se wahi $\dot{m}$ cross karta hai.

KYU HAI. Ek sealed nozzle ke andar koi gas create ya destroy nahi hoti, aur na hi koi pile up hoti — yahi steady assumption kaam karti hai. Jo bhi kilograms-per-second ek region mein enter karte hain, woh use exit bhi karte hain. Toh $\dot{m}$ poore channel mein ek constant hai.

PICTURE. Figure mein, wide station par mote slow arrows hain; narrow station par patli fast arrows hain. Same kilograms har second squeeze through hote hain, toh jahan darwaza narrow hai wahan gas zyada fast honi chahiye (ya zyada dense pack honi chahiye).

Yeh ek fact — narrow means fast — hi sab kuch ka beej hai. Sabse narrow point ko throat kehte hain, area $A^{\ast }$. Yahi bottleneck hai jis par hum apna sara dhyan lagaenge.

---

Step 3 — Gas physics ke do tools (aur exactly yahi kyun)

"Narrow means fast" ko ek number mein badalne ke liye, humein gas ke behave karne ke baare mein do facts chahiye.

Tool 1 — ideal gas law. Yeh teen slug-ingredients ko link karta hai: $p=\rho RT⟹\rho =\frac{p}{RT}$ Yahan $p$ = pressure ($\text{Pa}$), $T$ = temperature ($\text{K}$), aur $R$ = is gas ka gas constant ($\text{J/kg⋅K}$). YEH TOOL KYU: humara slug formula $\rho$ chahta hai, lekin density directly measure karna mushkil hai. Pressure aur temperature aasaan hain. Yeh law $\rho$ ko $p$ aur $T$ se trade karta hai.

Tool 2 — speed of sound. Gas mein chote pressure nudges kitni tezi se travel karte hain? $a=\sqrt{\gamma RT}$ $\gamma$ (Greek "gamma") = specific heats ka ratio, gas ke liye ek fixed number (lagbhag $1.2$–$1.4$). YEH TOOL KYU: choking ka poora drama speed of sound par hota hai. Humein gas speed $v$ ko $a$ se compare karna aana chahiye.

PICTURE. Figure ek gas ki bell dikhata hai; pressure ka ek chota "ping" speed $a$ par baahar failt hai. Agar gas khud $a$ se zyada fast flow kar rahi hai, toh ping kabhi upstream nahi ja sakti — yeh yaad rakho, yeh Step 7 ka punchline hai.

---

Step 4 — Chamber vs local: stagnation relations

KYA HAI. Combustion chamber ke andar gas (almost) rest par hoti hai — inhe stagnation conditions $p_0,T_0$ kaho (subscript zero = "rest par laya gaya"). Flow mein gas move kar rahi hai, local values $p,T$ ke saath. Humein unke beech ka bridge chahiye.

KYU HAI. Hum sirf chamber conditions $p_0,T_0$ ko measure aur control kar sakte hain. Throat ki local $p,T$ ek supersonic gas stream ke andar chupi hain. Toh final formula sirf $p_0,T_0$ mein bolna chahiye — lekin physics local $p,T$ par hoti hai. Humein translation chahiye.

PICTURE. Figure temperature ratio ko line-by-line walk karta hai, taaki tum exactly dekh sako ki $a^2=\gamma RT$ aur $M=v/a$ kahan enter karte hain. Neeche diye algebra ke saath saath ise follow karo.

Temperature ratio — energy se (adiabatic), line by line derive kiya. Jab walls se koi heat cross nahi hoti, tab har kilogram ki total energy chamber mein (sab thermal, gas at rest) wahi hoti hai jo stream mein (thermal $+$ kinetic): $\underset{\text{sab thermal, at rest}}{\underbrace{c_p{T}_0}}=\underset{\text{thermal}}{\underbrace{c_{p}T}}+\underset{\text{kinetic}}{\underbrace{\frac{1}{2}{v}^2}}$ Har term ko $c_{p}T$ se divide karo taaki yeh ek ratio ban jaaye: $\frac{T_0}{T}=1+\frac{v^2}{2c_{p}T}.$ Ab do tools feed karo. $c_p=\frac{\gamma R}{\gamma -1}$ replace karo (upar wale box se), aur note karo $v^2=M^2{a}^2=M^2(\gamma RT)$ using $M=v/a$ aur $a^2=\gamma RT$ Step 3 se: $\frac{v^2}{2c_{p}T}=\frac{M^2\gamma RT}{2\cdot \frac{\gamma R}{\gamma -1}\cdot T}=\frac{M^2\gamma RT(\gamma -1)}{2\gamma RT}=\frac{\gamma -1}{2}{M}^2.$ $\gamma RT$ upar aur neeche cancel ho jaata hai — yahi reason hai ki messy constants gayab ho jaate hain aur sirf $\frac{\gamma -1}{2}{M}^2$ bachta hai. Is step mein sirf adiabatic assumption chahiye.

Pressure ratio — isentropic se (adiabatic + frictionless + no shocks). Ek constant-entropy ideal gas ke liye, $p\propto T^{\gamma /(\gamma -1)}$. Poore temperature ratio ko us power tak raise karo: $\frac{p_0}{p}={\left(\frac{T_0}{T}\right)}^{\frac{\gamma }{\gamma -1}}={\left(1+\frac{\gamma -1}{2}{M}^2\right)}^{\frac{\gamma }{\gamma -1}}.$ Yahan hum poora isentropic assumption kharach karte hain — friction ya shock entropy raise kar deta aur $p\propto T^{\gamma /(\gamma -1)}$ tod deta.

$M=0$ par (gas already rest par) dono ratios $1$ hain — local chamber ke barabar, jaisa hona chahiye.

---

Step 5 — $\dot{m}$ ko Mach number ke function ke roop mein assemble karo

KYA HAI. Slug formula $\dot{m}=\rho Av$ lo aur har piece ko apne tools use karke rewrite karo, taaki kuch bhi chuapa na rahe.

KYU HAI. Hum chahte hain ki $\dot{m}$ purely un cheezein mein express ho jo hum jaante hain: area, chamber pressure/temperature, aur Mach number $M$. Phir hum woh $M$ dhundh sakte hain jo flow ko maximise karta hai.

PICTURE. Figure ek "cancellation map" hai: yeh assembled expression ko gas ke thinning ka picture ke upar rakhta hai jaise woh accelerate hoti hai, taaki tum $M$ ke do competing powers dekh sako — upar wala lone $M$ jo flow boost karta hai, aur neeche bada hota bracket jo use thin karta hai — jo Step 6 mein draw tak ladenge.

Density (Tool 1) aur velocity ($v=Ma=M\sqrt{\gamma RT}$, Step 3 se) swap karo:

= p\,A\,M\sqrt{\frac{\gamma}{R T}}$$ Ab local $p,T$ ko chamber $p_0,T_0$ se trade karo Step 4 use karke. $D = 1+\tfrac{\gamma-1}{2}M^2$ likho (shared bracket). Phir $p = p_0\,D^{-\gamma/(\gamma-1)}$ aur $T = T_0\,D^{-1}$, toh $\tfrac{1}{\sqrt T} = D^{1/2}/\sqrt{T_0}$. Substitute karo: $$\dot m = \underbrace{p_0\,D^{-\frac{\gamma}{\gamma-1}}}_{p}\;A\;M\;\sqrt{\frac{\gamma}{R}}\;\underbrace{\frac{D^{1/2}}{\sqrt{T_0}}}_{1/\sqrt T} = A\,p_0\,M\sqrt{\frac{\gamma}{R T_0}}\;\; D^{\,-\frac{\gamma}{\gamma-1}+\frac12}$$ **FINAL EXPONENT KYU.** $D$ ke do powers add karo jo humne abhi collect kiye — ek pressure se, ek $1/\sqrt T$ factor se: $$-\frac{\gamma}{\gamma-1}+\frac12 = \frac{-2\gamma+(\gamma-1)}{2(\gamma-1)} = \frac{-(\gamma+1)}{2(\gamma-1)}.$$ Upar wala negative power neeche positive power ke barabar hai, jo exactly wahi exponent hai jo box mein hai: $$\boxed{\;\dot m(M) = \frac{A\,p_0\,M\sqrt{\dfrac{\gamma}{R T_0}}}{\left(1+\dfrac{\gamma-1}{2}M^2\right)^{\frac{\gamma+1}{2(\gamma-1)}}}\;}$$ Right side par har symbol ya toh fixed hai ($\gamma, R$), controllable hai ($A, p_0, T_0$), ya woh dial hai jise hum ab ghoomayenge ($M$). --- ## Step 6 — Dial ghoomao: kaunsa $M$ sabse zyada mass push karta hai? **KYA HAI.** Area $A$ fixed rakho aur pucho: jaise hum $M$ ko $0$ se upar badhate hain, kya $\dot m$ hamesha badhta rahega? **KYU HAI.** Numerator $\propto M$ flow *grow* karta hai (faster gas). Lekin denominator bhi grow karta hai (gas speed badhne par thin hoti jaati hai). Upar growth, neeche growth — zaroor koi sweet spot hoga jahan dono balance karte hain aur unit area per $\dot m$ *sabse bada* hota hai. **PICTURE.** Figure *flux per unit area*, $\dot m/A$, ko $M$ ke against plot karta hai. Yeh rise hoti hai, peak karti hai, aur fall hoti hai. Peak exactly $M=1$ par hai (chalk-pink dot se mark kiya gaya). ![[deepdives/dd-physics-3.3.07-d2-s06.png]] > [!intuition] Peak exactly $M=1$ par kyun hai > $M=1$ se neeche: gas ko speed up karna thinning par jeet jaata hai — flux badhta hai. $M=1$ se upar: thinning jeet jaati hai — flux girta hai. Exact crossover, jahan dono effects cancel hote hain, sonic point $M=1$ hai. Yeh algebra ka coincidence nahi hai; yeh throat ki defining property hai. Toh ek *fixed* area se maximum flow tab hoti hai jab wahan gas exactly sonic ho. Nozzle ka narrowest point usi peak par baithne ke liye banaya jaata hai. --- ## Step 7 — Choke karo: throat par $M=1$ set karo **KYA HAI.** Throat par $A = A^\ast$ aur $M=1$ set karte hain. Step-5 ke boxed formula mein $M=1$ substitute karo. **KYU HAI.** Yahan do reasons milte hain. (1) Step 6: $M=1$ mass-flow peak hai, toh ek given throat itna hi pass kar sakta hai. (2) Step 3 ki picture: $M=1$ par pressure pings frozen hain — woh throat ke past *upstream travel nahi kar sakti*. Toh chamber literally downstream jo ho raha hai woh "sun" nahi sakta. Exit pressure jitni chahiye kam karo; message kabhi chamber nahi pahunchta, aur $\dot m$ **freeze** ho jaata hai. Yeh freezing ==choking== hai. **PICTURE.** Figure do pings dikhata hai: ek subsonic wala upstream crawl karta hua (chamber sun raha hai) aur ek sonic wala throat par pinned (chamber bahra hai). ![[deepdives/dd-physics-3.3.07-d2-s07.png]] > [!formula] Back-pressure kitna low hona chahiye $M=1$ *force* karne ke liye? > Step-4 pressure ratio mein $M=1$ set karo. Sonic flow ke liye zaruri throat pressure ==critical pressure== hai: > $$\frac{p^\ast}{p_0} = \left(1+\frac{\gamma-1}{2}\cdot 1^2\right)^{-\frac{\gamma}{\gamma-1}} = \left(\frac{2}{\gamma+1}\right)^{\frac{\gamma}{\gamma-1}}$$ > Throat tab choke hota hai ($M=1$) jab downstream (exit/back) pressure itna low ho: > $$\boxed{\;\frac{p_{\text{exit}}}{p_0} \le \left(\frac{2}{\gamma+1}\right)^{\frac{\gamma}{\gamma-1}}\;}$$ > $\gamma=1.2$ ke liye yeh threshold lagbhag $0.56$ hai; $\gamma=1.4$ ke liye, lagbhag $0.53$. Back-pressure ko chamber pressure ke lagbhag aadhe se neeche drop karo aur throat sonic par lock ho jaata hai. Ise us ratio se *upar* push karo aur throat subsonic rehta hai — yeh Step 8 ka "not-yet-choked" Case D hai. Yahi critical ratio [[Choked Flow and Sonic Conditions]] ko drive karta hai. $M=1$ rakhne se numerator ka $M$ 1 ban jaata hai aur bracket $(1+\tfrac{\gamma-1}{2})^{\ldots} = (\tfrac{\gamma+1}{2})^{\ldots}$ ban jaata hai, jo flip hokar deta hai: > [!formula] Choked (maximum) mass flow rate > $$\dot m = \underbrace{A^\ast}_{\text{darwaze ki width}}\;\underbrace{p_0}_{\text{kitna push karte hain}}\;\underbrace{\sqrt{\frac{\gamma}{R\,T_0}}}_{\text{thanda = densa}}\;\underbrace{\left(\frac{2}{\gamma+1}\right)^{\frac{\gamma+1}{2(\gamma-1)}}}_{\text{fixed choking coefficient}}$$ > $\dot m \propto$ ==$A^\ast$== (linear), $\propto$ ==$p_0$== (linear), $\propto$ ==$1/\sqrt{T_0}$==. Yeh [[Choked Flow and Sonic Conditions]] se aage connect hota hai, aur wahi $A^\ast$ downstream expansion ko [[Nozzle Area Ratio and Expansion]] mein drive karta hai. --- ## Step 8 — Edge aur degenerate cases (koi gap mat chhodo) **KYA HAI & KYU HAI.** Ek formula tabhi trustworthy hai jab woh apni extremes par sahi behave kare. Har dial ko uski limit tak push karte hain aur check karte hain ki picture abhi bhi sense banata hai. ![[deepdives/dd-physics-3.3.07-d2-s08.png]] > [!example] Case A — throat kuch bhi nahi rah jaata, $A^\ast \to 0$ > $\dot m \propto A^\ast \to 0$. Ek closed pinch koi gas pass nahi karti. ✔ (Figure ka sabse left wala bar.) > [!example] Case B — chamber pressure zero ho jaata hai, $p_0 \to 0$ > $\dot m \propto p_0 \to 0$. Koi push nahi, koi flow nahi, chahe kitna bhi garam ho. ✔ > [!example] Case C — chamber temperature bahut zyada badh jaata hai, $T_0 \to \infty$ > $\dot m \propto 1/\sqrt{T_0} \to 0$. Fixed pressure par, infinitely hot gas infinitely thin hoti hai ($\rho = p/RT \to 0$), toh lagbhag koi *mass* cross nahi karta chahe woh bahut tezi se bahar nikle. Heat [[Specific Impulse|exhaust speed]] mein madad karti hai, raw kilograms mein nahi. ✔ > [!example] Case D — abhi choke nahi hua ($M<1$, subsonic throat) > Back-pressure Step 7 ke critical ratio se *upar* hai, toh $M<1$: hum Step 6 mein peak ke *left* par hain, aur back-pressure kam karna abhi bhi $\dot m$ badhata hai. Sirf jab $p_{\text{exit}}/p_0$ girkar $\left(\tfrac{2}{\gamma+1}\right)^{\gamma/(\gamma-1)}$ ho jaata hai tab $M$ 1 tak pahunchta hai aur $\dot m$ lock ho jaata hai. Choked formula *ceiling* hai, sirf choking ke baad valid. ✔ > [!mistake] "Choking coefficient $\left(\tfrac{2}{\gamma+1}\right)^{\frac{\gamma+1}{2(\gamma-1)}}$ pressure ke saath badalta hai." > **Fix:** yeh **sirf** $\gamma$ par depend karta hai — gas ki ek pure property. $\gamma=1.2$ ke liye yeh $\approx 0.585$ hai; $\gamma=1.4$ ke liye $\approx 0.578$. Chamber pressure aur temperature *doosre* factors mein rehte hain. --- ## Ek picture ka summary ![[deepdives/dd-physics-3.3.07-d2-s09.png]] Poori derivation ek board par — har arrow apne step number ke saath tagged hai taaki tum ise trace kar sako: **S1** honest count $\dot m = \rho A v$, **S2** continuity ise constant fix karti hai, **S3** gas law aur sound speed, **S4** chamber conditions mein translate karo (isentropic), **S5** $\dot m(M)$ assemble karo, **S6** $M=1$ par flux peak dhundho, **S7** wahan freeze karo boxed choked-flow law paane ke liye, **S8** edge cases sab check out hote hain. Har arrow upar wala ek step hai. > [!recall]- Feynman retelling — walkthrough plain words mein > Main jaanna chahta tha ki engine se har second kitne kilograms garam gas nikalty hain. Maine khud se promise kiya ki flow tidy hogi: steady (time mein nahi badlegi), one-dimensional (har slice mein same), frictionless, aur walls se koi heat leak nahi — woh tidiness hi hai jo sare neat formulas ko sach banati hai. Toh maine pipe ke across ek line draw ki aur us crossing slug ko count kiya: isse mila $\dot m = \rho\,A\,v$ — kitna packed, kitna wide, kitna fast. Kyunki koi gas gayab nahi hoti, wahi count har line se cross hoti hai, toh jahan pipe pinch hoti hai wahan gas rush karni chahiye. "Rush" ko numbers mein badalne ke liye maine do facts udhare: ideal-gas law (pressure aur temperature mujhe packing batate hain) aur speed of sound (taaki main keh sakun kitna fast *sound ke compare mein*, Mach number). Main sirf calm chamber measure kar sakta tha, toh maine energy conservation use ki — ek kilogram ko ek kelvin heat karne mein $c_p$ joules lagte hain — temperature ratio paane ke liye, aur constant-entropy rule $p\propto T^{\gamma/(\gamma-1)}$ pressure ratio paane ke liye. Dono plug karne par, us shared bracket ki powers exactly $-(\gamma+1)/[2(\gamma-1)]$ add ho gayeen — wahan se woh funny exponent aata hai. Phir maine poocha: agar main pinch ko same width rakhu aur gas ko speed up karun, kya hamesha zyada mass through jaayegi? Nahi — speed up karna use thin kar deta hai, aur dono effects exactly speed of sound par tie karte hain. Toh pinch tab sabse zyada mass push karti hai jab wahan gas exactly sonic ho. Aur us speed par, pressure whispers upstream nahi ja sakti, toh chamber bahari duniya ke liye bahra ho jaata hai — lekin sirf tab jab bahari pressure chamber pressure ke lagbhag aadhe se neeche aa jaaye; wahi critical ratio hai jo throat ko lock *force* karta hai. Formula ko wahan freeze karo aur mujhe boxed law milta hai: mass flow pinch width aur push pressure ke saath badhti hai, aur temperature ke square root ke saath ghatti hai. Wahi poori kahani hai, aur yeh [[Thrust Equation and Effective Exhaust Velocity|thrust equation]] aur [[Tsiolkovsky Rocket Equation|Tsiolkovsky's rocket equation]] ko bhi underlie karta hai. --- > [!recall]- > Is walkthrough mein har relation ke neeche kaunse char assumptions hain? ::: Steady, one-dimensional, frictionless, aur adiabatic/no-shocks — saath milke flow ko isentropic banate hain. > Specifically kaunsa assumption pressure ratio $p_0/p=(1+\tfrac{\gamma-1}{2}M^2)^{\gamma/(\gamma-1)}$ ko justify karta hai? ::: Isentropic (constant entropy), $p\propto T^{\gamma/(\gamma-1)}$ ke through; friction ya shock ise tod deta. > $c_p$ kya hai aur uski units kya hain? ::: Specific heat at constant pressure — fixed pressure par ek kilogram ko ek kelvin raise karne ke liye joules, $\text{J/kg·K}$ mein; ideal gas ke liye $c_p=\gamma R/(\gamma-1)$. > Exponent $-(\gamma+1)/[2(\gamma-1)]$ kahan se aata hai? ::: Bracket ke powers ko pressure factor ($-\gamma/(\gamma-1)$) aur $1/\sqrt T$ factor ($+\tfrac12$) se add karne par. > Choked formula mein kaunsa factor sirf gas par depend karta hai, engine settings par nahi? ::: Choking coefficient $\left(\tfrac{2}{\gamma+1}\right)^{\frac{\gamma+1}{2(\gamma-1)}}$ — sirf $\gamma$ ka function. > Throat ko $M=1$ force karne ke liye back-pressure kitna low hona chahiye? ::: $p_{\text{exit}}/p_0 \le \left(\tfrac{2}{\gamma+1}\right)^{\gamma/(\gamma-1)}$ (lagbhag $0.5$–$0.56$). > Throat sonic ho jaane par $\dot m$ kyun freeze ho jaata hai? ::: $M=1$ par pressure signals upstream nahi ja sakti, toh chamber kabhi lower back-pressure nahi sun sakta. > Kaun se Mach number par unit area per flux maximum hoti hai? ::: $M=1$, throat/sonic point. > Jab $T_0 \to \infty$ fixed $p_0$ par, $\dot m$ ka kya hota hai aur kyun? ::: Yeh $\to 0$ ho jaata hai; gas infinitely thin ho jaati hai ($\rho = p/RT \to 0$).