Intuition The one-line ideac ∗ c^* c ∗ "How much pressure did my chamber build up for the amount of propellant I stuffed in?" It measures ONLY what happens inside the combustion chamber and at the throat — it is deliberately blind to the nozzle . If your fuel burns hot and produces low-molecular-weight gas, c ∗ c^* c ∗ combustion efficiency score.
Definition Characteristic velocityc ∗ = P c A ∗ m ˙ c^* = \frac{P_c\, A^*}{\dot m} c ∗ = m ˙ P c A ∗ P c P_c P c A ∗ A^* A ∗ m ˙ \dot m m ˙ velocity (m/s) even though it is not the speed of any physical gas particle — it is a characteristic (bookkeeping) velocity.
The three cousins in rocketry:
c ∗ c^* c ∗ combustion chamber + throat only.C F C_F C F nozzle expansion only.c c c c = c ∗ C F c = c^* \, C_F c = c ∗ C F
This clean split is the whole reason c ∗ c^* c ∗ : it isolates chemistry from geometry.
Intuition Why the throat controls the flowAt the throat the gas reaches Mach 1 . Once flow is choked (sonic at the throat), the mass flow rate m ˙ \dot m m ˙ locked — it depends only on chamber conditions ( P c , T c ) (P_c, T_c) ( P c , T c ) A ∗ A^* A ∗ m ˙ \dot m m ˙ P c P_c P c c ∗ c^* c ∗
We derive the choked mass-flow formula, then rearrange.
Step 1 — Mass flow at the throat.
m ˙ = ρ ∗ A ∗ v ∗ \dot m = \rho^* A^* v^* m ˙ = ρ ∗ A ∗ v ∗ Why this step? Mass conservation: mass flow = density × area × speed, evaluated at the throat where things are sonic.
Step 2 — At the throat the gas is sonic , v ∗ = a ∗ = γ R T ∗ v^* = a^* = \sqrt{\gamma R T^*} v ∗ = a ∗ = γ R T ∗ Why? By definition the throat is where Mach = 1 =1 = 1 γ R T \sqrt{\gamma R T} γ R T
Step 3 — Relate throat conditions to chamber (stagnation) conditions. Using isentropic relations with chamber values T c , P c T_c, P_c T c , P c M = 1 M=1 M = 1 T ∗ T c = 2 γ + 1 , P ∗ P c = ( 2 γ + 1 ) γ γ − 1 , ρ ∗ ρ c = ( 2 γ + 1 ) 1 γ − 1 \frac{T^*}{T_c} = \frac{2}{\gamma+1}, \qquad \frac{P^*}{P_c} = \left(\frac{2}{\gamma+1}\right)^{\frac{\gamma}{\gamma-1}}, \qquad \frac{\rho^*}{\rho_c} = \left(\frac{2}{\gamma+1}\right)^{\frac{1}{\gamma-1}} T c T ∗ = γ + 1 2 , P c P ∗ = ( γ + 1 2 ) γ − 1 γ , ρ c ρ ∗ = ( γ + 1 2 ) γ − 1 1 Why? Flow accelerates from ~rest in the chamber to Mach 1 at the throat isentropically; these are the standard isentropic ratios at M = 1 M=1 M = 1
Step 4 — Ideal gas in the chamber: ρ c = P c R T c \rho_c = \dfrac{P_c}{R T_c} ρ c = R T c P c
Step 5 — Assemble m ˙ \dot m m ˙ Substitute density and speed:
m ˙ = ρ c ( 2 γ + 1 ) 1 γ − 1 ⏟ ρ ∗ ⋅ A ∗ ⋅ γ R T c 2 γ + 1 ⏟ a ∗ \dot m = \underbrace{\rho_c\left(\tfrac{2}{\gamma+1}\right)^{\frac{1}{\gamma-1}}}_{\rho^*} \cdot A^* \cdot \underbrace{\sqrt{\gamma R T_c \tfrac{2}{\gamma+1}}}_{a^*} m ˙ = ρ ∗ ρ c ( γ + 1 2 ) γ − 1 1 ⋅ A ∗ ⋅ a ∗ γ R T c γ + 1 2 ρ c = P c / ( R T c ) \rho_c = P_c/(RT_c) ρ c = P c / ( R T c ) T c T_c T c R R R m ˙ = P c A ∗ R T c γ ( 2 γ + 1 ) γ + 1 2 ( γ − 1 ) \dot m = \frac{P_c A^*}{\sqrt{R T_c}}\,\sqrt{\gamma}\left(\frac{2}{\gamma+1}\right)^{\frac{\gamma+1}{2(\gamma-1)}} m ˙ = R T c P c A ∗ γ ( γ + 1 2 ) 2 ( γ − 1 ) γ + 1
Why the exponent γ + 1 2 ( γ − 1 ) \tfrac{\gamma+1}{2(\gamma-1)} 2 ( γ − 1 ) γ + 1 Because 1 γ − 1 \tfrac{1}{\gamma-1} γ − 1 1 1 2 \tfrac12 2 1 2 / ( γ + 1 ) \sqrt{2/(\gamma+1)} 2/ ( γ + 1 )
Step 6 — Solve for P c A ∗ / m ˙ P_c A^*/\dot m P c A ∗ / m ˙ is c ∗ c^* c ∗
c ∗ = P c A ∗ m ˙ = γ R T c γ ( γ + 1 2 ) γ + 1 2 ( γ − 1 ) \boxed{\;c^* = \frac{P_c A^*}{\dot m} = \frac{\sqrt{\gamma R T_c}}{\gamma}\left(\frac{\gamma+1}{2}\right)^{\frac{\gamma+1}{2(\gamma-1)}}\;} c ∗ = m ˙ P c A ∗ = γ γ R T c ( 2 γ + 1 ) 2 ( γ − 1 ) γ + 1
Intuition Read the boxed formula physicallyc ∗ ∝ T c / M c^* \propto \sqrt{T_c/\mathcal M} c ∗ ∝ T c / M R = R u / M R = R_u/\mathcal M R = R u / M High flame temperature T c T_c T c M \mathcal M M c ∗ c^* c ∗ This is the design mantra: burn hot, make light gases (why H 2 H_2 H 2
c ∗ c^* c ∗ η c ∗ = c measured ∗ c ideal ∗ = ( P c A ∗ / m ˙ ) actual c theoretical ∗ \eta_{c^*} = \frac{c^*_{\text{measured}}}{c^*_{\text{ideal}}} = \frac{(P_c A^*/\dot m)_{\text{actual}}}{c^*_{\text{theoretical}}} η c ∗ = c ideal ∗ c measured ∗ = c theoretical ∗ ( P c A ∗ / m ˙ ) actual clean measure of how completely and efficiently the propellant burned. Typical good engines: η c ∗ ≈ 0.92 \eta_{c^*} \approx 0.92 η c ∗ ≈ 0.92 0.99 0.99 0.99
Why is this the right efficiency metric? A low η c ∗ \eta_{c^*} η c ∗ incomplete combustion, heat loss to walls, or mixing losses . It cannot be blamed on a bad nozzle (that's C F C_F C F
Worked example Example 1 — Computec ∗ c^* c ∗ P c = 5.0 MPa P_c = 5.0\ \text{MPa} P c = 5.0 MPa A ∗ = 0.010 m 2 A^* = 0.010\ \text{m}^2 A ∗ = 0.010 m 2 m ˙ = 25 kg/s \dot m = 25\ \text{kg/s} m ˙ = 25 kg/s c ∗ = P c A ∗ m ˙ = ( 5.0 × 10 6 ) ( 0.010 ) 25 = 2000 m/s c^* = \frac{P_c A^*}{\dot m} = \frac{(5.0\times10^6)(0.010)}{25} = 2000\ \text{m/s} c ∗ = m ˙ P c A ∗ = 25 ( 5.0 × 1 0 6 ) ( 0.010 ) = 2000 m/s Why this step? Direct plug-in of the measured definition — no nozzle info needed, confirming c ∗ c^* c ∗ Pa ⋅ m 2 kg/s = N kg/s = m/s \frac{\text{Pa}\cdot\text{m}^2}{\text{kg/s}} = \frac{\text{N}}{\text{kg/s}} = \text{m/s} kg/s Pa ⋅ m 2 = kg/s N = m/s
Worked example Example 2 — Theoreticalc ∗ c^* c ∗ γ = 1.20 \gamma = 1.20 γ = 1.20 T c = 3200 K T_c = 3200\ \text{K} T c = 3200 K M = 22 g/mol \mathcal M = 22\ \text{g/mol} M = 22 g/mol R = 8314 / 22 = 378 J/kg⋅K R = 8314/22 = 378\ \text{J/kg·K} R = 8314/22 = 378 J/kg⋅K Γ \Gamma Γ = γ + 1 2 ( γ − 1 ) = 2.2 0.4 = 5.5 =\frac{\gamma+1}{2(\gamma-1)} = \frac{2.2}{0.4} = 5.5 = 2 ( γ − 1 ) γ + 1 = 0.4 2.2 = 5.5 ( 2 2.2 ) 5.5 = ( 0.909 ) 5.5 ≈ 0.588 \left(\frac{2}{2.2}\right)^{5.5} = (0.909)^{5.5}\approx0.588 ( 2.2 2 ) 5.5 = ( 0.909 ) 5.5 ≈ 0.588 Γ = 1.2 × 0.588 = 0.644 \Gamma = \sqrt{1.2}\times0.588 = 0.644 Γ = 1.2 × 0.588 = 0.644 c ideal ∗ = γ R T c γ Γ ... equivalently 1 Γ R T c / γ c^*_{\text{ideal}} = \frac{\sqrt{\gamma R T_c}}{\gamma\,\Gamma}\quad\text{... equivalently } \frac{1}{\Gamma}\sqrt{R T_c/\gamma} c ideal ∗ = γ Γ γ R T c ... equivalently Γ 1 R T c / γ R T c = 378 × 3200 = 1.21 × 10 6 = 1100 \sqrt{R T_c} = \sqrt{378\times3200} = \sqrt{1.21\times10^6}=1100 R T c = 378 × 3200 = 1.21 × 1 0 6 = 1100 R T c / γ = 1100 / 1.2 = 1004 \sqrt{R T_c/\gamma}=1100/\sqrt{1.2}=1004 R T c / γ = 1100/ 1.2 = 1004 Γ \Gamma Γ c ideal ∗ = 1004 / 0.644 ≈ 1559 m/s c^*_{\text{ideal}}=1004/0.644\approx 1559\ \text{m/s} c ideal ∗ = 1004/0.644 ≈ 1559 m/s Why this step? We built Γ \Gamma Γ γ \gamma γ R T c \sqrt{RT_c} R T c
Worked example Example 3 — EfficiencyIf the same propellant (ideal c ∗ = 1559 c^* = 1559 c ∗ = 1559 c ∗ = 1480 c^* = 1480 c ∗ = 1480 η c ∗ = 1480 1559 = 0.95 = 95 % \eta_{c^*} = \frac{1480}{1559} = 0.95 = 95\% η c ∗ = 1559 1480 = 0.95 = 95% Why this step? Ratio of measured to ideal → 5% of combustion potential lost to incomplete burning / wall heat loss.
c ∗ c^* c ∗ Why it feels right: it has units of m/s and lives in rocket equations next to real velocities. Fix: the physical exhaust speed is v e v_e v e c = c ∗ C F c=c^*C_F c = c ∗ C F c ∗ c^* c ∗ bookkeeping velocity for the chamber — no gas actually moves at c ∗ c^* c ∗
Common mistake "A bigger nozzle raisesc ∗ c^* c ∗ Why it feels right: bigger nozzles raise thrust and I s p I_{sp} I s p Fix: nozzle expansion only changes C F C_F C F c ∗ c^* c ∗ P c , A ∗ , m ˙ P_c, A^*, \dot m P c , A ∗ , m ˙ T c , M , γ T_c, \mathcal M, \gamma T c , M , γ exit does nothing to c ∗ c^* c ∗
Common mistake "Use exit areaA e A_e A e Why it feels right: you see "area" and grab the obvious big one. Fix: it MUST be the throat area A ∗ A^* A ∗ there ; that's the whole physical basis.
Common mistake "Forgetting flow is choked."Why it feels right: the derivation looks like generic gas dynamics. Fix: the P c A ∗ / m ˙ P_c A^*/\dot m P c A ∗ / m ˙ only holds when the throat is sonic (M = 1 M=1 M = 1 c ∗ c^* c ∗
Recall Feynman: explain to a 12-year-oldImagine you blow up a balloon and pinch the neck to a small hole. If you use warm, light air, the balloon gets nice and tight (high pressure) for the same amount of air you blew in. c ∗ c^* c ∗ shape of the nozzle after the pinch — that's a different score. Big c ∗ c^* c ∗
Mnemonic Remember the split"c c c
c ∗ = P c A ∗ / m ˙ c^* = P_c A^*/\dot m c ∗ = P c A ∗ / m ˙ P ressure, throat A rea, mass flow — the three chamber/throat quantities . Nozzle stuff lives in C F C_F C F "hot & light flies right" for c ∗ ∝ T c / M c^*\propto\sqrt{T_c/\mathcal M} c ∗ ∝ T c / M
Choked Flow and the de Laval Nozzle — why M = 1 M=1 M = 1 m ˙ \dot m m ˙ Thrust Coefficient C_F — the nozzle's score; c = c ∗ C F c = c^* C_F c = c ∗ C F Effective Exhaust Velocity and Specific Impulse — I s p = c / g 0 I_{sp} = c/g_0 I s p = c / g 0 Isentropic Flow Relations — source of the 2 / ( γ + 1 ) 2/(\gamma+1) 2/ ( γ + 1 ) Combustion Chamber Thermochemistry — where T c T_c T c M \mathcal M M γ \gamma γ Vandenkerckhove Function Γ — packaging of the γ \gamma γ
Define characteristic velocity by measurement. c ∗ = P c A ∗ / m ˙ c^* = P_c A^*/\dot m c ∗ = P c A ∗ / m ˙ (chamber pressure × throat area ÷ mass flow).
What part of the engine does c ∗ c^* c ∗ Only the combustion chamber + throat (combustion efficiency); it ignores the nozzle.
What part does C F C_F C F C F C_F C F measures nozzle expansion; effective exhaust velocity
c = c ∗ C F c = c^* C_F c = c ∗ C F .
Which area goes in the formula and why? Throat area
A ∗ A^* A ∗ , because the flow is choked (sonic,
M = 1 M=1 M = 1 ) there.
c ∗ c^* c ∗ c ∗ ∝ T c / M c^* \propto \sqrt{T_c/\mathcal M} c ∗ ∝ T c / M — hot flame, low molecular weight → high
c ∗ c^* c ∗ .
Give the theoretical formula for c ∗ c^* c ∗ c ∗ = γ R T c γ ( γ + 1 2 ) γ + 1 2 ( γ − 1 ) c^* = \dfrac{\sqrt{\gamma R T_c}}{\gamma}\left(\dfrac{\gamma+1}{2}\right)^{\frac{\gamma+1}{2(\gamma-1)}} c ∗ = γ γ R T c ( 2 γ + 1 ) 2 ( γ − 1 ) γ + 1 .
Define c ∗ c^* c ∗ η c ∗ = c measured ∗ / c ideal ∗ \eta_{c^*} = c^*_{\text{measured}}/c^*_{\text{ideal}} η c ∗ = c measured ∗ / c ideal ∗ , typically 0.92–0.99.
Why is c ∗ c^* c ∗ clean combustion metric? Both measured and ideal values exclude the nozzle, so a low value points to incomplete burn / heat loss, not nozzle design.
Is c ∗ c^* c ∗ No — it's a bookkeeping "characteristic" velocity; no particle moves at
c ∗ c^* c ∗ .
What is the Vandenkerckhove function Γ \Gamma Γ Γ = γ ( 2 γ + 1 ) γ + 1 2 ( γ − 1 ) \Gamma = \sqrt\gamma\left(\frac{2}{\gamma+1}\right)^{\frac{\gamma+1}{2(\gamma-1)}} Γ = γ ( γ + 1 2 ) 2 ( γ − 1 ) γ + 1 ; then
c ∗ = 1 Γ R T c / γ c^*=\frac1\Gamma\sqrt{RT_c/\gamma} c ∗ = Γ 1 R T c / γ ... equivalently
R T c / ( Γ γ ) \sqrt{RT_c}/(\Gamma\sqrt\gamma) R T c / ( Γ γ ) .
Units check of P c A ∗ / m ˙ P_cA^*/\dot m P c A ∗ / m ˙ Pa ⋅ m 2 kg/s = N kg/s = m/s \frac{\text{Pa}\cdot\text{m}^2}{\text{kg/s}} = \frac{\text{N}}{\text{kg/s}} = \text{m/s} kg/s Pa ⋅ m 2 = kg/s N = m/s .
Choked flow Mach 1 at throat
Combustion chamber + throat
Mass conservation rho A v
Ideal gas rho_c = Pc / R Tc
Effective exhaust velocity c = c-star C_F
Intuition Hinglish mein samjhoDekho, c ∗ c^* c ∗ c ∗ = P c A ∗ / m ˙ c^* = P_c A^*/\dot m c ∗ = P c A ∗ / m ˙ c ∗ c^* c ∗
Sabse important baat — c ∗ c^* c ∗ choked ho jaata hai (Mach 1, sonic). Ek baar sonic ho gaya, to m ˙ \dot m m ˙ P c , T c , A ∗ P_c, T_c, A^* P c , T c , A ∗ c ∗ c^* c ∗ C F C_F C F c = c ∗ × C F c = c^* \times C_F c = c ∗ × C F
Derivation ka core: mass flow m ˙ = ρ ∗ A ∗ v ∗ \dot m = \rho^* A^* v^* m ˙ = ρ ∗ A ∗ v ∗ v ∗ v^* v ∗ c ∗ = γ R T c γ ( γ + 1 2 ) γ + 1 2 ( γ − 1 ) c^* = \frac{\sqrt{\gamma R T_c}}{\gamma}\left(\frac{\gamma+1}{2}\right)^{\frac{\gamma+1}{2(\gamma-1)}} c ∗ = γ γ R T c ( 2 γ + 1 ) 2 ( γ − 1 ) γ + 1 T c / M \sqrt{T_c/\mathcal M} T c / M
Efficiency nikaalne ke liye: η c ∗ = c measured ∗ / c ideal ∗ \eta_{c^*} = c^*_{\text{measured}}/c^*_{\text{ideal}} η c ∗ = c measured ∗ / c ideal ∗ P c A ∗ / m ˙ P_c A^*/\dot m P c A ∗ / m ˙ c ∗ c^* c ∗
