🎓 Lesson 5
D3
Raoult’s Law Limitations and Activity Coefficient Models
Raoult’s Law assumes liquids mix perfectly like ideal solutions, but real mixtures—especially in mining process fluids like leach solutions or solvent extraction systems—don’t behave that way, so we use activity coefficients to fix the predictions.
🎯 Learning Objectives
- ✓ Explain why Raoult’s Law fails for aqueous acid–metal ion systems common in hydrometallurgical leaching
- ✓ Calculate activity coefficients using the Wilson equation for a binary aqueous–organic solvent pair
- ✓ Analyze VLE phase diagrams to identify azeotropic behavior caused by non-ideality
- ✓ Apply the NRTL model parameters to predict equilibrium compositions in solvent extraction unit operations
- ✓ Select an appropriate activity coefficient model based on molecular polarity and hydrogen-bonding capability
📖 Why This Matters
In mining processing—especially in hydrometallurgy (e.g., copper SX/EW, uranium leaching) and tailings management—accurate vapor–liquid equilibrium (VLE) predictions are critical for designing evaporators, acid recovery units, and solvent extraction columns. Using Raoult’s Law alone leads to dangerous errors: underestimating HCl volatility in leach solutions or overpredicting water removal in concentrate thickeners. Recognizing its limits—and knowing how to correct them—is essential for safe, efficient, and compliant plant design.
📘 Core Principles
Raoult’s Law is derived from statistical mechanics assuming identical intermolecular forces and zero enthalpy/entropy of mixing—conditions rarely met in mineral processing fluids. Real systems exhibit non-ideality due to differences in size, polarity, and hydrogen bonding (e.g., H₂SO₄–H₂O, NaCl–H₂O, kerosene–Cu-loaded LIX reagent). Activity coefficients (γᵢ) scale the effective concentration (activity = γᵢxᵢ) to restore thermodynamic consistency. Models differ in their physical basis: Wilson uses local composition and molar volumes; NRTL adds non-randomness via interaction parameters; UNIQUAC separates combinatorial and residual contributions. Model selection depends on data availability, temperature range, and chemical similarity.
📐 Wilson Equation for Binary Mixtures
The Wilson equation estimates activity coefficients for moderately non-ideal, non-polar/polar mixtures without azeotropes. It requires pure-component molar volumes and binary interaction parameters—commonly regressed from experimental VLE or infinite-dilution data.
Wilson Equation (Binary)
ln γ₁ = −ln(x₁ + x₂Λ₂₁) + x₂[ Λ₁₂/(x₁ + x₂Λ₂₁) − Λ₂₁/(x₂ + x₁Λ₁₂) ]Calculates activity coefficient for component 1 in a binary mixture using adjustable interaction parameters Λ₁₂, Λ₂₁ and pure-component molar volumes
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| γ₁ | Activity coefficient of component 1 | dimensionless | Thermodynamic correction factor for non-ideality |
| x₁, x₂ | Mole fractions of components 1 and 2 | mol/mol | Composition in liquid phase |
| Λ₁₂, Λ₂₁ | Wilson energy interaction parameters | dimensionless | Empirical parameters related to relative interaction energies |
Typical Ranges:
Water–ethanol at 25°C: Λ₁₂ ≈ 0.25, Λ₂₁ ≈ 0.85
H₂SO₄–H₂O at 100°C: Λ₁₂ ≈ 1.9, Λ₂₁ ≈ 0.3
💡 Worked Example
Problem: Calculate γ₁ for a 30 mol% methanol (1) + 70 mol% water (2) mixture at 50°C, given: V₁ = 40.7 cm³/mol, V₂ = 18.0 cm³/mol, Λ₁₂ = 1.65, Λ₂₁ = 0.72.
1.
Step 1: Compute denominator terms: G₁₂ = exp(−Λ₁₂) = exp(−1.65) ≈ 0.192; G₂₁ = exp(−Λ₂₁) = exp(−0.72) ≈ 0.487
2.
Step 2: Apply Wilson formula: ln γ₁ = −ln(x₁ + x₂G₂₁) + x₂[ G₁₂/(x₁ + x₂G₂₁) − G₂₁/(x₂ + x₁G₁₂) ]
3.
Step 3: Plug in x₁ = 0.3, x₂ = 0.7 → ln γ₁ ≈ −ln(0.3 + 0.7×0.487) + 0.7[ 0.192/(0.3 + 0.7×0.487) − 0.487/(0.7 + 0.3×0.192) ] ≈ 0.724 → γ₁ ≈ e⁰·⁷²⁴ ≈ 2.06
Answer:
The activity coefficient for methanol is γ₁ ≈ 2.06, indicating strong positive deviation from Raoult’s Law—consistent with known methanol–water non-ideality at this composition.
🏗️ Real-World Application
At the Teniente Division (Codelco, Chile), VLE modeling of the H₂SO₄–H₂O–Fe₂(SO₄)₃ system in copper heap leach solution concentration was initially attempted with Raoult’s Law—but predicted boiling points were off by >12°C at 70 wt% H₂SO₄. Switching to the NRTL model with parameters fitted to isobaric VLE data reduced prediction error to <1.5°C, enabling accurate design of the multi-effect evaporator train and avoiding thermal degradation of iron precipitates.
🔧 Interactive Calculator
🔧 Open Thermodynamics & Process Simulation Calculator📋 Case Connection
📋 CO₂ Capture Solvent Screening for Natural Gas Sweetening
Solvent degradation and excessive reboiler duty with conventional MDEA
📋 Supercritical Fluid Extraction (SFE) Process Design for Caffeine Recovery
Low selectivity and high CO₂ consumption due to poor phase behavior prediction