🎓 Lesson 15
D5
eNRTL Theory: Local Composition and Long-Range Electrostatic Terms
eNRTL is a thermodynamic model that predicts how ions and molecules interact in salty water solutions by combining local mixing effects with the long-distance pull of electric charges.
🎯 Learning Objectives
- ✓ Explain the physical origin and relative contribution of local-composition versus long-range electrostatic terms in eNRTL
- ✓ Calculate activity coefficients for ionic species in a mixed-electrolyte solution using eNRTL equations and published binary interaction parameters
- ✓ Analyze convergence behavior and parameter sensitivity in eNRTL-based process simulations of acid leaching or brine purification
- ✓ Apply eNRTL to design a solvent extraction feed conditioning step by predicting free acid and metal speciation at varying ionic strength
📖 Why This Matters
In mining hydrometallurgy—like copper SX-EW, lithium brine processing, or uranium heap leaching—accurate prediction of ion speciation, solubility limits, and phase equilibria is critical for reagent efficiency, precipitate control, and corrosion mitigation. Traditional models (e.g., UNIQUAC, ideal-solution) fail dramatically in high-ionic-strength aqueous electrolytes. eNRTL is the industry-standard thermodynamic framework embedded in commercial simulators (Aspen Plus®, HSC Chemistry®, OLI Systems) for modeling these systems reliably—making it essential for designing and optimizing real-world extractive processes.
📘 Core Principles
eNRTL builds on three foundational ideas: (1) Local composition: Molecules near an ion are arranged non-randomly due to polarity and hydration shells—captured via modified NRTL energy parameters (α_ij, τ_ij). (2) Long-range electrostatics: Ion–ion Coulombic forces dominate at larger distances and are modeled using a Debye–Hückel limiting law extension (with McMillan–Mayer expansion), scaled by ionic strength and dielectric constant. (3) Reference state consistency: Neutral molecules use symmetric (pure-component) reference states; ions use unsymmetric (infinite-dilution) reference states—requiring careful handling of charge balance and electroneutrality constraints. The total excess Gibbs energy is the sum of both contributions, ensuring thermodynamic consistency across wide concentration ranges (0–6 mol/kg H₂O).
📐 eNRTL Total Excess Gibbs Energy
The eNRTL model expresses the dimensionless excess Gibbs energy (gᴱ/RT) as the sum of local-composition and long-range electrostatic terms. The local term uses NRTL-style interaction parameters between species (molecular or ionic), while the long-range term depends on ionic strength, temperature, and solvent dielectric properties. This decomposition allows robust fitting and extrapolation where pure-electrolyte data is sparse.
💡 Worked Example
Problem: Calculate the long-range contribution (gᴱ,LR/RT) for a 2.5 mol/kg aqueous NaCl solution at 25°C (εᵣ = 78.4, ρ = 0.997 g/cm³). Assume Debye screening parameter A_Φ = 1.173 (unitless, for water at 25°C).
1.
Step 1: Compute ionic strength I = ½ Σ c_i z_i² = ½[(2.5)(+1)² + (2.5)(−1)²] = 2.5 mol/kg
2.
Step 2: Apply extended Debye–Hückel term: gᴱ,LR/RT = −A_Φ · I^(3/2) / (1 + B·I^(1/2)) — but for the simplified limiting-law form used in many eNRTL implementations: gᴱ,LR/RT ≈ −A_Φ · I^(1/2)
3.
Step 3: Evaluate: gᴱ,LR/RT ≈ −1.173 × √2.5 = −1.173 × 1.581 = −1.855
Answer:
The long-range electrostatic contribution is −1.855 (dimensionless). This dominates over the local term (~−0.3 to −0.8 for same system), confirming its critical role in concentrated electrolytes.
🏗️ Real-World Application
At the Salar de Atacama lithium brine operation (SQM, Chile), eNRTL is used in Aspen Plus® to simulate the multi-stage evaporation–crystallization train for LiCl and NaCl separation. The model accurately predicts K⁺/Mg²⁺/Li⁺/Cl⁻/SO₄²⁻ speciation, gypsum (CaSO₄·2H₂O) scaling onset, and lithium recovery >92%—replacing empirical charts and reducing pilot-plant iterations by 70%. Key eNRTL parameters were regressed from 350+ experimental osmotic coefficient and solubility points measured across 5–80°C and I = 0.5–5.2 mol/kg.
🔧 Interactive Calculator
🔧 Open Thermodynamics & Process Simulation Calculator📋 Case Connection
📋 Supercritical Fluid Extraction (SFE) Process Design for Caffeine Recovery
Low selectivity and high CO₂ consumption due to poor phase behavior prediction