Combustion of hypergolics — N₂O₄ + UDMH - MMH; ignition delay: Visual walkthrough

Step 1 — What does "ignition delay" even measure?

WHAT: we place two events on a timeline — contact and flame — and call the gap between them $τ_{i g n}$ .

WHY: before we can find a formula for a quantity, we must agree exactly what the quantity is. If you and I time from different starting whistles we get different numbers. So we pin down: start the clock at first liquid contact, stop it when pressure/light jumps.

PICTURE: the red bar is the delay. It is short (milliseconds) but it is the whole game — too long and fuel piles up.

Figure — Combustion of hypergolics — N₂O₄ + UDMH - MMH; ignition delay

Step 2 — A reaction has a speed, and speed has a name: $k$

Before "delay" we need "how fast does the pre-ignition chemistry go?"

WHAT: we introduce a single dial, $k$ , that turns the reaction's speed up or down.

WHY this tool and not another? We could track concentrations, pressures, dozens of things. But the one knob that ignition delay hangs on is the intrinsic speed of the rate-limiting step. Isolating it as $k$ lets us reason about time without carrying all the chemistry.

PICTURE: two beakers — a fast reaction empties quickly (tall red arrow), a slow one dawdles (short red arrow). Same reactants, different $k$ .

Step 3 — Time is one over speed: $τ \propto 1/ k$

WHAT: we claim $τ_{i g n}$ is proportional to $1/ k$ .

WHY: think of a car. Speed $= 60$ km/h covers a fixed distance in half the time of speed $= 30$ km/h. Time and speed are reciprocals for a fixed "distance to travel". Here the fixed "distance" is how much chemistry must happen before flame; the speed is $k$ .

PICTURE: the classic $y = 1/ x$ curve. As $k$ (horizontal) grows, $τ$ (vertical, red) plunges toward zero. As $k \to 0$ (a dead-cold, unreactive mix) $τ \to \infty$ — it never lights. That is the "cold engine won't start" case, captured for free.

Step 4 — Where does $k$ come from? The Arrhenius idea

Now we open the box marked $k$ . Why is one reaction fast and another slow?

WHAT: we say $k$ depends on two things — the height of the hill $E_{a}$ , and the temperature $T$ that decides how many molecules can climb it.

WHY the exponential? The fraction of molecules with energy above a barrier $E_{a}$ follows the Boltzmann factor $e^{- E_{a} / R T}$ — an exponential, because energetic molecules get exponentially rarer as the barrier rises. This is why we need the tool $e^{x}$ and not a simple ratio: nature counts rare energetic molecules exponentially. (See Arrhenius equation.)

PICTURE: the reaction coordinate — a hill of height $E_{a}$ . A red dashed line marks the barrier; only the molecules above it (shaded) make it over.

Step 5 — Substitute: the minus sign flips to a plus

WHAT: put the Arrhenius $k$ (Step 4) into the reciprocal law (Step 3).

WHY: we want $τ$ in terms of the physical knobs $E_{a}$ and $T$ , not the abstract $k$ . So we eliminate $k$ by substitution.

$τ_{i g n} = \frac{C}{k} = \frac{C}{A e ^{- E_{a} / R T}}$

Now the crucial algebra move: dividing by $e^{- x}$ is the same as multiplying by $e^{+ x}$ (a negative exponent in the denominator climbs to the numerator and flips sign).

$τ_{i g n} = \frac{C}{A} e^{+ E_{a} / R T}$

Fold the two constants into one, still call it $C$ :

$τ_{i g n} = C e^{E_{a} / R T}$

PICTURE: the sign flip shown as a mirror — Arrhenius $k$ falls with $1/ T$ (its curve drops); ignition delay $τ$ , being $1/ k$ , rises with $1/ T$ (red curve, the mirror image).

Step 6 — Take the log: bend the curve into a straight line

An exponential curve is hard to read by eye. Straight lines are easy. So we apply the tool that undoes exponentials: the natural logarithm.

WHAT: take $ln$ of both sides of the boxed law.

WHY: a product inside a log splits into a sum ( $ln (ab) = ln a + ln b$ ), and $ln e^{x} = x$ . Together these turn a curvy exponential into a straight line.

$ln τ_{i g n} = ln (C e^{E_{a} / R T}) = ln C + = E_{a} / R T ln e^{E_{a} / R T}$

$ln τ_{i g n} = ln C + \frac{E _{a}}{R} \cdot \frac{1}{T}$

Match it to the shape of a line, $y = c + m x$ :

PICTURE: the straight line on axes $y = ln τ$ versus $x = 1/ T$ . The red line has slope $E_{a} / R$ and intercept $ln C$ .

Step 7 — Two points give the slope: the worked extraction

WHAT: the lab measures $τ$ at just two temperatures and reads $E_{a}$ .

WHY subtract? Writing the line at $T_{1}$ and at $T_{2}$ and subtracting cancels the unknown intercept $ln C$ , leaving one equation in one unknown, $E_{a}$ .

$ln τ_{1} - ln τ_{2} = \frac{E _{a}}{R} (\frac{1}{T _{1}} - \frac{1}{T _{2}}) \Rightarrow ln \frac{τ _{1}}{τ _{2}} = \frac{E _{a}}{R} (\frac{1}{T _{1}} - \frac{1}{T _{2}})$

Plug in the parent's numbers: $τ_{1} = 8 ms$ at $T_{1} = 250 K$ , $τ_{2} = 2 ms$ at $T_{2} = 300 K$ .

$ln \frac{8}{2} = ln 4 = 1.386, \frac{1}{250} - \frac{1}{300} = 6.667 \times 1 0^{- 4} K^{- 1}$

$E_{a} = \frac{R l n 4}{6.667 \times 1 0 ^{- 4}} = \frac{8.314 \times 1.386}{6.667 \times 1 0 ^{- 4}} \approx 1.73 \times 1 0^{4} J/mol = 17.3 kJ/mol$

PICTURE: two red dots on the straight line; the slope drawn between them (rise $= ln 4$ , run $= 6.667 \times 1 0^{- 4}$ ) is literally $E_{a} / R$ .

Step 8 — Edge & degenerate cases (never leave the reader stranded)

PICTURE: all three limits on one plot — the horizontal floor at $τ = C$ (Case A/B) and a curving dataset (Case C) shown in red against the ideal straight line.

The one-picture summary

Everything above collapses into a single chain, and a single plot. Read left to right: an energy hill sets $k$ ; $k$ inverts to $τ$ ; the exponential straightens under $ln$ into a line whose slope hands you $E_{a}$ .

Recall Feynman retelling — the whole walkthrough in plain words

We started with a stopwatch: the ignition delay is just how long you wait from "the two liquids touch" to "flame". Then we asked what controls that wait. A reaction has a speed — call it $k$ ; go faster, finish sooner, so the wait is one over the speed. Next we opened up the speed: molecules have to climb an energy hill before they can react, and only the hottest ones make it over. That fraction shrinks exponentially as the hill gets taller or the temperature drops — that's the Arrhenius $e^{- E_{a} / R T}$ . Because the wait is one-over-speed, its exponent flips to plus: cold means a long, dangerous wait; hot means a quick light. Finally we tamed the curve with a logarithm, which turns the exponential into a dead-straight line — plot the log of the delay against one-over-temperature and the tilt of that line, times the gas constant, is the height of the hill. Two lab measurements are enough: subtract them, the messy constant cancels, and out pops $E_{a} \approx 17 kJ/mol$ — a low hill, which is the real reason these propellants light themselves even in the cold of space.

Recall Quick self-check

Why does $τ$ have $+ E_{a} / R T$ while $k$ has $- E_{a} / R T$ ? ::: Because $τ \propto 1/ k$ ; reciprocating an exponential flips the sign of its exponent. What does the slope of $ln τ$ vs $1/ T$ equal? ::: $E_{a} / R$ ; multiply the measured slope by $R$ to get the activation energy. As $T \to \infty$ , what does $τ$ approach and why? ::: The constant $C$ — the physical mixing floor; chemistry can vanish but you still must bring the liquids together.