This page is the worked-example workshop for Characteristic velocity $c^*$ . We do not re-derive the theory here — we use it, and we deliberately walk into every corner: normal numbers, zero inputs, limiting values of γ , a word problem, and an exam trap. If you have not seen the boxed formula yet, read the parent first.
Definition The chamber symbols
P c , T c , and Mach number M
P c = combustion chamber pressure (the "stagnation" pressure, where the gas is nearly at rest before it rushes out). Measured in pascals (Pa).
T c = combustion chamber temperature (the flame temperature), in kelvin (K).
M = Mach number = the gas's local speed divided by the local speed of sound. M = 1 means the gas is moving exactly at the speed of sound — this is the "choked" condition that happens at the throat (see Choked Flow and the de Laval Nozzle ).
γ — the ratio of specific heats
γ (Greek "gamma") is the ratio of a gas's specific heat at constant pressure to that at constant volume , γ = c p / c v . Physically it tells you how "springy" the gas is when squeezed. For rocket exhaust it sits between about 1.1 (heavy, many-atom molecules) and 1.67 (monatomic gases like helium). It is a pure number (no units) and it is the only gas property inside the Vandenkerckhove function Γ .
Everything below rests on the two faces of c ∗ :
Intuition Why the two theoretical forms are the same
They look different but are algebraically identical. Start from the right form Γ 1 γ R T c and pull Γ 1 apart. Flipping the fraction inside the Γ bracket flips the sign of its exponent, turning γ + 1 2 into 2 γ + 1 :
Γ 1 = γ 1 ( 2 γ + 1 ) 2 ( γ − 1 ) γ + 1 .
Now multiply by γ R T c :
Γ 1 γ R T c = γ 1 ( 2 γ + 1 ) 2 ( γ − 1 ) γ + 1 ⋅ γ R T c = ( 2 γ + 1 ) 2 ( γ − 1 ) γ + 1 R T c .
The two γ factors cancel cleanly (one from 1/Γ , one from inside γ R T c ), leaving ( 2 γ + 1 ) 2 ( γ − 1 ) γ + 1 R T c . To recover the left form γ γ R T c ( 2 γ + 1 ) ⋯ , note γ γ R T c = γ γ R T c = γ R T c — so the left form equals γ 1 ( 2 γ + 1 ) ⋯ R T c , i.e. the left form is γ 1 times the right-hand expression above. In other words the truly-equal pair is
c ideal ∗ = Γ 1 γ R T c = γ γ R T c ( 2 γ + 1 ) 2 ( γ − 1 ) γ + 1 .
No stray γ appears or disappears anywhere — every one is accounted for. Left form is handy when you already know Γ ; right form is handy when you don't.
Two symbols you must never mix up:
A ∗ = throat area (the pinch, where the flow is choked ). This is the "star" area.
m ˙ = mass flow rate in kg/s (the "m-dot" — the little dot means "per second").
Every question this topic can throw at you falls into one of these cells. Each example below is tagged with the cell it fills, and every cell carries its own unit sanity-check inside the worked example.
Cell
What makes it tricky
Example
A. Plain measured
just plug into P c A ∗ / m ˙
Ex 1
B. Plain theoretical
build Γ first, then γ R T c
Ex 2
C. Efficiency
ratio measured / ideal
Ex 3
D. Solve backwards
given c ∗ , find m ˙ or P c
Ex 4
E. Degenerate: m ˙ → 0
limiting behaviour, why it blows up
Ex 5
F. Limiting γ
γ → 1 + vs γ ≈ 1.67 — sanity of Γ
Ex 6
G. Design mantra: change M
swap propellant, why H 2 wins
Ex 7
H. Word problem
real test stand, unit soup
Ex 8
I. Exam trap
someone hands you A e or a nozzle detail
Ex 9
We hit every cell A–I. Units matter as much as the arithmetic — most exam blood is lost to units.
Worked example Example 1 — Cell A: plain measured
c ∗
Test stand reads P c = 6.0 MPa , A ∗ = 0.008 m 2 , m ˙ = 20 kg/s . Find c ∗ .
Forecast: Guess the order of magnitude — chamber pressures give c ∗ around 1500 –2500 m/s. Where will this land?
Write the measured definition: c ∗ = m ˙ P c A ∗ .
Why this step? No nozzle info was given, and none is needed — that is the whole point of c ∗ being a chamber quantity.
Convert MPa to Pa: 6.0 MPa = 6.0 × 1 0 6 Pa .
Why this step? SI units only, or the m/s at the end is meaningless.
Plug in: c ∗ = 20 ( 6.0 × 1 0 6 ) ( 0.008 ) = 20 48000 = 2400 m/s .
Why this step? Pure arithmetic once units are clean.
Verify: Units: kg/s Pa ⋅ m 2 = kg/s N/m 2 ⋅ m 2 = kg/s N = kg/s kg⋅m/s 2 = m/s . ✓ Lands in the expected band.
Worked example Example 2 — Cell B: theoretical
c ∗ from chemistry
γ = 1.20 , T c = 3200 K , M = 22 g/mol . Find c ideal ∗ .
Forecast: Theoretical c ∗ is usually a bit higher than what a real stand measures (real engines lose a few percent). Expect ~1500 –1600 m/s.
Gas constant: R = R u / M . Fix the units first: R u = 8314 J/(kmol⋅K) and M = 22 g/mol = 22 kg/kmol , so R = 8314/22 = 378 J/(kg⋅K) .
Why this step? The formula needs the specific gas constant R (per kg). Matching R u in J/(kmol·K) with M in kg/kmol keeps the prefixes consistent.
Exponent: 2 ( γ − 1 ) γ + 1 = 0.4 2.2 = 5.5 .
Why this step? This single power packages the density ratio (γ − 1 1 ) plus the sound-speed from the isentropic relations .
Build Γ = 1.2 ( 2.2 2 ) 5.5 = 1.0954 × ( 0.9091 ) 5.5 = 1.0954 × 0.5875 = 0.6436 .
Why this step? Γ collects all the γ -dependence into one factor, so the rest of the formula is just thermochemistry.
Thermochemical part: γ R T c = 1.2 × 378 × 3200 = 1.452 × 1 0 6 = 1205 m/s .
Why this step? This is the "burn hot, make light gas" core — it carries γ R T c so we can divide by Γ directly.
Divide: c ideal ∗ = Γ γ R T c = 0.6436 1205 = 1872 m/s .
Why this step? The clean identity c ideal ∗ = Γ 1 γ R T c assembles the two pieces in one division.
Verify (the other form, step by step): Use the left form c ideal ∗ = γ γ R T c ( 2 γ + 1 ) 5.5 .
γ γ R T c = 1.2 1205 = 1004 .
( 2 γ + 1 ) 5.5 = ( 1.1 ) 5.5 = 1.8642 .
Product: 1004 × 1.8642 = 1872 m/s . ✓
Both forms give 1872 m/s exactly (to rounding). Note this is a touch above the 1500 –1600 forecast because γ = 1.2 with this R is a slightly higher-energy case than the mental estimate — the forecast was a rough guess, and the math is the arbiter.
Worked example Example 3 — Cell C: combustion efficiency
Same propellant as Ex 2 (c ideal ∗ = 1872 m/s). The stand actually measures c ∗ = 1780 m/s. Find η c ∗ .
Forecast: Real engines run 0.92 –0.99 . Guess where 1780/1872 lands.
η c ∗ = c ideal ∗ c measured ∗ = 1872 1780 .
Why this step? Efficiency is always real-over-ideal; both numbers ignore the nozzle, so this cleanly scores the combustion.
= 0.9509 ≈ 95% .
Why this step? The missing 5% is combustion potential lost to incomplete burning, wall heat loss, or imperfect mixing — not a nozzle fault.
Verify: 0.9509 × 1872 = 1780 . ✓ Dimensionless (m/s over m/s). Sits inside the healthy 0.92 –0.99 band.
Worked example Example 4 — Cell D: solve backwards for
m ˙
A design targets c ∗ = 1700 m/s with P c = 4.0 MPa and A ∗ = 0.015 m 2 . What mass flow must the pumps deliver?
Forecast: Bigger throat + high pressure means a lot of propellant per second — tens of kg/s.
Rearrange the measured formula: m ˙ = c ∗ P c A ∗ .
Why this step? c ∗ is a fixed ratio for the chosen propellant, so it sets the flow needed to hold that pressure through that throat.
Plug in: m ˙ = 1700 ( 4.0 × 1 0 6 ) ( 0.015 ) = 1700 60000 = 35.3 kg/s .
Why this step? Straight algebra; this is how pump sizing actually starts.
Verify: Feed it back: 35.29 ( 4.0 × 1 0 6 ) ( 0.015 ) = 1700 m/s. ✓ Units: m/s Pa⋅m 2 = m/s N = kg/s . ✓
Worked example Example 5 — Cell E: degenerate input
m ˙ → 0
What does c ∗ = P c A ∗ / m ˙ do as you throttle the flow toward zero (with P c , A ∗ held)? Is that physical?
Forecast: The formula divides by m ˙ . Guess what dividing by "almost nothing" does.
Take the limit: as m ˙ → 0 + , c ∗ = m ˙ P c A ∗ → + ∞ .
Why this step? This is the mathematical behaviour of a 1/ m ˙ curve — it explodes.
Ask if the premise survives: to hold a fixed P c with no propellant flowing is impossible. Real chambers need m ˙ and the choking condition (M = 1 , i.e. gas moving at the speed of sound, at the throat — see the M definition above) for the tidy formula to hold at all.
Why this step? Cell E teaches that a formula's math and its physical validity are different questions — the limit is real math but forbidden physics.
The honest physical limit: for a fixed propellant, c ∗ is a constant (∝ T c / M stuff). If m ˙ genuinely drops, P c drops in lockstep so the ratio stays put. You cannot move one without the other while choked.
Verify: Numerically with the Ex-1 numbers, halving m ˙ from 20 to 10 while forcing P c fixed doubles the apparent c ∗ from 2400 to 4800 m/s — physically flagging that you have violated the choked coupling. ✓ The figure below makes this concrete.
Reading the figure: the horizontal axis is mass flow m ˙ (kg/s), the vertical axis is the apparent c ∗ (km/s), both explicitly labelled on the plot. The coral curve is c ∗ = P c A ∗ / m ˙ with P c artificially held fixed. Follow it leftward: as m ˙ shrinks into the buttery-shaded danger zone it rockets toward the sky — that vertical wall is the m ˙ → 0 blow-up of step 1, annotated on the figure. The mint dot (also labelled) is the only honest reading: the real operating point where P c and m ˙ are locked together by choking, giving the sensible 2400 m/s of Ex 1. The picture's lesson: the blow-up is an artefact of pretending P c stays put — nature never lets you sit on the left of that dot.
Worked example Example 6 — Cell F: limiting values of
γ
Probe the theoretical formula across the whole physical range of γ : the soft limit γ → 1 + (heavy, many-atom exhaust) and the hard limit γ ≈ 1.67 (monatomic gas like helium). Hold R T c = 1.21 × 1 0 6 m 2 / s 2 (the Ex-2 value) and use the left form c ideal ∗ = γ γ R T c ( 2 γ + 1 ) 2 ( γ − 1 ) γ + 1 = γ R T c ( 2 γ + 1 ) 2 ( γ − 1 ) γ + 1 .
Forecast: Which γ gives a bigger c ∗ — the soft low-γ gas or the stiff monatomic one? Guess before computing.
We evaluate c ideal ∗ = Γ 1 γ R T c = γ R T c ( 2 γ + 1 ) 2 ( γ − 1 ) γ + 1 at four points: γ = 1.10 , γ = 1.13 , γ = 1.30 , and the monatomic extreme γ = 1.67 . Here R T c = 1.21 × 1 0 6 = 1100 .
Why this step? Four points bracket the entire physical range and let us actually see how Γ behaves.
At γ = 1.10 : exponent = 0.2 2.1 = 10.5 ; Γ = 1.1 ( 2/2.1 ) 10.5 = 1.0488 × 0.6035 = 0.6329 ; then c ∗ = Γ γ R T c = 0.6329 1.1 ⋅ 1100 = 0.6329 1153.7 = 1823 m/s .
Why this step? The "soft gas" edge, computed explicitly.
At γ = 1.13 : exponent = 0.26 2.13 = 8.192 ; Γ = 1.13 ( 2/2.13 ) 8.192 = 1.0630 × 0.6062 = 0.6444 . At γ = 1.30 : exponent = 0.6 2.3 = 3.833 ; Γ = 1.3 ( 2/2.3 ) 3.833 = 1.1402 × 0.5254 = 0.5991 ; c ∗ = 0.5991 1.3 ⋅ 1100 = 0.5991 1254.2 = 2094 m/s .
Why this step? Γ rises from 0.6329 (γ = 1.10 ) to 0.6444 (γ = 1.13 ) then falls to 0.5991 (γ = 1.30 ) — a genuine hump. Note that c ∗ itself keeps rising because the γ factor grows; the hump is in Γ , not in c ∗ .
Monatomic extreme γ = 1.67 : exponent = 1.34 2.67 = 1.993 ; Γ = 1.67 ( 2/2.67 ) 1.993 = 1.2923 × 0.5619 = 0.7262 ; c ∗ = 0.7262 1.67 ⋅ 1100 = 0.7262 1421.5 = 1958 m/s .
Why this step? Closes the high-γ end: the stiff monatomic gas gives a lower c ∗ than the γ = 1.30 peak but still above the softest gas — the interplay of γ (rising) and the Γ -bracket (humped) produces a broad maximum near γ ≈ 1.3 .
Sanity of Γ itself across the whole sweep: Γ stays bounded in ≈ [ 0.60 , 0.73 ] and never touches 0 or blows up for any γ ∈ ( 1 , 1.67 ] .
Why this step? A formula is only trustworthy if it is well-behaved on its whole domain — checking that Γ has no poles or zeros confirms c ideal ∗ is finite and positive for every real propellant.
Verify: The four Γ values 0.6329 , 0.6444 , 0.5991 , 0.7262 (all checked in the appendix) confirm the hump near γ ≈ 1.13 in Γ and the bounded range. Resulting c ∗ : 1823 , 2094 , 1958 m/s for γ = 1.10 , 1.30 , 1.67 — a broad peak, not a monotone slide. ✓
Worked example Example 7 — Cell G: swap the propellant (design mantra)
Keep T c = 3200 K , γ = 1.20 but switch exhaust from M = 22 g/mol to a hydrogen-rich M = 10 g/mol. What happens to c ∗ ?
Forecast: The mantra says "make light gases." How much does more-than-halving M buy you?
c ideal ∗ ∝ γ R T c = γ R u T c / M , so c ∗ ∝ 1/ M at fixed T c , γ .
Why this step? Everything else (Γ , T c , γ ) is unchanged, so only the 1/ M factor moves.
Ratio: c 22 ∗ c 10 ∗ = 10 22 = 2.2 = 1.483 .
Why this step? Direct proportionality lets us skip re-computing Γ .
New value: c 10 ∗ = 1872 × 1.483 = 2776 m/s .
Why this step? A 48% jump in combustion score just by making the exhaust lighter — this is exactly why H 2 propellants dominate high-performance stages.
Verify: Recompute from scratch. R = 8314/10 = 831.4 ; γ R T c = 1.2 × 831.4 × 3200 = 3.193 × 1 0 6 = 1787 ; Γ = 0.6436 (unchanged); c ∗ = 1787/0.6436 = 2776 m/s . ✓ Matches the ratio method. Units of the ratio: dimensionless (g/mol over g/mol). ✓
Worked example Example 8 — Cell H: full word problem (real test stand)
An engineer logs a 12 -second burn. The chamber held P c = 3.5 MPa steady, throat diameter d ∗ = 90 mm , and the propellant tanks lost 216 kg of mass over the burn. The predicted (ideal) c ∗ for this propellant is 1720 m/s . Report c measured ∗ and the combustion efficiency.
Forecast: Efficiency should land in the healthy 0.9 –0.99 range if the engine is any good.
Mass flow: m ˙ = 12 s 216 kg = 18 kg/s .
Why this step? m ˙ is total mass burned over burn time — the tanks tell you this directly.
Throat area from diameter: A ∗ = π ( d ∗ /2 ) 2 = π ( 0.045 ) 2 = 6.362 × 1 0 − 3 m 2 .
Why this step? You were handed a diameter in mm — must convert to m and turn into area, or units collapse.
Measured c ∗ = m ˙ P c A ∗ = 18 ( 3.5 × 1 0 6 ) ( 6.362 × 1 0 − 3 ) = 18 22267 = 1237 m/s .
Why this step? The measured definition — every quantity here is chamber/throat, nothing from the nozzle.
Efficiency: η c ∗ = 1237/1720 = 0.719 = 72% .
Why this step? Real-over-ideal. This is low — a red flag.
Verify: 72% is far below the healthy 0.92 –0.99 band → this engine has serious incomplete combustion or a mis-measured throat/flow. The math checks (0.719 × 1720 = 1237 ✓), but the engineering verdict is "investigate." ✓ Units all resolve to m/s and dimensionless.
Worked example Example 9 — Cell I: the exam trap (wrong area)
Exam gives: P c = 5.0 MPa , throat area A ∗ = 0.010 m 2 , exit area A e = 0.080 m 2 , m ˙ = 25 kg/s . Compute c ∗ . What is the trap?
Forecast: Two areas are dangled. Which one belongs in c ∗ — and what wrong answer does the trap produce?
Identify the trap: c ∗ uses the throat area A ∗ , because the flow is choked there (M = 1 ). The exit area A e belongs to the nozzle and only affects the thrust coefficient $C_F$ .
Why this step? Grabbing the bigger "obvious" area is the single most common c ∗ error.
Correct: c ∗ = 25 ( 5.0 × 1 0 6 ) ( 0.010 ) = 2000 m/s .
Why this step? Same as parent Ex 1 — throat only.
The trap answer (using A e ): 25 ( 5.0 × 1 0 6 ) ( 0.080 ) = 16000 m/s — physically absurd (eight times too big; no chemistry gives that).
Why this step? Seeing the wrong number is what trains you to reject it on sight.
Verify: Correct c ∗ = 2000 m/s (matches parent). The trap gives 16000 m/s = 8 × too large, exactly the ratio A e / A ∗ = 8 . ✓ Sanity: any c ∗ above ~3000 m/s for chemical propellants should trigger suspicion.
Recall Which area, and why?
In c ∗ = P c A ∗ / m ˙ , why must the area be the throat, not the exit? ::: Because the flow is choked (sonic, M = 1 ) at the throat — that is where m ˙ locks to ( P c , T c , A ∗ ) . The exit area only governs C F (nozzle expansion).
Recall Degenerate
m ˙
As m ˙ → 0 the formula → ∞ — is that a real prediction? ::: No. Holding P c fixed while m ˙ → 0 violates the choked coupling; physically P c falls with m ˙ so c ∗ stays a constant set by T c , M , γ .
Recall Design mantra check
Halving exhaust molecular mass M multiplies c ∗ by what factor? ::: 2 ≈ 1.41 , since c ∗ ∝ 1/ M .
Recall Which
γ wins?
Across γ ∈ ( 1 , 1.67 ] at fixed R T c , is the peak c ∗ at the soft end, the stiff end, or in the middle? ::: In the middle — c ∗ rose to a broad maximum near γ ≈ 1.3 (2094 m/s) and was lower at both γ = 1.10 (1823 ) and γ = 1.67 (1958 ).
Mnemonic The whole page in one line
"Throat, not exit; hot and light; real over ideal." — the three things every c ∗ problem tests.