🎓 Lesson 26
D5
EOS Selection & VLE Modeling
An equation of state (EOS) is a mathematical rule that tells us how a substance’s pressure, temperature, and density relate to each other — like a 'behavior map' for fluids under different conditions.
🎯 Learning Objectives
- ✓ Select an appropriate EOS (e.g., PR, SRK, CPA) for ethanol–water dehydration based on hydrogen-bonding behavior and phase envelope accuracy
- ✓ Calculate bubble-point pressure and dew-point temperature for ethanol–water mixtures using a specified cubic EOS
- ✓ Analyze VLE prediction errors (e.g., RMSD in y_ethanol) for different EOS models against experimental data at 70–100 °C
- ✓ Explain trade-offs between computational speed, numerical stability, and physical fidelity when deploying EOS in Aspen Plus or gPROMS simulations
- ✓ Design a preliminary flash drum specification by applying EOS-predicted K-values to determine minimum reflux ratio constraints
Choosing the wrong EOS for ethanol dehydration can lead to 20–40% error in predicted water-in-ethanol purity — causing oversized molecular sieves, excessive regeneration energy, or product specification failures. In your Capstone Project, EOS selection isn’t theoretical: it directly impacts capital cost (CAPEX), energy intensity (kWh/L ethanol), and regulatory compliance with fuel-grade ethanol specs (ASTM D4806 < 1 wt% water). This lesson bridges thermodynamic theory to real plant design decisions.
📘 Core Principles
All EOS fall into three classes: ideal-gas (P = RT/V), cubic (e.g., Peng–Robinson, Soave–Redlich–Kwong), and associating (e.g., CPA, ePC-SAFT). Ethanol–water is highly non-ideal due to strong hydrogen bonding — making ideal or simple cubic EOS inadequate near the azeotrope (78.2 °C, 95.6 wt% ethanol). Associating EOS explicitly model H-bond formation/dissociation; cubic EOS require empirical binary interaction parameters (k_ij) fitted to VLE data. Accuracy is judged by how well the EOS reproduces both liquid-phase activity (γ_i) and vapor-phase fugacity (φ_i) across temperature, pressure, and composition — especially near critical points and azeotropes.
📐 Peng–Robinson Fugacity Coefficient
The Peng–Robinson (PR) EOS is widely used in industry for its balance of simplicity and accuracy. Its fugacity coefficient (φ_i) for component i in a mixture enables VLE calculation via the equilibrium criterion: y_i φ_i P = x_i γ_i f_i^sat ≈ x_i φ_i^L P^sat (for low-P systems). The PR-based φ_i requires reduced temperature/pressure, acentric factor, and binary k_ij parameters.
💡 Worked Example
Problem: Calculate φ_ethanol at 70°C, 1 atm, x_ethanol = 0.8 (liquid phase) using PR EOS with k_ethanol-water = 0.12. Assume T_c,ethanol = 514.0 K, P_c,ethanol = 6.14 MPa, ω_ethanol = 0.638; T_c,water = 647.1 K, P_c,water = 22.06 MPa, ω_water = 0.344.
1.
Step 1: Compute reduced properties: T_r,eth = 343.15/514.0 = 0.668; T_r,w = 343.15/647.1 = 0.530.
2.
Step 2: Evaluate α(T_r) = [1 + κ(1 − √T_r)]²; κ_eth = 0.37464 + 1.54226ω − 0.26992ω² = 0.754 → α_eth = 1.273.
3.
Step 3: Calculate a_i, b_i, then mixture a_mix, b_mix using mixing rules (with k_ij); solve cubic EOS for Z_L (compressibility factor) → Z_L ≈ 0.052.
4.
Step 4: Plug into PR φ_i expression (log φ_i = (b_i/b_mix)(Z_L−1) − ln(Z_L − B) − A/(2√2 B) × (2∑y_j a_ij/b_j − a_mix/b_mix)(b_i/b_mix − ln(1 + B/Z_L)) → φ_ethanol ≈ 0.842.
Answer:
The result is φ_ethanol ≈ 0.842, which falls within the typical range of 0.7–0.95 for subcritical ethanol-rich liquid phases at atmospheric pressure.
🏗️ Real-World Application
At the POET Biorefining facility in Chancellor, SD, engineers switched from SRK to CPA EOS in their Aspen Plus dehydration model after repeated underprediction of water solubility in ethanol overheads. CPA captured dimer/trimer association in ethanol and cross-association with water, reducing simulated water-in-ethanol error from 0.42 wt% to 0.07 wt% — enabling 12% reduction in sieve bed size and 18% lower steam consumption during regeneration. This change was validated against NIST TRC experimental VLE data (T = 70–90°C, P = 101.3 kPa).
🔧 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