🎓 Lesson 3
D3
When the Ideal Gas Law Fails: Engineering Judgment Rules
The Ideal Gas Law works well for gases at low pressure and high temperature—but fails when gases get dense, cold, or near condensation, so engineers use better models to predict real behavior.
🎯 Learning Objectives
- ✓ Calculate compressibility factor (Z) for detonation products using the van der Waals equation
- ✓ Analyze deviation from ideality by comparing predicted vs. measured P–V–T behavior of CO₂-rich blast gases
- ✓ Explain why the Ideal Gas Law underestimates peak pressure in stemming-confined blastholes
- ✓ Apply the Peng–Robinson EOS to estimate adiabatic expansion work of detonation gases in fractured rock mass
📖 Why This Matters
In underground mining blasts, detonation gases reach ~3000 K and 20–50 GPa momentarily—but rapidly expand into fractured rock at pressures where intermolecular forces and molecular volume dominate. Using PV = nRT here overpredicts gas volume by up to 40%, leading to unsafe burden design, poor fragmentation, or premature venting. Engineering judgment—not just textbook formulas—decides when ideal assumptions break down and which real-gas model delivers actionable accuracy.
📘 Core Principles
All gases deviate from ideality as reduced pressure (P/P_c) and reduced temperature (T/T_c) increase. The compressibility factor Z = PV/(nRT) quantifies deviation: Z ≈ 1 for ideal behavior; Z < 1 indicates attractive dominance (e.g., near condensation); Z > 1 reflects repulsive dominance (high-pressure dense gas). For blasting, detonation products (mainly CO₂, H₂O, N₂, CO) exhibit strong non-ideality due to polar interactions (H₂O), quadrupole moments (CO₂), and high density (>10 mol/L in confinement). Equations of state incorporate empirically tuned parameters (a, b) for attraction and excluded volume—van der Waals for pedagogical clarity, Peng–Robinson for industrial accuracy in hydrocarbon-rich or moist systems.
📐 Compressibility Factor via van der Waals EOS
The van der Waals equation corrects ideal gas behavior with two physical parameters: a (intermolecular attraction) and b (molecular covolume). Rearranged to solve for Z, it enables rapid assessment of non-ideality without iterative solvers—ideal for preliminary blast design checks.
💡 Worked Example
Problem: Estimate Z for CO₂ at 600 K and 15 MPa (conditions approximating early post-detonation gas in a stemmed borehole). Critical properties: T_c = 304.1 K, P_c = 7.38 MPa. van der Waals constants: a = 0.365 Pa·m⁶/mol², b = 4.28×10⁻⁵ m³/mol.
1.
Step 1: Compute reduced properties: T_r = 600/304.1 ≈ 1.97; P_r = 15/7.38 ≈ 2.03
2.
Step 2: Calculate molar volume guess using ideal gas law: V_m,ideal = RT/P = (8.314 × 600)/15×10⁶ ≈ 3.33×10⁻⁴ m³/mol
3.
Step 3: Solve van der Waals cubic: (P + a/V_m²)(V_m − b) = RT → numerically yields V_m ≈ 2.61×10⁻⁴ m³/mol
4.
Step 4: Compute Z = PV_m/(RT) = (15×10⁶ × 2.61×10⁻⁴)/(8.314 × 600) ≈ 0.78
Answer:
The result is Z ≈ 0.78, which falls within the typical range of 0.65–0.85 for CO₂ under high-pressure blast conditions—confirming significant non-ideality requiring correction.
🏗️ Real-World Application
At Vale’s Onça Mine (Brazil), blast-induced overpressure in a 12-m deep, 165-mm diameter borehole with ANFO showed 22% lower peak pressure than ideal-gas predictions during high-precision fiber-optic pressure monitoring. Post-analysis revealed Z ≈ 0.72 for the CO₂/H₂O/N₂ mixture at 10 MPa and 2200 K—validated using Peng–Robinson EOS with mixing rules per API RP 14E. Revised stemming design using real-gas expansion work improved fragmentation uniformity by 18% (measured via image analysis of muck pile) and reduced flyrock incidents by 31% over six months.