🎓 Lesson 22
D5
Validating Simulation Results Against Hand Calculations
Checking if a computer simulation gives the same answer as a quick hand calculation to make sure it’s working correctly.
🎯 Learning Objectives
- ✓ Calculate equilibrium phase compositions using Raoult’s law for binary mixtures and compare them to simulation outputs
- ✓ Analyze discrepancies between simulated and hand-calculated dew points to diagnose EOS selection errors
- ✓ Apply reduced-property correlations to estimate compressibility factors and validate simulation Z-factor outputs
- ✓ Explain how non-idealities (e.g., hydrogen bonding, polarity) necessitate departure from ideal-gas or Raoult’s-law assumptions
- ✓ Apply the Peng–Robinson EOS in spreadsheet form to compute vapor pressure and assess simulation convergence tolerance
📖 Why This Matters
In mining and blasting engineering, simulation tools like HYSYS, Aspen Plus, or custom Python-based EOS solvers are used to model explosive gas decomposition, venting dynamics, or thermal energy release in confined blasts. But if the simulation uses an inappropriate equation of state—or misconfigured component properties—it can predict unsafe overpressures or under-predict fragmentation energy by 20–40%. Hand calculations serve as a 'sanity check' anchor: they’re transparent, traceable, and rooted in first principles. Skipping validation is like calibrating a blast seismometer without a known reference signal—results look plausible, but may be dangerously wrong.
📘 Core Principles
Validation rests on three pillars: (1) Thermodynamic consistency — ensuring simulations obey mass/energy balance and phase equilibrium criteria (e.g., fugacity equality); (2) Model hierarchy — recognizing that ideal-gas or Raoult’s law approximations are valid only at low pressures and high temperatures, while real systems require cubic or advanced EOS; (3) Error attribution — distinguishing numerical convergence artifacts (e.g., poor initial guesses, tight tolerances) from fundamental model limitations (e.g., missing association terms in ANTOINE parameters). For blasting applications, key validations include adiabatic flame temperature (AFT) prediction, detonation product composition (e.g., CO/CO₂ ratio), and post-detonation gas-phase compressibility — all sensitive to EOS choice and mixing rules.
📐 Raoult’s Law Vapor Pressure Validation
Raoult’s law provides a first-principles benchmark for vapor–liquid equilibrium (VLE) in ideal mixtures. Though limited for high-pressure or polar systems (e.g., H₂O–NH₃ in ammonium nitrate emulsions), it’s indispensable for validating simulation VLE initialization and identifying gross EOS misconfigurations.
💡 Worked Example
Problem: Given: A 60 mol% ammonium nitrate (AN) and 40 mol% fuel oil (FO) mixture at 150°C and 1 atm. Assume FO behaves as n-decane (C₁₀H₂₂). Estimate mole fraction of n-decane in vapor phase using Raoult’s law. Antoine constants for n-decane: A=6.9572, B=1565.3, C=214.5 (P in mmHg, T in °C).
1.
Step 1: Calculate vapor pressure of pure n-decane at 150°C using Antoine equation: log₁₀(P) = A − B/(T + C) → log₁₀(P) = 6.9572 − 1565.3/(150 + 214.5) = 6.9572 − 4.287 ≈ 2.670 → P = 10^2.670 ≈ 468 mmHg = 62.4 kPa.
2.
Step 2: Apply Raoult’s law: y₁·P = x₁·P₁ˢᵃᵗ → y₁ = (x₁·P₁ˢᵃᵗ)/P = (0.40 × 62.4 kPa) / 101.3 kPa ≈ 0.246.
3.
Step 3: Compare to Aspen Plus simulation result (using NRTL-RK property package): y₁ = 0.238. Absolute error = 0.008 (3.3% relative). Within acceptable tolerance (<5%) — confirms simulation initialization is physically reasonable.
Answer:
The calculated vapor-phase mole fraction of n-decane is 0.246, which falls within the acceptable validation tolerance of ±5% relative to simulation output (0.238).
🏗️ Real-World Application
At the Bingham Canyon Mine (Rio Tinto), engineers validated ANFO detonation product composition using a hand-calculated oxygen balance and simplified equilibrium code (based on NASA CEA methodology) before deploying full Aspen HYSYS batch reactor simulations. Discrepancies >8% in predicted CO yield triggered re-evaluation of the mixing rule (van der Waals vs. Wong–Sandler) and revealed improper specification of the ‘solid AN’ pseudo-component enthalpy of formation. Correcting this reduced simulated overpressure prediction error from 17% to 2.1% — directly impacting blast design safety margins for adjacent infrastructure.
🔧 Interactive Calculator
🔧 Open Thermodynamics and Equations of State Calculator📋 Case Connection
📋 Geothermal Power Plant Working Fluid Selection
Low net power output (<12% thermal efficiency) using isobutane due to poor match with 140°C geofluid