17.3 Flammability, Explosions, and Toxic Dispersion Limits

17.3 Flammability, Explosions, and Toxic Dispersion Limits

Chemical process engineers routinely handle flammable solvents, liquefied fuel gases, combustible powders, and highly toxic reagents. On the NCEES PE Chemical Exam, questions in this domain require calculations of flammability envelopes, blast parameters, dust explosibility metrics, toxic occupational exposure standards, and Gaussian downwind vapor dispersion concentrations.

---

1. Flammability Fundamentals & Temperature Dependencies

Flammability Limits (LFL and UFL)

A vapor-air mixture will sustain combustion only within a specific concentration window:

• Lower Flammability Limit (LFL): The minimum volumetric concentration of combustible gas or vapor in air below which flame propagation does not occur ("too lean" to burn).
• Upper Flammability Limit (UFL): The maximum volumetric concentration of combustible gas or vapor in air above which flame propagation does not occur ("too rich" to burn due to oxygen starvation).

   Combustion Regimes in Fuel-Air Mixtures
   +------------------+-----------------------------+-------------------+
   | Lean (Non-Flam)  |      FLAMMABLE REGIME       | Rich (Non-Flam)   |
   | Insufficient     | Flame Propagates Upon       | Insufficient      |
   | Fuel Molecules   | Ignition: LFL to UFL        | Oxygen Molecules  |
   +------------------+-----------------------------+-------------------+
   0 vol%            LFL                           UFL               100 vol%

Effect of Temperature on Flammability Limits

As temperature increases, molecules possess higher thermal energy, reducing the chemical energy required from the combustion reaction to sustain the flame. Consequently, increasing temperature widens the flammable range (LFL decreases, while UFL increases). The modified Burgess-Wheeler correlations quantify this temperature effect:

$LFL(T) = LFL(25^\circ\text{C}) \left[ 1 - \frac{0.75 (T - 25)}{\Delta H_c} \right]$ $UFL(T) = UFL(25^\circ\text{C}) \left[ 1 + \frac{0.75 (T - 25)}{\Delta H_c} \right]$

Where $T$ is in $^\circ\text{C}$, and $\Delta H_c$ is the net heat of combustion (in $\text{kcal/mol}$, typically positive in this empirical formula).

Effect of Pressure on Flammability Limits

• Pressure effect on LFL: LFL is virtually independent of pressure between vacuum and moderate pressures ($P \le 50\text{ bar}$) because reaction kinetics at the lean limit are governed almost entirely by temperature and thermal quenching.
• Pressure effect on UFL: UFL increases significantly with increasing pressure (e.g., UFL of methane expands from $15\text{ vol}%$ at $1\text{ atm}$ to $> 40\text{ vol}%$ at $100\text{ bar}$) because elevated pressure accelerates multi-body radical recombination reactions.

Le Chatelier's Mixing Rule for Fuel Blends

For a mixture of $n$ combustible gases, Le Chatelier's rule predicts the mixture flammability limit:

$LFL_{\text{mix}} = \frac{1}{\sum_{i=1}^n \frac{y_i}{LFL_i}} = \frac{100}{\sum_{i=1}^n \frac{C_i}{LFL_i}}$

Where $y_i$ is the mole (or volume) fraction of combustible component $i$ on a combustible-only basis, and $C_i$ is the volume percent on a combustible-only basis.

[!IMPORTANT] Inert Gases in Le Chatelier's Rule: When inert gases (such as $\text{N}_2, \text{CO}2, \text{Ar}, \text{H}2\text{O}$) are present in the mixture, they MUST be excluded from the denominator summation. First normalize the combustible components so that $\sum y{i,\text{comb}} = 1.0$, compute $LFL{\text{comb}}$, and then divide by the total fraction of combustibles to determine the as-received mixture LFL: $LFL_{\text{mix,as-received}} = \frac{LFL_{\text{comb}}}{\sum_{i=1}^n y_{i,\text{feed}}}$

Flash Point ($T_f$), Autoignition Temperature (AIT), and MIE

• Flash Point ($T_f$): The lowest liquid temperature at which a liquid produces sufficient vapor to form an ignitable mixture with air near the liquid surface at $1.0\text{ atm}$. Thermodynamically, this is the temperature where the liquid's saturation vapor pressure equals the partial pressure corresponding to the LFL: $P^{\text{sat}}(T_f) = LFL \times P_{\text{total}}$ Testing methods: Closed-cup testers (Pensky-Martens, Tagliabue) prevent vapor loss and yield flash points $5-10^\circ\text{F}$ lower (more conservative) than open-cup testers (Cleveland Open Cup).
• Autoignition Temperature (AIT): The lowest temperature at which a vapor-air mixture spontaneously ignites without an external spark or open flame. Crucial chemical trend: Unlike flash point (which increases with molecular weight), AIT decreases as hydrocarbon molecular weight and chain length increase! Longer straight-chain alkanes produce stable free radicals at lower temperatures (e.g., Methane AIT = $580^\circ\text{C}$, Propane AIT = $450^\circ\text{C}$, n-Decane AIT = $210^\circ\text{C}$).
• Minimum Ignition Energy (MIE): The smallest electrical spark discharge energy required to ignite the most easily ignitable fuel-air mixture. Typical hydrocarbons (methane, propane, benzene) exhibit $\text{MIE} \approx 0.20 - 0.30\text{ mJ}$. In stark contrast, hydrogen has an extremely low $\text{MIE} \approx 0.017\text{ mJ}$, meaning a nearly imperceptible static discharge will ignite hydrogen.
• Limiting Oxygen Concentration (LOC): The minimum oxygen concentration (in volume percent) in a fuel-air-inert gas mixture below which combustion cannot propagate, regardless of fuel concentration. For hydrocarbon blanketing and inerting, safety guidelines mandate operating at $\text{LOC} - 2%$ (if continuously monitored) or $\text{LOC} - 4%$ (if unmonitored).

Flammability Reference Values (at $25^\circ\text{C}, 1\text{ atm}$ in Air)

---

2. Explosions and Combustible Dusts

Deflagration vs. Detonation

An explosion is a rapid expansion of gases producing a damaging pressure shock wave:

• Deflagration: A chemical reaction wave that propagates at subsonic velocity ($< 340\text{ m/s}$ in unburned gas). Heat transfer occurs via molecular conduction and mass diffusion. Peak overpressures in enclosed equipment reach $P_{\text{max}} \approx 7 - 10\text{ times initial absolute pressure}$ ($P_{\text{max}} \approx 8\text{ bar}$ from atmospheric).
• Detonation: A chemical reaction wave that propagates at supersonic velocity ($> 1500 - 3000\text{ m/s}$ relative to unburned gas). The reaction is initiated and driven by an intense, propagating shock compression wave. Peak overpressures reach $P_{\text{max}} \approx 15 - 25+\text{ times initial absolute pressure}$ ($P_{\text{max}} > 20\text{ bar}$), with destructive impulse and localized shattering.

The Cubic Law and Deflagration Index ($K_G$ and $K_{\text{St}}$)

Inside an enclosed vessel of volume $V$, the maximum rate of pressure rise during deflagration scales inversely with the vessel radius according to the Cubic Law:

$\left( \frac{dP}{dt} \right)_{\text{max}} V^{1/3} = K_G \quad (\text{Vapor/Gas}) \qquad \left( \frac{dP}{dt} \right)_{\text{max}} V^{1/3} = K_{\text{St}} \quad (\text{Combustible Dust})$

Where:

• $(dP/dt)_{\text{max}}$ = maximum rate of pressure rise ($\text{bar/s}$).
• $V$ = vessel volume ($\text{m}^3$).
• $K_G, K_{\text{St}}$ = deflagration index ($\text{bar}\cdot\text{m/s}$).

Combustible Dust Explosion Hazard Classes (NFPA 68)

Finely divided organic or metallic particulate suspended in air can explode violently. Dusts are categorized into St classes based on standardized 20-liter or $1\text{-m}^3$ spherical vessel tests:

[!NOTE] Primary vs. Secondary Dust Explosions: A small primary explosion inside process equipment (such as a bucket elevator or grinder) dislodges accumulated fugitive dust from overhead beams, piping runs, and light fixtures. When this large dispersed cloud ignites, it produces a devastating secondary dust explosion that frequently causes structural plant collapse.

---

3. Toxic Exposure Benchmarks & Standards

Chemical engineers design ventilation systems, scrubbers, and safe evacuation radii using established toxicological limits:

---

4. Toxic Vapor Dispersion: The Gaussian Plume Model

When a continuous airborne release of toxic vapor occurs from a chemical facility, the downwind atmospheric concentration field is modeled using the Gaussian Plume Dispersion Model.

                    Downwind Plume Coordinate Geometry
                                    Z (Vertical)
                                    ^
                                    |       ..---..
                                    |    .-'       '-.
      Wind Speed u ------------->   +---/--- Plume ---\---> Centerline (y=0, z=H)
                                    |    '-.       .-'
                                    |       ''---''
   Source (0, 0, H) *               |  
  ==================|=========================================> X (Downwind Distance)
  GRADE LEVEL       |                               /
                    |                              / 
                    +-----------------------------+---------> Y (Crosswind Distance)

General Gaussian Plume Equation

For a continuous point source emitting at mass rate $Q_m$ ($\text{g/s}$) at effective release height $H$ ($\text{m}$) into a steady wind speed $u$ ($\text{m/s}$):

$C(x,y,z) = \frac{Q_m}{2 \pi u \sigma_y \sigma_z} \exp\left( -\frac{y^2}{2 \sigma_y^2} \right) \left[ \exp\left( -\frac{(z - H)^2}{2 \sigma_z^2} \right) + \exp\left( -\frac{(z + H)^2}{2 \sigma_z^2} \right) \right]$

Where:

• $C(x,y,z)$ = concentration at coordinates $(x, y, z)$ ($\text{g/m}^3$).
• $x$ = downwind distance parallel to wind ($\text{m}$).
• $y$ = crosswind distance perpendicular to plume centerline ($\text{m}$).
• $z$ = vertical height above grade ($\text{m}$).
• $\sigma_y(x), \sigma_z(x)$ = Pasquill-Gifford dispersion coefficients (standard deviations of concentration spread, in $\text{m}$), which increase with downwind distance $x$ according to atmospheric stability classes (A: extremely unstable, through D: neutral, to F: moderately stable).

Ground-Level Centerline Concentration ($y = 0, z = 0$)

At ground level directly along the downwind plume centerline, the vertical terms combine due to ground reflection ($[\exp(-H^2/2\sigma_z^2) + \exp(-H^2/2\sigma_z^2)] = 2\exp(-H^2/2\sigma_z^2)$):

$C(x,0,0) = \frac{Q_m}{\pi u \sigma_y \sigma_z} \exp\left( -\frac{H^2}{2 \sigma_z^2} \right)$

For a ground-level release ($H = 0$):

$C(x,0,0) = \frac{Q_m}{\pi u \sigma_y \sigma_z}$

Conversion between Mass Concentration and Volumetric ppm

To compare dispersion model outputs ($\text{mg/m}^3$) against occupational limits ($\text{ppm}$), apply the ideal gas law at standard industrial conditions ($25^\circ\text{C}$ and $1.0\text{ atm}$, where molar volume $V_m = 24.45\text{ L/mol}$):

$\text{ppm} = \frac{C\;(\text{mg/m}^3) \times 24.45}{M}$

Where $M$ is the toxic gas molecular weight ($\text{g/mol}$).

---

5. Comprehensive Worked Numerical Example

Problem Statement

A natural gas processing facility evaluates two critical safety evaluations:

1. Fuel Blend Flammability via Le Chatelier's Rule: A process vent stream contains $50.0\text{ mol}%$ methane ($LFL_1 = 5.0\text{ vol}%$, $UFL_1 = 15.0\text{ vol}%$), $20.0\text{ mol}%$ ethane ($LFL_2 = 3.0\text{ vol}%$, $UFL_2 = 12.4\text{ vol}%$), $10.0\text{ mol}%$ propane ($LFL_3 = 2.1\text{ vol}%$, $UFL_3 = 9.5\text{ vol}%$), and $20.0\text{ mol}%$ inert nitrogen ($\text{N}2$). Calculate the combustible-basis $LFL{\text{comb}}$ and the actual as-received mixture $LFL_{\text{mix}}$.
2. Atmospheric Toxic Vapor Dispersion: An outdoor ground-level flange failure on a sulfur dioxide ($\text{SO}_2$, $M = 64.06\text{ g/mol}$) transfer header discharges vapor continuously at $Q_m = 2.50\text{ kg/s}$ ($2500\text{ g/s}$). Ambient wind speed is $u = 4.0\text{ m/s}$ under Pasquill Stability Class D (neutral). At a downwind receptor distance of $x = 1.0\text{ km}$ ($1000\text{ m}$), the dispersion coefficients are $\sigma_y = 68.0\text{ m}$ and $\sigma_z = 32.0\text{ m}$. Calculate the ground-level centerline $\text{SO}_2$ concentration in $\text{mg/m}^3$ and in $\text{ppm}$ ($25^\circ\text{C}, 1\text{ atm}$). Compare the result against $\text{ERPG-2} = 3.0\text{ ppm}$ and $\text{IDLH} = 100.0\text{ ppm}$.

---

Step 1: Multi-Component Le Chatelier Flammability Limit

First, isolate the combustible components from the inert nitrogen:

• Combustible fraction in feed: $\sum y_{\text{comb}} = 0.50 + 0.20 + 0.10 = 0.80$ ($80.0\text{ mol}%$).
• Normalize combustible mole fractions on an inert-free basis ($y_i'$): $y_{\text{CH}_4}' = \frac{0.50}{0.80} = 0.625$ $y_{\text{C}_2\text{H}_6}' = \frac{0.20}{0.80} = 0.250$ $y_{\text{C}_3\text{H}_8}' = \frac{0.10}{0.80} = 0.125$

Apply Le Chatelier's rule to the combustible fraction:

$\frac{1}{LFL_{\text{comb}}} = \sum_{i=1}^3 \frac{y_i'}{LFL_i} = \frac{0.625}{5.0\%} + \frac{0.250}{3.0\%} + \frac{0.125}{2.1\%}$ $\frac{1}{LFL_{\text{comb}}} = 0.12500 + 0.08333 + 0.05952 = 0.26785\text{ (vol\%)}^{-1}$ $LFL_{\text{comb}} = \frac{1}{0.26785} = \mathbf{3.733\text{ vol}\%}$

Now calculate the as-received mixture LFL accounting for the $20.0%$ inert $\text{N}_2$:

$LFL_{\text{mix,as-received}} = \frac{LFL_{\text{comb}}}{\sum y_{\text{comb}}} = \frac{3.733\text{ vol}\%}{0.80} = \mathbf{4.67\text{ vol}\%}$

(Interpretation: The mixture in air must reach $4.67\text{ vol}%$ total gas to contain the required $3.733\text{ vol}%$ of combustible molecules to propagate flame).

---

Step 2: Gaussian Plume Dispersion Evaluation

For a ground-level release ($H = 0$) directly along the plume centerline ($y = 0, z = 0$):

$C(1000, 0, 0) = \frac{Q_m}{\pi u \sigma_y \sigma_z}$ $\text{Denominator} = \pi \times (4.0\text{ m/s}) \times (68.0\text{ m}) \times (32.0\text{ m}) = \pi \times 8704 = 27{,}344.6\text{ m}^3/\text{s}$ $C = \frac{2500\text{ g/s}}{27{,}344.6\text{ m}^3/\text{s}} = 0.091426\text{ g/m}^3 = \mathbf{91.43\text{ mg/m}^3}$

Convert mass concentration to volumetric parts per million at $25^\circ\text{C}$ and $1.0\text{ atm}$:

$\text{ppm} = \frac{C\;(\text{mg/m}^3) \times 24.45}{M} = \frac{91.43 \times 24.45}{64.06} = \frac{2235.46}{64.06} = \mathbf{34.9\text{ ppm}}$

Evaluate against toxic criteria:

• $C = 34.9\text{ ppm} > \text{ERPG-2} = 3.0\text{ ppm}$ (Severe exposure! Individuals at $1.0\text{ km}$ will suffer irreversible respiratory impairment preventing self-rescue; evacuation zone must extend beyond $1\text{ km}$).
• $C = 34.9\text{ ppm} < \text{IDLH} = 100.0\text{ ppm}$ (Below the 30-minute immediate fatality/life-threatening ceiling).

---

6. Critical PE Exam Traps & Pitfalls

[!WARNING] Trap 1: Directly Averaging Flammability Limits Instead of Reciprocals
Never calculate a mixture LFL as a simple weighted arithmetic average ($LFL_{\text{mix}} \ne \sum y_i LFL_i$). Le Chatelier's rule is an inverse harmonic relation ($1 / \sum [y_i / LFL_i]$). An arithmetic average produces substantial non-conservative errors.

[!WARNING] Trap 2: Autoignition Temperature Trends vs. Flash Point
Candidates frequently assume that heavy oils are harder to ignite in every respect. While heavier hydrocarbons have higher flash points (lower volatility), they have significantly lower autoignition temperatures! Liquid diesel or lube oil contacting a hot uninsulated steam pipe ($T > 450^\circ\text{F}$) will spontaneously ignite, whereas light methane gas requires $> 1000^\circ\text{F}$.

[!WARNING] Trap 3: Ground Reflection Factor in Dispersion Modeling
In the Gaussian dispersion formula, ground reflection doubles the ground-level concentration, replacing the $2\pi$ in the denominator with $\pi$. Forgetting ground reflection cuts the predicted concentration in half, leading to hazardous under-prediction of toxic exposure zones.