🎓 Lesson 20
D5
Applying API RP 521 with EOS-Based Flash Calculations
API RP 521 tells engineers how to safely design pressure relief systems for equipment—like tanks or reactors—that might overpressurize, and EOS-based flash calculations predict how much liquid turns into vapor during sudden pressure drops, which is critical for sizing those relief systems.
🎯 Learning Objectives
- ✓ Calculate equilibrium vapor fraction and temperature drop during adiabatic two-phase flashing using a cubic EOS
- ✓ Apply API RP 521 methodology to select and justify the governing relief scenario for a hydrocarbon storage tank
- ✓ Design a pressure relief valve (PRV) set pressure and required orifice area using EOS-derived mass flux and thermodynamic properties
- ✓ Explain the impact of non-ideal behavior (e.g., acentric factor, binary interaction parameters) on flash composition predictions
- ✓ Analyze discrepancies between ideal-gas/liquid assumptions and EOS-based results for real mixtures (e.g., LPG, condensate)
📖 Why This Matters
In mining and mineral processing, pressure vessels—such as autoclaves for pressure oxidation, solvent extraction tanks, or LNG/LPG storage at remote sites—are subject to runaway reactions, fire exposure, or blocked outlets. Undersized relief systems can lead to catastrophic rupture; oversized ones waste capital and complicate control. API RP 521 mandates scenario-based relief design, but its accuracy hinges on thermodynamic realism—especially for flashing liquids where >90% of the relieving flow may be vapor generated *in situ*. Guessing vapor fractions using ideal models fails for heavy hydrocarbons or near-critical mixtures. That’s why EOS-based flash calculations aren’t optional—they’re the technical backbone of compliant, safe, and economical relief system design.
📘 Core Principles
The lesson builds from fundamental thermodynamics: (1) A flash calculation solves for phase equilibrium at specified T, P, and overall composition by enforcing equality of fugacities across phases—a condition only rigorously satisfied with an Equation of State. (2) API RP 521 defines ‘fire case’ relief as the dominant scenario for vessels containing volatile liquids, requiring calculation of wet-bulb–limited heat absorption and resulting two-phase mass flux. (3) The ‘adiabatic flash’ model assumes no heat transfer during rapid depressurization, making it suitable for initial sizing; it relies on isenthalpic expansion from operating conditions to relief set pressure. (4) Cubic EOS (e.g., PR, SRK) provide accurate fugacity coefficients for hydrocarbons and polar components when tuned with proper mixing rules and binary interaction parameters (kij), unlike Raoult’s law or Antoine-based approximations.
📐 Isenthalpic Flash Calculation (EOS-Based)
This formula computes vapor fraction (β) and phase compositions during adiabatic depressurization. It requires iterative solution of material and energy balances coupled with EOS-derived fugacity constraints. Widely implemented in process simulators (Aspen HYSYS, PRO/II), but conceptually grounded in the Rachford-Rice equation and enthalpy balance.
💡 Worked Example
Problem: A 50 m³ propane-butane (60 mol% C3H8, 40 mol% C4H10) storage tank operates at 35°C and 1.2 MPa. Fire exposure raises vessel pressure to 1.8 MPa, then PRV opens at 1.7 MPa set pressure. Using Peng-Robinson EOS, calculate vapor fraction β and relieving temperature after isenthalpic expansion to 1.7 MPa.
1.
Step 1: Obtain inlet (operating) enthalpy h_in = −245.3 kJ/mol (from PR EOS at 35°C, 1.2 MPa)
2.
Step 2: At outlet P = 1.7 MPa, solve for T and β such that h_out = h_in, using PR EOS to compute h_L and h_V at each trial T
3.
Step 3: Converge to T ≈ 41.2°C, β = 0.382 (i.e., 38.2% vapor by mole), with vapor phase enriched in propane (82.6 mol%)
4.
Step 4: Compute mass flux G = β·ρ_V·v_sonic (using EOS-derived ρ_V and sonic velocity) → G = 1,240 kg/m²·s
Answer:
The result is β = 0.382 and relieving temperature = 41.2°C, which falls within the typical fire-case vapor fraction range of 0.2–0.5 for LPG mixtures.
🏗️ Real-World Application
At the Cadia East copper-gold autoclave facility (New South Wales, Australia), a pressure oxidation reactor handling acidic sulfide slurry and oxygen experienced unexpected overpressure during a cooling water failure. Initial relief sizing using ideal-gas flash predicted β = 0.12; EOS-based flash (SRK with kij tuning) revealed β = 0.41 due to strong non-ideality in the H2O–H2SO4–O2 system. This 3.4× increase in vapor generation rate necessitated upsizing the pilot-operated relief valve from 25 mm to 50 mm orifice and revising the flare header hydraulics—preventing potential overpressure during future upset scenarios. The revision was documented per API RP 521 Section 4.3.2 and accepted by NSW Resources Regulator during safety audit.
🔧 Interactive Calculator
🔧 Open Thermodynamics and Equations of State Calculator📋 Case Connection
📋 Ethylene Oxide Absorption Column Design
Poor CO₂ co-absorption leading to catalyst poisoning and frequent shutdowns
📋 Supercritical CO₂ Extraction of Caffeine
Low caffeine yield and inconsistent selectivity due to inaccurate P–T–x phase diagrams