MaximaLabs
All unit operations

Chromatography capture (Protein A) simulation

feed -> eluate (product pool) + flowthrough

Governing equations

The exact equations the solver works for a chromatography capture (protein a) — the same math shown in the app's "Theory" panel, not a black box.

cc0=11+ek(qmaxz/uc0t)(Bohart-Adams)\frac{c}{c_0} = \frac{1}{1 + e^{k(q_{max} z/u - c_0 t)}}\quad(\text{Bohart-Adams})
DBC=qmaxzu/tb  basis,Vcol=mDBCfsafetyncycles\mathrm{DBC} = \frac{q_{max} z}{u}\Big/ t_{b}\;\text{basis},\quad V_{col} = \frac{m}{\mathrm{DBC}\cdot f_{safety}\cdot n_{cycles}}
c0c_0
feed titre [kg/m^3]
qmaxq_{max}
equilibrium resin capacity [kg/m^3]
kk
Bohart-Adams rate constant [m^3/(kg.s)]
zz
bed height [m]
uu
superficial velocity [m/s]
tbt_b
time to the breakthrough fraction [s]
DBC\mathrm{DBC}
dynamic binding capacity [kg/m^3] — the capacity actually usable
VcolV_{col}
column volume [m^3]
mm
product mass in the harvest [kg]
fsafetyf_{safety}
load safety factor
ncyclesn_{cycles}
cycles per batch

Parameters

product [component id], product_molar_mass [kg/mol], q_max [kg/m^3 of resin], rate_constant [m^3/(kg·s)] — the last two measured for your molecule on your resin. Optional: bed_height [m], superficial_velocity [m/s], breakthrough_fraction, max_cycles_per_batch, load_safety_factor, step_yield, batch_time [s], equilibration/wash/elution/strip_volumes [CV], bound_fractions [{component: 0-1} for anything else the resin binds], resin_price, resin_lifetime_cycles.

Example flowsheets that use it

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