CEA (Chemical Equilibrium with Applications) — using NASA-CEA tool to compute Isp, Tc, products: Foundations

This page assumes you have seen none of the notation in the parent note. We build each symbol from a plain-words meaning, a picture, and a reason the topic needs it — in an order where each rung of the ladder rests on the one below.

1. Moles and mole numbers — $n_{i}$

Figure — CEA (Chemical Equilibrium with Applications) — using NASA-CEA tool to compute Isp, Tc, products

Look at the figure: each coloured bin holds a pile of one kind of molecule. The height of each pile is its $n_{i}$ . Combustion is the act of shuffling atoms between these bins — the total number of atoms of each element never changes, but which bins they sit in does.

Why the topic needs it: every output of CEA (temperature, product fractions, $I_{s p}$ ) depends on how much of each species ends up in the mix. The $n_{i}$ are the unknowns CEA solves for.

2. Mole fraction — $X_{i}$ and total moles $n_{tot}$

The strange Greek-looking symbol $\sum$ (capital sigma) just means "add up all the terms". $\sum_{i} n_{i}$ = "add up $n_{i}$ over every species $i$ ." Nothing more.

Why the topic needs it: CEA reports products as mole fractions $X_{i}$ because they are size-independent — they describe the composition whether your chamber is a thimble or a Saturn V.

3. Pressure $p$ and partial pressure $p_{i}$

In the figure the total wall-push (black arrows) splits into coloured contributions, one per species, in proportion to each $X_{i}$ .

Why the topic needs it: a molecule's chemical "eagerness to react or escape" depends on how crowded it is — i.e. on its partial pressure $p_{i}$ . This is the quantity that walks into the log term in the next section.

4. Temperature $T$ and the chamber temperature $T_{c}$

Why the topic needs it: $T_{c}$ is one of CEA's three headline outputs, and it also controls the exhaust speed ( $I_{s p} \propto T_{c} / M$ ). Hotter chamber, faster exhaust.

5. The reference state — $T^{\circ}$ , $p^{\circ}$ , and the little circle $^{\circ}$

Why the topic needs it: CEA's thermodynamic tables list formation enthalpies and potentials at $T^{\circ}, p^{\circ}$ , then correct to the real chamber conditions. Without the baseline, none of the numbers could be added together.

6. Enthalpy $h$ , formation enthalpy $Δ_{f} h^{\circ}$ , heat capacity $c_{p}$

Definition The pieces

$h_{i} (T) = Δ_{f} h_{i}^{\circ} + \int_{T^{\circ}}^{T} c_{p, i} d T$

$Δ_{f} h_{i}^{\circ}$ = standard enthalpy of formation: the energy locked in by building one mole of species $i$ from its raw elements, measured at the standard state ( $T^{\circ} = 298.15$ K, $p^{\circ} = 1$ bar) — that is what the $^{\circ}$ superscript flags. Units: $J mol^{- 1}$ (same as $h$ ). This is the "chemical" part — see Hess's Law and Enthalpy of Formation.
$c_{p, i}$ = heat capacity at constant pressure: how many joules it takes to warm one mole by one Kelvin. Units: $J mol^{- 1} K^{- 1}$ .
The $\int_{T^{\circ}}^{T} \dots d T$ (an integral) is a running total: "sum up $c_{p} d T$ as we heat from the reference temperature $T^{\circ} = 298.15$ K up to the actual temperature $T$ ." Note the units multiply correctly: $c_{p} [J mol^{- 1} K^{- 1}]$ times $d T [K]$ gives $J mol^{- 1}$ , matching $h$ . It is the "thermal" (sensible-heat) part, added on top of the standard-state chemical energy.

Why the tool "integral"? Because $c_{p}$ itself changes with temperature. You cannot just multiply "heat capacity times temperature rise" with a single number; you must accumulate the little contributions $c_{p} d T$ across the whole climb. The integral is the machine that does that accumulation — starting exactly at the reference temperature $T^{\circ}$ where $Δ_{f} h^{\circ}$ was measured.

7. Gibbs free energy $G$ and chemical potential $μ_{i}$

Let us earn every symbol in that last line:

$μ_{i}^{\circ} (T)$ = the standard chemical potential: the value $μ_{i}$ takes at the reference pressure $p^{\circ}$ , as a function of temperature. Units $J mol^{- 1}$ . CEA reads it from fitted polynomial tables.
$p^{\circ}$ = the reference pressure, fixed at $p^{\circ} = 1 bar$ (the standard state from §5). The ratio $p_{i} / p^{\circ}$ is dimensionless, which is required before you can take its logarithm.
$R$ = the universal gas constant, $8.314 J mol^{- 1} K^{- 1}$ — the conversion factor between temperature and energy per mole. So $R T$ has units $J mol^{- 1}$ , matching $μ_{i}$ .
$ln$ = the natural logarithm, the function that answers " $e$ to what power gives this number?" ( $e \approx 2.718$ ). See Chemical Equilibrium and Kp for why logs appear.

Why the topic needs it: minimizing $G = \sum_{i} n_{i} μ_{i}$ is the equilibrium calculation. The $ln$ term is what drives dissociation: splitting one crowded molecule into two dilute ones can lower $G$ even though it costs bond energy. See Dissociation at High Temperature.

8. Element-conservation: $a_{ij}$ , $b_{j}$ , and the constraints

Why the topic needs it: minimizing $G$ with no rules would send every $n_{i} \to 0$ (empty is lowest energy). The atom-conservation constraints are the walls that keep the minimization honest.

9. Lagrange multipliers $λ_{j}$

Intuition Where the master condition comes from (derivation sketch)

We want the smallest $G$ , but we may only wiggle the $n_{i}$ in ways that keep every atom count $b_{j}$ fixed. The trick of Lagrange multipliers is to build a combined function $L = G i \sum n_{i} μ_{i} - \sum_{j} λ_{j} (= 0 at any valid mix i \sum a_{ij} n_{i} - b_{j}),$ where each $λ_{j}$ is an unknown "price tag" that penalises breaking constraint $j$ . Because the bracketed constraint is zero for any physically valid mixture, $L$ equals $G$ there — so minimizing $L$ freely (no leftover constraints) finds the same answer as minimizing $G$ under constraints. At a free minimum every partial derivative vanishes: $\frac{\partial L}{\partial n _{i}} = μ_{i} - \sum_{j} λ_{j} a_{ij} = 0 ⟹ μ_{i} = j \sum λ_{j} a_{ij}$ What each line means: the derivative of $\sum_{i} n_{i} μ_{i}$ with respect to $n_{i}$ gives $μ_{i}$ ; the derivative of the constraint term gives $\sum_{j} λ_{j} a_{ij}$ . Setting them equal says: at equilibrium, each species' chemical potential is exactly the sum of its atoms' elemental potentials.

This one boxed line secretly contains every $K_{p}$ mass-action relation at once (linked topic: Chemical Equilibrium and Kp): for any reaction that shuffles atoms, the $λ_{j}$ cancel and you recover $\sum_{i} ν_{i} μ_{i} = 0$ .

Why this tool? Lagrange multipliers are the standard machine for "minimize this, keep those fixed." CEA solves the resulting equations by Newton iteration on the $λ_{j}$ and $n_{i}$ .

10. Nozzle symbols: $v_{e}$ , $γ$ , $M$ , $g_{0}$ , $I_{s p}$

Definition The other exhaust quantities

$v_{e}$ = exhaust velocity, how fast gas leaves the nozzle (units $m s^{- 1}$ ). Faster exhaust = more push.
$M$ = average (mean) molar mass of the exhaust in $kg mol^{- 1}$ — how heavy the average escaping molecule is. Because it is an average over all species it is often written with an overbar, $\overset{ˉ}{M} = \sum_{i} X_{i} M_{i}$ ; the overbar simply means "averaged over the mixture." In this note $M$ and $\overset{ˉ}{M}$ mean the same thing. CEA prints this quantity in its output table under the heading M (its molecular-weight column).
$g_{0} = 9.81 m s^{- 2}$ = standard gravity, used only to define $I_{s p}$ in seconds.
$I_{s p}$ = specific impulse $= v_{e} / g_{0}$ , the thrust delivered per unit weight-flow of propellant — the efficiency score (units: seconds). Feeds Tsiolkovsky Rocket Equation and Rocket Nozzle and Thrust Equation.

Where the exponent $(γ - 1) / γ$ comes from

Why the square root? Kinetic energy $\frac{1}{2} v_{e}^{2}$ is proportional to the enthalpy released, so speed goes as the square root of that energy — doubling the available heat multiplies speed by $2$ , not 2.

Prerequisite map

Read top to bottom: counting molecules leads to composition, composition plus temperature and the reference state builds chemical potential, that feeds Gibbs energy, which — constrained by atom conservation via the multipliers — pins down the equilibrium products and $T_{c}$ . Those then drive the nozzle to give $I_{s p}$ .

Equipment checklist

Cover the right side and test yourself:

What does

n_{i}

count?

The number of moles (a counting unit) of species

i

in the mixture.

What is the mole fraction

X_{i}

in terms of

n_{i}

X_{i} = n_{i} / n_{tot}

, the share of the total soup that is species

i

How is partial pressure

p_{i}

related to

X_{i}

p_{i} = X_{i} p

— each species' slice of the total pressure (ideal-gas result).

Why is temperature always in Kelvin here?

T

is never negative, since it multiplies energies that must stay physical.

What do the reference values

T^{\circ}

and

p^{\circ}

equal?

T^{\circ} = 298.15

K and

p^{\circ} = 1

bar — the standard state the

^{\circ}

superscript flags.

What are the units of enthalpy

h

and heat capacity

c_{p}

h

is in

J mol^{- 1}

;

c_{p}

is in

J mol^{- 1} K^{- 1}

What two parts make up the enthalpy

h_{i} (T)

Standard-state chemical part

Δ_{f} h^{\circ}

(formation energy at

T^{\circ}, p^{\circ}

) plus thermal part

\int_{T^{\circ}}^{T} c_{p} d T

(sensible heat).

What quantity does a constant-

T

, constant-

p

system minimize, and in what units?

Its total Gibbs free energy

G = \sum_{i} n_{i} μ_{i}

, measured in joules.

Why does a logarithm appear in

μ_{i}

It encodes mixing entropy — being dilute lowers a species' chemical potential, driving dissociation.

What model underlies

μ_{i} = μ_{i}^{\circ} + R T ln (p_{i} / p^{\circ})

The ideal-gas model — no intermolecular forces, molecules of negligible volume.

How does the method of Lagrange multipliers give

μ_{i} = \sum_{j} λ_{j} a_{ij}

Set

\partial / \partial n_{i}

L = G - \sum_{j} λ_{j} (\sum_{i} a_{ij} n_{i} - b_{j})

to zero; the two derivatives

μ_{i}

and

\sum_{j} λ_{j} a_{ij}

must balance.

What is the difference between

c_{p}

and

c_{v}

c_{p}

heats at constant pressure (gas free to expand and do work);

c_{v}

at constant volume (boxed in). So

c_{p} > c_{v}

and

γ = c_{p} / c_{v}

Where does the exponent

(γ - 1) / γ

come from?

The adiabatic ideal-gas law

p V^{γ} =

const combined with

p V = n R T

, giving

T / T_{c} = (p / p_{c})^{(γ - 1) / γ}

What do

a_{ij}

and

b_{j}

enforce?

Conservation of atoms:

\sum_{i} a_{ij} n_{i} = b_{j}

for each element

j

What is a Lagrange multiplier

λ_{j}

physically?

The "elemental potential" of element

j

— the cost of supplying one more mole of it.

What is

I_{s p}

and its formula?

Specific impulse

= v_{e} / g_{0}

, thrust per unit weight-flow, measured in seconds.

Why does

I_{s p} \propto T_{c} / M

favour

H_{2} / O_{2}

Hydrogen-rich exhaust is very light (small

\overset{ˉ}{M}

), so exhaust speed and thus

I_{s p}

are large.