Equilibrium Relationships in Distillation: Raoult’s Law, Henry’s Law, and Activity Coefficients

⚠️ Why It Matters

1

Incorrect equilibrium model selection

↓

2

Inaccurate VLE prediction

↓

3

Poor column stage/energy requirement estimation

↓

4

Overdesign or underdesign of distillation columns

↓

5

Increased capital and operating costs

↓

6

Process safety or product purity failure

📘 Definition

Equilibrium relationships in distillation describe the composition relationship between coexisting liquid and vapor phases at equilibrium. Raoult’s Law applies to ideal solutions where partial vapor pressure equals mole fraction times pure-component saturation pressure; Henry’s Law governs dilute solutes with linear proportionality to concentration; activity coefficients quantify deviations from ideality via γᵢ = (yᵢP)/(xᵢPᵢ^sat), enabling rigorous phase-equilibrium calculations for non-ideal mixtures.

🎨 Concept Diagram

AI-generated illustration for visual understanding

💡 Engineering Insight

Never assume ideality—even for hydrocarbons, γ₁₂ > 1.02 at 50°C for C₇/C₁₀ mixtures alters minimum reflux by 8%. Always cross-check Raoult’s Law predictions against published VLE isotherms before finalizing tray spacing or packing height. When γᵢ deviates >5% from unity, the shortcut Fenske-Underwood-Gilliland method loses >15% accuracy in stage count.

📖 Detailed Explanation

At its core, distillation relies on differences in volatility—how readily molecules escape liquid into vapor. Raoult’s Law captures this simply: for an ideal mixture, each component’s vapor pressure drops linearly with its liquid concentration. This works well for chemically similar substances like benzene–toluene, where intermolecular forces are nearly identical.

But real mixtures rarely behave ideally. Ethanol–water strongly attracts via hydrogen bonding, suppressing ethanol’s volatility—its activity coefficient γ_ethanol exceeds 6 near x=0.5 at 1 atm. Henry’s Law steps in when one component is so dilute it ‘feels’ only the solvent’s environment, not other solutes—like oxygen dissolving in cooling water. Here, proportionality holds, but the constant Hᵢ reflects solute–solvent affinity, not pure-component vapor pressure.

Advanced practice demands thermodynamic consistency: activity coefficient models (NRTL, Wilson) must satisfy the Gibbs–Duhem equation across composition; EOS-based methods (Peng–Robinson) become mandatory above 500 kPa or near critical regions. Modern design also accounts for temperature-dependent non-ideality—e.g., γ_acetone in chloroform drops from 1.8 at 30°C to 1.3 at 60°C—requiring polynomial τᵢⱼ(T) fits, not single-point regressions.

🔄 Engineering Workflow

Step 1

Step 1: Identify components, operating pressure range, and purity targets

→

Step 2

Step 2: Screen for non-ideality (polarity, H-bonding, size disparity) using UNIFAC group contribution or Margules constants

→

Step 3

Step 3: Select equilibrium model (Raoult’s → Henry’s → γ-model → EOS) based on concentration, pressure, and component classes

→

Step 4

Step 4: Regress or source reliable parameters (Pᵢ^sat, Hᵢ, τᵢⱼ, kᵢⱼ) from NIST, DIPPR, or plant-specific VLE data

→

Step 5

Step 5: Validate model against experimental bubble/dew point or azeotrope data (±0.5°C, ±2 kPa tolerance)

→

Step 6

Step 6: Integrate into column simulation (Aspen Plus, CHEMCAD) for stage count, reflux ratio, and heat duty

→

Step 7

Step 7: Perform sensitivity analysis on γᵢ uncertainty (±10%) to assess design margin for feed composition drift

📋 Decision Guide

Rock/Field Condition	Recommended Design Action
Binary mixture with ΔT_b.p. > 50°C and low polarity (e.g., n-hexane/n-octane)	Apply Raoult’s Law with Antoine equation for Pᵢ^sat; neglect activity coefficients (γᵢ ≈ 1.0).
Polar/nonpolar mixture (e.g., ethanol/water) near azeotrope or at high purity (>99.5 mol%)	Use NRTL or UNIQUAC with regressed binary parameters; validate with isobaric VLE data at 1 atm and 10–90°C.
Trace VOC (<0.01 mol%) in aqueous stream (e.g., benzene in groundwater)	Apply Henry’s Law for solute + Raoult’s Law for solvent; use tabulated Hᵢ values from NIST Chemistry WebBook.
High-pressure hydrocarbon mixture (>1000 kPa) with supercritical components (e.g., CH₄/C₂H₆/C₃H₈)	Use cubic EOS (PR or SRK) with mixing rules—Raoult’s/Henry’s laws invalid above ~500 kPa for light hydrocarbons.

📊 Key Properties & Parameters

Activity Coefficient (γᵢ)

0.1–10.0 (unitless) for most hydrocarbon and polar systems

Dimensionless factor correcting Raoult’s Law for non-ideal interactions, defined as γᵢ = (yᵢP)/(xᵢPᵢ^sat).

⚡ Engineering Impact:

Values >2 or <0.5 indicate strong positive/negative deviations—dictating need for advanced thermodynamic models (NRTL, UNIQUAC) and affecting minimum reflux ratio.

Relative Volatility (αᵢⱼ)

1.05–25 (unitless) for separable binaries; <1.05 implies high energy demand or azeotropy

Ratio of effective volatilities of two components: αᵢⱼ = (yᵢ/xᵢ)/(yⱼ/xⱼ) = (γᵢPᵢ^sat)/(γⱼPⱼ^sat).

⚡ Engineering Impact:

Directly determines minimum theoretical stages (via Fenske equation) and feasibility of separation—α < 1.1 often requires extractive or reactive distillation.

Henry’s Constant (Hᵢ)

1–10⁵ kPa for gases in water (e.g., CO₂: ~1.6×10³ kPa; O₂: ~4.3×10⁴ kPa at 25°C)

Proportionality constant relating vapor-phase partial pressure of a dilute solute to its liquid-phase mole fraction: Pᵢ = Hᵢ·xᵢ.

⚡ Engineering Impact:

Critical for modeling trace volatile organics or dissolved gases in wastewater stripping or air pollution control absorbers—errors >20% cause >30% sizing error in packed towers.

Bubble Point Pressure (P_bub)

10–3000 kPa (e.g., atmospheric crude stabilization: ~100–200 kPa; high-pressure ethylene recovery: ~2500 kPa)

Total pressure at which the first vapor forms from a subcooled liquid mixture at fixed T and x.

⚡ Engineering Impact:

Used to set reboiler duty and column pressure profile—underestimation leads to excessive vaporization and tray flooding.

📐 Key Formulas

Raoult’s Law

y_i P = x_i P_i^{sat}(T)

Relates vapor-phase mole fraction y_i to liquid-phase mole fraction x_i for ideal mixtures.

Variables:

Symbol	Name	Unit	Description
y_i	vapor-phase mole fraction of component i	dimensionless	mole fraction of component i in the vapor phase
P	total system pressure	Pa	absolute pressure of the vapor-liquid equilibrium system
x_i	liquid-phase mole fraction of component i	dimensionless	mole fraction of component i in the liquid phase
P_i^{sat}(T)	saturation vapor pressure of pure component i	Pa	vapor pressure of pure component i at temperature T

Typical Ranges:

Benzene–toluene at 1 atm

0.1–0.9 mol/mol

Light naphtha fractionation

0.05–0.95 mol/mol

⚠️ Valid only if γ_i ∈ [0.95, 1.05] and P < 500 kPa

🔧 Open Calculator

Modified Raoult’s Law

y_i P = x_i \gamma_i P_i^{sat}(T)

Corrects Raoult’s Law for non-ideal liquid-phase behavior using activity coefficient γ_i.

Variables:

Symbol	Name	Unit	Description
y_i	mole fraction of component i in vapor phase	dimensionless	fraction of component i in the vapor phase
P	total system pressure	Pa	pressure of the vapor-liquid equilibrium system
x_i	mole fraction of component i in liquid phase	dimensionless	fraction of component i in the liquid phase
γ_i	activity coefficient of component i	dimensionless	dimensionless correction factor for non-ideal behavior in the liquid phase
P_i^{sat}(T)	saturation vapor pressure of pure component i	Pa	vapor pressure of pure component i at temperature T

Typical Ranges:

Ethanol–water at x_ethanol=0.5

γ_ethanol = 6.2–7.1

Acetone–chloroform at x_acetone=0.3

γ_acetone = 1.4–1.7

⚠️ γ_i must satisfy Gibbs–Duhem; avoid extrapolation beyond fitted T/x range

🔧 Open Calculator

Henry’s Law

y_i P = x_i H_i(T)

Models vapor pressure of very dilute solutes where Raoult’s Law fails.

Variables:

Symbol	Name	Unit	Description
y_i	mole fraction of component i in vapor phase	-	dimensionless mole fraction
P	total pressure	Pa	system total pressure
x_i	mole fraction of component i in liquid phase	-	dimensionless mole fraction
H_i(T)	Henry's law constant for component i	Pa	temperature-dependent Henry's law constant