Absorption Column Design Tool Back
Chemical Engineering · Separation

Absorption Column Design Tool — NTU / HTU

Visualize gas absorption operating line and equilibrium curve on a y-x diagram. Calculate NTU graphical integration, packed column height, and absorption factor in real time.

Flow & Composition Parameters
Gas flow rate G (kmol/h)
Liquid flow rate L (kmol/h)
Inlet gas y₁ (mol fr.)
Outlet gas target y₂
Henry's constant m
Column diameter D (m)
⚠ Absorption factor A ≤ 1: absorption infeasible. Increase liquid flow rate.
Results
NTU
HTU
Column height Z
Absorption factor A
Yx

Blue: operating line  |  Green: equilibrium curve y*=mx  |  Orange dashed: NTU integration path

Colschematic
Theory & Key Formulas
Operating line: $L(x-x_0) = G(y-y_2)$
Equilibrium: $y^* = mx$
NTU: $N_{OG}= \int_{y_2}^{y_1}\dfrac{dy}{y-y^*}$
Absorption factor: $A = \dfrac{L}{mG}$

What is the NTU/HTU Method?

🙋
What exactly is the NTU/HTU method for designing an absorption column? I see the terms "Number of Transfer Units" and "Height of a Transfer Unit" but I'm not sure what they physically represent.
🎓
Basically, it's a clever way to separate the design problem into two parts. Think of NTU as the "difficulty" of the separation—how many theoretical steps are needed to scrub the gas clean. HTU is the "effectiveness" of the packing material—how tall each of those steps is in a real column. The total height is just Z = NTU × HTU. In this simulator, you can see how changing the gas flow rate (G) or liquid flow rate (L) directly changes that difficulty.
🙋
Wait, really? So the NTU is calculated from an integral? That seems complicated. How do I see it in practice here?
🎓
That's the beauty of this tool! The integral $N_{OG}= \int \frac{dy}{y-y^*}$ is solved graphically. You see the blue operating line and the orange equilibrium line on the plot? The vertical distance between them at any point is $(y-y^ )$. The NTU is the area under the curve of $1/(y-y^ )$. Try it: slide the "Outlet gas target y₂" to make it lower. See how the lines get closer at the bottom, making that distance tiny? That makes the $1/(y-y^*)$ huge, and the NTU integral—the required number of steps—skyrockets!
🙋
Okay, that makes sense visually. What about this "Absorption Factor" A = L/(mG)? What's its role?
🎓
Great question. The absorption factor is the key driver of your column's economics. It's the ratio of liquid "scrubbing power" (L) to gas "load" (mG). In practice, if A is much greater than 1, you have plenty of liquid, the operating line is steep, and separation is easy (low NTU). If A is less than 1, you're trying to do a lot with a little liquid, the lines pinch, and NTU goes to infinity! Play with the sliders for L and G. Watch how the slope of the blue operating line changes and see the direct impact on the calculated NTU and final column height.

Physical Model & Key Equations

The design is based on a mass balance across the column, which gives the Operating Line. It relates the compositions of the gas (y) and liquid (x) at any height in the column.

$$L(x - x_0) = G(y - y_2)$$

Where L is liquid molar flow rate (kmol/h), G is gas molar flow rate (kmol/h), x is liquid-phase mole fraction of solute, y is gas-phase mole fraction of solute. Subscripts 0 and 2 refer to the bottom (inlet) and top (outlet) of the column, respectively. This line is plotted in blue on the simulator.

The driving force for mass transfer is the difference between the actual gas composition (y) and the composition that would be in equilibrium with the liquid (y*). For dilute systems, equilibrium is given by Henry's Law.

$$y^* = m x$$

Where m is Henry's constant (dimensionless). This is the orange Equilibrium Line on the plot. The Number of Transfer Units (NTU) is then found by integrating the reciprocal of the driving force over the column's concentration range.

$$N_{OG}= \int_{y_2}^{y_1}\dfrac{dy}{y - y^*}$$

The physical height of the column is this "number of steps" multiplied by the "height per step": Z = NOG × HOG, where HOG is the Height of a Transfer Unit, a property of the packing and system.

Frequently Asked Questions

The intersection of the operating line and the equilibrium line indicates that the gas and liquid compositions are in equilibrium at that point. If this point does not match the conditions at the top or bottom of the column, it means that composition cannot be reached in actual operation, serving as a reference for checking design constraints or limitations.
On the y-x diagram, measure the vertical distance (driving force) between the operating line and the equilibrium line. Divide the y-axis into intervals, calculate 1/(y - y*) for each interval, and integrate. The tool performs this integration automatically and displays the NTU value in real time.
The packed column height is obtained by multiplying the NTU (Number of Transfer Units) by the HTU (Height of a Transfer Unit). Since HTU depends on the type of packing material, gas and liquid flow rates, and physical properties (such as diffusion coefficient and viscosity), input of experimental or literature values is necessary.
When the absorption factor A = L/(mG) is less than 1, the slope of the operating line becomes gentler than that of the equilibrium line, reducing the driving force within the column. In this case, achieving the required absorption rate may require an extremely tall column, so consider increasing the liquid flow rate or selecting an absorbent with a smaller m value.

Real-World Applications

Natural Gas Sweetening: Raw natural gas straight from the well contains corrosive hydrogen sulfide (H₂S) and carbon dioxide (CO₂). Absorption columns using amine solvents (like MEA) are designed with the NTU/HTU method to calculate the required packed bed height to "sweeten" the gas to pipeline specifications before transport.

Ammonia Production Scrubbing: In the Haber process, the product gas stream contains unreacted ammonia. Absorption columns using water as the solvent are designed to recover this valuable ammonia. Engineers use this method to balance the column height (capital cost) against water circulation rate (operating cost).

Removing CO₂ from Flue Gas: In carbon capture and storage (CCS) systems for power plants, absorption columns are the first major unit operation. The NTU/HTU method is critical for designing an efficient scrubber that can handle huge volumes of flue gas to meet specific CO₂ removal targets, directly impacting the feasibility and cost of the whole CCS project.

Odor Control in Wastewater Treatment: Air stripped from wastewater treatment processes often contains foul-smelling compounds like hydrogen sulfide. Small-scale packed tower scrubbers, designed using these principles, use chemicals like sodium hydroxide to absorb and neutralize the odors before the air is released to the environment.

Common Misunderstandings and Points to Note

When you start using this tool, there are a few common pitfalls to watch out for. First, understand that the equilibrium constant m is not a fixed value. The tool treats it as constant for simplicity, but in actual absorption processes, m changes with temperature and concentration. For example, in ammonia absorption into water, a 10°C increase in liquid temperature can cause m to increase by a factor of 1.5. This can lead to a higher actual NTU than the NTU assumed during design, potentially resulting in insufficient separation performance. In practice, remember that detailed calculations, such as using different m values at the top and bottom of the column to account for the heat of dissolution generated inside, are often necessary.

Next, a common misunderstanding is that "HTU depends solely on flow rate." Looking at the tool's formula $H_{OG}= G / (K_y a A)$, it seems HTU is proportional to gas flow rate G, right? However, the mass transfer capacity coefficient $K_y a$ itself also changes with flow rate (more precisely, with gas-liquid shear forces). For typical Raschig ring packing, $K_y a$ tends to increase as G increases. Therefore, HTU does not simply scale proportionally; in reality, it often has a minimum value (an optimum point) at a certain flow rate . When using catalog values for $K_y a$, always check the experimental conditions (at which flow rate it was measured).

Finally, a practical pitfall: "You can't just calculate the column diameter D." While the tool calculates column height Z when you change the diameter, in the field, pressure drop, liquid holdup, and the quality of liquid distribution are critical. For instance, arbitrarily increasing liquid flow in a 1m diameter column can cause "liquid channeling," preventing you from achieving the calculated performance. Generally, you need to separately satisfy hydrodynamic constraints, such as ensuring the liquid spray density does not exceed $5 \,\mathrm{m^3/(m^2 \cdot h)}$.

How to Use

  1. Enter gas inlet mole fraction (y1, e.g., 0.08 for 8% CO2) and outlet mole fraction (y2, e.g., 0.01)
  2. Input gas mass flow rate (e.g., 500 kg/h) and liquid flow rate (e.g., 1200 kg/h) in consistent units
  3. The simulator plots the operating line and equilibrium curve; graphical integration over the y-difference zone yields NTU
  4. Multiply NTU by HTU (height of transfer unit, derived from mass transfer coefficients) to obtain packed column height Z

Worked Example

Design an SO2 absorption column using 10% NaOH solution. Gas inlet y1 = 0.06 mol SO2/mol, outlet y2 = 0.002, gas flow = 800 kg/h, liquid flow = 2000 kg/h. Equilibrium relation: y* = 0.15x. Graphical integration yields NTU = 3.2. With HTU = 0.45 m (based on 25 mm Raschig rings, gas-phase mass transfer coefficient kg = 0.18 kmol/m²·s·Pa, molar density = 40 kmol/m³), column height Z = 3.2 × 0.45 = 1.44 m. Absorption factor A = (L/G)m_eq = (2000/800) × 0.15 = 0.375, confirming absorption feasibility (A < 1).

Practical Notes

  1. For dilute gases (y < 0.15), assume constant molar overflow; check Henry's law validity against actual equilibrium data from literature (e.g., NIST or DDBST)
  2. HTU increases with gas velocity; use flooding correlation (e.g., Sherwood-Holloway for random packing) to verify 60–70% of flood velocity
  3. Absorption factor A near 0.5 requires tall columns; consider staged contactors (bubble caps, sieve trays) or reduce inlet loading for cost-effective design
  4. Temperature sensitivity: equilibrium slope m changes ~3–5% per 10 K; run sensitivity analysis for ±5 K inlet temperature variance