Mass Transfer Coefficients: kₗ, k₉, Kₗ, K₉ and Their Correlations

⚠️ Why It Matters

1

Inaccurate kₗ estimation in CO₂ absorption

↓

2

Underdesigned column height and packing volume

↓

3

Failure to meet 90% removal target

↓

4

Regulatory noncompliance (e.g., EPA 40 CFR Part 63)

↓

5

Plant shutdown or costly retrofit

📘 Definition

Mass transfer coefficients (kₗ, k₉, Kₗ, K₉) quantify the rate of interphase transport per unit driving force (concentration difference) under convective-diffusive conditions. Local coefficients (kₗ, k₉) describe transfer at the interface using film theory, while overall coefficients (Kₗ, K₉) account for resistances across both phases and are referenced to bulk-phase driving forces. They are defined via flux equations: N_A = kₗ(c* − cₗ) = k₉(p₉ − p*) = Kₗ(cₗ,b − cₗ*) = K₉(p₉,b − p₉*).

🎨 Concept Diagram

_A = k₉(p₉,b − p*)N_A = kₗ(cₗ* − cₗ,b)

AI-generated illustration for visual understanding

💡 Engineering Insight

Never assume kₗ and k₉ scale linearly with flow — at low Re, k ∝ Re^0.33 (laminar film), but above Re > 1000 in packed beds, k ∝ Re^0.67 (turbulent dispersion dominates). Always check if your correlation’s Re range matches your design point; extrapolation beyond ±20% introduces >40% error in HTU.

📖 Detailed Explanation

Mass transfer coefficients originate from two limiting models: the stagnant film theory (Whitman, 1923), where transport occurs across a hypothetical immobile layer, and the penetration theory (Higbie, 1935), which accounts for transient diffusion into fluid elements exposed intermittently at the interface. These yield k ∝ D^(2/3)·u^(1/3) and k ∝ D^(1/2)·f^(1/2), respectively — foundational for all empirical correlations.

Modern correlations (e.g., Onda, Bravo–Rocha–Fair) embed these scaling laws within dimensionless groups: Sh = f(Re, Sc), where Sherwood number (Sh = k·L/D) represents dimensionless k, Reynolds (Re = u·L/ν) captures flow regime, and Schmidt (Sc = ν/D) encodes fluid properties. Critical to application is selecting the correct characteristic length L (e.g., column diameter for trays, equivalent sphere diameter for packings, hydraulic diameter for channels).

Advanced treatment requires recognizing that kₗ and k₉ are not intrinsic fluid properties but *system-specific* responses to hydrodynamics, interfacial area (a), and local turbulence. In reactive systems, the enhancement factor E = kₗ,eff / kₗ,phys depends on the Hatta number (Ha = δ·√(k₂·C_B)/Dₗ), linking kinetics to mass transfer. When Ha > 3, reaction occurs entirely within the liquid film — enabling dramatic kₗ amplification without increasing energy input.

🔄 Engineering Workflow

Step 1

Step 1: Identify controlling phase using Henry’s law constant and diffusivities

→

Step 2

Step 2: Select appropriate correlation (e.g., Onda for packed beds, Chilton–Colburn for trays)

→

Step 3

Step 3: Calculate Reynolds, Schmidt, and Sherwood numbers from operating conditions

→

Step 4

Step 4: Compute kₗ and k₉ using correlation; verify resistance distribution (1/k₉ vs. H/kₗ)

→

Step 5

Step 5: Derive Kₗ or K₉ using overall resistance summation (1/Kₗ = 1/kₗ + H/k₉)

→

Step 6

Step 6: Size equipment via HTU–NTU or HETP methods using calculated K-values

→

Step 7

Step 7: Validate with pilot-scale data or tracer tests (e.g., dynamic response of pH or conductivity)

📋 Decision Guide

Rock/Field Condition	Recommended Design Action
Gas-film controlled system (H > 10⁶ Pa·m³/mol, e.g., O₂ in water)	Optimize gas turbulence (e.g., increase superficial velocity, use high-efficiency distributors); liquid flow rate has minimal effect on Kₗ.
Liquid-film controlled system (H < 10⁴ Pa·m³/mol, e.g., SO₂ in water)	Enhance liquid-side mixing (e.g., high-velocity spray, rotating packed beds, or chemical enhancement with NaOH).
Chemically reactive absorption (e.g., CO₂ + MEA)	Use enhancement factor (E) correlations (e.g., van Krevelen & Hoftijzer) to scale kₗ; design for 3–5× higher effective Kₗ than physical absorption.

📊 Key Properties & Parameters

kₗ (liquid-phase coefficient)

1×10⁻⁵ – 5×10⁻³ m/s

Local mass transfer coefficient for the liquid phase, defined as flux divided by the liquid-side concentration driving force (mol/m²·s per mol/m³).

⚡ Engineering Impact:

Directly governs required contact time and column diameter in packed absorption towers.

k₉ (gas-phase coefficient)

1×10⁻⁴ – 2×10⁻² m/s

Local mass transfer coefficient for the gas phase, defined as flux divided by the gas-side partial pressure driving force (mol/m²·s per Pa).

⚡ Engineering Impact:

Controls minimum gas velocity to avoid flooding and determines pressure drop across trays or packings.

Kₗ (overall liquid-phase coefficient)

5×10⁻⁶ – 1×10⁻³ m/s

Overall coefficient referenced to the liquid bulk concentration driving force, incorporating both gas- and liquid-film resistances.

⚡ Engineering Impact:

Used with HTU–NTU method to size absorbers; low Kₗ demands taller columns or higher solvent flow rates.

H (Henry’s law constant)

10³ – 10⁸ Pa·m³/mol (for dilute aqueous systems)

Equilibrium constant relating solute partial pressure in gas to its mole fraction or concentration in liquid (p = H·x or p = H·c).

⚡ Engineering Impact:

Dominates resistance distribution: high H implies gas-film control (k₉ governs); low H implies liquid-film control (kₗ governs).

📐 Key Formulas

Onda Correlation (packed beds)

Sh = 0.023·Re^0.83·Sc^0.45·(μ/μ_w)^0.14

Predicts kₗ and k₉ for random packings (e.g., Raschig rings, Pall rings) under turbulent flow

Variables:

Symbol	Name	Unit	Description
Sh	Sherwood number	dimensionless	Dimensionless mass transfer coefficient
Re	Reynolds number	dimensionless	Ratio of inertial to viscous forces
Sc	Schmidt number	dimensionless	Ratio of momentum diffusivity to mass diffusivity
μ	dynamic viscosity of fluid	Pa·s	Viscosity of the bulk fluid
μ_w	dynamic viscosity at wall temperature	Pa·s	Viscosity of fluid at the wall (or surface) temperature

Typical Ranges:

Ceramic Raschig rings, Re = 40–2000

Sh = 10–120

Metal Pall rings, Re = 500–5000

Sh = 50–350

⚠️ Re < 40 invalidates correlation (laminar film regime); Sc > 10⁴ requires correction for high-viscosity liquids

🔧 Open Calculator

Overall coefficient (liquid-referenced)

1/Kₗ = 1/kₗ + H/k₉

Converts local coefficients into overall Kₗ using Henry’s law constant H (Pa·m³/mol)

Variables:

Symbol	Name	Unit	Description
Kₗ	Overall liquid-phase mass transfer coefficient	m/s	Overall coefficient referenced to the liquid phase
kₗ	Local liquid-phase mass transfer coefficient	m/s	Local coefficient at the liquid side
k₉	Local gas-phase mass transfer coefficient	m/s	Local coefficient at the gas side
H	Henry's law constant	Pa·m³/mol	Ratio of partial pressure to concentration, relating gas and liquid phases