Fugacity Coefficients from EOS for VLE Calculations - Chemical Engineering Knowledge Platform

⚠️ Why It Matters

1

Inaccurate φᵢ estimation

↓

2

Incorrect K-value prediction

↓

3

Poor column stage/energy targeting

↓

4

Excessive reboiler duty or reflux ratio

↓

5

Higher CAPEX/OPEX for separation units

↓

6

Failure to meet product purity specs

📘 Definition

The fugacity coefficient (φᵢ) is the ratio of the fugacity of component i in a mixture to its partial pressure, quantifying non-ideal behavior in vapor-phase equilibria. It is derived from an equation of state (EOS) by integrating the residual Helmholtz or Gibbs energy, and serves as the fundamental correction factor linking real-gas phase compositions to chemical potential equality at equilibrium. For vapor–liquid equilibrium (VLE), φᵢ enables rigorous calculation of K-values (Kᵢ = yᵢ/xᵢ) when activity coefficients alone are insufficient (e.g., high-pressure hydrocarbon systems).

🎨 Concept Diagram

AI-generated illustration for visual understanding

💡 Engineering Insight

Never treat φᵢ as a 'black-box output' — always back-calculate it from first-principles EOS integrals (not lookup tables) and verify sign consistency: φᵢ < 1 indicates attractive dominance (common at low Tᵣ), while φᵢ > 1 signals repulsive dominance (high Pᵣ, light gases). In sour gas dehydration, a φ_H₂O > 1.5 at 8 MPa means water is *less* likely to partition into vapor than ideal-gas models suggest — a counterintuitive but critical insight for glycol regenerator design.

📖 Detailed Explanation

Fugacity coefficient originates from the need to extend the ideal-gas law to real fluids under process conditions where intermolecular forces and molecular volume matter. For a pure component, φᵢ = exp[(1/RT) ∫₀^P (Z−1) dP], where Z is the compressibility factor — this integral captures how the EOS deviates from ideality across pressure. Engineers use this to replace partial pressure (yᵢP) with fugacity (fᵢ = φᵢ yᵢ P) in the equilibrium condition μᵢ^L = μᵢ^V.

For mixtures, φᵢ depends not only on overall P, T, and yᵢ, but also on composition-dependent EOS parameters (a_mix, b_mix) and mixing rules. The most widely used approach employs quadratic mixing for a and linear for b, with kᵢⱼ correcting for unlike-pair interactions. Because φᵢ is defined at constant T and composition, its numerical evaluation requires solving the cubic EOS for vapor-phase roots and computing partial derivatives — a task embedded in commercial simulators (Aspen HYSYS, CHEMCAD) but often misused when default kᵢⱼ are unverified.

At advanced levels, φᵢ must be reconciled with quantum-chemical insights: for hydrogen-bonding systems (e.g., MEA–CO₂), the cubic EOS fails unless augmented with association terms (CPA EOS), and φᵢ then reflects both physical and chemical contributions to nonideality. Similarly, near-critical regions (|T − T_c| < 5 K), φᵢ exhibits singular behavior requiring crossover EOS or renormalized perturbation theory — making experimental VLE data indispensable for kᵢⱼ regression and uncertainty quantification per ISO 5167 and API RP 14E guidelines.

🔄 Engineering Workflow

Step 1

Step 1: Identify system pressure, temperature, and composition bounds (feed, top, bottom streams)

→

Step 2

Step 2: Select appropriate EOS based on P, T, and component polarity (PR, SRK, CPA, GERG)

→

Step 3

Step 3: Assign critical properties (T_c, P_c, ω) and binary interaction parameters (kᵢⱼ) from trusted databases (DIPPR, NIST, GPSA)

→

Step 4

Step 4: Solve EOS for vapor-phase molar volume and compute φᵢ via departure function integration

→

Step 5

Step 5: Validate φᵢ against experimental K-data or phase envelope (e.g., bubble-point pressure residuals < 2%)

→

Step 6

Step 6: Embed φᵢ in flash calculation (Rachford–Rice + Newton–Raphson) for rigorous VLE convergence

→

Step 7

Step 7: Sensitivity-test φᵢ to kᵢⱼ and ω uncertainties; document uncertainty bands for safety-critical separations

📋 Decision Guide

Rock/Field Condition	Recommended Design Action
P < 0.5 MPa and T > T_c of all components	Use ideal-gas assumption (φᵢ ≈ 1); skip EOS-based φᵢ calculation
1.0 MPa < P < 8.0 MPa, hydrocarbon mix, no polar components	Apply Peng–Robinson EOS with literature kᵢⱼ values; validate φᵢ against Debye–Hückel–type VLE correlations
P > 10 MPa or contains H₂O, CO₂, H₂S, or alcohols	Use GERG-2008 or CPA EOS; regress kᵢⱼ using high-pressure VLE data (e.g., DDBST) before φᵢ calculation

📊 Key Properties & Parameters

Fugacity Coefficient (φᵢ)

0.1 – 5.0 (unitless) for hydrocarbons at 1–15 MPa

Dimensionless ratio of component fugacity to its partial pressure; measures deviation from ideal-gas behavior.

⚡ Engineering Impact:

Directly determines vapor-phase nonideality in K-value calculations—errors >10% cause >20% stage misestimation in distillation.

Acentric Factor (ω)

−0.3 to 0.4 (unitless) for common process fluids (e.g., methane: 0.011, water: 0.344)

Empirical parameter characterizing molecular nonsphericity and polarity, derived from vapor pressure data at Tᵣ = 0.7.

⚡ Engineering Impact:

Critical input for cubic EOS (e.g., PR, SRK); omission or misassignment causes systematic φᵢ bias >15% in polar or heavy components.

Reduced Pressure (Pᵣ)

0.1 – 3.0 (unitless) for refinery and petrochemical VLE design

Ratio of system pressure to critical pressure of the pure component (P/P_c).

⚡ Engineering Impact:

Controls EOS convergence behavior; φᵢ becomes highly sensitive above Pᵣ > 2.0—requiring EOS validation against experimental VLE data.

Binary Interaction Parameter (kᵢⱼ)

−0.2 to +0.3 (unitless) for hydrocarbon pairs; up to ±0.5 for CO₂–amine or H₂O–hydrocarbon systems

Empirical correction term in mixing rules that adjusts cross-term attraction between components i and j.

⚡ Engineering Impact:

Uncalibrated kᵢⱼ introduces >30% error in φᵢ for trace polar components—leading to false predictions of water dew point or acid gas solubility.

📐 Key Formulas

Fugacity Coefficient (Cubic EOS, component i)

ln φᵢ = (bᵢ/b)(Z−1) − ln(Z−B) − (A/B)(2∑ⱼyⱼαᵢⱼ/b − bᵢ/b)(1/B)ln((Z+B)/Z)

Peng–Robinson EOS expression for ln φᵢ in vapor phase; requires mixture compressibility root Z.

Variables:

Symbol	Name	Unit	Description
φᵢ	Fugacity coefficient of component i	dimensionless	Ratio of fugacity of component i in mixture to its fugacity in ideal gas state at same T, P
bᵢ	Cubic EOS parameter b for component i	m3/mol	Covolume parameter for component i in Peng-Robinson EOS
b	Mixture covolume parameter	m3/mol	Mixture-averaged covolume parameter, b = ∑ⱼ yⱼ bⱼ
Z	Compressibility factor	dimensionless	Real gas compressibility factor, root of cubic EOS for vapor phase
B	Reduced mixture covolume	dimensionless	B = bP/(RT), dimensionless mixture covolume
A	Reduced mixture attraction parameter	dimensionless	A = aP/(R²T²), dimensionless mixture attraction parameter
yⱼ	Mole fraction of component j in vapor phase	dimensionless	Vapor-phase mole fraction of component j
αᵢⱼ	Binary interaction parameter	dimensionless	Empirical binary interaction parameter between components i and j
R	Universal gas constant	J/(mol·K)	Ideal gas constant
T	Temperature	K	System temperature
P	Pressure	Pa	System pressure

Typical Ranges:

Natural gas processing (P = 2–8 MPa)

ln φᵢ = −1.2 to +0.8

Ammonia synthesis loop (P = 15 MPa)

ln φᵢ = −0.5 to +0.3

⚠️ |ln φᵢ| > 1.5 warrants EOS revalidation or CPA usage

🔧 Open Calculator

Binary Interaction Parameter Estimation (Lee–Kesler)

kᵢⱼ ≈ 0.0116(1 − √(T_{c,i}T_{c,j})/T_{c,mix})(1 − √(P_{c,i}P_{c,j})/P_{c,mix})

Empirical estimate for kᵢⱼ when experimental data unavailable; uses critical property geometric means.

Variables:

Symbol	Name	Unit	Description
k_ij	Binary Interaction Parameter		Empirical parameter for phase equilibrium calculations in mixtures
T_c,i	Critical Temperature of Component i	K	Critical temperature of pure component i
T_c,j	Critical Temperature of Component j	K	Critical temperature of pure component j
T_c,mix	Critical Temperature of Mixture	K	Mixture critical temperature, often estimated as mole-fraction-weighted average or using mixing rules
P_c,i	Critical Pressure of Component i	Pa	Critical pressure of pure component i
P_c,j	Critical Pressure of Component j	Pa	Critical pressure of pure component j
P_c,mix	Critical Pressure of Mixture	Pa	Mixture critical pressure, often estimated using mixing rules

Typical Ranges:

Alkane–alkane pairs

0.00 to 0.05

CO₂–n-decane

0.12 to 0.18

⚠️ Do not use for polar–nonpolar or hydrogen-bonding systems; max error >40%

🔧 Open Calculator

🏭 Engineering Example

QatarEnergy LNG Train 7 (Ras Laffan)

N/A

φ_C1

0.79

System

C₁–C₄–N₂–H₂S feed to demethanizer

φ_H2S

0.41

k_C1-H2S

0.127 (regressed from DDBST VLE data)

P_operating

6.8 MPa

T_operating

245 K

🏗️ Applications

Natural gas sweetening & dew-point control
LNG liquefaction train design
Refinery fractionation (CDU, VDU, FCC)
Hydrogen purification membranes

🔧 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]

GPSA Engineering Data Book — Gas Processors Suppliers Association

[2]

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

[3]

Thermodynamic Models: Fundamentals & Computational Aspects — Brule, A., 2021

[4]

DDBST GmbH DDBST Database — Dortmund Data Bank Software & Separation Technology

🔗 Related Knowledge

Phase Envelope Construction Using PR-SRK Hybrid Models Residual Property Calculations (Hᴿ, Sᴿ, Gᴿ) Thermodynamic Consistency Testing for Experimental VLE Data