Back to Publications
Process Safety8 min read

Advanced Reactive System Screening Tool (ARSST) in Hazard Identification: Low-Phi Calorimetry, Gas Evolution & DIERS Emergency Vent Sizing

Kiran SeepanaSeptember 9, 202611 Views
Executive Summary & Scope

Detailed process safety guide on the Advanced Reactive System Screening Tool (ARSST). Covers low thermal inertia (Phi ~ 1.05) calorimetry, polynomial heating mode, P-T flow regime diagnostics (Tempered, Gassy, Hybrid), and DIERS Leung emergency relief vent sizing.

# Advanced Reactive System Screening Tool (ARSST) in Hazard Identification: Low-Phi Calorimetry, Gas Evolution & DIERS Emergency Vent Sizing

The Advanced Reactive System Screening Tool (ARSST)—developed by Fauske & Associates under the Design Institute for Emergency Relief Systems (DIERS)—is a specialized benchtop calorimeter designed for rapid thermal hazard screening and quantitative emergency pressure relief sizing (PSV and Rupture Disk sizing).

Unlike conventional Accelerating Rate Calorimetry (ARC), which requires heavy metal bombs resulting in high thermal inertia (Φ1.53.0\Phi \approx 1.5 - 3.0), the ARSST utilizes an open glass test cell suspended in a pressurized vessel, achieving an extraordinarily low thermal inertia factor (Φ1.041.10\Phi \approx 1.04 - 1.10). This low-Φ\Phi factor allows direct, uncorrected measurement of runaway self-heating rates (dT/dt)(dT/dt), gas generation rates (QgasQ_{gas}), and pressure rise rates (dP/dt)(dP/dt) required for DIERS Leung vent sizing calculations.

ARSST Low-Phi Thermal Trajectory and DIERS Relief Parameters
ARSST Low-Phi Thermal Trajectory and DIERS Relief Parameters


# 1. Operating Principles & Instrument Architecture

The ARSST system consists of a 10 mL10\text{ mL} spherical open glass test cell containing the liquid sample (510 g5 - 10\text{ g}), surrounded by a lightweight internal electric heater belt. The glass cell sits inside a 350 mL350\text{ mL} stainless steel high-pressure containment vessel rated to 70 bar70\text{ bar} (1000 psig1000\text{ psig}).

# Key Features of ARSST Operation

  1. Pressure Equalization System: The glass cell is open to the internal containment vessel volume. Nitrogen backpressure (1030 bar10 - 30\text{ bar}) is applied before testing to prevent low-temperature boiling of volatile solvents, suppressing evaporative masking.
  2. Polynomial Heating Mode: The internal heater applies a constant or polynomial power input to heat the sample at a fixed baseline rate (0.52.0 C/min0.5 - 2.0\text{ }^\circ\text{C/min}). When an exothermic reaction occurs, the self-heating rate superimposes on the background ramp.
  3. Low Phi-Factor (Φ1.05\Phi \approx 1.05): Because the glass cell wall is thin (mcell1.5 gm_{cell} \approx 1.5\text{ g}) and insulation is minimal, almost all chemical heat generation goes directly into raising the sample temperature:
Φ=1+mglassCp,glassmsampleCp,sample1.041.10\Phi = 1 + \frac{m_{glass} \cdot C_{p,glass}}{m_{sample} \cdot C_{p,sample}} \approx 1.04 - 1.10

# 2. Governing Equations & Mathematical Data Integration

# A. Power & Energy Balance in ARSST

The total non-isothermal heat balance for the ARSST sample cell during a thermal run is:

(msCp,s+mgCp,g)dTdt=qrxn(t)+qheater(t)qloss(t)\left( m_s C_{p,s} + m_g C_{p,g} \right) \frac{dT}{dt} = q_{rxn}(t) + q_{heater}(t) - q_{loss}(t)

Because the background heater power qheater(t)q_{heater}(t) is calibrated against thermal losses qloss(t)q_{loss}(t) prior to exotherm onset:

qheater(t)qloss(t)q_{heater}(t) \approx q_{loss}(t)

The net measured self-heating rate (dT/dt)meas(dT/dt)_{meas} directly mirrors the true chemical self-heating rate:

qrxn(t)(msCp,s)(dTdt)measq_{rxn}(t) \approx \left( m_s C_{p,s} \right) \cdot \left( \frac{dT}{dt} \right)_{meas}

# B. Gas Evolution Rate Integration (QgasQ_{gas})

For non-condensable gas-generating decomposition reactions, the total moles of gas evolved Ngas(t)N_{gas}(t) and volumetric gas generation rate Qgas(t)Q_{gas}(t) (m3/s\text{m}^3/\text{s} or L/min\text{L/min}) are derived from the containment vessel pressure rise rate dP/dtdP/dt:

Ngas(t)=VcontainmentRTgas[P(t)P0]N_{gas}(t) = \frac{V_{containment}}{R \cdot T_{gas}} \cdot \left[ P(t) - P_0 \right]
Qgas(t)=dNgasdtRTstdPstd=VcontainmentPstd(TstdTgas)(dPdt)Q_{gas}(t) = \frac{dN_{gas}}{dt} \cdot \frac{R \cdot T_{std}}{P_{std}} = \frac{V_{containment}}{P_{std}} \cdot \left( \frac{T_{std}}{T_{gas}} \right) \cdot \left( \frac{dP}{dt} \right)

where VcontainmentV_{containment} is the free gas headspace volume inside the containment vessel (350 mLVsample350\text{ mL} - V_{sample}).


# 3. Characterization of Runaway Flow Regimes: PTP-T Diagnostics

To apply DIERS emergency relief sizing equations correctly, the runaway mechanism must be categorized into one of three distinct thermodynamic regimes using the ARSST Pressure versus Temperature (PTP-T) plot.

ARSST P-T Diagnostic Curves for Tempered, Gassy, and Hybrid Regimes
ARSST P-T Diagnostic Curves for Tempered, Gassy, and Hybrid Regimes

# Comparison of DIERS Flow Regimes

Flow RegimePhysical MechanismPTP-T Curve CharacteristicsRelief Sizing Philosophy
Tempered SystemReaction involves a volatile solvent/reagent. Vaporization tempers runaway temperature.Follows Antoine vapor pressure curve Psat(T)P_{sat}(T). Pressure drops back to P0P_0 upon cooling.Sized using Leung Tempered Two-Phase Flow Model. Sizing governed by latent heat ΔHvap\Delta H_{vap}.
Gassy SystemDecomposition generates non-condensable gases (N2,CO2,O2\text{N}_2, \text{CO}_2, \text{O}_2) with negligible solvent volatility.Linear or step pressure increase. Permanent pressure rise ΔP\Delta P remains after cooling.Sized using DIERS Gassy Model. Sizing governed by max gas generation rate (dQgas/dt)max(dQ_{gas}/dt)_{max}.
Hybrid SystemCombined non-condensable gas generation AND volatile solvent vapor pressure.Initial gassy pressure step followed by steep exponential vapor pressure surge.Sized using DIERS Hybrid Model (combined gas + vapor mass flow).

# 4. DIERS Emergency Relief Vent Sizing & Pressure Trajectory

When an emergency pressure relief valve or rupture disk opens, the pressure inside the reactor is governed by the dynamic discharge capacity through the vent pipe.

DIERS Emergency Relief Venting Pressure Trajectory
DIERS Emergency Relief Venting Pressure Trajectory

# A. Leung vs. Fauske Equations for Tempered Two-Phase Venting

In DIERS methodology, two primary formulations are used to calculate the required vent area for tempered (vapor-controlled) two-phase runaway reactions depending on whether exact Clausius-Clapeyron thermodynamics or allowable overpressure approximations are used:

# Formulation 1: Leung Differential Equation (Exact Clausius-Clapeyron Form)

For a tempered runaway reaction, the required vent area per unit vessel volume A/VA/V (m1\text{m}^{-1}) evaluated at the relief pressure setpoint PreliefP_{relief} is expressed as:

AV=ρmG[(dTdt)PreliefTrelief(Cp,mΔvv)]\frac{A}{V} = \frac{\rho_m}{G} \left[ \frac{\left(\frac{dT}{dt}\right)_{P_{relief}}}{T_{relief}} \left( \frac{C_{p,m}}{\Delta v_v} \right) \right]

where:

  • ρm\rho_m is the homogeneous two-phase mixture density inside the vessel at relief setpoint (kg/m3\text{kg/m}^3).
  • GG is the flashing two-phase critical mass flux through the vent pipe (kg/(m2s)\text{kg/(m}^2\cdot\text{s)}).
  • (dT/dt)Prelief(dT/dt)_{P_{relief}} is the ARSST self-heating rate evaluated at the temperature corresponding to relief setpoint pressure PreliefP_{relief} (C/s^\circ\text{C/s}).
  • TreliefT_{relief} is the absolute temperature at relief setpressure (K\text{K}).
  • Cp,mC_{p,m} is the specific heat capacity of the liquid reaction mass (J/(kgK)\text{J/(kg}\cdot\text{K)}).
  • Δvv=vgvlRTreliefPreliefMwt\Delta v_v = v_g - v_l \approx \frac{R \cdot T_{relief}}{P_{relief} \cdot M_{wt}} is the specific volume change upon vaporization (m3/kg\text{m}^3/\text{kg}).

# Formulation 2: Fauske Integrated Nomograph Equation (Allowable Overpressure Form)

When an allowable overpressure ΔPover=PmaxPset\Delta P_{over} = P_{max} - P_{set} (typically 10%30%10\% - 30\% above set pressure) is specified for hand calculation:

Av=m0Cp(dTdt)set2Gchfg0(ΔPoverPset)A_v = \frac{m_0 \cdot C_p \cdot \left(\frac{dT}{dt}\right)_{set}}{2 \cdot G_c \cdot h_{fg0} \cdot \left( \frac{\Delta P_{over}}{P_{set}} \right)}

where:

  • m0=ρmVm_0 = \rho_m \cdot V is the total mass of the reaction mixture in the vessel (kg\text{kg}).
  • hfg0h_{fg0} is the latent heat of vaporization of the volatile solvent (J/kg\text{J/kg}).
  • (ΔPoverPset)\left( \frac{\Delta P_{over}}{P_{set}} \right) is the fractional overpressure above set pressure (e.g.,0.10 for 10% overpressuree.g., 0.10\text{ for }10\%\text{ overpressure}).
ℹ️ Note
Mathematical Equivalence: By Clausius-Clapeyron, hfg=TΔvv(dPdT)h_{fg} = T \cdot \Delta v_v \cdot \left(\frac{dP}{dT}\right). In Fauske's overpressure integration across a small interval ΔPover\Delta P_{over}, TsetPset(dPdT)2(ΔPoverPset)\frac{T_{set}}{P_{set}} \left(\frac{dP}{dT}\right) \approx 2 \cdot \left(\frac{\Delta P_{over}}{P_{set}}\right). Substituting this into Formulation 1 yields Formulation 2 directly!

# B. DIERS Vent Area Equation for Gassy Systems

For non-tempered gassy decompositions, tempering does not occur, and relief must discharge gas fast enough to prevent vessel overpressurization above maximum allowable working pressure (MAWPMAWP):

AV=1Ggas(msVbatch)(Qgas,maxms)\frac{A}{V} = \frac{1}{G_{gas}} \cdot \left( \frac{m_s}{V_{batch}} \right) \cdot \left( \frac{Q_{gas,max}}{m_s} \right)

where:

  • GgasG_{gas} is the ideal gas sonic choke flux (kg/(m2s)\text{kg/(m}^2\cdot\text{s)}):
Ggas=PreliefMwtγRTrelief(2γ+1)γ+1γ1G_{gas} = P_{relief} \sqrt{ \frac{M_{wt} \cdot \gamma}{R \cdot T_{relief}} \left( \frac{2}{\gamma + 1} \right)^{\frac{\gamma + 1}{\gamma - 1}} }
  • (Qgas,max/ms)(Q_{gas,max}/m_s) is the maximum mass-specific gas evolution rate (m3/(kgs)\text{m}^3/(\text{kg}\cdot\text{s})) extracted directly from ARSST pressure data.

# 5. Sample Analysis Illustration & Step-by-Step Result Derivation

To demonstrate how ARSST data directly yields an industrial emergency pressure relief vent sizing result, consider a case study of an Aqueous Diazonium Salt Decomposition in a 2,500 L2,500\text{ L} (2.5 m32.5\text{ m}^3) batch vessel.

# A. ARSST Experimental Parameters

  • Sample Mass (msm_s): 6.80 g6.80\text{ g} (0.00680 kg0.00680\text{ kg}).
  • Containment Free Volume (VcontV_{cont}): 350 mL350\text{ mL} (3.50×104 m33.50 \times 10^{-4}\text{ m}^3).
  • Pad Pressure (P0P_0): 15.0 psig15.0\text{ psig} (2.04 bar2.04\text{ bar}).
  • Measured Phi-Factor (Φ\Phi): 1.061.06 (negligible thermal dampening).
  • Target PSV Relief Setpoint (PreliefP_{relief}): 3.50 bar3.50\text{ bar} (50.7 psia50.7\text{ psia}).
  • MAWP of Plant Reactor: 4.50 bar4.50\text{ bar} (65.2 psia65.2\text{ psia}).
+---------------------------------------------------------------------------------------+
|                    ARSST SIGNAL INTEGRATION & VENT SIZING FLOW                        |
+---------------------------------------------------------------------------------------+
| 1. Log P-T Diagnostic: P increases by 18.2 bar and STAYS elevated upon cooling.       |
|    --> FLOW REGIME DIAGNOSIS = GASSY SYSTEM (Non-Condensable N2 Gas Generation).      |
| 2. Read ARSST Data at Prelief = 3.50 bar:                                             |
|    --> Temp at Relief Setpoint: Trelief = 92.0 °C (365.15 K)                          |
|    --> Pressure Rise Rate:      (dP/dt)max = 4.20 bar/min (7,000 Pa/s)                 |
+---------------------------------------------------------------------------------------+
                                           |
                                           v
+---------------------------------------------------------------------------------------+
| 3. Specific Gas Rate: Qgas/ms = (Vcont / Pstd) × (dP/dt) / ms                         |
|                       Qgas/ms = 1.85 × 10⁻⁴ m³/(kg·s)                                 |
| 4. Sonic Choke Flux:  Ggas = 185.4 kg/(m²·s)                                          |
| 5. Vent Area Ratio:   A/V = (1 / 185.4) × 1020 kg/m³ × (1.85 × 10⁻⁴) = 0.001018 m²/m³   |
+---------------------------------------------------------------------------------------+

# B. Step-by-Step Emergency Vent Sizing Result Derivation Table

Derivation StepParameterMathematical ValuePhysical & Safety Interpretation
Step 1Flow Regime DiagnosisGASSY SYSTEMPermanent non-condensable N2\text{N}_2 gas. Evaporative cooling cannot stop pressure rise.
Step 2Relief State PressurePrelief=3.50 barP_{relief} = 3.50\text{ bar} (350,000 Pa350,000\text{ Pa})Setpoint of emergency rupture disk device.
Step 3Relief State TemperatureTrelief=92.0 CT_{relief} = 92.0\text{ }^\circ\text{C} (365.15 K365.15\text{ K})Temperature of reaction mass when relief setpoint is reached.
Step 4Gas Evolution Rate(dQgasdt)max=1.85×104 m3/(kgs)\left(\frac{dQ_{gas}}{dt}\right)_{max} = 1.85 \times 10^{-4}\text{ m}^3/(\text{kg}\cdot\text{s})Maximum volumetric gas generation per unit mass derived from ARSST dP/dtdP/dt.
Step 5Sonic Gas Vent FluxGgas=185.4 kg/(m2s)G_{gas} = 185.4\text{ kg/(m}^2\cdot\text{s)}Maximum discharge mass flow capacity through relief nozzle.
Step 6Required Area Ratio (A/VA/V)AV=0.001018 m2/m3\frac{A}{V} = 0.001018\text{ m}^2/\text{m}^3 (0.0040 in2/gal0.0040\text{ in}^2/\text{gal})DIERS required relief vent area per cubic meter of batch liquid.
Step 7Final Vent Diameter (2,500 L2,500\text{ L})Avent=25.45 cm2    dvent=5.69 cmA_{vent} = 25.45\text{ cm}^2 \implies d_{vent} = 5.69\text{ cm}FINAL PLANT DECISION: Install a 2.5 inch2.5\text{ inch} (65 mm65\text{ mm} NB) rupture disk on the reactor.

# 6. Where to Use ARSST in Process Safety Engineering

  1. DIERS Emergency Relief Valve & Rupture Disk Sizing: Provide primary experimental inputs for relief valve nozzle and rupture disk orifice sizing for reactors, distillation columns, and storage vessels.
  2. Screening Reactive Chemical Incompatibilities: Rapidly evaluate binary compatibility matrices (e.g., acid + base, oxidizer + solvent, catalyst + impurity) within 34 hours3 - 4\text{ hours}.
  3. Effluent Handling System Sizing: Sizing quench tanks, knockout drums, and flare stacks based on two-phase vent discharge rates.
💡 Pro Tip
Summary Checklist for ARSST Vent Sizing: - [x] Verify pad pressure P0P_0 is set high enough to suppress solvent boiling prior to exotherm onset. - [x] Check PTP-T plot to determine if system is Tempered, Gassy, or Hybrid. - [x] Read self-heating rate (dT/dt)Prelief(dT/dt)_{P_{relief}} at the exact temperature corresponding to relief setpoint pressure. - [x] Apply Leung Omega method to calculate two-phase flux GG and specify required vent area AventA_{vent}.
ARSSTLow Phi CalorimetryDIERS Vent SizingEmergency Relief SystemFlow Regime DiagnosticsTempered Gassy HybridLeung EquationProcess Safety
Comments (0)

Discussion

Please Log In to participate in the technical discussion.

No comments posted yet. Be the first to share your input!