🎓 Lesson 6
D4
Vapor–Liquid Equilibrium: Raoult’s Law Limitations and Activity Models
Raoult’s Law assumes liquids mix perfectly like ideal solutions, but real chemical mixtures—especially in distillation of solvents used in mining reagents or explosives manufacturing—don’t behave that way, so we need corrections called 'activity models'.
🎯 Learning Objectives
- ✓ Calculate activity coefficients using the Wilson equation for binary solvent systems relevant to hydrometallurgical reagent recovery
- ✓ Explain why Raoult’s Law fails for strongly polar or hydrogen-bonding mixtures encountered in acid regeneration or solvent extraction circuits
- ✓ Apply NRTL parameters to predict bubble-point temperature in a ternary mixture of water–HNO₃–NH₄NO₃ used in ANFO stabilization
- ✓ Analyze VLE phase diagrams to identify azeotropes and assess feasibility of separation in solvent recovery units
📖 Why This Matters
In mining and explosives engineering, distillation is critical for recovering and purifying solvents (e.g., kerosene in froth flotation, ethanol in emulsion explosives, or nitric acid in ammonium nitrate production). Using Raoult’s Law alone leads to dangerous design errors—like undersized condensers or incorrect reflux ratios—because it ignores molecular interactions in real mixtures. Understanding its limitations and applying activity models ensures safe, energy-efficient, and compliant separation system design per ISO 2631 and IChemE Safety Guidelines.
📘 Core Principles
Raoult’s Law assumes identical intermolecular forces between all components (A–A ≈ A–B ≈ B–B), which rarely holds in practice. Deviations arise from differences in polarity, hydrogen bonding, size, and association—e.g., water–ethanol forms a minimum-boiling azeotrope at 78.2°C (95.6% ethanol), making complete separation impossible via simple distillation. Activity models introduce γᵢ—the activity coefficient—to scale Raoult’s Law: yᵢP = γᵢxᵢPᵢˢᵃᵗ. The Wilson model uses local composition concepts with temperature-dependent energy parameters; NRTL adds non-randomness via αᵢⱼ, making it superior for highly asymmetric or associating systems like HNO₃–H₂O common in explosive precursor processing.
📐 Wilson Equation for Binary Mixtures
The Wilson equation estimates activity coefficients (γ₁, γ₂) for binary mixtures using molar volumes and interaction parameters derived from VLE data. It’s widely used in process simulators (Aspen Plus®, CHEMCAD®) for preliminary design of solvent recovery columns in leach solution purification or explosive formulation plants.
💡 Worked Example
Problem: Given: binary mixture of water (1) and methanol (2) at 60°C; x₁ = 0.3, x₂ = 0.7. Pure-component molar volumes: V₁ = 18.07 cm³/mol, V₂ = 40.73 cm³/mol. Interaction parameters: Λ₁₂ = 0.296, Λ₂₁ = 1.352 (from literature at 60°C). Calculate γ₁ and γ₂.
1.
Step 1: Compute G₁₂ = exp(−Λ₁₂) = exp(−0.296) = 0.744; G₂₁ = exp(−Λ₂₁) = exp(−1.352) = 0.259
2.
Step 2: Compute denominator D = x₁ + x₂G₂₁ = 0.3 + 0.7×0.259 = 0.481; and x₂ + x₁G₁₂ = 0.7 + 0.3×0.744 = 0.923
3.
Step 3: Apply Wilson equations: lnγ₁ = −ln(x₁ + x₂G₂₁) + x₂[ G₂₁/(x₁ + x₂G₂₁) − G₁₂/(x₂ + x₁G₁₂) ] = −ln(0.481) + 0.7[0.259/0.481 − 0.744/0.923] = 0.732 + 0.7[0.538 − 0.806] = 0.732 − 0.188 = 0.544 → γ₁ = e⁰·⁵⁴⁴ = 1.72. Similarly, lnγ₂ = −ln(x₂ + x₁G₁₂) + x₁[ G₁₂/(x₂ + x₁G₁₂) − G₂₁/(x₁ + x₂G₂₁) ] = −ln(0.923) + 0.3[0.744/0.923 − 0.259/0.481] = 0.080 + 0.3[0.806 − 0.538] = 0.080 + 0.080 = 0.160 → γ₂ = e⁰·¹⁶⁰ = 1.17.
Answer:
γ₁ = 1.72, γ₂ = 1.17 — indicating strong positive deviation from Raoult’s Law for water, consistent with known hydrophilic–hydrophobic interactions in solvent recovery streams.
🏗️ Real-World Application
At a Chilean copper SX–EW facility, spent kerosene-based extractant (containing ~5 wt% Cu-loaded LIX®984N) requires thermal regeneration via vacuum distillation. Initial Raoult’s Law-based column design predicted 12 theoretical stages—but actual operation showed flooding and poor separation due to strong dipole–induced interactions between diluent and complexant. Switching to NRTL-based simulation (using parameters regressed from pilot-scale VLE data) reduced required stages to 8 and increased kerosene recovery from 82% to 96.5%, meeting ISO 14001 waste minimization targets.