A solid rocket is basically a controlled fire in a can . The shape of the burning surface (the "grain") sets how much propellant burns each second; the chemistry and nozzle set how fast the gas leaves. Specific impulse I s p I_{sp} I s p I s p I_{sp} I s p thrust level and burn time , but not the efficiency of each kilogram. Understanding this split is the 80/20 of solid motor design.
Definition Specific impulseI s p I_{sp} I s p total impulse delivered per unit weight of propellant consumed :
I s p = ∫ F d t g 0 ∫ m ˙ d t I_{sp} = \frac{\int F\,dt}{g_0 \int \dot m\,dt} I s p = g 0 ∫ m ˙ d t ∫ F d t F F F m ˙ \dot m m ˙ g 0 = 9.81 m/s 2 g_0 = 9.81\ \text{m/s}^2 g 0 = 9.81 m/s 2 bookkeeping constant , not local gravity).
The grain is the shaped block of solid propellant. Its exposed burning surface area A b A_b A b burn rate r r r
Why start here? Thrust needs mass flow, and mass flow is manufactured by the burning grain. This is the one place geometry enters.
The grain recedes normal to its surface at rate r r r d t dt d t r d t r\,dt r d t A b A_b A b
d m = ρ p A b ( r d t ) dm = \rho_p \, A_b \, (r\,dt) d m = ρ p A b ( r d t )
so the mass generation rate is
m ˙ = ρ p A b r \boxed{\dot m = \rho_p A_b r} m ˙ = ρ p A b r
ρ p \rho_p ρ p A b A_b A b r r r
Why this step? Volume burned per second = area × recession speed; multiply by density to get kg/s. Pure geometry + kinematics.
Empirically the burn rate rises with chamber pressure p c p_c p c r = a p c n r = a\, p_c^{\,n} r = a p c n a a a n n n pressure exponent (typically 0.2 0.2 0.2 0.5 0.5 0.5 < 1 <1 < 1
Thrust from a rocket nozzle (derived from momentum conservation on the control volume):
F = m ˙ v e + ( p e − p a ) A e F = \dot m\, v_e + (p_e - p_a)A_e F = m ˙ v e + ( p e − p a ) A e
The exhaust velocity for isentropic expansion of a hot gas from chamber temperature T c T_c T c h c = h e + 1 2 v e 2 h_c = h_e + \tfrac12 v_e^2 h c = h e + 2 1 v e 2 h = c p T h = c_p T h = c p T
v e = 2 γ γ − 1 R T c M [ 1 − ( p e p c ) γ − 1 γ ] v_e = \sqrt{\frac{2\gamma}{\gamma-1}\frac{R T_c}{M}\left[1-\left(\frac{p_e}{p_c}\right)^{\frac{\gamma-1}{\gamma}}\right]} v e = γ − 1 2 γ M R T c [ 1 − ( p c p e ) γ γ − 1 ]
Why this form? All the enthalpy c p T c c_pT_c c p T c p c p_c p c p e p_e p e R R R M M M γ = c p / c v \gamma = c_p/c_v γ = c p / c v
For steady operation, define the effective exhaust velocity c = F / m ˙ c = F/\dot m c = F / m ˙
I s p = c g 0 = F m ˙ g 0 \boxed{I_{sp} = \frac{c}{g_0} = \frac{F}{\dot m\, g_0}} I s p = g 0 c = m ˙ g 0 F
Substitute F = m ˙ v e + ( p e − p a ) A e F = \dot m v_e + (p_e-p_a)A_e F = m ˙ v e + ( p e − p a ) A e
I s p = 1 g 0 [ v e + ( p e − p a ) A e m ˙ ] I_{sp} = \frac{1}{g_0}\left[v_e + \frac{(p_e-p_a)A_e}{\dot m}\right] I s p = g 0 1 [ v e + m ˙ ( p e − p a ) A e ]
For an ideal (optimum-expanded, p e = p a p_e = p_a p e = p a motor the pressure term vanishes:
I s p = v e g 0 = 1 g 0 2 γ γ − 1 R T c M [ 1 − ( p e p c ) γ − 1 γ ] \boxed{I_{sp} = \frac{v_e}{g_0} = \frac{1}{g_0}\sqrt{\frac{2\gamma}{\gamma-1}\frac{R T_c}{M}\left[1-\left(\frac{p_e}{p_c}\right)^{\frac{\gamma-1}{\gamma}}\right]}} I s p = g 0 v e = g 0 1 γ − 1 2 γ M R T c [ 1 − ( p c p e ) γ γ − 1 ]
Intuition THE KEY INSIGHT (80/20)Look at that final formula: it contains γ , T c , M , p e / p c \gamma, T_c, M, p_e/p_c γ , T c , M , p e / p c all chemistry and nozzle quantities . The grain properties ρ p , A b , a , n \rho_p, A_b, a, n ρ p , A b , a , n cancelled out of I s p I_{sp} I s p m ˙ \dot m m ˙
Want more I s p I_{sp} I s p T c T_c T c M M M H 2 O H_2O H 2 O H 2 H_2 H 2 p e / p c → 0 p_e/p_c \to 0 p e / p c → 0
Want more thrust? → increase A b A_b A b
Grain geometry does set chamber pressure , via the throat. Mass in = mass out through throat A t A_t A t ρ p A b r = p c A t c ∗ , c ∗ = 1 Γ R T c M , Γ = γ ( 2 γ + 1 ) γ + 1 2 ( γ − 1 ) \rho_p A_b r = \frac{p_c A_t}{c^*}, \qquad c^* = \frac{1}{\Gamma}\sqrt{\frac{RT_c}{M}},\quad \Gamma=\sqrt{\gamma}\left(\tfrac{2}{\gamma+1}\right)^{\frac{\gamma+1}{2(\gamma-1)}} ρ p A b r = c ∗ p c A t , c ∗ = Γ 1 M R T c , Γ = γ ( γ + 1 2 ) 2 ( γ − 1 ) γ + 1
Using r = a p c n r = a p_c^n r = a p c n p c = ( ρ p a c ∗ A b A t ) 1 1 − n p_c = \left(\frac{\rho_p a\, c^* A_b}{A_t}\right)^{\frac{1}{1-n}} p c = ( A t ρ p a c ∗ A b ) 1 − n 1
Why n < 1 n<1 n < 1 the exponent 1 1 − n \tfrac{1}{1-n} 1 − n 1 n → 1 n\to1 n → 1 n ≥ 1 n\ge1 n ≥ 1 steel-man of "why not use fast-burning high-n n n
Worked example Example 1 — IdealI s p I_{sp} I s p γ = 1.2 \gamma=1.2 γ = 1.2 T c = 3000 K T_c=3000\ \text{K} T c = 3000 K M = 0.025 kg/mol M=0.025\ \text{kg/mol} M = 0.025 kg/mol p e / p c = 0.01 p_e/p_c = 0.01 p e / p c = 0.01 R = 8.314 R=8.314 R = 8.314 Step 1: 2 γ γ − 1 = 2.4 0.2 = 12 \frac{2\gamma}{\gamma-1} = \frac{2.4}{0.2}=12 γ − 1 2 γ = 0.2 2.4 = 12 Why? Prefactor from enthalpy conversion.
Step 2: R T c M = 8.314 ⋅ 3000 0.025 = 9.98 × 10 5 J/kg \frac{RT_c}{M}=\frac{8.314\cdot3000}{0.025}=9.98\times10^{5}\ \text{J/kg} M R T c = 0.025 8.314 ⋅ 3000 = 9.98 × 1 0 5 J/kg Why? Specific gas energy scale.
Step 3: γ − 1 γ = 0.2 1.2 = 0.1667 \frac{\gamma-1}{\gamma}=\frac{0.2}{1.2}=0.1667 γ γ − 1 = 1.2 0.2 = 0.1667 ( 0.01 ) 0.1667 = 0.464 (0.01)^{0.1667}=0.464 ( 0.01 ) 0.1667 = 0.464 = 1 − 0.464 = 0.536 =1-0.464=0.536 = 1 − 0.464 = 0.536 Step 4: v e = 12 ⋅ 9.98 × 10 5 ⋅ 0.536 = 6.42 × 10 6 = 2534 m/s v_e=\sqrt{12\cdot9.98\times10^5\cdot0.536}=\sqrt{6.42\times10^6}=2534\ \text{m/s} v e = 12 ⋅ 9.98 × 1 0 5 ⋅ 0.536 = 6.42 × 1 0 6 = 2534 m/s Step 5: I s p = 2534 / 9.81 = 258 s I_{sp}=2534/9.81 = \mathbf{258\ s} I s p = 2534/9.81 = 258 s Why divide by g 0 g_0 g 0 Converts velocity to "seconds" definition. Note: no grain data used!
Worked example Example 2 — Grain setsm ˙ \dot m m ˙ I s p I_{sp} I s p ρ p = 1800 kg/m 3 \rho_p=1800\ \text{kg/m}^3 ρ p = 1800 kg/m 3 A b = 0.5 m 2 A_b=0.5\ \text{m}^2 A b = 0.5 m 2 r = 0.008 m/s r=0.008\ \text{m/s} r = 0.008 m/s Step 1: m ˙ = ρ p A b r = 1800 ⋅ 0.5 ⋅ 0.008 = 7.2 kg/s \dot m=\rho_p A_b r=1800\cdot0.5\cdot0.008=7.2\ \text{kg/s} m ˙ = ρ p A b r = 1800 ⋅ 0.5 ⋅ 0.008 = 7.2 kg/s Why? Mass generation law.
Step 2: With I s p = 258 I_{sp}=258 I s p = 258 c = g 0 I s p = 2531 m/s c=g_0 I_{sp}=2531\ \text{m/s} c = g 0 I s p = 2531 m/s Step 3: F = m ˙ c = 7.2 ⋅ 2531 = 18 200 N F=\dot m c = 7.2\cdot2531 = \mathbf{18\,200\ N} F = m ˙ c = 7.2 ⋅ 2531 = 18 200 N Why? Thrust = mass flow × effective exhaust speed. Double A b A_b A b same I s p I_{sp} I s p .
Worked example Example 3 — Equilibrium chamber pressurea = 5 × 10 − 5 a=5\times10^{-5} a = 5 × 1 0 − 5 n = 0.35 n=0.35 n = 0.35 c ∗ = 1500 m/s c^*=1500\ \text{m/s} c ∗ = 1500 m/s ρ p = 1800 \rho_p=1800 ρ p = 1800 A b / A t = 200 A_b/A_t=200 A b / A t = 200 Step 1: base = ρ p a c ∗ ( A b / A t ) = 1800 ⋅ 5 × 10 − 5 ⋅ 1500 ⋅ 200 = 2.7 × 10 4 =\rho_p a c^* (A_b/A_t)=1800\cdot5\times10^{-5}\cdot1500\cdot200 = 2.7\times10^{4} = ρ p a c ∗ ( A b / A t ) = 1800 ⋅ 5 × 1 0 − 5 ⋅ 1500 ⋅ 200 = 2.7 × 1 0 4 Step 2: exponent 1 1 − n = 1 0.65 = 1.538 \tfrac{1}{1-n}=\tfrac{1}{0.65}=1.538 1 − n 1 = 0.65 1 = 1.538 Step 3: p c = ( 2.7 × 10 4 ) 1.538 ≈ 8.3 MPa p_c=(2.7\times10^4)^{1.538}\approx 8.3\ \text{MPa} p c = ( 2.7 × 1 0 4 ) 1.538 ≈ 8.3 MPa Why? Balance of gas made vs gas vented.
Common mistake Common mistakes (Steel-manned)① "Bigger grain surface → higher I s p I_{sp} I s p
Why it feels right: more surface = more thrust, and thrust feels like performance. Fix: A b A_b A b m ˙ \dot m m ˙ F F F I s p = F / ( m ˙ g 0 ) I_{sp}=F/(\dot m g_0) I s p = F / ( m ˙ g 0 ) A b A_b A b
② "Use g 0 g_0 g 0
Why it feels right: it's called "gravity." Fix: g 0 g_0 g 0 9.81 9.81 9.81 I s p I_{sp} I s p
③ "High pressure exponent n n n
Why it feels right: higher r r r Fix: if n ≥ 1 n\ge1 n ≥ 1 p c ∝ ( ⋅ ) 1 / ( 1 − n ) p_c\propto(\cdot)^{1/(1-n)} p c ∝ ( ⋅ ) 1/ ( 1 − n ) n < 1 n<1 n < 1
④ Forgetting the pressure-thrust term when nozzle isn't matched. At altitude p a p_a p a p e ≠ p a p_e\neq p_a p e = p a ( p e − p a ) A e (p_e-p_a)A_e ( p e − p a ) A e
Recall Feynman: explain to a 12-year-oldImagine a firework that burns from the inside out. How big the burning wall is decides how much smoke it makes every second — that's the push (thrust). But how fast the smoke shoots out depends on how hot and how light the smoke is — that's the "quality" (specific impulse). A giant firework and a tiny one made of the same stuff shoot smoke at the same speed ; the big one just makes more of it. I s p I_{sp} I s p
"GRAIN gives GRAMS, CHEM gives GO."
Grain (area/density/burn rate) → mass flow (grams per second). Chemistry (T c T_c T c M M M γ \gamma γ I s p I_{sp} I s p
Recall Forecast-then-VerifyBefore checking: if I halve exhaust molar mass M M M I s p I_{sp} I s p 2 \sqrt2 2 (Predict, then see VERIFY.) Because v e ∝ 1 / M v_e\propto\sqrt{1/M} v e ∝ 1/ M 2 ≈ 1.41 \sqrt2 \approx 1.41 2 ≈ 1.41
What quantity does grain burning-surface area A b A_b A b The mass flow rate
m ˙ = ρ p A b r \dot m=\rho_p A_b r m ˙ = ρ p A b r (hence thrust & burn time), NOT
I s p I_{sp} I s p .
Write the mass generation law of a solid grain. m ˙ = ρ p A b r \dot m = \rho_p A_b r m ˙ = ρ p A b r (density × burning area × burn rate).
Saint-Robert burn-rate law? r = a p c n r = a\,p_c^{\,n} r = a p c n , with pressure exponent
n < 1 n<1 n < 1 for stability.
Definition of specific impulse? I s p = ∫ F d t g 0 ∫ m ˙ d t = c g 0 I_{sp}=\dfrac{\int F\,dt}{g_0\int\dot m\,dt}=\dfrac{c}{g_0} I s p = g 0 ∫ m ˙ d t ∫ F d t = g 0 c , effective exhaust velocity over
g 0 g_0 g 0 .
Ideal-nozzle I s p I_{sp} I s p I s p = 1 g 0 2 γ γ − 1 R T c M [ 1 − ( p e / p c ) ( γ − 1 ) / γ ] I_{sp}=\dfrac{1}{g_0}\sqrt{\dfrac{2\gamma}{\gamma-1}\dfrac{RT_c}{M}\left[1-(p_e/p_c)^{(\gamma-1)/\gamma}\right]} I s p = g 0 1 γ − 1 2 γ M R T c [ 1 − ( p e / p c ) ( γ − 1 ) / γ ] .
Why don't grain properties appear in ideal I s p I_{sp} I s p They cancel:
I s p = F / ( m ˙ g 0 ) I_{sp}=F/(\dot m g_0) I s p = F / ( m ˙ g 0 ) and
m ˙ = ρ p A b r \dot m=\rho_p A_b r m ˙ = ρ p A b r divides out; only chemistry/nozzle remain.
Three ways to raise I s p I_{sp} I s p Increase
T c T_c T c , decrease exhaust molar mass
M M M , expand more (
p e / p c → 0 p_e/p_c\to0 p e / p c → 0 ).
Why must n < 1 n<1 n < 1 Equilibrium
p c ∝ ( ⋯ ) 1 / ( 1 − n ) p_c\propto(\cdots)^{1/(1-n)} p c ∝ ( ⋯ ) 1/ ( 1 − n ) ;
n ≥ 1 n\ge1 n ≥ 1 gives runaway pressure → explosion.
Equilibrium chamber pressure expression? p c = ( ρ p a c ∗ A b A t ) 1 / ( 1 − n ) p_c=\left(\dfrac{\rho_p a c^* A_b}{A_t}\right)^{1/(1-n)} p c = ( A t ρ p a c ∗ A b ) 1/ ( 1 − n ) .
What is g 0 g_0 g 0 I s p I_{sp} I s p Fixed constant
9.81 m/s 2 9.81\,\text{m/s}^2 9.81 m/s 2 used only for unit conversion (gives seconds), not local gravity.
Saint-Robert r equals a p_c^n
F equals m-dot v_e plus pressure term
sets thrust and burn time, not efficiency
Nozzle expansion p_e / p_c
Effective exhaust velocity c
Intuition Hinglish mein samjhoDekho, solid rocket basically ek "can ke andar controlled aag" hai. Propellant ka jo shape hota hai usko grain kehte hain, aur jo surface jal rahi hoti hai uska area A b A_b A b m ˙ = ρ p A b r \dot m = \rho_p A_b r m ˙ = ρ p A b r r = a p c n r = a\,p_c^n r = a p c n kitna gas banana — matlab thrust aur burn time.
Ab asli twist: specific impulse I s p I_{sp} I s p matlab "per kilogram kitni kick". Jab hum I s p = F / ( m ˙ g 0 ) I_{sp} = F/(\dot m\, g_0) I s p = F / ( m ˙ g 0 ) F = m ˙ v e F=\dot m v_e F = m ˙ v e m ˙ \dot m m ˙ v e v_e v e T c T_c T c M M M p e / p c p_e/p_c p e / p c I s p I_{sp} I s p
Isliye agar zyada I s p I_{sp} I s p T c T_c T c n n n 1 1 1 p c p_c p c ( ⋯ ) 1 / ( 1 − n ) (\cdots)^{1/(1-n)} ( ⋯ ) 1/ ( 1 − n ) Grain = grams, Chemistry = go.
