BESS Container Deflagration: Explosion Vent Area Estimate

The calculation implements the NFPA 68:2018 low-pressure gas route, §§7.2.1, 7.2.6 and 7.4, at 100 kPa initial absolute pressure. It retains the general mixture coefficient, not the hydrocarbon shortcut $C=0.0223\lambda S_u$.

For mole fractions $x_i$, oxygen demand is $\nu_{O_2}=0.5(x_{H_2}+x_{CO})+2x_{CH_4}+3x_{C_2H_4}$ and $x_{st}=[1+\nu_{O_2}/0.2095]^{-1}$. Ideal gas density and sound speed use the gas–air mixture at the measured maximum-velocity equivalence ratio. Sonic mass flux is calculated with ideal isentropic flow.

With $p_0=100\,\mathrm{kPa}$ and $p_{max}$ a gauge pressure, the full low-inertia equation is

The solver updates the square-panel hydraulic diameter $D_v=\sqrt{A_{v0}/n}$, flame/vent Reynolds enhancement, obstacle multiplier and length-to-diameter correction until relative area change is below $10^{-11}$ (32 iterations maximum). Panel threshold and translating/hinged corrections follow. The same sizing kernel is used for target sizing, pressure sweeps and the NIST worked reference. Pressure inversion searches each connected valid interval separately, using 45 bisections. A continuous root must close the proposed-area relation; a skipped release interval or a downward inertia jump returns the first feasible pressure as an explicitly identified upper bound, with its actual nonzero area residual.

Henriksen et al. (2021), Tables 2, 5 and 7 (attached s6 pp.2,5) provide the two paired compositions and measurements. Sizing uses tabulated maxima; $K_G=(dp/dt)_{max}V_{test}^{1/3}$, with $V_{test}=0.020\,\mathrm{m^3}$. PCHIP interpolation draws the measured velocity curves and does not invent a kinetic mixing law. Molkov–Bragin (2015) is hydrogen-specific and is used to assess applicability, not to transfer hydrogen coefficients to mixed battery gas.

Verification includes NIST TN 2183 (2021), Appendix A3: the production solver gives 0.921 m² against its rounded 0.9 m² hinged-panel result, with an active inertia multiplier of 1.588. Its installed 1 m² wall also satisfies the published 6.9 kPa pressure limit through the production inverse. Nominal transport properties omitted by that appendix are declared in the reference case; this checks the published rounded design example, not a measured battery-container event.

• One isolated rectangular enclosure, initially at 100 kPa absolute; published gas cases are at 300 K. Custom measured properties must match the entered initial temperature and this pressure.
• Full enclosure volume participates; no credit for displaced rack volume, limited gas inventory, HVAC dilution or a partial cloud.
• Equal square clear openings on the selected roof or single end wall, unrestricted outdoor discharge, $C_d=0.70$, no duct or flame arrester and no obstructions to panel movement. A row-packing check includes the entered edge clearance between openings and at the face perimeter; reserve frames, ribs and other unusable surface with that clearance. This is a preliminary footprint check, not panel detailing.
• Flame-path assessment is separate from mounting. Full enclosure length is the conservative default even for roof mounting; a shorter served section requires a documented assessment of vent distribution. The schematic identifies the selected mounting face without prescribing individual locations.
• Initial air velocity is at most 5 m/s; no ignition by a flame jet from another enclosure. Obstacle surface area is a screening surrogate, not a resolved rack layout.
• Specific heat ratios and viscosity are explicit screening inputs; replace them with mixture transport/thermochemical data. Molar masses use standard atomic weights. Heat release chemistry and laminar flame speed are not inferred by linear blending.

Method basis reviewed 2 November 2026. Estimate only, using the named 2018 method rather than claiming compliance with a later edition. Check against the standard in force. Limit to $L/D\le5$, $S_u<3\,\mathrm{m/s}$, closed-vessel overpressure below 10 bar, panel mass at most 40 kg/m² and enclosure volume above 1 m³. This screening tool also limits obstacle surface area to twice the enclosure surface; denser congestion needs a validated analysis.

Release limits are applied at every reduced pressure: at or below 10 kPa, activation must be at least 2.4 kPa lower; above 10 kPa, activation must be strictly less than 75% of reduced pressure. This can leave disconnected valid intervals. Invalid targets have no area result, and sweeps omit invalid pressures. The upper limit is 50 kPa. Outside the valid range, pressure and its inversion residual are unavailable.

The empirical panel-inertia threshold produces a discontinuity. If the proposed area falls within that jump, no exact inverse exists: the pressure output is a conservative upper bound at the first feasible branch, and the assessment names the inertia threshold. A release-pressure gap similarly yields a labelled upper bound on the valid upper branch. A bound's area residual is retained and need not be zero. No interpolation or bisection root is claimed inside a gap; these bounds describe this empirical sizing relation, not a simulated panel-opening trajectory.

Published $P_{ex}$ tables do not explicitly label their pressure datum. They are conservatively treated as gauge pressures here; if they are absolute, this interpretation increases the estimate. Confirm the datum for a real gas report. The source's prose also differs from its tabulated velocity/rate maxima; the tables are retained. The NFPA Annex I example's later obstruction increase appears to omit the stated −0.2 term; that transfer check therefore uses its consistent first iteration. Final-result verification instead uses the independent NIST hinged-panel example.

Pressure sweeps use 121 solved points per connected interval up to the entered enclosure resistance, or 81 points on each side of an inertia threshold. Plotted lines break at release gaps and inertia jumps. Flame curves use 141 PCHIP points over measured fuel–air ratios. No CFD, transient panel-opening dynamics, external pressure peak, detonation, structural resistance calculation or partial-volume model is implied. Use qualified engineering and representative tests for congested BESS equipment, structural response, safe discharge, restraints and system integration.