Absorption Tower Sizing Using Kremser Method - Chemical Engineering Knowledge Platform

Typical Scale

Packed heights: 3–15 m; diameters: 1–6 m; capacities: 10–500 kmol/h gas

Key Standards

GPSA Engineering Data Book (14th ed.), Perry’s Chemical Engineers’ Handbook (9th ed.)

Industry Applications

Flue gas CO₂ capture, natural gas sweetening (H₂S/CO₂ removal), VOC abatement, HF alkylation off-gas scrubbing

⚠️ Why It Matters

1

Incorrect tower height estimation

↓

2

Insufficient solute removal

↓

3

Non-compliant emissions or product purity

↓

4

Retrofitting or re-rating required later

↓

5

Increased capital and operating cost

↓

6

Process safety or environmental risk

📘 Definition

The Kremser method is an analytical solution to the material balance equations for countercurrent multistage absorption (or stripping), assuming constant molar flow rates, linear equilibrium (y = mx), and no chemical reaction. It expresses fractional solute removal as a function of stage efficiency, absorption factor (A = L/mV), and number of theoretical stages (N). Though approximate, it remains widely used in preliminary design due to its algebraic simplicity and physical transparency.

🎨 Concept Diagram

AI-generated illustration for visual understanding

💡 Engineering Insight

The Kremser method shines not as a final design tool—but as a diagnostic lens: if your calculated Nₜ exceeds 15–20, the process is likely thermodynamically constrained—not equipment-limited. That’s your signal to revisit solvent selection, operating pressure, or consider reactive absorption instead of chasing taller towers.

📖 Detailed Explanation

The Kremser method begins with the intuitive idea that gas absorption happens in discrete, idealized stages where gas and liquid reach equilibrium before moving to the next stage. By writing overall and stage-wise material balances and assuming linear equilibrium (y = mx), the method derives a closed-form expression linking fractional removal (φ), absorption factor (A), and number of theoretical stages (N): φ = (1 − A⁻ᴺ)/(1 − A⁻ᴺ⁺¹). This avoids iterative McCabe–Thiele construction while preserving physical insight.

Its power lies in revealing the nonlinear penalty of high recovery: going from φ = 0.95 to φ = 0.99 requires ~2.5× more stages when A = 1.2—but only ~1.6× when A = 1.8. This explains why industrial absorbers targeting >99% removal almost always operate at A > 1.5, even though higher A increases solvent pumping cost. The trade-off is explicit and quantifiable.

Advanced use extends beyond the textbook form: for non-constant flows (e.g., significant solvent absorption of carrier gas), the method can be adapted using pseudo-stages or numerical integration of the Kremser differential form. When chemical reaction enhances absorption (e.g., amine-CO₂), the apparent ‘m’ becomes composition- and pH-dependent—requiring effective mₑff estimation from reaction kinetics or pilot data. In such cases, Kremser serves best as a bounding check—not a standalone solution.

🔄 Engineering Workflow

Step 1

Step 1: Define inlet gas composition, flow rate, target outlet concentration, and operating T/P

→

Step 2

Step 2: Select solvent and determine equilibrium relationship (m) via Henry’s law, NRTL, or experimental data

→

Step 3

Step 3: Estimate minimum solvent flow (Lₘᵢₙ) and choose practical L/V ratio to fix absorption factor A

→

Step 4

Step 4: Apply Kremser equation to calculate required theoretical stages (Nₜ) for specified φ

→

Step 5

Step 5: Convert Nₜ to actual height using vendor-supplied HETP or correlation-based Eₘ (e.g., Onda correlation)

→

Step 6

Step 6: Perform hydraulic check (capacity, pressure drop, flooding margin) and adjust diameter/height iteratively

→

Step 7

Step 7: Validate against pilot data or process simulator (Aspen Plus, CHEMCAD) before FEED approval

📋 Decision Guide

Rock/Field Condition	Recommended Design Action
High-purity requirement (φ ≥ 0.995) with low-m system (m < 0.01, e.g., CO₂ in water)	Use A ≈ 1.4–1.8; verify with rigorous simulation; consider chemically enhanced solvent or staged towers
High-m system (m > 2, e.g., SO₂ in NaOH) with moderate φ (0.90–0.95)	Operate at A = 0.8–1.0 to minimize solvent circulation; single packed bed < 6 m tall usually sufficient
Viscous solvent (μ > 20 cP) or foaming tendency observed in lab tests	Reduce design A by 15–20%; increase HETP estimate by 30%; specify low-pressure-drop structured packing (e.g., Mellapak 250.Y)

📊 Key Properties & Parameters

Absorption Factor (A)

0.8–2.5 (dimensionless)

Ratio of liquid-to-gas molar flow rate divided by equilibrium slope (A = L/(m·V)) — determines driving force and feasibility of separation.

⚡ Engineering Impact:

A < 0.8 leads to impractical stage requirements; A > 2.5 wastes solvent flow and increases pumping energy.

Equilibrium Slope (m)

0.001–10 (dimensionless, highly system-dependent)

Henry’s law constant expressed as m = y*/x (ratio of solute mole fraction in gas phase to liquid phase at equilibrium).

⚡ Engineering Impact:

Low m (e.g., CO₂ in water) implies high resistance to absorption and demands taller towers; high m (e.g., NH₃ in water) enables efficient removal with fewer stages.

Fractional Removal (φ)

0.90–0.999 (90%–99.9%)

Proportion of solute removed from inlet gas stream: φ = (y₁ − y₂)/y₁.

⚡ Engineering Impact:

Each additional decimal (e.g., 99% → 99.9%) typically doubles required theoretical stages—and thus tower height—under fixed A.

Stage Efficiency (Eₘ)

0.3–0.8 m per theoretical stage (for structured packings); 0.5–1.5 m for random packings

Ratio of actual height to theoretical stage height (HETP = Hₜ/Nₜ), representing mass-transfer effectiveness per unit packing height.

⚡ Engineering Impact:

Underestimating Eₘ overpredicts performance and risks underdesign; typical values depend strongly on fluid properties, loading, and packing type.

📐 Key Formulas

Kremser Equation (Absorption)

φ = (1 − A⁻ᴺ) / (1 − A⁻ᴺ⁺¹)

Calculates fractional solute removal given absorption factor A and theoretical stages N

Variables:

Symbol	Name	Description
φ	Fractional solute removal	Dimensionless fraction of solute removed in absorption process
A	Absorption factor	Ratio of liquid to gas flow rates times equilibrium constant, dimensionless
N	Number of theoretical stages	Integer count of ideal equilibrium stages in absorber

Typical Ranges:

Natural gas sweetening (H₂S)

A = 1.2–1.8, N = 4–10, φ = 0.95–0.999

Post-combustion CO₂ capture

A = 1.4–2.0, N = 8–14, φ = 0.90–0.99

⚠️ A must be >0.5 for convergence; avoid A < 0.8 unless φ < 0.9

🔧 Open Calculator

Minimum Solvent Flow

Lₘᵢₙ/V = m·y₁/x₁*

Determines lowest feasible liquid flow to achieve equilibrium at top stage

Variables:

Symbol	Name	Unit	Description
Lₘᵢₙ	Minimum Solvent Flow Rate	mol/s or kg/s	Lowest feasible liquid (solvent) molar or mass flow rate to achieve equilibrium at the top stage
V	Inert Gas Flow Rate	mol/s or kg/s	Molar or mass flow rate of inert carrier gas (e.g., air, nitrogen)
m	Distribution Coefficient	dimensionless	Equilibrium slope (y = mx) in solute distribution between gas and liquid phases
y₁	Solute Mole (or Mass) Fraction in Inlet Gas	dimensionless	Concentration of solute in entering gas stream
x₁*	Equilibrium Solute Mole (or Mass) Fraction in Exit Liquid	dimensionless	Solute concentration in liquid phase in equilibrium with inlet gas composition y₁