Introduction & Context

Depth filtration is a critical unit operation in process engineering, widely utilized in water treatment, chemical processing, and environmental remediation. Unlike surface filtration, where particles are captured on a membrane or screen, depth filtration involves the removal of suspended solids as they pass through a porous medium (e.g., sand, anthracite, or granular activated carbon). The particles are trapped within the interstitial spaces of the bed through mechanisms such as Brownian diffusion, interception, and sedimentation. This calculation is essential for sizing filter beds using a depth filter scale-up calculation to meet specific effluent quality requirements while ensuring that hydraulic pressure drops remain within the operational limits of the system.

Methodology & Formulas

The kinetics of depth filtration are governed by the log-penetration law, which describes the spatial decay of particle concentration as a function of bed depth. The filter coefficient, λ, represents the efficiency of the media in removing particles per unit depth.

The fundamental relationship for particle removal is defined as:

\[ \ln\left(\frac{C}{C_{0}}\right) = -\lambda \cdot Z \]

To determine the filter coefficient λ, we utilize the correlation for single collector efficiency, which accounts for the physical properties of the media and the transport mechanisms of the particles:

\[ \lambda = \frac{3}{2} \cdot \frac{(1 - \varepsilon) \cdot \alpha \cdot \eta}{d_{g}} \]

The hydraulic feasibility of the design is verified using the Ergun equation to calculate the pressure drop per unit length of the bed:

\[ \frac{\Delta P}{L} = \frac{150 \cdot \mu \cdot v \cdot (1 - \varepsilon)^{2}}{d_{g}^{2} \cdot \varepsilon^{3}} + \frac{1.75 \cdot \rho \cdot v^{2} \cdot (1 - \varepsilon)}{d_{g} \cdot \varepsilon^{3}} \]

Where the variables are defined as follows:

C: Outlet concentration
C₀: Inlet concentration
Z: Bed depth (m)
λ: Filter coefficient (m^-1)
ε: Bed porosity (void fraction)
α: Attachment efficiency
η: Single collector efficiency
d_g: Media grain diameter (m)
ΔP/L: Pressure drop per unit length (Pa/m)
μ: Dynamic viscosity (Pa·s)
v: Superficial velocity (m/s)
ρ: Fluid density (kg/m³)

Parameter	Symbol	Validity Constraint
Interception Parameter	N_R	N_R ≤ 0.1
Particle Reynolds Number	Re_p	Re_p ≤ 1.0
Peclet Number	Pe	Pe ≤ 10⁶
Gravity Number	N_G	N_G ≤ 0.01

To determine the optimal flux rate, you must balance throughput against the pressure differential across the filter media. Follow these steps:

Perform a constant flux trial at multiple flow rates to establish the pressure_rise_rate.
Calculate the specific_throughput by dividing the total volume processed by the filter surface area.
Identify the flux point where the pressure_differential exceeds the manufacturer's maximum recommended limit before the target volume is reached.
Select a flux rate that maintains a linear pressure increase to ensure process stability.

Accurate kinetic modeling requires the integration of several key parameters to predict filter life and capacity. The essential variables include:

feed_concentration of the particulate matter.
filter_pore_size_distribution and media_porosity.
fluid_viscosity and temperature, which influence the drag force.
particle_size_distribution relative to the internal void structure of the depth matrix.

The Hermia model is used to characterize the mechanism of pore blockage during the filtration process. For depth filters, the internal_pore_constriction model is typically most applicable. It assumes that:

Particles deposit within the internal structure of the filter media rather than forming a cake on the surface.
The effective pore diameter decreases over time as a function of the volume filtered.
The pressure_drop increases exponentially as the internal void volume is consumed by captured particles.

Worked Example: Depth Filtration Kinetics Calculation

Scenario: A granular deep bed filter is used to remove 10 µm particles from water. The goal is to determine the clean bed depth required to achieve 99 % removal (i.e., outlet concentration ratio C/C₀ = 0.01). The filter operates under favorable attachment conditions.

Known Parameters:

Target removal ratio: C/C₀ = 0.01
Media grain diameter: d_g = 0.5 mm
Bed porosity: ε = 0.4
Fluid viscosity: μ = 1.0 cP
Fluid density: ρ = 998.0 kg/m³
Filtration rate: v = 10.0 m/h
Attachment efficiency: α = 0.1
Single collector efficiency: η = 0.004
Particle diameter: d_p = 10.0 µm

Step-by-Step Calculation:

Convert units to SI base:
d_g = 0.5 mm / 1000 = 0.0005 m
v = 10.0 m/h / 3600 = 0.002777777777777778 m/s
μ = 1.0 cP / 1000 = 0.001 Pa·s
d_p = 10.0 µm / 1×10⁶ = 1×10^-5 m
Compute the filter coefficient λ:
\[ \lambda = \frac{3}{2} \cdot \frac{(1-\varepsilon) \cdot \alpha \cdot \eta}{d_g} \] Substituting the known values: \[ \lambda = \frac{1.5 \cdot (1-0.4) \cdot 0.1 \cdot 0.004}{0.0005} = 0.72 m^-1 \]
Compute the required bed depth Z:
From the depth removal law: ln(C/C₀) = -λZ
\[ Z = -\frac{\ln(C/C_0)}{\lambda} \] With C/C₀ = 0.01: \[ Z = -\frac{\ln(0.01)}{0.72} = \frac{4.605170185988091}{0.72} = 6.396 m \]
Compute pressure drop gradient ΔP/L (Ergun equation):
ΔP/L = 150 μv ^(1-ε)² / _{d_g²ε³} + 1.75 ρv² ^(1-ε) / _{d_gε³}
First term: \[ \text{term1} = \frac{150 \cdot 0.001 \cdot 0.002777777777777778 \cdot (1-0.4)^2}{(0.0005)^2 \cdot (0.4)^3} = 9374.999999999998 Pa/m \]
Second term: \[ \text{term2} = \frac{1.75 \cdot 998.0 \cdot (0.002777777777777778)^2 \cdot (1-0.4)}{0.0005 \cdot (0.4)^3} = 252.67650462962962 Pa/m \]
Total: \[ \Delta P/L = 9374.999999999998 + 252.67650462962962 = 9627.677 Pa/m \]
Validity checks (empirical bounds):
Interception parameter: N_R = d_p / d_g = 1×10^-5 / 0.0005 = 0.02
N_R = 0.02 is less than the upper bound of 0.1 → valid.
Particle Reynolds number: Re_p = ρvd_p / μ = 998.0 · 0.002777777777777778 · 1×10^-5 / 0.001 = 0.028
Re_p = 0.028 is less than 1.0 → valid.

Final Answer:

A clean bed depth of 6.396 m is required to achieve 99 % removal of 10 µm particles under the given conditions. The corresponding pressure drop gradient is 9627.677 Pa/m.

"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

Depth Filtration Kinetics Calculation

Introduction & Context

Methodology & Formulas

Worked Example: Depth Filtration Kinetics Calculation

🚀 Skip the Manual Math!