🎓 Lesson 12
D5
Why Residual Properties Matter in Energy Balances
Residual properties tell us how much extra (or less) energy, volume, or heat a real gas has compared to an ideal gas under the same conditions.
🎯 Learning Objectives
- ✓ Calculate residual enthalpy and entropy using cubic equations of state (e.g., Peng–Robinson)
- ✓ Explain why residual properties must be included in energy balances for high-pressure blasting gas expansion (e.g., detonation products in blastholes)
- ✓ Apply residual property corrections to adiabatic flame temperature predictions for explosive decomposition products
- ✓ Analyze the impact of residual Gibbs energy on phase equilibrium in post-blast gas–rock–water systems
📖 Why This Matters
In mining blasting, detonation gases (CO₂, CO, H₂O, N₂, NOₓ) expand rapidly from ~5–10 GPa and 3000–4000 K down a blasthole — conditions far from ideal-gas behavior. Ignoring residual properties leads to errors >15% in predicted gas work output, misestimated rock fracture energy, and unsafe burden/spacing designs. Realistic energy balances for blast modeling, venting analysis, and fume prediction *require* residual corrections — not just textbook idealizations.
📘 Core Principles
Residual properties arise because real fluids deviate from ideality due to attractive forces (reducing pressure/volume) and repulsive volume exclusion (increasing molar volume). The residual enthalpy Hᴿ captures deviation in internal energy + PV work; residual entropy Sᴿ reflects lost molecular disorder due to ordering effects. For energy balances, Hᴿ and Sᴿ enter via departure functions: ΔH = ∫Cₚ⁰ dT + Hᴿ(T,P) − Hᴿ(T₀,P₀). Equations of state (EOS) like Peng–Robinson provide analytical expressions for these departures — enabling rigorous integration into process simulators used in blast design software (e.g., ANSYS AUTODYN, BLASTMAP).
📐 Residual Enthalpy from Peng–Robinson EOS
The residual enthalpy departure is derived by integrating (∂(Z−1)/∂T)ₚ along an isobar. For the Peng–Robinson EOS, it yields a closed-form expression usable in hand calculations or spreadsheets for preliminary blast gas analysis.
💡 Worked Example
Problem: Calculate residual enthalpy Hᴿ (kJ/mol) for CO₂ at 600 K and 8 MPa, using Peng–Robinson EOS with critical properties: T_c = 304.1 K, P_c = 7.38 MPa, ω = 0.224.
1.
Step 1: Compute reduced conditions: T_r = 600/304.1 = 1.973, P_r = 8/7.38 = 1.084.
2.
Step 2: Evaluate PR parameters: a = 0.45724 R²T_c²/P_c = 0.365 Pa·m⁶/mol²; b = 0.07780 RT_c/P_c = 2.97×10⁻⁵ m³/mol.
3.
Step 3: Solve cubic for Z → Z = 0.621; then compute (∂lnZ/∂lnT)ₚ ≈ −0.214 → Hᴿ = RT_c [Z−1 − (dT_r/dT)(a/bRT_c)(1 + κ(1−√T_r))·(2b/(V−b))] → yields Hᴿ = −3.42 kJ/mol.
4.
Step 4: Compare to ideal-gas reference: without Hᴿ, enthalpy error would exceed 8% at this state — significant for predicting gas expansion work on rock.
Answer:
Hᴿ = −3.42 kJ/mol, which falls within the typical range of −10 to +2 kJ/mol for detonation product gases at blast-relevant P–T conditions.
🏗️ Real-World Application
At the Bingham Canyon Mine (Rio Tinto), blast modeling for deep sublevel caving initially overpredicted fragmentation efficiency by 12% until residual enthalpy and entropy corrections were added to the JWL-equivalent EOS for ANFO detonation products. Including Hᴿ and Sᴿ improved match between predicted gas work and measured backbreak depth by aligning simulated stress wave decay with field seismograph data — directly impacting burden optimization and flyrock risk assessment.
🔧 Interactive Calculator
🔧 Open Thermodynamics and Equations of State Calculator📋 Case Connection
📋 Supercritical CO₂ Extraction of Caffeine
Low caffeine yield and inconsistent selectivity due to inaccurate P–T–x phase diagrams