MaximaLabs
All unit operations

Extraction column simulation

feed + solvent -> extract + raffinate (countercurrent LLE)

Governing equations

The exact equations the solver works for a extraction column — the same math shown in the app's "Theory" panel, not a black box.

Ei=KD,iSR(extraction factor; KD=extract/raffinate solute ratio)E_i = \frac{K_{D,i}\,S}{R}\quad(\text{extraction factor};\ K_D = \text{extract/raffinate solute ratio})
rj1+Erj+1=(1+E)rj(per-stage solute balancetridiagonal cascade)r_{j-1} + E\,r_{j+1} = (1+E)\,r_j\quad(\text{per-stage solute balance} \Rightarrow \text{tridiagonal cascade})
rNr0=E1EN+11(Kremser limit at constant KD; 1N+1 at E=1)\frac{r_N}{r_0} = \frac{E-1}{E^{N+1}-1}\quad(\text{Kremser limit at constant } K_D;\ \tfrac{1}{N+1} \text{ at } E=1)
KDK_D
distribution coefficient — solute concentration in the extract over the raffinate; the driving force of the whole separation
SS
solvent (extract-phase) molar flow [mol/s]
RR
carrier (raffinate-phase) molar flow [mol/s]
EE
extraction factor K_D·S/R — the analog of the absorber's A; >1 favours the solvent
NN
number of equilibrium stages
rjr_j
solute molar flow leaving stage j in the raffinate [mol/s]
r0r_0
solute fed to the cascade [mol/s]
rNr_N
solute left unrecovered in the final raffinate [mol/s]

Parameters

n_stages [int >=1], distribution_coefficients {solute: K_D = extract/raffinate ratio}, pressure/temperature [optional] — two inlets (feed, solvent) -> extract + raffinate. Countercurrent LLE cascade; recovery rises with stages (Kremser limit at constant K_D)

Example flowsheets that use it

Stop fighting legacy software. Build your first flowsheet in 60 seconds.