Level 4 — ApplicationCombustion Chemistry (Propulsion Bridge)

Combustion Chemistry (Propulsion Bridge)

60 minutes60 marksprintable — key stays hidden on paper

Level 4 — Application (novel/unseen problems, no hints) Time limit: 60 minutes Total marks: 60

Useful data (use where relevant; assume ideal gas, R=8.314 J mol1K1R = 8.314\ \text{J mol}^{-1}\text{K}^{-1}):

Species ΔHf\Delta H_f^\circ (kJ/mol, 298 K) cˉp\bar{c}_p (J mol⁻¹K⁻¹, mean)
CH₄(g) −74.8
H₂O(g) −241.8 52
CO₂(g) −393.5 62
N₂(g) 0 35
O₂(g) 0 38
H₂(g) 0
CO(g) −110.5
N₂O₄(l) −19.6
MMH CH₃NHNH₂(l) +54.2

Air is 21% O₂ / 79% N₂ by mole (N₂/O₂ = 3.76).


Q1. (12 marks) — Equivalence ratio & product streams

A gas turbine burns methane with air at an equivalence ratio ϕ=0.80\phi = 0.80.

(a) Write the stoichiometric combustion equation for CH₄ in air (complete combustion). (2)

(b) For ϕ=0.80\phi = 0.80, determine the actual moles of air supplied per mole of CH₄, and state whether the mixture is fuel-lean or fuel-rich. (3)

(c) Write the balanced product equation (assume complete combustion of all fuel, excess O₂ appears in products). (4)

(d) Compute the mole fraction of O₂ in the wet exhaust. (3)


Q2. (14 marks) — Adiabatic flame temperature

Methane burns with the stoichiometric amount of air. Reactants enter at 298 K. Products are CO₂, H₂O(g) and N₂ only (no dissociation).

(a) Using enthalpies of formation, compute the heat released per mole of CH₄ at 298 K (products as gases). (4)

(b) Using the mean cˉp\bar{c}_p values in the table, estimate the adiabatic flame temperature TadT_{ad} (constant pressure). (6)

(c) A real measurement gives ~2230 K. Give two distinct physical reasons your value is an overestimate, and explain how each lowers TadT_{ad}. (4)


Q3. (12 marks) — High-temperature dissociation equilibrium

At the flame temperature, consider the dissociation CO2CO+12O2\text{CO}_2 \rightleftharpoons \text{CO} + \tfrac{1}{2}\text{O}_2 with Kp=0.010K_p = 0.010 (bar-based) at that temperature. Start with 1 mol pure CO₂ at total pressure 1 bar, and let α\alpha be the fraction of CO₂ dissociated at equilibrium.

(a) Express the mole fractions of all species and total moles in terms of α\alpha. (4)

(b) Write KpK_p in terms of α\alpha and total pressure PP (in bar), then solve for α\alpha (small-α\alpha approximation acceptable; state your assumption). (5)

(c) State qualitatively (with reasoning) how α\alpha changes if the chamber pressure is raised to 20 bar, and why this matters for computing TadT_{ad}. (3)


Q4. (12 marks) — Detonation vs deflagration; solid propellant burn rate

(a) In a table, contrast a deflagration and a Chapman–Jouguet detonation on: propagation mechanism, wave speed relative to sound, and pressure change across the wave. (4)

(b) A composite solid propellant (AP/HTPB/Al) follows Vieille's law r=aPnr = a\,P^n. In a static test the burn rate is 6.0 mm/s6.0\ \text{mm/s} at 5.0 MPa5.0\ \text{MPa} and 9.0 mm/s9.0\ \text{mm/s} at 10.0 MPa10.0\ \text{MPa}. (i) Determine the pressure exponent nn and the coefficient aa (with units, PP in MPa, rr in mm/s). (4) (ii) A stable motor requires n<1n < 1. State whether this propellant is stable, and explain physically what n>1n>1 would imply for chamber-pressure control. (2)

(c) State the role of aluminium in the AP/HTPB/Al formulation and one drawback it introduces in the exhaust. (2)


Q5. (10 marks) — Hypergolics, flames & pollutants

(a) N₂O₄ oxidises MMH (CH₃NHNH₂). Assuming products are N₂, CO₂ and H₂O(g), balance the reaction and state the oxidiser-to-fuel mole ratio. (4)

(b) Define ignition delay for a hypergolic pair and explain why a short ignition delay is desirable for engine restart/throttling. (2)

(c) A hydrogen–oxygen rocket produces essentially no CO₂ or soot, yet a kerosene (RP-1) engine produces both plus NOₓ. For each of soot and NOₓ, name the combustion condition that promotes it and one mitigation strategy. (4)

Answer keyMark scheme & solutions

Q1 (12)

(a) Stoichiometric complete combustion: CH4+2O2+2(3.76)N2CO2+2H2O+7.52N2\text{CH}_4 + 2\text{O}_2 + 2(3.76)\text{N}_2 \rightarrow \text{CO}_2 + 2\text{H}_2\text{O} + 7.52\,\text{N}_2 Correct O₂ (2) and N₂ (7.52) — (2). Why: 2 mol O₂ needed per CH₄; N₂ carried along at 3.76× O₂.

(b) ϕ=(F/A)actual/(F/A)stoich\phi = (\text{F/A})_{actual}/(\text{F/A})_{stoich}, so air 1/ϕ\propto 1/\phi. Stoich air = 2 mol O₂ = 9.52 mol air. Actual air =9.52/0.80=11.9= 9.52/0.80 = 11.9 mol air per mol CH₄, i.e. 2.5 mol O₂ + 9.40 mol N₂. (2) ϕ<1\phi<1 \Rightarrow excess air \Rightarrow fuel-lean. (1)

(c) With 2.5 mol O₂ supplied, 2 consumed, 0.5 excess: CH4+2.5O2+9.40N2CO2+2H2O+0.5O2+9.40N2\text{CH}_4 + 2.5\text{O}_2 + 9.40\text{N}_2 \rightarrow \text{CO}_2 + 2\text{H}_2\text{O} + 0.5\text{O}_2 + 9.40\text{N}_2 Balanced C, H, O, N — (4) (1 each element balance / correct excess O₂).

(d) Total product moles =1+2+0.5+9.40=12.90= 1 + 2 + 0.5 + 9.40 = 12.90. xO2=0.5/12.90=0.03883.9%x_{O_2} = 0.5/12.90 = 0.0388 \approx \mathbf{3.9\%} (3) (total moles 1, formula 1, value 1).

Q2 (14)

(a) Reaction: CH₄ + 2O₂ → CO₂ + 2H₂O(g). ΔHrxn=[(393.5)+2(241.8)][(74.8)+0]=877.1+74.8=802.3 kJ/mol\Delta H_{rxn} = [(-393.5) + 2(-241.8)] - [(-74.8) + 0] = -877.1 + 74.8 = -802.3\ \text{kJ/mol} Heat released 802.3 kJ\approx \mathbf{802.3\ kJ} per mol CH₄. (4)

(b) Products (stoich air): 1 CO₂, 2 H₂O, 7.52 N₂. All heat goes to raise products from 298 K. ncˉp=1(62)+2(52)+7.52(35)=62+104+263.2=429.2 J/K\sum n\bar c_p = 1(62) + 2(52) + 7.52(35) = 62 + 104 + 263.2 = 429.2\ \text{J/K} ΔT=802300429.2=1869 K\Delta T = \frac{802\,300}{429.2} = 1869\ \text{K} Tad=298+1869=2167 K (2170 K)T_{ad} = 298 + 1869 = \mathbf{2167\ K \ (\approx 2170\ K)} (6) (heat balance setup 2, ncp\sum nc_p 2, TadT_{ad} 2).

(c) Any two, each 2 marks:

  • Dissociation (CO₂→CO+½O₂, H₂O→OH+H) is endothermic, absorbing energy → lower T.
  • cpc_p rises with temperature; using low mean values underestimates heat capacity, so real ΔT\Delta T is smaller.
  • Heat loss to walls / radiation / incomplete combustion — real process not perfectly adiabatic.

Q3 (12)

(a) CO₂ → CO + ½O₂. Start 1 mol CO₂; dissociate α\alpha:

  • CO₂: 1α1-\alpha
  • CO: α\alpha
  • O₂: α/2\alpha/2

Total nt=1+α/2n_t = 1 + \alpha/2. (4) (species 3, total 1) Mole fractions: xCO2=(1α)/(1+α/2)x_{CO_2}=(1-\alpha)/(1+\alpha/2), xCO=α/(1+α/2)x_{CO}=\alpha/(1+\alpha/2), xO2=(α/2)/(1+α/2)x_{O_2}=(\alpha/2)/(1+\alpha/2).

(b) Kp=xCOxO21/2xCO2P1/2K_p = \frac{x_{CO}\,x_{O_2}^{1/2}}{x_{CO_2}}\,P^{1/2} =α1α(α/21+α/2)1/2P1/2= \frac{\alpha}{1-\alpha}\left(\frac{\alpha/2}{1+\alpha/2}\right)^{1/2} P^{1/2} Small-α (so 1α11-\alpha\approx1, 1+α/211+\alpha/2\approx1): Kpα(α2)1/2P1/2=α3/22P1/2K_p \approx \alpha\left(\frac{\alpha}{2}\right)^{1/2}P^{1/2} = \frac{\alpha^{3/2}}{\sqrt2}P^{1/2} With P=1P=1: α3/2=2Kp=1.4142(0.010)=0.014142\alpha^{3/2} = \sqrt2\,K_p = 1.4142(0.010) = 0.014142. α=(0.014142)2/3=0.0585 (5.9%)\alpha = (0.014142)^{2/3} = \mathbf{0.0585 \ (\approx 5.9\%)} (5) (Kp expression 2, substitution/approx 1, solve 2). (Assumption stated: α small.)

(c) Raising PP increases the P1/2P^{1/2} factor, so to hold KpK_p fixed α\alpha must decrease (Le Chatelier: dissociation increases moles, so high pressure suppresses it). (2) Less dissociation → less endothermic absorption → higher TadT_{ad}, so pressure raises flame temperature and must be included in CEA-type calculations. (1)

Q4 (12)

(a) (4) — 1 mark per correct contrast row:

Property Deflagration CJ Detonation
Mechanism subsonic flame, heat/species diffusion shock-coupled reaction
Speed vs sound subsonic (m/s) supersonic (km/s), = sound speed of burnt gas (CJ point)
Pressure ~constant / slight drop large pressure rise

(b)(i) r=aPnr = aP^n. Ratio: 9.0/6.0=(10/5)n1.5=2n9.0/6.0 = (10/5)^n \Rightarrow 1.5 = 2^n. n=ln1.5ln2=0.585n = \frac{\ln 1.5}{\ln 2} = 0.585 a=r/Pn=6.0/50.585=6.0/2.560=2.34a = r/P^n = 6.0/5^{0.585} = 6.0/2.560 = \mathbf{2.34} (mm/s·MPa⁻ⁿ). (4) (n 2, a 2)

(ii) n=0.585<1n = 0.585 < 1stable. (1) If n>1n>1, a small pressure rise raises burn rate faster than the nozzle can vent mass → runaway pressure (uncontrolled), risking rupture. (1)

(c) Al raises flame temperature / specific impulse and adds energy density (combustion enthalpy); it also damps combustion instabilities. (1) Drawback: solid Al₂O₃ particles cause two-phase flow losses / smoky exhaust. (1)

Q5 (10)

(a) MMH = CH₃NHNH₂ = CH₆N₂. Oxidise with N₂O₄: Balance C: 1 CO₂; H: 6 H → 3 H₂O; O: CO₂+H₂O need 2+3=5 O → 2.5 N₂O₄ (5 O); N: fuel 2 + oxidiser 5 = 7 N → 3.5 N₂. CH6N2+2.5N2O4CO2+3H2O+3.5N2\text{CH}_6\text{N}_2 + 2.5\,\text{N}_2\text{O}_4 \rightarrow \text{CO}_2 + 3\text{H}_2\text{O} + 3.5\,\text{N}_2 (×2 for integers: 2MMH+5N2O42CO2+6H2O+7N22\text{MMH} + 5\text{N}_2\text{O}_4 \to 2\text{CO}_2 + 6\text{H}_2\text{O} + 7\text{N}_2.) O/F mole ratio = 2.5 (N₂O₄ per MMH). (4) (balance 3, ratio 1)

(b) Ignition delay = time between contact/mixing of oxidiser and fuel and the onset of sustained combustion (pressure rise). (1) Short delay → prompt, repeatable ignition on valve opening; avoids accumulation of unburned propellant that can cause a hard-start/pressure spike; enables reliable restart & throttling. (1)

(c) Soot: promoted by fuel-rich / diffusion-flame (locally O₂-deficient) conditions; mitigate by better air–fuel mixing / premixing / higher O/F. (2) NOₓ: promoted by high flame temperature (thermal/Zeldovich NO); mitigate by lowering peak T — lean/staged combustion, dilution, or shorter residence time. (2) (H₂/O₂ has no carbon → no CO₂/soot; still can form NOₓ only if N present.)

[
  {"claim":"Q2 heat of reaction is -802.3 kJ/mol",
   "code":"dH=(-393.5+2*(-241.8))-(-74.8+0); result=abs(dH+802.3)<0.1"},
  {"claim":"Q2 adiabatic T approx 2167 K",
   "code":"ncp=1*62+2*52+7.52*35; dT=802300/ncp; Tad=298+dT; result=abs(Tad-2167)<5"},
  {"claim":"Q3 dissociation fraction alpha approx 0.0585",
   "code":"import sympy as sp; Kp=sp.Rational(1,100); a=(sp.sqrt(2)*Kp)**(sp.Rational(2,3)); result=abs(float(a)-0.0585)<0.001"},
  {"claim":"Q4 Vieille exponent n approx 0.585 and a approx 2.34",
   "code":"import math; n=math.log(1.5)/math.log(2); a=6.0/5**n; result=abs(n-0.585)<0.005 and abs(a-2.34)<0.02"},
  {"claim":"Q1 O2 mole fraction in exhaust approx 0.0388",
   "code":"tot=1+2+0.5+9.40; x=0.5/tot; result=abs(x-0.0388)<0.001"}
]