Combustion Chemistry: The Science of Flames and Fire

3D Simulations

Chemistry & Physics

Chemistry · Thermodynamics · ⌛ ~15 min read · Last updated: 22 June 2026

Combustion Chemistry: The Science of Flames and Fire

Fire appears simple — fuel meets oxygen and heat is released — yet the reality is a cascade of hundreds of elementary reactions sustained by reactive radicals that branch, propagate, and terminate faster than a microsecond. Understanding combustion chemistry means grasping why a candle glows yellow while a Bunsen burner burns blue, how the same reaction that powers a gas turbine can detonate an engine, and why the equivalence ratio controls everything from flame temperature to pollutant formation.

Written by MySimulator Team · Reviewed by MySimulator Editorial Review

1. Fuel Oxidation and the Combustion Reaction

At its most fundamental, combustion is rapid oxidation: a fuel molecule donates electrons (is oxidised) while oxygen accepts them (is reduced). The global stoichiometry for a general hydrocarbon C_xH_y burning in air is:

C_xH_y + (x + y/4) O2 \to x CO2 + (y/2) H2O Methane (CH4): CH4 + 2 O2 \to CO2 + 2 H2O ΔH° = -890 kJ/mol Octane (C8H18): C8H18 + 12.5 O2 \to 8 CO2 + 9 H2O ΔH° = -5,471 kJ/mol Hydrogen (H2): H2 + ½ O2 \to H2O ΔH° = -286 kJ/mol Stoichiometric air-fuel ratio (AFR) by mass: Methane: AFR_s = 17.2 g air / g fuel Octane: AFR_s = 15.1 g air / g fuel Hydrogen: AFR_s = 34.3 g air / g fuel Equivalence ratio: Φ = (fuel/air)_actual / (fuel/air)_stoichiometric Φ = 1.0 \to stoichiometric Φ < 1.0 \to lean (excess air) Φ > 1.0 \to rich (excess fuel)

These global equations are thermodynamic summaries, not mechanistic descriptions. In reality, no molecule of methane simply collides with two oxygen molecules and produces CO₂ and water in a single step. Instead, thousands of intermediate species and hundreds of elementary reactions occur simultaneously within the thin reaction zone we call a flame. The global heat release is the net result of breaking C–H and C–C bonds (endothermic) and forming C=O and O–H bonds (strongly exothermic).

The lower heating value (LHV) counts only the sensible heat released, assuming water leaves as vapour; the higher heating value (HHV) includes the latent heat of condensation. For natural gas, LHV ≈ 50 MJ/kg and HHV ≈ 55 MJ/kg. Power plant efficiencies are always quoted against LHV; domestic boilers often use HHV because condensing boilers recover that latent heat.

2. Radical Chain Reactions and Flame Propagation

Flames are sustained by a branching chain reaction. The chain carriers are radicals — species with unpaired electrons — chiefly H•, O•, OH•, and HO₂•. The hydrogen–oxygen sub-mechanism underpins all hydrocarbon oxidation:

Initiation (requires thermal energy): H2 + O2 \to 2 OH\cdot (slow, high temperature) RH + O2 \to R\cdot + HO2\cdot (fuel radical formation) Chain branching (multiplies radicals — key to rapid combustion): H\cdot + O2 \to OH\cdot + O\cdot k_b \approx 2\times10¹⁴ exp(-16,800/T) cm³/mol\cdots O\cdot + H2 \to OH\cdot + H\cdot Chain propagation (converts fuel, transfers radical): OH\cdot + CH4 \to CH3\cdot + H2O (H-abstraction from fuel) CO + OH\cdot \to CO2 + H\cdot (most CO oxidised this way) Chain termination (radicals destroyed at walls or by 3-body collisions): H\cdot + OH\cdot + M \to H2O + M H\cdot + O2 + M \to HO2\cdot + M (dominant below 1,000 K — inhibits branching) Laminar burning velocity S_L (m/s): S_L = \sqrt(2α\cdotω_T) where α = thermal diffusivity, ω_T = heat release rate CH4/air at Φ=1.0: S_L \approx 0.40 m/s H2/air at Φ=1.0: S_L \approx 2.65 m/s Acetylene/air: S_L \approx 1.58 m/s

The branching step H• + O₂ → OH• + O• is the most important reaction in combustion chemistry. It converts one radical into two, exponentially accelerating reaction rates. Its competition with the three-body recombination H• + O₂ + M → HO₂• + M determines the crossover temperature (~1,000 K) below which flames cannot self-sustain. Above the crossover, branching dominates and the mixture ignites; below it, oxidation proceeds only as a slow, cool reaction.

Laminar flame thickness δ ≈ α/S_L is typically 0.1–1 mm for hydrocarbon–air flames at atmospheric pressure. Within this thin zone, temperature rises from the unburnt gas temperature (~300 K) to the adiabatic flame temperature in a distance comparable to the mean free path multiplied by the Damköhler number. Turbulence wrinkles and stretches the flame front, increasing its surface area and the global burning rate by factors of 5–50 in practical combustors.

3. Adiabatic Flame Temperature and Enthalpy Balance

The adiabatic flame temperature (AFT) is the upper bound on combustion temperature: the temperature attained when all heat of reaction raises the product gases and no energy escapes. It is found by an enthalpy balance at constant pressure:

Enthalpy balance (adiabatic, constant pressure): H_reactants(T_in) = H_products(T_ad) ΔH_comb(T_ref) + Σ nᵢ Cp,i(T_in - T_ref) = Σ nⱼ Cp,j(T_ad - T_ref) Solving for T_ad (iterative, as Cp = f(T)): T_ad = T_ref + ΔH_comb / [Σ nⱼ Cp,j(T_ad)] Adiabatic flame temperatures at Φ=1.0, 25°C initial: CH4 / air: T_ad \approx 2,230 K C8H18 / air: T_ad \approx 2,276 K H2 / air: T_ad \approx 2,480 K CH4 / O2: T_ad \approx 3,054 K (pure oxygen, no N2 diluent) C2H2 / O2: T_ad \approx 3,430 K (oxy-acetylene torch) Effect of equivalence ratio on T_ad (methane/air): Φ = 0.6: T_ad \approx 1,820 K (excess air absorbs heat) Φ = 0.8: T_ad \approx 2,050 K Φ = 1.0: T_ad \approx 2,230 K (maximum) Φ = 1.2: T_ad \approx 2,170 K (unburnt fuel absorbs heat) Φ = 1.5: T_ad \approx 1,900 K Real flames are cooler than T_ad due to: - Radiation losses (10-30% in diffusion flames) - Dissociation: CO2 \to CO + ½O2 above ~2,500 K absorbs energy - Heat transfer to burner walls

At temperatures above about 2,500 K, thermal dissociation of CO₂ and H₂O becomes significant; real products are a mixture of CO₂, CO, H₂O, H₂, OH, and atomic species. This dissociation absorbs energy and depresses the actual temperature below the ideal AFT. Oxy-fuel combustion (burning in pure oxygen rather than air) eliminates the massive heat sink of nitrogen (~79% of air), raising AFT by 600–800 K and enabling cutting torches and specialty glass-making.

4. Flame Structure: Premixed versus Diffusion Flames

Combustion engineers classify flames by how fuel and oxidiser are brought together before the reaction zone.

Premixed flames

Fuel and air are mixed upstream of the reaction zone. The flame front propagates back through the mixture at the laminar burning velocity S_L. Bunsen burners, spark-ignition (petrol) engines, and lean-premixed gas turbine combustors all use premixed flames. Structure from unburnt to burnt side:

Unburnt zone: T = T_u, [fuel] = initial, radicals \approx 0 Preheat zone: T rises from T_u to ~0.8 T_ad; fuel begins decomposing Reaction zone: δ \approx 0.1-0.5 mm; heat release rate peaks; radical pool maximum Post-flame zone: T \approx T_ad; equilibration; slow CO \to CO2 burnout

Diffusion (non-premixed) flames

Fuel emerges from a jet or surface and mixes with surrounding air by molecular and turbulent diffusion. The flame sheet forms where the mixture fraction Z equals the stoichiometric value Z_st:

Mixture fraction: Z = (β − β_ox) / (β_fuel − β_ox) where β is a conserved scalar (Bilger's formula) Z_st = 1 / (1 + s·Φ), s = stoichiometric oxygen-fuel mass ratio Burke–Schumann limit (infinitely fast chemistry): T(Z) = T_u + Z(T_f − T_u) / Z_st for Z ≤ Z_st (oxidiser side) T(Z) = T_f + (Z − Z_st)(T_f − T_u)/(Z_st − 1) for Z > Z_st (fuel side) T_f ≈ T_ad at Z = Z_st Candle flame: Flame height H ∝ Q_fuel (volumetric flow) — linear in laminar regime H = 4 Q_fuel / (π D_m) where D_m = fuel mass diffusivity

Diffusion flames are inherently rich on the fuel side and lean on the oxidiser side, passing through stoichiometry only at the flame sheet. They cannot flash back (no premixed fuel ahead of the front) but they cannot be as controlled or as clean as premixed flames, producing more soot and CO.

5. Pollutant Formation: Soot, NO_x, and CO

Combustion chemistry determines not only energy release but also the formation of pollutants with serious health and environmental consequences.

Thermal NOx (Zeldovich mechanism, dominant above ~1,800 K): O\cdot + N2 \to NO + N\cdot Ea \approx 314 kJ/mol (very temperature-sensitive) N\cdot + O2 \to NO + O\cdot N\cdot + OH\cdot \to NO + H\cdot Rate: d[NO]/dt = 2k₁[O][N2], k₁ = 1.8\times10⁸ exp(-38,370/T) cm³/mol\cdots Doubling temperature from 1,800 to 2,200 K increases NOx rate ~100-fold Soot formation pathway: Fuel pyrolysis \to C2H2 (acetylene) \to C6H6 (benzene) \to PAH growth PAH + C2H2 \to larger ring systems (HACA mechanism) PAH stacking \to nascent soot particle (1-5 nm diameter) Surface growth: C2H2 adds to particle surface \to 10-100 nm mature soot Oxidation: C(soot) + OH\cdot \to CO + H; C + O2 \to CO2 Soot yield: Φ > ~1.8 for most hydrocarbons CO formation and burnout: CO produced: RCO\cdot \to R + CO (fuel fragments) CO consumed: CO + OH\cdot \to CO2 + H\cdot (rate-limiting at burnout) CO/CO2 equilibrium at 2,000 K: CO \approx 1-5% (significant) At 1,000 K: CO < 10 ppm (equilibrium nearly complete) High-temperature quench (rapid cooling) freezes CO at high levels

The trade-off between NOx and soot is fundamental: fuel-rich conditions suppress NOx (lower temperature, less O available) but generate soot; fuel-lean conditions eliminate soot but raise temperature and NOx. This “NOx–soot trade-off” drove the development of staged combustion in diesel engines, where a rich primary zone burns fuel to CO and soot, then a lean secondary zone oxidises these at lower peak temperature with reduced NOx.

6. Real-World Applications

Combustion chemistry underlies several major technologies, each exploiting different aspects of flame physics.

Gas turbines: Lean premixed combustion at Φ ≈ 0.5–0.6 reduces peak temperatures to ~1,800 K, cutting thermal NOx by 90% compared with diffusion flames. Compressor outlet pressures of 30–50 bar require careful lean blowout margin management.
Petrol engines: Spark ignition initiates a premixed flame kernel that propagates at ~25 m/s across the cylinder in ~2 ms. Engine knock occurs when end-gas autoignites ahead of the flame front, producing a pressure wave that can destroy pistons. Octane rating quantifies resistance to autoignition via the research octane number (RON) and motor octane number (MON) tests.
Rocket propulsion: Liquid oxygen / liquid hydrogen combustion in a Vulcain 2 engine reaches chamber temperatures of ~3,500 K at 11 MPa. Injector design controls mixing at microscale; combustion instability — acoustic coupling with heat release — can destroy engines within milliseconds if not suppressed.
Fire safety: Halogenated suppression agents interrupt radical chains; water mist cools below ignition temperature and displaces oxygen. Intumescent coatings expand under heat to form an insulating char, slowing heat transfer to structural steel.
Industrial furnaces: Regenerative burners pre-heat combustion air to >1,000°C using waste heat, cutting fuel consumption by 30–50%. Oxy-fuel furnaces use pure oxygen to raise temperatures for glass-melting and steelmaking slag processing.

Understanding the coupling between chemistry, heat transfer, and fluid mechanics in reacting flows remains one of the most computationally demanding problems in engineering. Direct numerical simulation (DNS) of a turbulent methane–air flame requires resolving scales from the Kolmogorov length (<0.1 mm) to the combustor diameter (>100 mm), spanning four orders of magnitude simultaneously.

Frequently Asked Questions

What makes a flame self-sustaining?

A flame is self-sustaining when the heat released by the exothermic reactions is sufficient to preheat fresh reactants to the ignition temperature and maintain the radical pool. The Damköhler number Da = t_flow / t_chem must exceed roughly unity: chemical reactions must complete faster than the residence time in the reaction zone. Blow-off occurs when Da drops below ~0.5 because the flow sweeps radicals away faster than they are replenished.

Why are hydrogen flames considered more hazardous than natural gas flames?

Hydrogen has a very wide flammability range in air (4–75% by volume versus 5–15% for methane), a very low minimum ignition energy (~0.017 mJ versus ~0.28 mJ for methane), and a laminar burning velocity roughly six times that of methane. Hydrogen flames are nearly invisible (weak UV, no soot), making detection difficult. Its small molecular size also makes it prone to leaking through seals and diffusing upwards rapidly, creating ignitable pockets in enclosed spaces.

How does a catalytic converter use combustion chemistry?

A three-way catalytic converter (TWC) uses platinum, palladium, and rhodium to simultaneously oxidise CO and unburnt hydrocarbons to CO₂ and H₂O, and reduce NO to N₂. Light-off temperature is ~250°C; below this, the catalyst is inactive and emissions are highest. The oxygen storage component (cerium oxide) buffers brief rich excursions, maintaining near-stoichiometric conditions essential for all three reactions to operate simultaneously.