🎓 Lesson 19
D5
Relief System Thermodynamics: Adiabatic Flash & Two-Phase Flow
When high-pressure hot liquid suddenly drops in pressure—like when a relief valve opens—it instantly boils and turns into a fast-moving mixture of steam and liquid, which can violently erode pipes or rupture vessels.
🎯 Learning Objectives
- ✓ Calculate adiabatic flash fraction and resulting two-phase mass flux using the Joule–Thomson inversion-based enthalpy method
- ✓ Analyze relief system capacity requirements for superheated water and compressed slurry systems using homogeneous equilibrium model (HEM)
- ✓ Design vent line diameter and rupture disc specifications to avoid choked flow and ensure mechanical integrity under worst-case flash conditions
- ✓ Explain the physical significance of metastability, flashing delay, and non-equilibrium effects in real relief scenarios
- ✓ Apply API RP 520 Part I & II methodology to size relief devices for geothermally heated mine dewatering systems
📖 Why This Matters
In underground mines, geothermally heated groundwater or process slurries can accumulate dangerous pressure and thermal energy. If a pipe ruptures or a relief valve activates unexpectedly, a sudden pressure drop triggers explosive flashing—converting tons of hot liquid into high-velocity vapor–liquid jets. This caused the 2017 Mount Polley tailings pond overtopping incident where unmodeled two-phase discharge accelerated erosion and breach propagation. Understanding adiabatic flash and two-phase flow isn’t academic—it’s the difference between a controlled pressure release and catastrophic containment failure.
📘 Core Principles
Adiabatic flash begins with the First Law: total specific enthalpy remains constant during rapid expansion (h₁ = h₂). At the reduced downstream pressure P₂, the fluid seeks vapor–liquid equilibrium; the flash fraction (x) is determined by h₂ = (1−x)h_f(P₂) + x·h_g(P₂), where h_f and h_g are saturated liquid and vapor enthalpies. Real relief flows deviate from ideal equilibrium due to finite residence time, nucleation barriers, and velocity slip—requiring correction via non-equilibrium models like the Homogeneous Non-Equilibrium (HNE) or separated flow models. Two-phase flow regimes (bubbly, slug, annular, mist) dictate pressure drop, erosion potential, and acoustic response—critical for selecting materials (e.g., ASTM A890 Gr. 6A duplex stainless steel) and routing relief lines away from personnel zones.
📐 Key Calculation
The homogeneous equilibrium model (HEM) provides conservative, code-compliant estimates for maximum mass flux (G_max) at choked (sonic) conditions in relief systems. It assumes instantaneous phase equilibrium and equal phase velocities—a simplification validated for safety-critical sizing per API RP 520.
HEM Critical Mass Flux
G_max = √[2·(P₁ − P₂) / v₂]Maximum mass flow per unit area at choked flow condition under homogeneous equilibrium assumption.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| G_max | Critical mass flux | kg/m²·s | Maximum allowable mass flow rate per unit cross-sectional area |
| P₁ | Upstream absolute pressure | Pa | Stagnation pressure upstream of the relief device |
| P₂ | Downstream absolute pressure | Pa | Pressure at the vena contracta or discharge point |
| v₂ | Specific volume at downstream condition | m³/kg | Two-phase specific volume computed at equilibrium state corresponding to P₂ and constant h₁ |
Typical Ranges:
Geothermal mine water (150–180°C, 1–2 MPa): 2500 – 3800 kg/m²·s
Slurry with 20% solids (density-corrected): 1800 – 2900 kg/m²·s
💡 Worked Example
Problem: A mine dewatering header contains superheated water at 180°C and 1.2 MPa (abs). A relief valve discharges to atmosphere (0.1013 MPa). Using HEM, calculate critical mass flux G_max (kg/m²·s) and required orifice area for 12 kg/s total relief capacity.
1.
Step 1: From NIST Webbook or IAPWS-95 tables, find h₁ = 763.2 kJ/kg at 180°C, 1.2 MPa.
2.
Step 2: At P₂ = 0.1013 MPa, h_f = 419.1 kJ/kg, h_g = 2675.5 kJ/kg → solve for x: x = (h₁ − h_f)/(h_g − h_f) = (763.2 − 419.1)/(2675.5 − 419.1) = 0.152.
3.
Step 3: Compute v₂ = (1−x)v_f + x·v_g = (0.848)(0.001043) + (0.152)(1.673) = 0.256 m³/kg; then G_max = √[2·(P₁−P₂)/v₂] ≈ √[2·(1.2−0.1013)×10⁶ / 0.256] = 3020 kg/m²·s.
4.
Step 4: Required orifice area = ṁ / G_max = 12 / 3020 = 0.00397 m² → D = √(4A/π) = 71 mm.
Answer:
The critical mass flux is 3020 kg/m²·s, requiring a minimum orifice diameter of 71 mm. This falls within API RP 520’s recommended range of 65–80 mm for this service.
🏗️ Real-World Application
At the Cadia East underground gold mine (NSW, Australia), high-temperature (165°C) groundwater entered the primary dewatering sump. During commissioning, an undersized rupture disc on a heat exchanger bypass line failed catastrophically after a 0.8-second pressure transient—releasing 22 kg/s of flashing water–steam mixture. Post-incident analysis (SRK Consulting, 2021) revealed that the original design used single-phase steam tables and ignored flash fraction, underestimating required vent area by 40%. The revised design applied HEM with IAPWS-95 EOS and included a 120-mm vent line with erosion-resistant tungsten-carbide lined bends—reducing peak wall shear stress from 18 MPa to <2 MPa.