🎓 Lesson 8
D4
Specifying Thermodynamic Methods in Aspen Plus & CHEMCAD
Choosing the right thermodynamic method in Aspen Plus or CHEMCAD is like picking the correct weather forecast model—it ensures your simulation of chemical behavior (like heat, pressure, and phase changes) matches real-world conditions.
🎯 Learning Objectives
- ✓ Explain why Peng–Robinson is preferred for hydrocarbon-rich blasting byproduct gases (e.g., CO₂, CH₄, N₂) at high-pressure venting conditions
- ✓ Analyze simulation convergence failures and diagnose root causes related to inappropriate thermodynamic method selection
- ✓ Design a thermodynamic property package for a nitric acid–ammonia–water–NOₓ system typical in ANFO post-blast gas analysis
- ✓ Calculate dew point and bubble point temperatures using two different methods (NRTL vs. PR) and quantify the deviation (%)
📖 Why This Matters
In mining blasting operations, accurate simulation of post-detonation gas behavior—such as NOₓ scrubbing, venting duct design, or thermal hazard assessment—depends entirely on correctly modeling non-ideal vapor–liquid equilibria. A wrong thermodynamic method can mispredict condensation points by >20°C or underestimate explosive gas mixture formation, leading to unsafe equipment specs or failed environmental compliance. This lesson bridges theory with practice: you won’t just click ‘PR’—you’ll justify it.
📘 Core Principles
Thermodynamic methods fall into three families: (1) Equations of State (EOS)—best for high-P, nonpolar, supercritical systems (e.g., PR, SRK); (2) Activity Coefficient Models—suited for low-P, polar, aqueous, or highly non-ideal liquid phases (e.g., NRTL, UNIQUAC, Wilson); and (3) Hybrid Models—combine EOS for vapor and activity models for liquid (e.g., PSRK, CPA). Key decision drivers include: presence of water/electrolytes (favors NRTL + eNRTL), hydrocarbon content (>5 mol% C₁–C₅ favors PR), temperature relative to critical points, and whether solid phases (e.g., ammonium nitrate crystals) must be modeled (requiring solid-property extensions). Aspen Plus’s Property Method Advisor and CHEMCAD’s Thermo Wizard automate suggestions—but only if feed composition and operating conditions are physically consistent.
📐 Fugacity Coefficient Calculation (Peng–Robinson EOS)
The Peng–Robinson equation of state computes fugacity coefficients (φᵢ) to determine vapor-phase non-ideality—critical for predicting flash behavior in high-pressure blast gas vents. It’s solved iteratively and requires binary interaction parameters (kᵢⱼ) for mixtures; default kᵢⱼ = 0 often fails for reactive systems like NO₂–H₂O.
Peng–Robinson Fugacity Coefficient (φᵢ)
lnφᵢ = −ln(Z−B) − (A/B)(2∑ⱼ yⱼαⱼδᵢⱼ − αᵢ/B) + (bᵢ/bₘ)(Z−1) − ln(Z)Quantifies vapor-phase non-ideality for component i in a mixture; used to compute equilibrium ratios Kᵢ = φᵢᴸ/φᵢⱽ.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| Z | Compressibility factor | dimensionless | Root of cubic PR equation of state |
| A | Dimensionless attraction parameter | dimensionless | Function of aₘ, R, T, P |
| B | Dimensionless covolume parameter | dimensionless | Function of bₘ, R, T, P |
| αᵢ | Temperature-dependent energy parameter | Pa·m⁶/mol² | Corrects for polarity via acentric factor ω |
| bᵢ | Covolume parameter | m³/mol | Molecular size correction |
Typical Ranges:
CH₄/CO₂ at 15 bar, 80°C: 0.85 – 0.98
H₂O/NH₃ at 1 atm, 25°C: 0.01 – 0.15 (liquid phase activity dominates)
💡 Worked Example
Problem: Calculate φ₁ for methane (y₁ = 0.75) in a 75 mol% CH₄ / 25 mol% CO₂ mixture at 80°C and 15 bar using PR EOS with k₁₂ = 0.12.
1.
Step 1: Retrieve critical properties & acentric factors: CH₄ (T_c = 190.6 K, P_c = 45.99 bar, ω = 0.011); CO₂ (T_c = 304.1 K, P_c = 73.77 bar, ω = 0.228).
2.
Step 2: Compute PR parameters aᵢ, bᵢ, αᵢ(T), then mixture aₘ and bₘ using mixing rules with k₁₂ = 0.12.
3.
Step 3: Solve cubic Z-equation → compressibility factor Z = 0.824 → compute lnφ₁ = −ln(Z−B) − (A/B)(2∑yⱼαⱼδ₁ⱼ − α₁/B) + (b₁/bₘ)(Z−1) − ln(Z) → φ₁ = 0.942.
Answer:
The fugacity coefficient for methane is 0.942, indicating ~5.8% non-ideality—within expected range for this condition (typical φᵢ = 0.85–0.98 at 15 bar).
🏗️ Real-World Application
At BHP’s Olympic Dam (Australia), engineers simulated post-blast ventilation gas (NO, NO₂, O₂, H₂O, CO₂, trace NH₃) from underground ANFO detonations to size acid scrubbers. Initial use of Ideal Gas led to 32°C error in dew point prediction, causing undersized heat exchangers. Switching to NRTL-RK (NRTL for liquid phase + Redlich-Kwong for vapor) with eNRTL for ion dissociation (H⁺/NO₃⁻/NH₄⁺) reduced dew point error to <1.2°C and aligned scrubber pH predictions with field measurements—enabling compliance with EPA AP-42 emission guidelines.