Interphase Mass Transfer: Two-Film Theory and Penetration Model

Typical Scale

Industrial absorbers: 1–4 m diameter, 10–40 m tall

Key Standard

AIChE Design Guidelines (2021), ISO 20815:2018 (Petroleum industry mass transfer)

Industry Impact

Accounts for ~65% of capital cost in post-combustion CO₂ capture plants

Computational Tool

Aspen Plus RADFRAC + rate-based models (e.g., EMPIRICAL, EQUILIBRIUM, or custom Fortran subroutines)

⚠️ Why It Matters

1

Inaccurate interphase mass transfer prediction

↓

2

Underdesigned absorber column height

↓

3

Insufficient CO₂ removal in amine scrubbers

↓

4

Exceedance of emissions limits

↓

5

Regulatory noncompliance and operational shutdown

📘 Definition

Interphase mass transfer describes the net movement of chemical species across the physical boundary between two immiscible phases (e.g., gas–liquid, liquid–liquid), governed by concentration gradients and interfacial resistance. The Two-Film Theory models this as diffusion through stagnant, laminar films on either side of the interface, while the Penetration Model treats the interface as transient and dynamic, where fluid elements are exposed to bulk phase for finite time intervals before being replaced.

🎨 Concept Diagram

AI-generated illustration for visual understanding

💡 Engineering Insight

Never assume the liquid film is the sole resistance — in modern amine-based CO₂ capture, gas-film resistance dominates above 0.1 MPa due to high partial pressure and low k_G. Always compute both k_G and k_L from first-principles correlations tied to your actual operating Re and geometry, not textbook averages.

📖 Detailed Explanation

At its core, interphase mass transfer arises because molecules naturally move from regions of high concentration to low concentration — a process called diffusion. When two phases meet (like air and water), molecules don’t instantly mix across the boundary; instead, they must cross thin, relatively still layers adjacent to the interface — these are the 'films' in the Two-Film Theory. Engineers use this model because it simplifies complex turbulent flow into tractable resistances that add in series, much like electrical resistors.

The Penetration Model improves realism by recognizing that fluid elements at the interface aren’t static — they’re constantly exchanged due to turbulence. In this view, mass transfer occurs during brief exposure times (τ), after which fresh fluid replaces the element. This leads to a time-dependent solution where average flux scales with √(D/τ), making it especially useful for agitated vessels and spray towers where residence time distributions matter more than steady-state films.

Advanced treatment merges both models via surface renewal theory (Danckwerts), where penetration time τ is linked to turbulent eddy frequency — yielding k_L ∝ √(D_L·E), with E as surface renewal rate. Modern CFD-coupled population balance models (PBM-CFD) now resolve droplet/bubble size evolution *and* local film dynamics simultaneously, enabling predictive design of next-generation contactors like falling-film microstructured reactors used in pharmaceutical API purification.

🔄 Engineering Workflow

Step 1

Step 1: Identify limiting phase (gas-side vs. liquid-side resistance via k_G/k_L ratio and Hatta analysis)

→

Step 2

Step 2: Measure or estimate equilibrium solubility (Henry’s law constant H or distribution coefficient K_D)

→

Step 3

Step 3: Characterize hydrodynamics (Re, Fr, Weber numbers) to select appropriate film or penetration time scale

→

Step 4

Step 4: Calculate local mass transfer coefficients using correlations (e.g., Onda, Lockett, or Calderbank)

→

Step 5

Step 5: Integrate across equipment geometry using NTU–HTU or stage-by-stage methods

→

Step 6

Step 6: Validate with pilot-scale mass balance and tracer studies (e.g., CO₂ absorption with pH titration)

→

Step 7

Step 7: Iterate design for pressure drop, entrainment, and fouling margins

📋 Decision Guide

Rock/Field Condition	Recommended Design Action
Gas–liquid system with irreversible fast reaction (e.g., SO₂ in NaOH)	Use packed column with high-specific-area structured packing (e.g., Sulzer BX); optimize liquid distribution to maximize interfacial renewal.
Liquid–liquid extraction with low interfacial tension (<5 mN/m) and high density difference	Select mixer-settler over centrifugal extractor; control agitation intensity to avoid stable emulsion formation.
High-viscosity solvent (μ > 50 cP) and low diffusivity (D < 1×10⁻¹⁰ m²/s)	Prefer rotating disc contactor (RDC) or pulsed sieve-plate column to enhance film renewal; avoid packed beds.

📊 Key Properties & Parameters

Overall Mass Transfer Coefficient (K_La)

0.02–0.5 s⁻¹ for packed absorption columns; 0.1–2.5 s⁻¹ for stirred-tank extractors

Volumetric mass transfer coefficient quantifying the combined resistance to transport across both phases, defined as K_La = 1/(1/k_L + H/k_G) for gas–liquid systems.

⚡ Engineering Impact:

Directly determines required equipment volume and energy consumption — low K_La forces larger, costlier units.

Film Thickness (δ_G, δ_L)

δ_G: 0.05–0.3 mm; δ_L: 0.1–1.0 mm (in turbulent flow, varies with Re and geometry)

Effective hydrodynamic boundary layer thickness in gas and liquid phases where molecular diffusion dominates transport.

⚡ Engineering Impact:

Controls local driving force and rate-limiting phase — thinner films increase flux but require higher turbulence energy input.

Hatta Number (Ha)

Ha < 0.3 (slow reaction), 0.3–3 (moderate), >3 (fast, diffusion-limited)

Dimensionless ratio comparing reaction rate to diffusion rate in the liquid film: Ha = √(k₂·D_L)/k_L.

⚡ Engineering Impact:

Determines whether chemical reaction occurs predominantly in bulk or within the film — dictates reactor configuration (e.g., spray tower vs. membrane contactor).

Interfacial Area (a)

10–100 m²/m³ for plate columns; 200–1200 m²/m³ for high-efficiency packed beds; 500–5000 m²/m³ for microdispersed emulsions

Specific contact area per unit volume of contacting equipment (m²/m³), determined by dispersion mechanism and phase properties.

⚡ Engineering Impact:

Primary lever for intensifying mass transfer — increasing 'a' reduces equipment size but raises pressure drop and entrainment risk.

📐 Key Formulas

Two-Film Overall Coefficient (Gas-Liquid)

1/K_G = 1/k_G + H/k_L

Relates overall gas-phase mass transfer coefficient to individual film resistances and equilibrium constant.

Variables:

Symbol	Name	Unit	Description
K_G	Overall gas-phase mass transfer coefficient	mol/(m^2·s·Pa)	Overall mass transfer coefficient based on gas-phase driving force
k_G	Gas-film mass transfer coefficient	mol/(m^2·s·Pa)	Mass transfer coefficient in the gas film
k_L	Liquid-film mass transfer coefficient	m/s	Mass transfer coefficient in the liquid film
H	Henry's law constant	Pa·m^3/mol	Equilibrium constant relating gas-phase partial pressure to liquid-phase concentration

Typical Ranges:

MEA absorption at 40°C

0.08–0.35 mol/m²·s·Pa

Water–air oxygenation

1.2–5.0 × 10⁻⁴ m/s

⚠️ K_G > 0.1 mol/m²·s·Pa required for commercial CO₂ capture targets

🔧 Open Calculator

Penetration Model Flux

N_A = 2·C*·√(D_A/πτ)

Instantaneous mass flux for a semi-infinite slab exposed for time τ, assuming surface concentration C* and diffusivity D_A.

Variables:

Symbol	Name	Unit	Description
N_A	Instantaneous mass flux	mol/(m²·s)	Molar flux of species A at the surface
C*	Surface concentration	mol/m³	Concentration of species A at the surface
D_A	Diffusivity of species A	m²/s	Diffusion coefficient of species A in the medium
τ	Exposure time	s	Time duration of exposure