The Arrhenius Driver: Why Thermal Runaway Accelerates
Stage-by-Stage Electrochemistry
Detection Architecture for Advanced BMS
Key Takeaways
Stage 1: SEI Decomposition (80–120°C)
Stage 2: Electrolyte Decomposition and Separator Failure (120–160°C)
Stage 3: Cathode Decomposition — Point of No Return
Stage 4: Full Thermal Runaway (>200°C)
Stage 5: Propagation

The Arrhenius relationship between temperature and reaction rate is the fundamental reason thermal runaway is self-accelerating — each 10°C rise doubles reaction rate, and the heat that temperature rise releases raises the temperature further.

TL;DRKey Points

Thermal runaway proceeds through five distinct electrochemical stages, each with a specific temperature window, gas signature, and intervention opportunity — by Stage 3 the process is irreversible.
The Arrhenius relationship causes exponential acceleration: each 10°C temperature rise roughly doubles the reaction rate, creating a positive feedback loop that cannot be stopped by cooling once the Frank-Kamenetskii criticality threshold is crossed.
LFP cathodes require ~90°C more heat to reach the Stage 3 equivalent and release far less oxygen than NMC 811, making LFP pack fires significantly less energetic and easier to contain.
Multi-sensor BMS detection (temperature rate-of-rise + CO gas + HF gas) provides up to 90 seconds more warning than temperature-threshold-only detection — critical for AIS-156 Phase 2 compliance.
Separator design (PE vs PP vs ceramic-coated) determines the Stage 2 intervention window; ceramic-coated separators extend usable response time by ~50°C.

This article describes thermal runaway at the mechanistic level — the specific electrochemical reactions, temperature windows, kinetic drivers, and gas compositions that define each stage. The focus is on what distinguishes LFP from NMC at the molecular level, what the intervention opportunities actually are in engineering terms, and how the Arrhenius kinetics explain why thermal runaway accelerates so catastrophically in the final stages.

Arrhenius temperature coefficient for SEI decomposition | 60–80 | kJ/mol activation energy

SEI exothermic heat generation rate at 120°C (NMC) | ~5–30 | W/kg

Cathode decomposition heat generation (NMC 811, 100% SOC) | ~300–700 | J/g

LFP cathode decomposition heat generation | ~200–400 | J/g (less intense, no O₂ release)

Separator integrity loss temperature (PE) | ~130–150 | °C

Peak ARC cell temperature (NMC 811, 100% SOC) | 700–950 | °C

The Arrhenius Driver: Why Thermal Runaway Accelerates

The fundamental driving force of thermal runaway is the Arrhenius relationship between temperature and reaction rate. For an exothermic reaction with activation energy Ea:

Heat generation rate Q_gen = ΔH × A × c × exp(-Ea/RT)

PYTHON

# Adiabatic rate calorimetry (ARC) data analysis
import numpy as np

def detect_arc_onset(temps_degC: list, times_min: list,
                     threshold_dT_dt: float = 0.02):
    """
    Find onset temperature in ARC experiment.
    ARC operates in heat-wait-search (HWS) mode.
    Onset = first temp where self-heating rate > threshold [deg C/min].
    """
    dT = np.diff(temps_degC)
    dt = np.diff(times_min)
    rates = dT / dt   # deg C/min

    for i, rate in enumerate(rates):
        if rate >= threshold_dT_dt:
            return temps_degC[i]   # onset temperature [deg C]
    return None

# Example ARC trace -- NMC811 charged to 4.3 V
# Typical onset: 185 deg C for NMC811 at full charge
arc_temps = [25, 60, 80, 100, 120, 140, 160, 180, 185, 190, 210, 280, 500]
arc_times  = [0,  5,  10, 15,  20,  25,  30,  35,  36,  36.5, 37,  37.5, 38]

T_onset = detect_arc_onset(arc_temps, arc_times)
print(f"ARC self-heating onset: {T_onset} deg C")
print(f"Safety margin vs operating limit (60 deg C): {T_onset - 60:.0f} deg C")

where ΔH is the reaction enthalpy, A is the pre-exponential factor, c is reactant concentration, R is the gas constant, and T is absolute temperature.

For most of the decomposition reactions in a lithium-ion cell (SEI, electrolyte, cathode), Ea is in the range of 60–120 kJ/mol. At these activation energies, a 10°C rise in temperature increases the reaction rate by approximately 2–4×. More critically, the heat generation rate increases exponentially while the heat dissipation rate (dominated by thermal conduction and convection) increases only linearly with temperature.

There is a critical temperature — the Frank-Kamenetskii criticality point — above which heat generation permanently exceeds dissipation. Below this point, the system can reach a stable (though elevated) temperature. Above it, temperature rises indefinitely.

◆

Key Insight

The Frank-Kamenetskii parameter δ = (Ea × Q_gen × V) / (RT² × λ) must exceed ~0.88 for a spherical geometry to initiate thermal runaway. For a 21700 cell at typical conditions, this criticality is reached around 80–120°C depending on the specific decomposition reactions active. This is why the Stage 1 window is the critical intervention zone — at this temperature, the system is still below or near criticality and cooling intervention can prevent progression.

QWhat does the Frank-Kamenetskii criticality condition mean in practical BMS terms?

It defines the temperature threshold above which heat generation from exothermic cell reactions permanently exceeds heat dissipation. For a typical 21700 cylindrical cell, this criticality is reached around 80–120°C depending on active reactions. Below it, cooling can stabilise the cell; above it, thermal runaway is inevitable regardless of cooling intervention. The BMS Stage 1 detection window — approximately 80–120°C — is the only period where the system is still below or near criticality and active cooling can prevent progression to Stage 2.

Stage-by-Stage Electrochemistry

Stage 1: SEI Decomposition (80–120°C)

The SEI (Solid Electrolyte Interphase) on the graphite anode consists of metastable organic compounds — primarily Li-alkyl carbonates (ROCO₂Li) and inorganic Li₂CO₃, LiF. The organic components have decomposition onset temperatures around 80–100°C:

ROCO₂Li → Li₂CO₃ + CO₂ + hydrocarbons

This reaction is exothermic: ~200–500 J/g of SEI material. The SEI mass is small relative to the cell (< 1% of total), so the total heat is modest. However, in a poorly cooled pack, even this small heat can raise cell temperature by 5–10°C — enough to accelerate into Stage 2.

Products: CO₂ (inert, non-flammable), hydrocarbons (trace), water vapour.

Detection: Temperature rate-of-rise above expected (1–3°C/min above baseline). CO₂ sensor signal if baseline is established. This stage is manageable — cooling intervention prevents progression.

Stage 2: Electrolyte Decomposition and Separator Failure (120–160°C)

At ~120°C, the liquid electrolyte solvents begin to react exothermically. For EC-based electrolytes:

EC + PF₅ → CO₂ + gaseous products DEC oxidation → CO + CO₂ + hydrocarbons

The heat generation rate is ~10–100 W/kg of electrolyte. Simultaneously, LiPF₆ salt decomposes: LiPF₆ → LiF + PF₅ PF₅ + H₂O → POF₃ + 2HF

HF generation is particularly significant — HF is acutely toxic (IDLH: 30 ppm), and corrosive to the aluminium current collectors and housing. Gas sensors in the pack headspace can detect HF as a Stage 2 indicator.

As the temperature approaches the separator's softening point (PE: ~130–135°C), the separator pores close (thermal shutdown feature designed in), briefly stopping ion transport — a self-protective mechanism. But if temperature continues rising, the separator physically deforms and can no longer maintain electrode separation. Local hotspots where the separator has melted create internal short circuits.

Event	Temperature	Gas Produced	Heat Rate	BMS Detectability
SEI organic decomp.	80–100°C	CO₂, trace HC	5–30 W/kg	Temperature rate-of-rise
Electrolyte oxidation	120–140°C	CO, CO₂, HC	10–100 W/kg	Gas sensor (CO, HF)
Salt decomp. (LiPF₆)	100–130°C	HF, POF₃	Moderate	Gas sensor (HF specific)
Separator softening	PE: 130°C / PP: 165°C	None	N/A	Voltage deviation (short)
Separator melting	PE: 150°C / PP: 175°C	None — cell shorts	Rapid	Voltage collapse

QWhy is HF the most dangerous gas produced during thermal runaway and how should a BMS respond to it?

HF (hydrogen fluoride) is produced by LiPF₆ salt decomposition beginning around 100–130°C. Its IDLH (Immediately Dangerous to Life and Health) is only 30 ppm — making it the most acutely toxic species in vent gas. Unlike CO, which requires thousands of ppm to be immediately lethal, a few breaths above 100 ppm HF can cause fatal pulmonary oedema with 24–48 hour delayed onset. A BMS detecting HF above 5 ppm should immediately trigger a Level 3 emergency response: open all contactors, activate emergency venting, alert the driver, and transmit a telematics emergency signal if available.

Stage 3: Cathode Decomposition — Point of No Return

For NMC, once fully delithiated above approximately 160–180°C:

Li_xNi₁₋y₋zMn_yCo_zO₂ → Li_xNi₁₋y₋zMn_yCo_zO₂₋δ + δ/2 O₂

The oxygen released reacts immediately with the electrolyte decomposition products (CO, H₂, hydrocarbons) in an oxidation reaction that generates intense heat — 300–700 J/g of cathode material. This reaction is exothermic, oxygen-producing, and self-sustaining. The external oxygen supply is irrelevant once it starts.

LFP cathode at Stage 3 equivalent temperatures:

LiFePO₄ → FeO + Li₃PO₄ + O₂ (this reaction requires ~270°C)

But the key difference: the FeO and Li₃PO₄ products are thermally stable, and the oxygen yield per gram is much lower. LFP gas is predominantly CO₂, not O₂-laden mixtures. The heat generation is ~200–400 J/g — lower than NMC, and without the self-sustaining oxidation mechanism.

Warning

The 90°C gap between NMC 811 and LFP Stage 3 onset is the most important safety parameter for Indian commercial EV pack design. In a 45°C ambient with poor cooling, an NMC 811 cell that is overcharged to 4.25 V per cell is approximately 135°C below its Stage 3 onset — marginally sufficient for most designs. The equivalent LFP cell at 3.70 V has 225°C margin. For summer fleet operation, the LFP margin provides substantially greater safety headroom for thermal management failures.

Stage 4: Full Thermal Runaway (>200°C)

By Stage 4, the cell is at 200–500°C and rising at 100–1,000°C/minute. All organic materials within the cell combust. The cell case (typically steel or aluminium) may rupture if venting is inadequate. Rupture ejects burning electrolyte droplets and solid cathode fragments at high velocity — a significant ignition hazard for adjacent cells and the pack structure.

Stage 5: Propagation

PYTHON

# Vent gas composition model for Li-ion cells
# Gases released during thermal runaway -- critical for safety room design

vent_gas_composition = {
    "LFP (at 200 deg C)": {
        "CO2":   45.0,   # % by volume
        "CO":    25.0,
        "H2":    15.0,
        "CH4":    8.0,
        "C2H4":  5.0,
        "HF":     2.0,   # extremely toxic
    },
    "NMC 811 (at 250 deg C)": {
        "CO2":   38.0,
        "CO":    30.0,
        "H2":    18.0,
        "CH4":    6.0,
        "C2H4":  4.0,
        "HF":     4.0,   # higher with high-voltage electrolyte
    },
}

def propagation_time_s(cell_spacing_mm: float, thermal_conductivity_W_mK: float, specific_heat_J_gK: float, cell_mass_g: float, onset_temp_C: float, ambient_temp_C: float = 25.0) -> float: """ Estimate time for thermal runaway to propagate to adjacent cell. Simplified 1D heat conduction model. """ delta_T = onset_temp_C - ambient_temp_C

Heat flux from runaway cell: ~1000-5000 W/cell peak

Q_peak_W = 2000.0 t_propagate = (cell_mass_g specific_heat_J_gK delta_T) / Q_peak_W return t_propagate

Cylindrical cells in module, 3mm air gap

t = propagation_time_s( cell_spacing_mm=3.0, thermal_conductivity_W_mK=0.026, # air specific_heat_J_gK=1.0, cell_mass_g=70.0, onset_temp_C=185.0 ) print(f"Propagation time (air gap, 3mm): {t:.0f} s ~= {t/60:.1f} min") print("-> BMS must detect and disconnect within this window")

TEXT




# Lethal concentration data (IDLH = Immediately Dangerous to Life/Health)
idlh_ppm = {"CO": 1200, "H2": 50000, "HF": 30, "CH4": 50000}

print("Vent gas composition and toxicity:")
for chemistry, gases in vent_gas_composition.items():
    print(f"\n{chemistry}:")
    for gas, pct in gases.items():
        ppm_in_1m3_vent = pct * 10000   # rough estimate: 1% ~= 10,000 ppm
        idlh = idlh_ppm.get(gas, None)
        danger = ""
        if idlh and ppm_in_1m3_vent > idlh:
            danger = f" !! EXCEEDS IDLH ({idlh} ppm)"
        print(f"  {gas:<8} {pct:.1f}%{danger}")

Heat transferred from a Stage 4 cell to its neighbours arrives via three mechanisms:

Conduction: Direct contact heat transfer — dominant for cells in tight contact
Radiation: Significant above 300°C — can heat adjacent cells through air gaps
Hot gas convection: Vented hot gas from failed cell heats adjacent cells

Propagation time (time from first cell failure to second cell entering Stage 3) depends on thermal barrier quality, cell-to-cell spacing, and cooling system response. Well-designed LFP packs achieve >15 minutes between first and second cell failure. NMC 811 packs without adequate barriers can propagate in <60 seconds.

QHow does cell-to-cell propagation occur in Stage 5 and what determines its speed?

Propagation occurs via three heat transfer mechanisms: conduction through direct contact (dominant in tightly packed modules), thermal radiation above 300°C (significant across air gaps), and hot gas convection from vented cell gases. The propagation time — from first cell failure to second cell entering Stage 3 — depends on cell spacing, inter-cell thermal barrier conductivity, and the cooling system's response speed. Well-designed LFP packs with intumescent barriers achieve >15 minutes cell-to-cell propagation time. NMC 811 packs without adequate barriers can propagate in under 60 seconds, which is why AIS-156 Phase 2 compliance requires validated thermal barriers, not just cell-level safety data.

Detection Architecture for Advanced BMS

Baseline characterisation

Run controlled TR tests on representative cell population at 25°C and 45°C. Record temperature profiles, gas signatures, and voltage behaviour at each stage. These define sensor thresholds.

Multi-sensor array design

Minimum: per-cell temperature + CO sensor in pack headspace. Enhanced: per-cell voltage deviation monitoring + HF sensor + pack pressure sensor.

Rate-of-rise algorithms

Implement exponential filter on dT/dt per cell. Threshold: 1.5×(expected dT/dt at current load) → Level 1 alert. 5× → emergency response.

Voltage deviation detection

Compare cell voltage to EKF-predicted voltage. Persistent negative deviation >50 mV at <1C load → possible internal short → Level 2 alert.

Gas sensor integration

CO baseline monitoring with drift compensation. >50 ppm above baseline → Level 2. HF detection >5 ppm → immediate Level 3 (emergency).

Response cascade

Level 1: alert VCU, increase cooling. Level 2: isolate module, maximum cooling, alert driver. Level 3: open all contactors, activate emergency venting, send emergency signal (if telematics present).

Note

The gas sensor (specifically CO and HF) provides the earliest definitive thermal runaway signal, typically 30–90 seconds before temperature sensors reach Level 2 thresholds. For new BMS designs targeting AIS-156 Phase 2 compliance, gas sensing is not mandatory but significantly improves the detection margin. The incremental BOM cost of a CO + HF sensor array is ₹500–1,500 per pack — minor relative to the pack cost and the regulatory liability of a propagation event.

Key Takeaways

✦Key Takeaways

The Arrhenius exponential relationship between temperature and reaction rate is the physical mechanism that makes thermal runaway self-accelerating. Once above the Frank-Kamenetskii criticality temperature (~80–120°C for a typical cell), cooling intervention cannot prevent progression.
The five stages have distinct gas signatures: CO₂ in Stage 1, CO + HF in Stage 2, O₂ + CO in Stage 3 (NMC only), full combustion products in Stage 4. Multi-gas sensor arrays can detect Stages 1–2, providing intervention windows that temperature sensors miss.
LFP's Stage 3 equivalent requires ~90°C more heating than NMC 811 and does not release cathode oxygen. This double advantage — higher onset plus no oxidiser supply — fundamentally changes pack propagation behaviour and reduces AIS-156 barrier material requirements.
Separator design (PE vs PP vs ceramic-coated) determines the Stage 2 temperature window. Ceramic-coated separators extend the intervention window by ~50°C, potentially providing additional BMS response time.
Production-grade thermal runaway detection requires multi-signal fusion: temperature rate-of-rise, voltage deviation, and gas sensing. Single-signal (temperature threshold) detection is late and provides minimal intervention time for NMC packs in Indian summer conditions.

Part of the thermal Series

← PreviousThermal Runaway — The 5 Stages Your BMS Is Racing Against Next →Thermal Runaway — Kinetics, ARC Testing, Propagation Modelling, and Pack Certification

Written by

Sai Chaitanya Dasari

Battery Systems Engineer · Volvo Eicher Commercial Vehicles

3+ years in commercial EV pack development. Writing about real battery engineering from the bench.

Frequently Asked Questions

What triggers thermal runaway in a lithium-ion battery?

Thermal runaway is triggered by any event that causes a cell to generate heat faster than it can dissipate it — including internal short circuits, external short circuits, mechanical abuse (crush, nail penetration), overcharge, or sustained overtemperature.

What is the Arrhenius relationship and how does it drive thermal runaway?

The Arrhenius equation describes how reaction rate increases exponentially with temperature: k = A × exp(-Ea/RT). For the exothermic decomposition reactions inside a battery (SEI, electrolyte, cathode), this means each ~10°C temperature rise roughly doubles the reaction rate and thus the heat generation rate. When heat generation rate exceeds heat dissipation rate (also temperature-dependent but more linearly), the system enters thermal runaway — a positive feedback loop where rising temperature accelerates reactions, which further raise temperature. The Frank-Kamenetskii criticality condition defines the threshold: when δ = (Ea × Q × A × exp(-Ea/RT) × V) / (RT² × λ) > 0.88 for a sphere.

How does the gas composition during thermal runaway differ between LFP and NMC?

NMC thermal runaway gas is dominated by CO₂, CO, and O₂ (from cathode decomposition), along with HF (from LiPF₆ salt), and various hydrocarbons. The oxygen content makes the gas mixture combustible and explosive. LFP thermal runaway gas is primarily CO₂ and CO, with much lower O₂ content because the olivine cathode does not release oxygen. LFP vent gases are still toxic and combustible but significantly less energetically than NMC gases. This is one mechanism (in addition to higher onset temperature) by which LFP pack fires are less intense.

What is the role of the separator in thermal runaway progression?

The separator is the physical barrier preventing electrode contact and thus internal short circuits. It is typically a microporous polyolefin film (PE, PP, or trilayer PP-PE-PP). PE starts softening at ~130°C and loses mechanical integrity above ~150°C. PP melts at ~165°C. Once the separator softens sufficiently, electrode materials can contact through the pores or through physical contact at defect sites, creating an internal short. The short converts stored electrochemical energy to heat instantaneously, driving the temperature to Stage 3 onset within seconds. Ceramic-coated separators maintain integrity to higher temperatures (~200°C), providing additional margin in the Stage 2 window.

Can accelerating rate calorimetry (ARC) predict thermal runaway in a new cell design?

ARC measures the self-heating rate of a cell or material sample under adiabatic conditions — no heat exchange with the environment. The measurement tracks when self-heating exceeds the instrument's sensitivity threshold (~0.02°C/min), maps the rate-of-rise as a function of temperature, and identifies onset temperature (T_onset), self-accelerating decomposition temperature (SADT), and the Arrhenius parameters (Ea, A) for each reaction stage. These parameters feed directly into pack-level thermal runaway propagation models (CFD or analytical). ARC is the standard characterisation tool used to parameterise thermal models for AIS-156 compliance simulations.

References

The Arrhenius Driver: Why Thermal Runaway Accelerates
Stage-by-Stage Electrochemistry
Detection Architecture for Advanced BMS
Key Takeaways
Stage 1: SEI Decomposition (80–120°C)
Stage 2: Electrolyte Decomposition and Separator Failure (120–160°C)
Stage 3: Cathode Decomposition — Point of No Return
Stage 4: Full Thermal Runaway (>200°C)
Stage 5: Propagation

TL;DRKey Points

Thermal runaway proceeds through five distinct electrochemical stages, each with a specific temperature window, gas signature, and intervention opportunity — by Stage 3 the process is irreversible.
The Arrhenius relationship causes exponential acceleration: each 10°C temperature rise roughly doubles the reaction rate, creating a positive feedback loop that cannot be stopped by cooling once the Frank-Kamenetskii criticality threshold is crossed.
LFP cathodes require ~90°C more heat to reach the Stage 3 equivalent and release far less oxygen than NMC 811, making LFP pack fires significantly less energetic and easier to contain.
Multi-sensor BMS detection (temperature rate-of-rise + CO gas + HF gas) provides up to 90 seconds more warning than temperature-threshold-only detection — critical for AIS-156 Phase 2 compliance.
Separator design (PE vs PP vs ceramic-coated) determines the Stage 2 intervention window; ceramic-coated separators extend usable response time by ~50°C.

Arrhenius temperature coefficient for SEI decomposition | 60–80 | kJ/mol activation energy

SEI exothermic heat generation rate at 120°C (NMC) | ~5–30 | W/kg

Cathode decomposition heat generation (NMC 811, 100% SOC) | ~300–700 | J/g

LFP cathode decomposition heat generation | ~200–400 | J/g (less intense, no O₂ release)

Separator integrity loss temperature (PE) | ~130–150 | °C

Peak ARC cell temperature (NMC 811, 100% SOC) | 700–950 | °C

The Arrhenius Driver: Why Thermal Runaway Accelerates

The fundamental driving force of thermal runaway is the Arrhenius relationship between temperature and reaction rate. For an exothermic reaction with activation energy Ea:

Heat generation rate Q_gen = ΔH × A × c × exp(-Ea/RT)

PYTHON

# Adiabatic rate calorimetry (ARC) data analysis
import numpy as np

def detect_arc_onset(temps_degC: list, times_min: list,
                     threshold_dT_dt: float = 0.02):
    """
    Find onset temperature in ARC experiment.
    ARC operates in heat-wait-search (HWS) mode.
    Onset = first temp where self-heating rate > threshold [deg C/min].
    """
    dT = np.diff(temps_degC)
    dt = np.diff(times_min)
    rates = dT / dt   # deg C/min

    for i, rate in enumerate(rates):
        if rate >= threshold_dT_dt:
            return temps_degC[i]   # onset temperature [deg C]
    return None

# Example ARC trace -- NMC811 charged to 4.3 V
# Typical onset: 185 deg C for NMC811 at full charge
arc_temps = [25, 60, 80, 100, 120, 140, 160, 180, 185, 190, 210, 280, 500]
arc_times  = [0,  5,  10, 15,  20,  25,  30,  35,  36,  36.5, 37,  37.5, 38]

T_onset = detect_arc_onset(arc_temps, arc_times)
print(f"ARC self-heating onset: {T_onset} deg C")
print(f"Safety margin vs operating limit (60 deg C): {T_onset - 60:.0f} deg C")

where ΔH is the reaction enthalpy, A is the pre-exponential factor, c is reactant concentration, R is the gas constant, and T is absolute temperature.

◆

Key Insight

QWhat does the Frank-Kamenetskii criticality condition mean in practical BMS terms?

Stage-by-Stage Electrochemistry

Stage 1: SEI Decomposition (80–120°C)

ROCO₂Li → Li₂CO₃ + CO₂ + hydrocarbons

Products: CO₂ (inert, non-flammable), hydrocarbons (trace), water vapour.

Stage 2: Electrolyte Decomposition and Separator Failure (120–160°C)

At ~120°C, the liquid electrolyte solvents begin to react exothermically. For EC-based electrolytes:

EC + PF₅ → CO₂ + gaseous products DEC oxidation → CO + CO₂ + hydrocarbons

The heat generation rate is ~10–100 W/kg of electrolyte. Simultaneously, LiPF₆ salt decomposes: LiPF₆ → LiF + PF₅ PF₅ + H₂O → POF₃ + 2HF

Event	Temperature	Gas Produced	Heat Rate	BMS Detectability
SEI organic decomp.	80–100°C	CO₂, trace HC	5–30 W/kg	Temperature rate-of-rise
Electrolyte oxidation	120–140°C	CO, CO₂, HC	10–100 W/kg	Gas sensor (CO, HF)
Salt decomp. (LiPF₆)	100–130°C	HF, POF₃	Moderate	Gas sensor (HF specific)
Separator softening	PE: 130°C / PP: 165°C	None	N/A	Voltage deviation (short)
Separator melting	PE: 150°C / PP: 175°C	None — cell shorts	Rapid	Voltage collapse

QWhy is HF the most dangerous gas produced during thermal runaway and how should a BMS respond to it?

Stage 3: Cathode Decomposition — Point of No Return

For NMC, once fully delithiated above approximately 160–180°C:

Li_xNi₁₋y₋zMn_yCo_zO₂ → Li_xNi₁₋y₋zMn_yCo_zO₂₋δ + δ/2 O₂

LFP cathode at Stage 3 equivalent temperatures:

LiFePO₄ → FeO + Li₃PO₄ + O₂ (this reaction requires ~270°C)

Warning

Stage 4: Full Thermal Runaway (>200°C)

Stage 5: Propagation

PYTHON

# Vent gas composition model for Li-ion cells
# Gases released during thermal runaway -- critical for safety room design

vent_gas_composition = {
    "LFP (at 200 deg C)": {
        "CO2":   45.0,   # % by volume
        "CO":    25.0,
        "H2":    15.0,
        "CH4":    8.0,
        "C2H4":  5.0,
        "HF":     2.0,   # extremely toxic
    },
    "NMC 811 (at 250 deg C)": {
        "CO2":   38.0,
        "CO":    30.0,
        "H2":    18.0,
        "CH4":    6.0,
        "C2H4":  4.0,
        "HF":     4.0,   # higher with high-voltage electrolyte
    },
}

Heat flux from runaway cell: ~1000-5000 W/cell peak

Q_peak_W = 2000.0 t_propagate = (cell_mass_g specific_heat_J_gK delta_T) / Q_peak_W return t_propagate

Cylindrical cells in module, 3mm air gap

TEXT




# Lethal concentration data (IDLH = Immediately Dangerous to Life/Health)
idlh_ppm = {"CO": 1200, "H2": 50000, "HF": 30, "CH4": 50000}

print("Vent gas composition and toxicity:")
for chemistry, gases in vent_gas_composition.items():
    print(f"\n{chemistry}:")
    for gas, pct in gases.items():
        ppm_in_1m3_vent = pct * 10000   # rough estimate: 1% ~= 10,000 ppm
        idlh = idlh_ppm.get(gas, None)
        danger = ""
        if idlh and ppm_in_1m3_vent > idlh:
            danger = f" !! EXCEEDS IDLH ({idlh} ppm)"
        print(f"  {gas:<8} {pct:.1f}%{danger}")

Heat transferred from a Stage 4 cell to its neighbours arrives via three mechanisms:

Conduction: Direct contact heat transfer — dominant for cells in tight contact
Radiation: Significant above 300°C — can heat adjacent cells through air gaps
Hot gas convection: Vented hot gas from failed cell heats adjacent cells

QHow does cell-to-cell propagation occur in Stage 5 and what determines its speed?

Detection Architecture for Advanced BMS

Baseline characterisation

Run controlled TR tests on representative cell population at 25°C and 45°C. Record temperature profiles, gas signatures, and voltage behaviour at each stage. These define sensor thresholds.

Multi-sensor array design

Minimum: per-cell temperature + CO sensor in pack headspace. Enhanced: per-cell voltage deviation monitoring + HF sensor + pack pressure sensor.

Rate-of-rise algorithms

Implement exponential filter on dT/dt per cell. Threshold: 1.5×(expected dT/dt at current load) → Level 1 alert. 5× → emergency response.

Voltage deviation detection

Compare cell voltage to EKF-predicted voltage. Persistent negative deviation >50 mV at <1C load → possible internal short → Level 2 alert.

Gas sensor integration

CO baseline monitoring with drift compensation. >50 ppm above baseline → Level 2. HF detection >5 ppm → immediate Level 3 (emergency).

Response cascade

Level 1: alert VCU, increase cooling. Level 2: isolate module, maximum cooling, alert driver. Level 3: open all contactors, activate emergency venting, send emergency signal (if telematics present).

Note

Key Takeaways

✦Key Takeaways

The Arrhenius exponential relationship between temperature and reaction rate is the physical mechanism that makes thermal runaway self-accelerating. Once above the Frank-Kamenetskii criticality temperature (~80–120°C for a typical cell), cooling intervention cannot prevent progression.
The five stages have distinct gas signatures: CO₂ in Stage 1, CO + HF in Stage 2, O₂ + CO in Stage 3 (NMC only), full combustion products in Stage 4. Multi-gas sensor arrays can detect Stages 1–2, providing intervention windows that temperature sensors miss.
LFP's Stage 3 equivalent requires ~90°C more heating than NMC 811 and does not release cathode oxygen. This double advantage — higher onset plus no oxidiser supply — fundamentally changes pack propagation behaviour and reduces AIS-156 barrier material requirements.
Separator design (PE vs PP vs ceramic-coated) determines the Stage 2 temperature window. Ceramic-coated separators extend the intervention window by ~50°C, potentially providing additional BMS response time.
Production-grade thermal runaway detection requires multi-signal fusion: temperature rate-of-rise, voltage deviation, and gas sensing. Single-signal (temperature threshold) detection is late and provides minimal intervention time for NMC packs in Indian summer conditions.

Tagsthermal runaway battery safety BMS LFP NMC SEI expert degradation cell failure electrochemistry

Part of the thermal Series

← PreviousThermal Runaway — The 5 Stages Your BMS Is Racing Against Next →Thermal Runaway — Kinetics, ARC Testing, Propagation Modelling, and Pack Certification

Written by

Sai Chaitanya Dasari

Battery Systems Engineer · Volvo Eicher Commercial Vehicles

3+ years in commercial EV pack development. Writing about real battery engineering from the bench.

Frequently Asked Questions

What triggers thermal runaway in a lithium-ion battery?

What is the Arrhenius relationship and how does it drive thermal runaway?

How does the gas composition during thermal runaway differ between LFP and NMC?

What is the role of the separator in thermal runaway progression?

Can accelerating rate calorimetry (ARC) predict thermal runaway in a new cell design?

Thermal Runaway — What Actually Happens Inside a Cell Before It Catches Fire

The Arrhenius Driver: Why Thermal Runaway Accelerates

Stage-by-Stage Electrochemistry

Stage 1: SEI Decomposition (80–120°C)

Stage 2: Electrolyte Decomposition and Separator Failure (120–160°C)

Stage 3: Cathode Decomposition — Point of No Return

Stage 4: Full Thermal Runaway (>200°C)

Stage 5: Propagation

Heat flux from runaway cell: ~1000-5000 W/cell peak

Cylindrical cells in module, 3mm air gap

Detection Architecture for Advanced BMS

Key Takeaways

Part of the thermal Series

Sai Chaitanya Dasari

Frequently Asked Questions

References

Thermal Runaway — What Actually Happens Inside a Cell Before It Catches Fire

The Arrhenius Driver: Why Thermal Runaway Accelerates

Stage-by-Stage Electrochemistry

Stage 1: SEI Decomposition (80–120°C)

Stage 2: Electrolyte Decomposition and Separator Failure (120–160°C)

Stage 3: Cathode Decomposition — Point of No Return

Stage 4: Full Thermal Runaway (>200°C)

Stage 5: Propagation

Heat flux from runaway cell: ~1000-5000 W/cell peak

Cylindrical cells in module, 3mm air gap

Detection Architecture for Advanced BMS

Key Takeaways

Part of the thermal Series

Sai Chaitanya Dasari

Frequently Asked Questions

References