This page is the exhaustive drill-book for the thrust-coefficient derivation . We take the master formula and fire it through every case class: perfectly matched nozzles, over-expanded ones, under-expanded ones, the vacuum limit, the infinite-expansion ceiling, degenerate zero-pressure inputs, a word problem, and an exam twist that hides a trap.
Before any numbers, meet the one formula every example below leans on.
Intuition Read the pressure term like a see-saw
The pressure term carries a sign . Picture a balance beam with P e on one pan and P a on the other (see the sketch below). Look at ( P e − P a ) :
P e = P a → beam level → term is zero → perfectly matched (see Over- and Under-Expanded Nozzles ).
P e > P a → exit pan sinks → term is positive → under-expanded, exit still pushing → bonus thrust.
P e < P a → ambient pan sinks → term is negative → over-expanded, atmosphere pushes back → thrust penalty.
Every "quadrant" of this topic is really a sign of this one bracket. The figure shows the nozzle cross-section and the three see-saw states side by side.
In the figure above, the left panel is the nozzle: the throat A ∗ is the pinch-point (violet), the exit A e is the wide mouth (orange). The right panel stacks the three see-saw states — read the sign of ( P e − P a ) straight off which pan is lower.
Here is every distinct case class this topic can throw at you. Each worked example below is tagged with the cell it fills.
#
Case class
Distinguishing condition
Filled by
A
Perfectly matched
P e = P a → pressure term = 0
Ex 1
B
Under-expanded
P e > P a → pressure term > 0
Ex 2
C
Over-expanded
P e < P a → pressure term < 0
Ex 3
D
Vacuum limit
P a = 0 → pressure term always > 0
Ex 4
E
Infinite-expansion ceiling
P e → 0 → C F → C F , m a x
Ex 5
F
Back-out from data
given F , find C F
Ex 6
G
Same engine, two altitudes
sign of pressure term flips
Ex 7
H
Word problem (design)
choose A e / A ∗ for a target
Ex 8
I
Exam twist (hidden trap)
student uses T c — it must vanish
Ex 9
Throughout we use γ = 1.2 (typical hot rocket gas) unless stated. Let us first nail the two building-block numbers we will reuse.
Definition The two reusable constants (for
γ = 1.2 )
K ≡ γ − 1 2 γ 2 ( γ + 1 2 ) γ − 1 γ + 1 , n ≡ γ γ − 1
For γ = 1.2 : γ − 1 2 γ 2 = 0.2 2 ( 1.44 ) = 14.4 , and ( 2.2 2 ) 2.2/0.2 = ( 0.9091 ) 11 = 0.3515 , so K = 14.4 × 0.3515 = 5.06 . And n = 0.2/1.2 = 0.1667 .
Then C F , mom = K [ 1 − ( P e / P c ) n ] .
Worked example Ex 1 — matched sea-level nozzle
γ = 1.2 , P c = 7 MPa, P e = P a = 0.1 MPa, A e / A ∗ = 25 . Find C F and thrust with A ∗ = 0.01 m².
Forecast: the pressure term is ( P e − P a ) = 0 , so C F should equal the momentum term alone — guess it lands near 1.5 –1.6 .
Pressure ratio power. ( P c P e ) n = ( 0.1/7 ) 0.1667 = ( 0.01429 ) 0.1667 = 0.4933 .
Why this step? This bracket measures how far the gas expanded; smaller P e / P c means more expansion.
Momentum term. C F , mom = 5.06 ( 1 − 0.4933 ) = 5.06 × 0.5067 = 2.564 = 1.601 .
Why this step? This is the directed-jet contribution set only by γ and P e / P c .
Pressure term. P c P e − P a A ∗ A e = 7 0 × 25 = 0 .
Why this step? Matched nozzle → the exit gas and atmosphere are in balance, nothing extra to add.
Total. C F = 1.601 + 0 = 1.601 . Thrust F = 1.601 × 7 × 1 0 6 × 0.01 = 112 , 070 N ≈ 112 kN.
Verify: units — C F is dimensionless; P c A ∗ = Pa·m² = N. 1.601 × 7 MPa × 0.01 m² gives ∼ 1.1 × 1 0 5 N ✓. Sanity: C F sits inside the typical 1.3 –2.0 band ✓.
Worked example Ex 2 — high-altitude nozzle, exit still over-pressured
Same engine: γ = 1.2 , P c = 7 MPa, P e = 0.1 MPa, A e / A ∗ = 25 , but now flying where P a = 0.04 MPa.
Forecast: P e > P a , so pressure term is positive → C F should beat the matched value 1.601 .
Momentum term unchanged. P e / P c is the same as Ex 1, so C F , mom = 1.601 .
Why this step? The momentum term never sees P a — only the exit-to-chamber ratio.
Pressure term. 7 0.1 − 0.04 × 25 = 7 0.06 × 25 = 0.2143 .
Why this step? The exit gas is still pushing harder than ambient over area A e → extra forward force.
Total. C F = 1.601 + 0.214 = 1.815 .
Verify: the pressure term is positive as predicted, and 1.815 > 1.601 ✓. Because P a dropped from 0.1 to 0.04 MPa, we gained thrust with the identical nozzle — exactly what the Thrust Equation promises for rising altitude.
Worked example Ex 3 — same nozzle forced down to sea level, but over-expanded
A big vacuum-style nozzle: γ = 1.2 , P c = 6 MPa, P e = 0.05 MPa, A e / A ∗ = 30 , tested at sea level P a = 0.1 MPa.
Forecast: P e < P a — atmosphere pushes back — so the pressure term goes negative and drags C F below its momentum term.
Momentum term. ( 0.05/6 ) 0.1667 = ( 0.008333 ) 0.1667 = 0.4500 . C F , mom = 5.06 ( 1 − 0.4500 ) = 5.06 × 0.5500 = 2.783 = 1.668 .
Why this step? Deeper expansion (smaller P e / P c ) than Ex 1, so a slightly bigger momentum term.
Pressure term. 6 0.05 − 0.1 × 30 = 6 − 0.05 × 30 = − 0.25 .
Why this step? Negative bracket = ambient pressure wins over the weak exit pressure across the big exit area → a real thrust loss . See Over- and Under-Expanded Nozzles .
Total. C F = 1.668 − 0.25 = 1.418 .
Verify: sign of pressure term is negative, and 1.418 < 1.668 ✓. This is the classic over-expansion penalty — the giant vacuum bell is wasteful at sea level.
Worked example Ex 4 — the same big nozzle now in space
Take the Ex 3 engine into vacuum: γ = 1.2 , P c = 6 MPa, P e = 0.05 MPa, A e / A ∗ = 30 , P a = 0 .
Forecast: with P a = 0 the pressure term is P c P e A ∗ A e , always positive — this should be the best C F the nozzle can reach at this P e .
Momentum term. Recompute (same P e / P c as Ex 3): ( 0.05/6 ) 0.1667 = 0.4500 , so C F , mom = 5.06 ( 1 − 0.4500 ) = 1.668 .
Why this step? Momentum term ignores P a ; only P e / P c matters, so it matches Ex 3.
Pressure term. 6 0.05 − 0 × 30 = 6 0.05 × 30 = 0.25 .
Why this step? No atmosphere pushes back, so the entire exit pressure over A e becomes pure bonus thrust.
Total. C F , vac = 1.668 + 0.25 = 1.918 .
Verify: compare with Ex 3 — the same nozzle went from 1.418 (sea level) to 1.918 (vacuum), a jump of exactly 0.50 , which equals P c P a A ∗ A e = 6 0.1 × 30 = 0.50 ✓. That is why space stages carry huge expansion ratios .
Worked example Ex 5 — the hard maximum set by
γ
Push expansion to the theoretical limit: P e → 0 , still γ = 1.2 . Find C F , m a x .
Forecast: as P e → 0 , the bracket [ 1 − ( P e / P c ) n ] → 1 , so C F , mom hits its ceiling K . Guess 5.06 ≈ 2.25 .
Bracket limit. ( P e / P c ) n → 0 0.1667 = 0 , so [ 1 − 0 ] = 1 .
Why this step? Infinite expansion extracts all available enthalpy — nothing left in the pressure ratio.
Ceiling value. C F , m a x = K = 5.06 = 2.249 .
Why this step? This depends only on γ — a fundamental limit no nozzle can beat, tied to Isentropic Nozzle Flow .
Interpret. Real engines with P a = 0 approach but never reach 2.25 because a finite nozzle leaves P e > 0 .
Verify: plug the Ex-4 numbers — its C F , vac = 1.918 is safely below 2.249 ✓. Take a bigger A e / A ∗ (smaller P e ) and C F creeps toward 2.25 but never crosses it, exactly as the ceiling demands.
The figure below shows how C F climbs with expansion ratio in vacuum, flattening toward this ceiling.
In that plot the horizontal axis is the expansion ratio A e / A ∗ and the vertical axis is C F . The magenta curve is the full vacuum C F , the dashed violet curve is the momentum term alone, and the dotted orange line is the ceiling K = 2.25 — notice the magenta curve bends toward it but never touches.
Worked example Ex 6 — grade an engine from a test-stand reading
A sea-level test reads F = 90 kN, P c = 5 MPa, A ∗ = 0.008 m². Find C F .
Forecast: just invert the definition C F = F / ( P c A ∗ ) ; no γ needed. Guess "a bit above 2 ".
Reference force. P c A ∗ = 5 × 1 0 6 × 0.008 = 40 , 000 N.
Why this step? This is the denominator of the coefficient — the "force yardstick".
Divide. C F = 40 , 000 90 , 000 = 2.25 .
Why this step? C F converts a raw thrust reading into a size-independent grade for the nozzle. See Characteristic Velocity c-star for the companion combustion grade.
Verify: units — N/(Pa·m²) = N/N = dimensionless ✓. A value of 2.25 is suspiciously high for sea level — it flags a very well-expanded, mildly under-pressured engine (or a measurement to double-check).
Worked example Ex 7 — watch the pressure term flip sign
Engine: γ = 1.2 , P c = 6 MPa, P e = 0.05 MPa fixed, A e / A ∗ = 30 . Compute C F at sea level and in vacuum, and the change.
Forecast: the momentum term is identical at both altitudes, so the whole difference is the pressure term. From sea level (P a = 0.1 MPa) to vacuum (P a = 0 ), P a drops by 0.1 MPa; the swing should equal P c Δ P a A ∗ A e = 6 0.1 × 30 = 0.50 . Guess a rise of exactly 0.50 .
Momentum term. Recompute: ( 0.05/6 ) 0.1667 = 0.4500 , so C F , mom = 5.06 ( 1 − 0.4500 ) = 1.668 .
Why this step? Same P e / P c as Ex 3/Ex 4 → same momentum term; it is altitude-blind.
Sea level (P a = 0.1 MPa). Pressure term = 6 0.05 − 0.1 × 30 = − 0.25 → C F = 1.668 − 0.25 = 1.418 .
Why this step? Over-expanded on the ground — atmosphere pushes back.
Vacuum (P a = 0 ). Pressure term = 6 0.05 × 30 = + 0.25 → C F = 1.668 + 0.25 = 1.918 .
Why this step? No back-pressure — full exit pressure becomes thrust.
Swing. Δ C F = 1.918 − 1.418 = 0.50 .
Why this step? The momentum term cancels, so the entire difference is P c Δ P a A ∗ A e .
Verify: 6 0.1 × 30 = 0.50 matches the swing exactly ✓. This is the same engine as Ex 3/Ex 4, so the numbers cross-check.
Worked example Ex 8 — size a vacuum nozzle to a thrust target
Mission spec: deliver F = 150 kN in vacuum (P a = 0 ) with P c = 7 MPa, A ∗ = 0.01 m². Ground rules give γ = 1.2 and a matched-at-P e = 0.1 MPa design so C F , mom = 1.601 (from Ex 1). What expansion ratio A e / A ∗ do we need?
Forecast: we need C F = F / ( P c A ∗ ) ; then subtract the momentum term to isolate what the pressure term (hence area ratio) must supply.
Required C F . C F = 7 × 1 0 6 × 0.01 150 , 000 = 70 , 000 150 , 000 = 2.143 .
Why this step? The definition sets the total coefficient we must hit.
Required pressure term. C F − C F , mom = 2.143 − 1.601 = 0.542 .
Why this step? Whatever the momentum term can't provide, the vacuum pressure term must.
Solve for area ratio. Vacuum: P c P e A ∗ A e = 0.542 , so A ∗ A e = 0.1/7 0.542 = 0.542 × 70 = 37.9 .
Why this step? P e / P c = 0.1/7 is fixed by the matched design, so only A e / A ∗ is free.
Verify: rebuild C F = 1.601 + 7 0.1 × 37.9 = 1.601 + 0.542 = 2.143 ✓, and F = 2.143 × 70 , 000 = 150 , 000 N ✓. An expansion ratio near 38 is realistic for an upper stage — see Specific Impulse Isp for why bigger ratios also raise I s p in vacuum.
Worked example Ex 9 — the deliberately useless temperature
Exam gives: γ = 1.2 , P c = 7 MPa, P e = P a = 0.1 MPa, A e / A ∗ = 25 , and T c = 3200 K, R = 350 J/(kg·K). Here T c is the chamber temperature and R is the specific gas constant — the gas constant per kilogram of the exhaust, with units J/(kg·K). Find C F . (Same physical case as Ex 1, but stuffed with extra numbers.)
Forecast: the trap is that T c and R tempt you to compute exit velocity. But C F is independent of T c and R — they cancel in the derivation. Answer should equal Ex 1's 1.601 .
Spot the cancellation. In the parent's Step 4, m ˙ ∝ R T c 1 and u e ∝ R T c , so m ˙ u e / ( P c A ∗ ) has R T c cancel.
Why this step? C F measures nozzle geometry and ratios only — never absolute temperature or the specific gas constant.
Ignore T c , R ; use the master formula. Identical to Ex 1: ( 0.1/7 ) 0.1667 = 0.4933 , so C F , mom = 5.06 ( 1 − 0.4933 ) = 1.601 , and the pressure term = 0 (matched).
Why this step? Same γ , P e / P c , A e / A ∗ , and matched exit → same coefficient.
Total. C F = 1.601 + 0 = 1.601 .
Why this step? The extra data T c and R never entered — proving they are decoration for a C F question.
Verify: unchanged from Ex 1 despite the added T c = 3200 K and R = 350 ✓. If you had wrongly "used" them you would compute an exit velocity u e ≈ γ − 1 2 γ R T c [ 1 − ( P e / P c ) n ] = 0.2 2 ( 1.2 ) ( 350 ) ( 3200 ) ( 0.5067 ) ≈ 2860 m/s — a real, meaningful number, but it belongs to c ∗ and I s p , not C F .
Common mistake The single trap behind Ex 9
C F never depends on T c or R . Those set the exhaust speed and mass flow (hence c ∗ ), but they cancel out of the coefficient. If a problem hands you T c or R for a pure C F question, they are decoration — plug only γ , P e / P c , P a / P c , and A e / A ∗ into the master formula.
Recall Quick self-test
Matched nozzle pressure term equals what? ::: Zero (since P e = P a ).
Over-expanded means which inequality, and what sign of pressure term? ::: P e < P a , negative term (thrust penalty).
Going from sea level to vacuum, the C F gain equals? ::: P c Δ P a A ∗ A e .
The infinite-expansion ceiling C F , m a x depends only on? ::: γ (for γ = 1.2 it is ≈ 2.25 ).
Does C F depend on chamber temperature T c or specific gas constant R ? ::: No — both cancel; they live in m ˙ and c ∗ .