Solid rocket Isp derivation from grain properties: Visual walkthrough

Step 0 — The words we are allowed to use

Before any algebra, let us agree on the physical pictures. Everything below is a thing you could point at.

Definition The starting objects

Grain — the solid block of fuel. Think of a candle, but shaped.
Burning surface $A_{b}$ — the exposed area that is on fire right now, measured in square metres ( $m^{2}$ ). It is the only place where solid turns to gas.
Burn rate $r$ — how fast the fire eats into the solid, measured in metres per second ( $m/s$ ). The surface marches inward at this speed.
Density $ρ_{p}$ — how many kilograms of solid sit in one cubic metre ( $kg/m^{3}$ ).
Chamber — the sealed space where the gas piles up hot and pressurised.
Chamber pressure $p_{c}$ — the gas pressure inside the chamber (pascals, Pa).
Nozzle exit — the wide open mouth at the very back where gas finally escapes.
Exit pressure $p_{e}$ — the gas pressure at that mouth, after the gas has expanded and cooled (Pa). Always lower than $p_{c}$ .
Molar mass $M$ — the mass of one mole of the exhaust gas ( $kg/mol$ ); "how heavy each gas molecule is." Light molecules (like water vapour, $H_{2}$ ) shoot out faster.
Exhaust — the hot gas that leaves out the back. Its speed is what we ultimately care about.

Prerequisite pictures live in Grain Geometry and Thrust Profiles and Saint-Robert Burn Rate Law; we rebuild the pieces we need here.

Step 1 — Watch a shell of solid turn into gas

WHAT. We ask: in one tiny slice of time $d t$ , how many kilograms of solid become gas?

WHY. Thrust is "mass thrown backward per second times how fast." So the very first thing we need is how much mass the grain manufactures per second. Nothing about rockets works until we have this number.

PICTURE. Look at the figure: the burning surface (cyan) is peeled back by a paper-thin shell of thickness $r d t$ (amber). That shell is the propellant that vaporised during $d t$ .

Figure — Solid rocket Isp derivation from grain properties

The volume of that thin shell is area × thickness:

$d V = how wide the fire is A_{b} \times how deep it ate in time d t (r d t)$

Multiply volume by density to get mass, then divide by $d t$ to get a rate:

$d m = ρ_{p} A_{b} (r d t) ⟹ \overset{m}{˙} = ρ_{p} A_{b} r$

$\overset{m}{˙}$ (kg/s) — the "dot" means per second: mass produced every second.
$ρ_{p} A_{b} r$ — density $\times$ area $\times$ speed. Units check: $\frac{kg}{m ^{3}} \cdot m^{2} \cdot \frac{m}{s} = \frac{kg}{s}$ . ✓

Step 2 — What sets the burn rate? Pressure pushes the fire

WHAT. We open up the $r$ inside $\overset{m}{˙}$ : burn rate is not a fixed number, it rises with chamber pressure $p_{c}$ (defined in Step 0).

WHY. A fire under pressure burns hotter and faster against the surface. Rocket engineers measured this and fit a curve — the Saint-Robert Burn Rate Law.

PICTURE. The figure plots $r$ against $p_{c}$ for a few exponents $n$ . Because $n < 1$ , the curve is concave down — it keeps rising but ever more gently, flattening as pressure climbs. That flattening is exactly what keeps the motor from exploding (a concave-up, steepening curve would be the danger case).

$r = a p_{c}^{n}$

$a$ — a coefficient baked in by the chemistry and temperature (the height of the curve).
$p_{c}$ — chamber pressure (pascals).
$n$ — the pressure exponent, the "steepness" of the response, typically $0.2$ – $0.5$ .

Step 3 — How fast does the gas leave? Trade heat for speed

WHAT. We now switch from "how much gas" to "how fast gas." The nozzle converts the heat energy of the chamber gas into directed speed $v_{e}$ .

WHY. Speed of exhaust is what actually pushes the rocket. The picture to hold: a hot, slow, crowded gas in the chamber squeezes through a throat and fans out, coming out cold and fast.

PICTURE. The figure shows the chamber (hot, high $T_{c}$ , high $p_{c}$ , slow) → throat → nozzle exit (cold, low $p_{e}$ , fast $v_{e}$ ). Energy is conserved: the thermometer falls as the speedometer rises.

Before the algebra, two new symbols. When you heat one kilogram of gas by one degree, the energy it takes depends on whether you let it expand or not:

Now the energy bookkeeping per kilogram of gas ("enthalpy in = enthalpy out + kinetic energy gained"):

$heat energy in chamber c_{p} T_{c} = leftover heat c_{p} T_{e} + kinetic energy of jet \frac{1}{2} v_{e}^{2}$

Solving this alone gives $v_{e} = 2 c_{p} (T_{c} - T_{e})$ — but we cannot measure $T_{e}$ directly. We only know the pressures $p_{c}$ and $p_{e}$ . So we need a bridge from temperatures to pressures.

The bridge — the isentropic relation. For a gas expanding smoothly with no heat leaking in or out (called isentropic), temperature and pressure are locked together:

$\frac{T _{e}}{T _{c}} = (\frac{p _{e}}{p _{c}})^{\frac{γ - 1}{γ}}$

$γ = c_{p} / c_{v}$ — the heat-capacity ratio just defined.
Why this shape? When pressure drops, the gas does work pushing itself outward, and that work is paid for out of its own heat — so temperature must fall alongside pressure. The exponent $\frac{γ - 1}{γ}$ is the exact exchange rate, derived in Nozzle Isentropic Expansion.

Substitute $T_{e} = T_{c} (\frac{p _{e}}{p _{c}})^{(γ - 1) / γ}$ into $v_{e} = 2 c_{p} (T_{c} - T_{e})$ , then use $c_{p} = \frac{γ}{γ - 1} \frac{R}{M}$ (the specific heat written in terms of the gas constant and molar mass). This yields:

$v_{e} = \frac{2 γ}{γ - 1} \frac{R T _{c}}{M} [1 - (\frac{p _{e}}{p _{c}})^{\frac{γ - 1}{γ}}]$

Read it term by term:

$\frac{2 γ}{γ - 1}$ — a pure number from the gas type.
$\frac{R}{M}$ — this is the gas's specific gas constant (units $J/(kg\cdotK)$ ): the universal constant $R = 8.314 J/(mol\cdotK)$ divided by the molar mass $M$ turns "per mole" into "per kilogram." So $\frac{R T _{c}}{M}$ is the energy per kilogram scale: hotter chamber ( $T_{c}$ ↑) or lighter molecules ( $M$ ↓, defined in Step 0) means more energy per kilogram to convert.
The bracket $[1 - (\frac{p _{e}}{p _{c}})^{(γ - 1) / γ}]$ — the fraction of energy actually cashed in. Expand more (drop $p_{e}$ far below $p_{c}$ ) → bracket → 1 → maximum speed.

Step 4 — Assemble the thrust from a control volume

WHAT. Combine "mass per second" (Step 1) with "speed" (Step 3) into the force $F$ .

WHY. Force is the rate of momentum leaving. Each second, $\overset{m}{˙}$ kilograms depart at speed $v_{e}$ , carrying momentum $\overset{m}{˙} v_{e}$ backward — the rocket feels the equal-and-opposite forward kick.

PICTURE. The figure draws a dashed box (control volume) around the motor. Gas momentum flows out the right; a leftover pressure mismatch at the exit adds a second, smaller push.

$F = momentum thrust \overset{m}{˙} v_{e} + pressure thrust (p_{e} - p_{a}) A_{e}$

$\overset{m}{˙} v_{e}$ — the main event: mass flow $\times$ exhaust speed.
$p_{e}$ — exit pressure (Step 0); $p_{a}$ — outside air pressure; $A_{e}$ — exit area.
$(p_{e} - p_{a}) A_{e}$ — a bonus (or penalty) if the jet leaves at a different pressure than the surroundings.

Step 5 — Define $I_{s p}$ and watch the grain cancel

WHAT. Form specific impulse — "kick per kilogram of weight burned."

WHY. We want a fairness number that compares a tiny motor and a giant motor made of the same stuff. Dividing thrust by mass flow removes the "how big" and keeps the "how good."

PICTURE. The figure shows two motors — one huge $A_{b}$ , one small — spitting gas at the same speed. Big one makes more gas (taller arrow bundle), but each kilogram is equally good.

Define the effective exhaust velocity $c = F / \overset{m}{˙}$ , then

$I_{s p} = \frac{c}{g _{0}} = \frac{F}{m ˙ g _{0}} = \frac{1}{g _{0}} [v_{e} + \frac{( p _{e} - p _{a} ) A _{e}}{m ˙}]$

$g_{0} = 9.81 m/s^{2}$ — a bookkeeping constant, so the answer comes out in seconds. Not local gravity.
Notice $\overset{m}{˙}$ in the denominator: whatever $A_{b}$ boosted in the numerator's $F$ is divided right back out.

For a matched nozzle ( $p_{e} = p_{a}$ ) the pressure term dies:

$I_{s p} = \frac{v _{e}}{g _{0}} = \frac{1}{g _{0}} \frac{2 γ}{γ - 1} \frac{R T _{c}}{M} [1 - (\frac{p _{e}}{p _{c}})^{\frac{γ - 1}{γ}}]$

Every symbol is chemistry or nozzle ( $γ, T_{c}, M, p_{e} / p_{c}$ ). The grain quantities $ρ_{p}, A_{b}, a, n$ are gone. This is the same conclusion as the Tsiolkovsky Rocket Equation view: efficiency lives in exhaust speed, not in how much mass you carry.

Step 6 — Where grain DOES matter: chamber pressure

WHAT. Grain geometry does control one thing that circles back: the equilibrium chamber pressure $p_{c}$ .

WHY. Gas made by the grain must equal gas vented through the throat, or the pressure changes. Setting made = vented pins $p_{c}$ .

Deriving the "gas vented" side. The narrow throat can only pass gas so fast. Choked-flow theory (from Characteristic Velocity c-star) says the mass leaving through a throat of area $A_{t}$ is set by the chamber pressure and a chemistry-only speed constant $c^{*}$ :

$c^{*} \equiv \frac{p _{c} A _{t}}{m ˙ _{vent}} ⟹ \overset{m}{˙}_{vent} = \frac{p _{c} A _{t}}{c ^{*}}$

$c^{*}$ — the characteristic velocity, defined as chamber pressure times throat area divided by the mass flow it lets through. It measures "how much pressure the chamber builds per unit of vented mass," and depends only on $T_{c}$ , $M$ , $γ$ .
Bigger $p_{c}$ or bigger throat $A_{t}$ → more gas vented. Makes sense.

Set made = vented. Gas made is $\overset{m}{˙} = ρ_{p} A_{b} r$ (Step 1); gas vented is $p_{c} A_{t} / c^{*}$ . In steady operation they are equal:

$ρ_{p} A_{b} r = \frac{p _{c} A _{t}}{c ^{*}}$

PICTURE. The figure is a balance scale: on the left "gas made" $= ρ_{p} A_{b} r$ ; on the right "gas vented" $= p_{c} A_{t} / c^{*}$ . They settle at one pressure.

Insert $r = a p_{c}^{n}$ (Step 2) and solve for $p_{c}$ :

$ρ_{p} A_{b} a p_{c}^{n} = \frac{p _{c} A _{t}}{c ^{*}} ⟹ p_{c}^{1 - n} = \frac{ρ _{p} a c ^{*} A _{b}}{A _{t}} ⟹ p_{c} = (\frac{ρ _{p} a c ^{*} A _{b}}{A _{t}})^{\frac{1}{1 - n}}$

Here grain terms ( $ρ_{p}, a, A_{b}$ ) finally appear — but only for pressure, which feeds thrust and burn time, not the $I_{s p}$ formula.

Step 7 — The degenerate case: $n \to 1$ blows up

WHAT. Examine what happens to $p_{c}$ as the exponent $n$ approaches (and passes) 1.

WHY. This is the edge case the parent flagged. We must show the reader exactly where stability dies.

PICTURE. The figure plots the exponent $\frac{1}{1 - n}$ against $n$ : flat near $n = 0.3$ , then shooting to infinity as $n \to 1$ . A vertical amber asymptote marks the cliff.

$p_{c} \propto (\dots)^{\frac{1}{1 - n}}$

At $n = 0.35$ : exponent $= 1/0.65 = 1.54$ — mild, self-correcting.
As $n \to 1^{-}$ : exponent $\to + \infty$ — the tiniest change in the base explodes $p_{c}$ .
At $n \geq 1$ : no stable equilibrium — a pressure bump feeds a bigger bump. Thermal runaway.

The one-picture summary

This single diagram compresses the whole walkthrough: the grain feeds mass flow (left, geometry lives here), the chemistry+nozzle set exhaust speed (right), and $I_{s p} = F / (\overset{m}{˙} g_{0})$ divides the geometry back out — leaving efficiency as a pure-chemistry number.

Recall Feynman retelling — the whole page in plain words

Picture a shaped candle burning inside a can. The wall of flame (its area) decides how much smoke pours off each second — that's Step 1. Squeeze the can harder and the flame eats faster, but only a little faster, which is what keeps it from blowing up — Steps 2 and 7. All that smoke rushes out the back nozzle, trading its heat for speed: hotter and lighter smoke leaves faster — Step 3. The forward kick is smoke-per-second times smoke-speed — Step 4. Now here's the magic trick: to judge quality we divide the kick by the smoke-per-second, and the "how big is the flame" part cancels itself out — Step 5. So a giant candle and a tiny candle of the same wax shoot smoke at the same speed and have the same quality score; the giant just makes more smoke. The flame size only comes back to set the pressure inside — Step 6 — never the quality. Quality is about smoke speed, and smoke speed is chemistry.

Quick self-check

\overset{m}{˙}

formula and what each factor means

\overset{m}{˙} = ρ_{p} A_{b} r

: density × burning area × burn rate = kg/s produced.

Why does

A_{b}

vanish from

I_{s p}

A_{b}

multiplies

F

and

\overset{m}{˙}

equally;

I_{s p} = F / (\overset{m}{˙} g_{0})

cancels it.

What is the one thing grain geometry sets that circles back?

Equilibrium chamber pressure

p_{c} = (ρ_{p} a c^{*} A_{b} / A_{t})^{1/ (1 - n)}

Why must

n < 1

The exponent

1/ (1 - n)

diverges as

n \to 1

, giving thermal runaway.