Van der Waals Equation: Derivation and Limitations - Chemical Engineering Knowledge Platform

⚠️ Why It Matters

1

Inaccurate EOS selection in process simulation

↓

2

Erroneous phase envelope prediction

↓

3

Misestimated dew/bubble points

↓

4

Poor column design (e.g., reflux ratio, tray count)

↓

5

Thermodynamic inconsistency in energy balances

↓

6

Plant startup delays or off-spec product

📘 Definition

The Van der Waals equation is an empirical cubic equation of state: (P + a(n/V)²)(V − nb) = nRT, where P is pressure, V is volume, T is temperature, n is moles, R is the universal gas constant, and a and b are substance-specific constants representing intermolecular attraction and excluded molar volume, respectively. It extends kinetic theory by introducing finite molecular size and pairwise attractive forces, enabling qualitative prediction of vapor–liquid coexistence and critical phenomena.

🎨 Concept Diagram

AI-generated illustration for visual understanding

💡 Engineering Insight

Van der Waals is never used for final process design—but it *is* the diagnostic lens that reveals when a more complex EOS is truly needed. If your Peng–Robinson model gives Z = 0.82 at 50 bar and 300 K, but Van der Waals gives Z = 0.71, that 14% divergence flags significant non-idealities—likely requiring binary interaction parameters or polar corrections. Never skip this sanity check.

📖 Detailed Explanation

The Van der Waals equation begins by correcting two ideal gas assumptions: that molecules have zero volume and exert no forces on each other. The term 'nb' subtracts the volume inaccessible due to finite molecular size, while 'a(n/V)²' adds internal pressure arising from net intermolecular attraction—derived from mean-field integration over a spherical potential well. These corrections introduce nonlinearity and enable multiple real roots in the P–V isotherm below the critical temperature.

This mathematical structure yields three important physical features absent in the ideal gas law: a critical point (where (∂P/∂V)_T = (∂²P/∂V²)_T = 0), metastable regions (local max/min in isotherms), and qualitatively correct vapor–liquid equilibrium via Maxwell’s equal-area rule. Though quantitatively inaccurate, it correctly predicts that all fluids obey corresponding states when reduced by their critical properties (P_r, V_r, T_r).

Advanced treatment reveals its limitations stem from oversimplified physics: b assumes rigid-sphere repulsion only, ignoring temperature-dependent packing; a assumes isotropic, pair-wise, short-range attraction—failing for hydrogen bonding or quadrupolar interactions. Modern cubic EOS (Soave–Redlich–Kwong, Peng–Robinson) retain Van der Waals’ functional form but replace a(T) with temperature-dependent alpha functions and introduce mixing rules for mixtures—making Van der Waals the indispensable conceptual ancestor, not a competitor.

🔄 Engineering Workflow

Step 1

Step 1: Identify fluid system and operating domain (P, T, phase regime)

→

Step 2

Step 2: Retrieve pure-component a and b from NIST Chemistry WebBook or DIPPR database

→

Step 3

Step 3: Evaluate Van der Waals predictions for P–V–T and saturation properties

→

Step 4

Step 4: Compare against experimental data or high-fidelity EOS (e.g., REFPROP) to quantify error

→

Step 5

Step 5: Assess suitability for target application (e.g., compressor discharge temp vs. flash drum split)

→

Step 6

Step 6: If error exceeds ±5% in key property (e.g., Z, h_vap), escalate to advanced EOS

→

Step 7

Step 7: Document EOS selection rationale in Process Safety Information (PSI) package

📋 Decision Guide

Rock/Field Condition	Recommended Design Action
Light hydrocarbons (C₁–C₄) at low-to-moderate P (< 10 bar), T > 1.2T_c	Van der Waals may be used for preliminary screening but prefer Peng–Robinson for rigorous simulation.
Polar or associating fluids (e.g., H₂O, NH₃, alcohols) near saturation	Avoid Van der Waals entirely—use SAFT or NRTL-RK with activity coefficient coupling.
Pedagogical context or first-principles derivation of EOS concepts	Use Van der Waals as foundational teaching tool—emphasize its role in motivating modern cubic EOS.

📊 Key Properties & Parameters

a (attraction parameter)

0.14 – 5.6 Pa·m⁶/mol² (e.g., CO₂: 0.365, propane: 0.902)

Empirical constant quantifying strength of intermolecular attractive forces; proportional to square of cohesive pressure.

⚡ Engineering Impact:

Directly affects predicted vapor pressure and saturation temperature—underestimation causes overprediction of volatility and poor condenser duty estimates.

b (excluded volume parameter)

2.6 × 10⁻⁵ – 2.0 × 10⁻⁴ m³/mol (e.g., methane: 4.28×10⁻⁵, water: 3.05×10⁻⁵)

Empirical constant representing effective molar volume occupied by gas molecules themselves.

⚡ Engineering Impact:

Controls compressibility factor Z at high pressure; underestimation leads to unsafe vessel sizing and erroneous relief valve capacity calculations.

Critical Compressibility Factor (Z_c)

0.23 – 0.31 for most real fluids (e.g., n-butane: 0.274, ammonia: 0.242)

Ratio of critical pressure × critical molar volume to R × critical temperature; theoretical value for Van der Waals fluid is 0.375.

⚡ Engineering Impact:

A low Z_c deviation signals poor EOS suitability for near-critical operations like supercritical extraction or LNG liquefaction.

📐 Key Formulas

Van der Waals Equation

(P + a(n/V)²)(V − nb) = nRT

Cubic equation of state for real gases

Typical Ranges:

Ethane at 100 bar, 300 K

P = 95–105 bar (predicted)

CO₂ at critical point (304.1 K, 73.8 bar)

a = 0.365 Pa·m⁶/mol², b = 4.27×10⁻⁵ m³/mol

⚠️ Not recommended for design if |Z_error| > 3% vs. REFPROP/NIST data

🔧 Open Calculator

Critical Constants from a and b

P_c = a/(27b²), T_c = 8a/(27Rb), V_c = 3b

Derives critical properties from Van der Waals parameters

Variables:

Symbol	Name	Unit	Description
P_c	critical pressure	Pa	Pressure at the critical point
T_c	critical temperature	K	Temperature at the critical point
V_c	critical molar volume	m³/mol	Molar volume at the critical point
a	Van der Waals constant a	Pa·m⁶/mol²	Measure of attractive forces between molecules
b	Van der Waals constant b	m³/mol	Excluded molar volume per mole of molecules
R	universal gas constant	J/(mol·K)	Fundamental physical constant relating energy and temperature

Typical Ranges:

Methane

P_c = 45.4 bar (actual: 45.9 bar), T_c = 190.4 K (actual: 190.6 K)

⚠️ Predicted T_c and P_c should fall within ±2% of measured values for acceptable use

🔧 Open Calculator

🏭 Engineering Example

BASF Ludwigshafen Olefins Complex (Germany)

N/A — fluid system: ethylene/ethane refrigeration loop

Error

−5.5%

Operating Pressure

24.5 bar

REFPROP Z (NIST 10)

0.679

Operating Temperature

-28 °C

Ethylene Mole Fraction

0.72

Van der Waals Z Prediction

0.642

🏗️ Applications

Preliminary compressor sizing
Teaching thermodynamic foundations
Benchmarking new EOS implementations
Critical point estimation for safety relief studies

🔧 Try It: Interactive Calculator

🔧

Thermodynamics and Equations of State Calculator

Multi-tab calculator for Thermodynamics and Equati...

Open →

📋 Real Project Case

Liquefied Natural Gas (LNG) Train Optimization

QatarEnergy North Field Expansion – 8 MTPA LNG train

Challenge: Excessive compressor power consumption and suboptimal refrigerant blend performance

Read full case study →

Ethylene Oxide Absorption Column Design Supercritical CO₂ Extraction of Caffeine Ammonia Synthesis Loop Purge Optimization Geothermal Power Plant Working Fluid Selection

🎨 Technical Diagrams

📚 References

[1]

Perry's Chemical Engineers' Handbook, 9th Edition — McGraw-Hill Education

[2]

NIST Chemistry WebBook (Standard Reference Database 69) — National Institute of Standards and Technology

[3]

DIPPR® Project 801: Physical Property Data — AIChE Design Institute for Physical Properties

🔗 Related Knowledge

Peng–Robinson EOS for Hydrocarbon Mixtures Soave–Redlich–Kwong EOS in Refinery Simulation Cubic EOS Parameter Estimation (a, b, α)