First-Principles Mathematical Derivation of Countercurrent Extraction
Liquid-liquid extraction columns operate under countercurrent flow where solute transfers across liquid-liquid phase boundaries driven by thermodynamic chemical activity differences.
1. Overall Solute Mass Balance
Solute entering in the feed ($F \cdot x_{in}$) and solvent ($S \cdot y_{in}$) equals solute exiting in raffinate ($R \cdot x_{out}$) and extract ($E_{str} \cdot y_{out}$):
F \cdot x_{in} + S \cdot y_{in} = F \cdot x_{out} + S \cdot y_{out}\implies y_{out} = y_{in} + \frac{F}{S} (x_{in} - x_{out})
2. The Kremser Equation for Theoretical Stages
The extraction factor $E = \frac{K_D \cdot S}{F}$ represents the ratio of the equilibrium line slope to the operating line slope. The exact analytical Kremser formula gives:
N_{theor} = \frac{\ln \left[ \left( \frac{x_{in} - y_{in} / K_D}{x_{out} - y_{in} / K_D} \right) \left( 1 - \frac{1}{E} \right) + \frac{1}{E} \right]}{\ln(E)}
3. Column Hydraulic Sizing & Packed Height
Total packed height accounts for stage efficiency or HETS:
H_{pack} = N_{theor} \cdot HETS\quad\text{and}\quad N_{actual} = \frac{N_{theor}}{E_o / 100}
5 Fatal Traps & Engineering Pitfalls in Extraction Column Design
1. Sub-Critical Extraction Factor ($E < 1.0$)
If $E = K_D \cdot S / F < 1.0$, the solvent capacity is thermodynamically insufficient to extract the solute. An infinite number of stages ($N \to \infty$) will still fail to reach the target raffinate purity, leading to permanent off-spec waste discharge.
2. Emulsification & Phase Inversion Flooding
Excessive rotor agitation speed in stirred columns (Karr, Scheibel, RDC) shatters droplets below $0.5\,\text{mm}$. Micro-droplets lack buoyant velocity to overcome continuous phase drag, causing stable emulsion inversion and spewing solvent directly into the raffinate discharge.
3. Assuming Constant $K_D$ Across Broad Concentration Ranges
In concentrated systems, high solute levels increase mutual solvent-water solubility (solutropy). The distribution coefficient $K_D$ drops significantly near the column inlet. Designing columns using dilute Henry's law $K_D$ underestimates required stage count by 40%.
4. Axial Backmixing in Wide Columns ($D_c > 1.5\,\text{m}$)
Large cross-sectional areas promote Taylor circulation cells where continuous phase fluid recirculates upstream. Axial dispersion degrades interstage concentration driving forces, requiring 25% to 50% more height than pilot-scale test columns.
5. Internal Surface Wettability Inversion
Internal packings or sieve trays must be preferentially wetted by the continuous phase. If the dispersed phase wets the packing material, droplets coalesce into thick liquid films that channel along the steel, collapsing specific interfacial area by 85%.
Frequently Asked Questions: Liquid-Liquid Extraction Columns
What is Liquid-Liquid Extraction (LLE) and when is it preferred over distillation? +
Liquid-Liquid Extraction (LLE), or solvent extraction, is a separation unit operation that separates components of a liquid solution by contacting it with an immiscible or partially miscible solvent that selectively dissolves one or more solutes. LLE is preferred when distillation is impossible or uneconomical, such as: separating close-boiling mixtures, recovering high-boiling solutes from dilute aqueous streams (e.g. phenol recovery), separating azeotropes, and isolating thermally labile biopharmaceuticals, vitamins, or antibiotics that degrade at distillation temperatures.
What is the Extraction Factor (E) and why must E > 1 for high recovery? +
The extraction factor ($E$) represents the ratio of the solute's equilibrium capacity in the solvent stream to that in the feed stream: $E = \frac{K_D \cdot S}{F}$, where $K_D$ is the distribution coefficient ($y^* / x$), $S$ is solvent mass flow rate, and $F$ is feed mass flow rate. If $E < 1$, the operating line slope is less than the equilibrium line slope, making it mathematically impossible to achieve high solute recovery regardless of how many theoretical stages are added ($N \to \infty$). For economical industrial column design, engineers target $E = 1.3\text{--}2.0$.
How does the Kremser equation determine the number of theoretical stages? +
For dilute, immiscible systems with linear equilibrium ($y^* = K_D \cdot x$), the Kremser shortcut equation analytically calculates the required number of theoretical stages $N_{theor}$: $N_{theor} = \frac{\ln\left[ \left(\frac{x_{in} - y_{in}/K_D}{x_{out} - y_{in}/K_D}\right) \left(1 - \frac{1}{E}\right) + \frac{1}{E} \right]}{\ln(E)}$. If $E = 1$, the equation simplifies to $N_{theor} = \frac{x_{in} - x_{out}}{x_{out} - y_{in}/K_D}$.
What is HETS and how does it determine column packed height? +
HETS stands for Height Equivalent to a Theoretical Stage. In packed extraction columns, mass transfer occurs continuously rather than in discrete equilibrium stages. The total required height of packing is $H_{pack} = N_{theor} \cdot HETS$. For standard structured packings or random dump packings (such as ceramic saddles or metal Pall rings), HETS typically ranges from $0.6$ to $1.2\,\text{m}$ ($2\text{--}4\,\text{ft}$), depending on interfacial tension, phase viscosity, and droplet Sauter mean diameter.
What causes flooding and emulsification in extraction columns? +
Column flooding occurs when the counter-current flow of the continuous and dispersed phases exceeds hydraulic capacity. Droplets coalesce into large stagnant liquid slugs that exit through the wrong outlet, causing carryover. Emulsification is triggered by intense mechanical shear (excessive rotor RPM) or trace surface-active contaminants that reduce interfacial tension below $2\,\text{mN/m}$, creating micro-droplets that refuse to settle.