Back to Publications
Process Engineering7 min read

Liquid-Liquid Extraction (LLE) & Mixer-Settler Sizing in API Downstream Isolation

Kiran SeepanaOctober 1, 20266 Views
Executive Summary & Scope

Master industrial Liquid-Liquid Extraction (LLE) for pharmaceutical API downstream workups. Learn partition coefficients, Kremser stage calculations, mixer power numbers, and settler hydraulic design.

Peer-Reviewed & PE Verified

ASME VIII • NFPA 68/69 • TEMA • ISO 9001 Alignment

This technical publication and associated design calculations have been reviewed for engineering consistency, unit integrity, and alignment with standard process design practices (Process Engineering).

# Liquid-Liquid Extraction (LLE) & Mixer-Settler Sizing in API Downstream Isolation

# Advanced Partition Thermodynamics, Hydraulic Dispersion & Counter-Current Stage Sizing for Pharmaceutical Workups

In active pharmaceutical ingredient (API) synthesis, natural product extraction, and fermentation broth recovery, Liquid-Liquid Extraction (LLE) is the quintessential unit operation for separating heat-labile molecules from aqueous reaction matrices, unreacted reagents, and polar inorganic salts.

Unlike thermal distillation, which risks thermal degradation, racemization, or product charring, liquid extraction operates at mild or ambient temperatures by exploiting differences in chemical potential and solubility between two immiscible liquid phases.


Counter-Current Mixer-Settler Battery Flowsheet
Counter-Current Mixer-Settler Battery Flowsheet


# 1. Thermodynamic Fundamentals: Distribution Coefficient & Selectivity

The core thermodynamic driving force of liquid-liquid extraction is the Distribution Coefficient (Partition Ratio), KDK_D, defined as the equilibrium ratio of solute mass fraction in the extract phase (yi∗y_i^*) to the raffinate phase (xix_i):

KD=yi∗xi=γi,R∞γi,E∞⋅MRMEK_D = \frac{y_i^*}{x_i} = \frac{\gamma_{i,R}^\infty}{\gamma_{i,E}^\infty} \cdot \frac{M_{R}}{M_{E}}

Where:

  • γi,R∞\gamma_{i,R}^\infty: Activity coefficient of solute in the raffinate phase at infinite dilution.
  • γi,E∞\gamma_{i,E}^\infty: Activity coefficient of solute in the extract phase at infinite dilution.
  • MR,MEM_R, M_E: Average molecular weight of raffinate and extract phases.

# 1.1. Separation Factor (Selectivity β\beta)

When separating a target API (AA) from an undesirable process impurity (BB), the solvent selectivity βA/B\beta_{A/B} dictates thermodynamic feasibility:

βA/B=KD,AKD,B=yA/xAyB/xB\beta_{A/B} = \frac{K_{D,A}}{K_{D,B}} = \frac{y_A / x_A}{y_B / x_B}
ℹ️ Note
If βA/B=1.0\beta_{A/B} = 1.0, separation by extraction is impossible regardless of stage count. For economic pharmaceutical manufacturing, βA/B≥3.5\beta_{A/B} \ge 3.5 is recommended.

# 1.2. The Extraction Factor (EE)

To evaluate equipment sizing feasibility, process engineers evaluate the dimensionless Extraction Factor (EE):

E=KD⋅(SF)E = K_D \cdot \left(\frac{S}{F}\right)

Where:

  • SS: Solvent mass flow rate (kg/h).
  • FF: Feed mass flow rate (kg/h).
                         THE EXTRACTION FACTOR REGIME SPECTRUM
   E < 1.0                     1.3 <= E <= 2.5                     E > 3.5
  ◄───────────────────────────┼─────────────────────────────────┼──────────────────────────►
   Thermodynamically Limited    Optimal Industrial Sizing Zone    Solvent Inefficiency Zone
   Incomplete recovery even     Balanced stage count & solvent    High recovery per stage,
   with infinite stages.        recovery utility costs.           but extreme downstream evap.

# 2. Multi-Stage Counter-Current Cascade Sizing: The Kremser Model

In multi-stage counter-current extraction batteries, fresh solvent enters at the opposite end of the feed stream, maintaining a nearly uniform chemical potential gradient across every stage.

flowchart LR
    F["Aqueous Feed (F, x_f)"] --> M1["Stage 1 (Mixer-Settler)"]
    M1 --> M2["Stage 2 (Mixer-Settler)"]
    M2 --> Mn["Stage N (Mixer-Settler)"]
    Mn --> R["Exhausted Raffinate (R, x_n)"]

    S["Fresh Solvent (S, y_s)"] --> Mn
    Mn --> M2
    M2 --> M1
    M1 --> E["Rich Extract (E, y_1)"]

    style F fill:#e0f2fe,stroke:#0284c7
    style S fill:#fef3c7,stroke:#d97706
    style E fill:#dcfce7,stroke:#16a34a
    style R fill:#fee2e2,stroke:#dc2626

# 2.1. Analytical Kremser Equation

For linear equilibrium operating lines (y=KD⋅xy = K_D \cdot x), the required number of theoretical stages (NtheorN_{theor}) is derived via the Kremser-Brown-Souders equation:

Ntheor=ln⁡[(xf−ysKDxn−ysKD)(1−1E)+1E]ln⁡(E)N_{theor} = \frac{\ln\left[ \left( \frac{x_f - \frac{y_s}{K_D}}{x_n - \frac{y_s}{K_D}} \right) \left(1 - \frac{1}{E}\right) + \frac{1}{E} \right]}{\ln(E)}

# 2.2. Murphree Stage Efficiency (EME_M)

Because liquid-liquid diffusion is significantly slower than vapor-liquid diffusion, practical stage efficiency in mixer-settlers is between 75% to 90%75\% \text{ to } 90\%:

Nactual=NtheorEMN_{actual} = \frac{N_{theor}}{E_M}

# 3. Mixer Hydrodynamics, Power Input & Droplet Breakup

The mixer tank must disperse the solvent into droplets fine enough to maximize specific interfacial area (aa), but coarse enough to prevent terminal emulsification.

                  MIXER HYDRODYNAMICS & AGITATOR SCHEMATIC
               ┌──────────────────────────────────────────────┐
               │              Drive Motor & VFD               │
               └──────────────────────┬───────────────────────┘
                                      │  Shaft
                                  ┌───┴───┐
                                  │       │
             Solvent In ──► ──────┤       ├────── ◄── Aqueous In
                                  │   ▲   │
                                  │  / \  │
                                  │ /   \ │
                                  │       │
                                  └───┬───┘
                                 ┌────┴────┐
                                 │ Impeller│ (Curved Pumping Turbine)
                                 └────┬────┘
                                      ▼
                           Discharge to Settler Chute

# 3.1. Specific Interfacial Mass Transfer Area (aa)

The interfacial area per unit mixer volume is governed by the dispersed phase hold-up (ϕ\phi) and Sauter mean diameter (d32d_{32}):

a=6⋅ϕd32a = \frac{6 \cdot \phi}{d_{32}}

# 3.2. Sauter Mean Diameter (d32d_{32}) via Hinze-Kolmogorov Theory

In fully turbulent liquid dispersion, droplet diameter is determined by the balance of turbulent kinetic energy dissipation (ε\varepsilon) and interfacial tension (σ\sigma):

d32=0.058⋅Dimp⋅(1+5.4ϕ)⋅We−0.6d_{32} = 0.058 \cdot D_{imp} \cdot \left(1 + 5.4\phi\right) \cdot \text{We}^{-0.6}

Where the Impeller Weber number is:

We=ρc⋅N2⋅Dimp3σ\text{We} = \frac{\rho_c \cdot N^2 \cdot D_{imp}^3}{\sigma}

# 3.3. Pumping Turbine Power & Tip Speed Criteria

To pump and disperse without external inter-stage transfer pumps, plants employ curved-blade bottom-suction turbines:

P=Np⋅ρmix⋅N3⋅Dimp5P = N_p \cdot \rho_{mix} \cdot N^3 \cdot D_{imp}^5
  • Power Density Target: 0.8 to 1.5 kW/m30.8 \text{ to } 1.5\text{ kW/m}^3 of mixer volume.
  • Impeller Tip Speed (vtip=πDimpNv_{tip} = \pi D_{imp} N): Kept strictly within 2.5 to 3.8 m/s2.5 \text{ to } 3.8\text{ m/s}. Speeds >4.2 m/s>4.2\text{ m/s} shatter droplets below 10 μm10\,\mu\text{m}, creating stubborn emulsions.

# 4. Gravity Settler Sizing, Coalescence & Hydraulic Balance

The settler must provide laminar horizontal flow allowing droplets to rise or settle to the interface.

                      SETTLER COALESCENCE WEDGE PROFILE
 ┌────────────────────────────────────────────────────────────────────────┐
 │                      Light Phase Layer (Extract)                       │
 │~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~│
 │ ◄── Coalescence Wedge (Dispersed Band) ──►                             │
 │........................................................................│
 │                      Heavy Phase Layer (Raffinate)                     │
 └────────────────────────────────────────────────────────────────────────┘

# 4.1. Stokes' Law Terminal Velocity

For rigid dispersed spherical droplets in the continuum:

vt=g⋅∣ρc−ρd∣⋅ddrop218⋅μcv_t = \frac{g \cdot |\rho_c - \rho_d| \cdot d_{drop}^2}{18 \cdot \mu_c}

# 4.2. Settler Area Sizing via Specific Settling Flux

In industrial batteries, settlers are sized using the Total Specific Volumetric Flux (Φsettler\Phi_{settler}):

Asettler=Qlight+QheavyΦsettlerA_{settler} = \frac{Q_{light} + Q_{heavy}}{\Phi_{settler}}
Phase SystemDensity Difference (Δρ\Delta \rho)Operating Flux (Φsettler\Phi_{settler})Settler Aspect Ratio (L/WL/W)
Water / Toluene130 kg/m3130\text{ kg/m}^325−35 m3/(m2⋅h)25 - 35\text{ m}^3/(\text{m}^2\cdot\text{h})3:1−4:13:1 - 4:1
Water / Ethyl Acetate98 kg/m398\text{ kg/m}^318−25 m3/(m2⋅h)18 - 25\text{ m}^3/(\text{m}^2\cdot\text{h})4:1−5:14:1 - 5:1
Water / DCM (Heavy)326 kg/m3326\text{ kg/m}^335−50 m3/(m2⋅h)35 - 50\text{ m}^3/(\text{m}^2\cdot\text{h})3:1−4:13:1 - 4:1
Water / MIBK200 kg/m3200\text{ kg/m}^328−40 m3/(m2⋅h)28 - 40\text{ m}^3/(\text{m}^2\cdot\text{h})3:1−4:13:1 - 4:1

# 4.3. Hydraulic Underflow / Overflow Weir Equation

To maintain a stable liquid-liquid interface inside the settler without active electronics, an external adjustable inverted loop (heavy phase underflow leg) is balanced by hydrostatic head:

Hlight⋅ρlight+Hheavy⋅ρheavy=Hweir⋅ρheavyH_{light} \cdot \rho_{light} + H_{heavy} \cdot \rho_{heavy} = H_{weir} \cdot \rho_{heavy}
Hinterface=Hweir⋅ρheavy−Htotal⋅ρlightρheavy−ρlightH_{interface} = \frac{H_{weir} \cdot \rho_{heavy} - H_{total} \cdot \rho_{light}}{\rho_{heavy} - \rho_{light}}

# 5. Comprehensive Worked Industrial Case Study: 1,500 kg/h API Extraction

# Problem Statement:

An antibiotic intermediate (MW=385 g/molMW = 385\text{ g/mol}) is synthesized in an aqueous broth at Cf=4.0 wt%C_f = 4.0\text{ wt}\% (40 g/kg40\text{ g/kg}).

  • Feed rate: F=1,500 kg/hF = 1,500\text{ kg/h} (ρ=1,020 kg/m3\rho = 1,020\text{ kg/m}^3).
  • Extractant: Pure Ethyl Acetate (S/F=0.8S/F = 0.8, ρ=900 kg/m3\rho = 900\text{ kg/m}^3).
  • Partition coefficient: KD=6.2K_D = 6.2.
  • Target extraction recovery: ≥99.0%\ge 99.0\% (xn≤0.04 wt%x_n \le 0.04\text{ wt}\%).

# Step 1: Calculate Extraction Factor (EE)

E=KD⋅(SF)=6.2⋅0.8=4.96E = K_D \cdot \left(\frac{S}{F}\right) = 6.2 \cdot 0.8 = 4.96

# Step 2: Determine Theoretical Stages via Kremser Formula

Ntheor=ln⁡[(0.040−00.0004−0)(1−14.96)+14.96]ln⁡(4.96)=ln⁡[100⋅0.7984+0.2016]ln⁡(4.96)=ln⁡(80.04)1.6014=4.3821.6014=2.74N_{theor} = \frac{\ln\left[ \left( \frac{0.040 - 0}{0.0004 - 0} \right) \left(1 - \frac{1}{4.96}\right) + \frac{1}{4.96} \right]}{\ln(4.96)} = \frac{\ln[100 \cdot 0.7984 + 0.2016]}{\ln(4.96)} = \frac{\ln(80.04)}{1.6014} = \frac{4.382}{1.6014} = 2.74

With an 85% Murphree stage efficiency (EM=0.85E_M = 0.85):

Nactual=2.740.85=3.22  ⟹  4 Counter-Current Stages RequiredN_{actual} = \frac{2.74}{0.85} = 3.22 \implies \mathbf{4\text{ Counter-Current Stages Required}}

# Step 3: Mixer Dimensions & Agitation Power

  • Total volumetric flow:
Qtot=1,5001,020+1,200900=1.47+1.33=2.80 m3/h=0.000778 m3/sQ_{tot} = \frac{1,500}{1,020} + \frac{1,200}{900} = 1.47 + 1.33 = 2.80\text{ m}^3/\text{h} = 0.000778\text{ m}^3/\text{s}
  • Target mixer residence time τ=90 seconds\tau = 90\text{ seconds}:
Vmixer=Qtot⋅90 s=0.070 m3=70 LitersV_{mixer} = Q_{tot} \cdot 90\text{ s} = 0.070\text{ m}^3 = 70\text{ Liters}
  • Assuming square geometry (H=DtankH = D_{tank}):
Dtank=(4⋅0.070π)1/3=0.446 m≈450 mmD_{tank} = \left(\frac{4 \cdot 0.070}{\pi}\right)^{1/3} = 0.446\text{ m} \approx 450\text{ mm}
  • Impeller diameter (Dimp=Dtank/3=150 mmD_{imp} = D_{tank}/3 = 150\text{ mm}):
    • Operating speed N=350 rpm=5.83 rpsN = 350\text{ rpm} = 5.83\text{ rps}.
    • Tip speed vtip=π⋅0.150⋅5.83=2.75 m/sv_{tip} = \pi \cdot 0.150 \cdot 5.83 = 2.75\text{ m/s} (well within the safe non-emulsifying 2.5−3.8 m/s2.5 - 3.8\text{ m/s} window).
  • Motor Power (Np=2.2N_p = 2.2, ρmix≈960 kg/m3\rho_{mix} \approx 960\text{ kg/m}^3):
P=2.2⋅960⋅(5.83)3⋅(0.150)5=2,112⋅198.1⋅0.0000759=31.8 WattsP = 2.2 \cdot 960 \cdot (5.83)^3 \cdot (0.150)^5 = 2,112 \cdot 198.1 \cdot 0.0000759 = 31.8\text{ Watts}

With motor losses and viscous margins, install a 0.25 kW0.25\text{ kW} motor with VFD.

# Step 4: Settler Area & Dimensions

  • Using design flux Φsettler=20 m3/(m2⋅h)\Phi_{settler} = 20\text{ m}^3/(\text{m}^2\cdot\text{h}):
Asettler=2.80 m3/h20 m3/(m2⋅h)=0.140 m2A_{settler} = \frac{2.80\text{ m}^3/\text{h}}{20\text{ m}^3/(\text{m}^2\cdot\text{h})} = 0.140\text{ m}^2
  • Applying 3.5:13.5:1 length-to-width aspect ratio:
W=0.1403.5=0.20 m=200 mmW = \sqrt{\frac{0.140}{3.5}} = 0.20\text{ m} = 200\text{ mm}
L=3.5⋅0.20 m=0.70 m=700 mmL = 3.5 \cdot 0.20\text{ m} = 0.70\text{ m} = 700\text{ mm}

Depth H=300 mmH = 300\text{ mm} (liquid depth 250 mm250\text{ mm}).


# 6. Industrial Troubleshooting & Operational Mitigations

Failure ModePhysical CauseImmediate Root CauseEngineered Corrective Action
Rag Layer / Emulsion BandingInterfacial accumulation of denatured protein, fine cell debris, or colloidal silicaFine sub-micron particulates stabilize droplets (Pickering Emulsion)1. Add 1 μm1\,\mu\text{m} bag pre-filter before Stage 1.
2. Install structured fluoropolymer (PVDF) coalescence plates in the front third of the settler.
3. Periodic hot water flush (60∘C60^\circ\text{C}) to dissolve interfacial scum.
Phase Inversion (Flooding)Dispersed solvent suddenly becomes continuous phaseDispersed volume fraction ϕ\phi exceeded critical inversion threshold (>45%>45\%)1. Implement automatic feed ratio interlock on DCS.
2. Install conductivity probe in mixer: immediately alert operators if continuous phase flips from polar to non-polar.
Solvent Loss in RaffinateExcessive entrainment of tiny solvent droplets in aqueous underflowSettler velocity exceeds droplet terminal Stokes velocity1. Lengthen settler by adding downstream calming zone.
2. Lower impeller tip speed by 15%15\%.
3. Install an inline coalescing cartridge downstream of raffinate discharge.

# Applicable Engineering Standards & Codes Used

  • ISO 5167: Measurement of fluid flow in closed conduits.
  • ASME BPE: Bioprocessing Equipment cGMP design for pharmaceutical extraction skids.
  • API RP 520 / 521: Sizing pressure relief valves for solvent extraction systems.
  • ASTM D971: Standard Test Method for Interfacial Tension of Liquid against Water.
Process EngineeringSeparationLLEMixer-SettlerScale-UpDownstream Processing
Comments (0)

Discussion

Please Log In to participate in the technical discussion.

No comments posted yet. Be the first to share your input!