Adsorption Column Capacity Estimation

Reference ID: MET-D849 | Process Engineering Reference Sheets Calculation Guide

Introduction & Context

The calculation determines the adsorption capacity of a fixed-bed, down-flow liquid-phase adsorber. It quantifies how much contaminant is retained per unit mass of adsorbent both at full exhaustion (saturation capacity, \(q_{\text{sat}}\)) and at a predefined breakthrough point (working capacity, \(q_{\text{work}}\)). This metric is essential for:

Sizing adsorption columns for water-treatment, air-cleaning, or solvent recovery processes.
Evaluating column utilization and the extent of the mass-transfer zone (MTZ).
Comparing different adsorbents or operating conditions without resorting to detailed kinetic models.

The method is applicable whenever a breakthrough curve (effluent concentration versus time) is available and the system operates isothermally with constant feed concentration and flow rate.

Methodology & Formulas

1. Fundamental mass balance

\[ q = \frac{Q \displaystyle\int_{0}^{T}\!\bigl(C_{0}-C(t)\bigr)\,dt}{M} \]

where
\(Q\) – volumetric flow rate \([\mathrm{m^{3}\,h^{-1}}]\)
\(C_{0}\) – influent concentration \([\mathrm{mg\,L^{-1}}]\)
\(C(t)\) – effluent concentration as a function of time \([\mathrm{mg\,L^{-1}}]\)
\(T\) – integration horizon (either saturation time \(T_{\text{sat}}\) or breakthrough time \(t_{\text{break}}\))
\(M\) – total mass of adsorbent \([\mathrm{kg}]\).

2. Adsorbent mass

\[ M = \rho_{\text{bulk}}\,V_{\text{bed}} \]

\(\rho_{\text{bulk}}\) – dry-packed bulk density \([\mathrm{kg\,m^{-3}}]\)
\(V_{\text{bed}}\) – column bed volume \([\mathrm{m^{3}}]\).

3. Empty-bed contact time (EBCT)

\[ \tau_{\text{EBCT}} = \frac{V_{\text{bed}}}{Q}\quad[\text{h}] \qquad \tau_{\text{EBCT,min}} = 60\,\tau_{\text{EBCT}}\quad[\text{min}] \]

4. Validity checks (empirical ranges)

Parameter	Acceptable Range	Check Expression
EBCT (min)	10–60 min	\(10 \le \tau_{\text{EBCT,min}} \le 60\)
Bulk density \(\rho_{\text{bulk}}\) (kg m\(^{-3}\))	350–500	\(350 \le \rho_{\text{bulk}} \le 500\)
Adsorbent mass \(M\) (kg)	\(M>0\)	\(M>0\)

5. Breakthrough time by linear interpolation

\[ t_{\text{break}} = t_{i} + \frac{C_{\text{b}} - C_{i}}{C_{i+1} - C_{i}}\,(t_{i+1} - t_{i}) \]

where \((t_{i},C_{i})\) and \((t_{i+1},C_{i+1})\) are the two consecutive data points that bracket the breakthrough concentration \(C_{\text{b}}\).

6. Trapezoidal approximation of the integral

\[ \int_{t_{a}}^{t_{b}}\!\bigl(C_{0} - C(t)\bigr)\,dt \;\approx\; \sum_{k}\Bigl[C_{0} - \tfrac{C_{k} + C_{k+1}}{2}\Bigr]\;\Delta t_{k} \]

with \(\Delta t_{k}=t_{k+1}-t_{k}\) and the summation performed over all segments that lie within \([t_{a},t_{b}]\).

7. Unit conversion (mg L\(^{-1}\) → kg m\(^{-3}\))

\[ 1\;\text{mg L}^{-1}=10^{-3}\;\text{kg m}^{-3} \qquad\Longrightarrow\qquad \frac{\displaystyle\int (C_{0}-C)\,dt\;[\text{mg·h L}^{-1}]} {1000} = \int (C_{0}-C)\,dt\;[\text{kg·h m}^{-3}] \]

8. Saturation and working capacities

\[ q_{\text{sat}} = \frac{Q\,\displaystyle\int_{0}^{T_{\text{sat}}}(C_{0}-C)\,dt}{M} \qquad q_{\text{work}} = \frac{Q\,\displaystyle\int_{0}^{t_{\text{break}}}(C_{0}-C)\,dt}{M} \]

9. Column utilization factor

\[ \eta = \frac{q_{\text{work}}}{q_{\text{sat}}}\times 100\;\% \]

Summary Reference Table

Symbol	Description	Unit
\(V_{\text{bed}}\)	Bed volume	m³
\(\rho_{\text{bulk}}\)	Bulk density of adsorbent	kg m\(^{-3}\)
\(M\)	Adsorbent mass	kg
\(Q\)	Volumetric flow rate	m³ h\(^{-1}\)
\(C_{0}\)	Influent concentration	mg L\(^{-1}\)
\(C_{\text{b}}\)	Breakthrough concentration limit	mg L\(^{-1}\)
\(\tau_{\text{EBCT}}\)	Empty-bed contact time	h (or min)
\(q_{\text{sat}}\)	Saturation capacity	kg kg\(^{-1}\) (or g kg\(^{-1}\))
\(q_{\text{work}}\)	Working capacity at breakthrough	kg kg\(^{-1}\)
\(\eta\)	Column utilization (%)	%

By applying the above algebraic framework to any measured breakthrough dataset, engineers can rapidly obtain both the theoretical maximum loading and the practically usable loading of an adsorption column, assess compliance with empirical design windows, and make informed scale-up decisions.

To estimate the breakthrough time from a measured breakthrough curve, apply linear interpolation between the two data points that bracket the allowable effluent concentration limit. For the method presented in this sheet:

Identify the time interval \((t_i, t_{i+1})\) where the effluent concentration passes through \(C_{\text{b}}\).
Calculate \(t_{\text{break}}\) using \( t_{\text{break}} = t_i + \frac{C_{\text{b}}-C_i}{C_{i+1}-C_i}(t_{i+1}-t_i) \).
Alternatively, more detailed approaches such as the bed depth service time (BDST) model or numerical integration of the mass-transfer zone can be used when isotherm parameters are known.

Several operational and physical variables dictate the performance of the column. Key factors include:

Adsorbent particle size and pore structure.
Operating temperature and pressure conditions.
The presence of competing species in the feed stream.
The flow rate, which influences the residence time and the shape of the mass transfer zone.

The selection of an isotherm model, such as Langmuir or Freundlich, defines the equilibrium relationship between the fluid phase and the solid phase. Using an inappropriate model can lead to:

Underestimation of the required bed volume for high-concentration feeds.
Inaccurate predictions of the adsorption capacity at low partial pressures.
Failure to account for surface saturation limits inherent in the Langmuir model.

Worked Example: Adsorption Column Capacity Estimation

A fixed-bed adsorber with granular activated carbon (GAC) is used to remove phenol from water. The column operates isothermally at a constant feed concentration and volumetric flow rate. Breakthrough curve data (effluent concentration vs. time) are collected to determine the saturation and working capacities of the adsorbent.

Knowns:

Bed volume, \( V_{\text{bed}} = 0.1\ \text{m}^{3} \)
Bulk density of GAC, \( \rho_{\text{bulk}} = 400.0\ \text{kg/m}^{3} \)
Mass of adsorbent, \( M = 40.0\ \text{kg} \)
Volumetric flow rate, \( Q = 0.5\ \text{m}^{3}/\text{h} \)
Influent concentration, \( C_{0} = 100.0\ \text{mg/L} \)
Breakthrough concentration limit, \( C_{\text{b}} = 10.0\ \text{mg/L} \)
Empty-bed contact time, \( \tau_{\text{EBCT,min}} = 12.0\ \text{min} \)
Breakthrough curve data points:
- Time (h): 0, 2, 4, 6, 8, 10, 12, 14
- Concentration (mg/L): 0, 0, 5, 20, 50, 80, 95, 100

Step-by-step calculation:

Determine breakthrough time by linear interpolation.
Between \( t = 4\ \text{h} \) and \( t = 6\ \text{h} \), concentration rises from 5 mg/L to 20 mg/L. The breakthrough time at \( C_{\text{b}} = 10\ \text{mg/L} \) is found by interpolation:
\[ t_{\text{break}} = 4 + \frac{10 - 5}{20 - 5} \times (6 - 4) = 4.667\ \text{h} \]
Compute total area above the breakthrough curve (from \( t = 0 \) to \( t = 14\ \text{h} \)).
Using trapezoidal integration of \( (C_{0} - C_{\text{out}}) \) over each time segment:
\[ I_{\text{total}} = 800.0\ \text{mg}\cdot\text{h/L} \]
(This is the cumulative difference in mass concentration units.)
Compute working area above the curve up to breakthrough.
Integrating from \( t = 0 \) to \( t = t_{\text{break}} = 4.667\ \text{h} \):
\[ I_{\text{work}} = 456.667\ \text{mg}\cdot\text{h/L} \]
Convert integrals to consistent SI units.
Since \( 1\ \text{mg/L} = 0.001\ \text{kg/m}^{3} \):
\[ I_{\text{total,SI}} = \frac{800.0}{1000} = 0.8\ \text{kg}\cdot\text{h/m}^{3} \]
\[ I_{\text{work,SI}} = \frac{456.667}{1000} = 0.456667\ \text{kg}\cdot\text{h/m}^{3} \]
Calculate saturation capacity.
The mass of adsorbate retained at full exhaustion is:
\[ q_{\text{sat}} = \frac{Q \times I_{\text{total,SI}}}{M} = \frac{0.5 \times 0.8}{40.0} = 0.01\ \text{kg/kg} = 10.0\ \text{g/kg} \]
Calculate working capacity.
Using the integral up to breakthrough:
\[ q_{\text{work}} = \frac{Q \times I_{\text{work,SI}}}{M} = \frac{0.5 \times 0.456667}{40.0} = 0.005708\ \text{kg/kg} = 5.708\ \text{g/kg} \]
Determine column utilization.
The ratio of working to saturation capacity gives bed utilization:
\[ \eta = \frac{q_{\text{work}}}{q_{\text{sat}}} \times 100 = \frac{5.708}{10.0} \times 100 = 57.08\% \]

Final answer:

Saturation capacity, \( q_{\text{sat}} = 10.0\ \text{g/kg} \)
Working capacity, \( q_{\text{work}} = 5.708\ \text{g/kg} \)
Column utilization, \( \eta = 57.08\% \)

"Un projet n'est jamais trop grand s'il est bien conçu." — André Citroën

"La difficulté attire l'homme de caractère, car c'est en l'étreignant qu'il se réalise." — Charles de Gaulle