MaximaLabs
All unit operations

Crystallizer simulation

1 feed -> solid crystals + mother liquor

Governing equations

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

B=kbΔCb,    G=kgΔCg,    ΔC=CCsatB = k_b\,\Delta C^{b},\;\; G = k_g\,\Delta C^{g},\;\; \Delta C = C-C_{sat}
m0=Bτ,  m1=m0Gτ,  m2=2m1Gτ,  m3=3m2Gτ,Lˉ=m1/m0=Gτm_0=B\tau,\;m_1=m_0 G\tau,\;m_2=2m_1 G\tau,\;m_3=3m_2 G\tau,\quad \bar L=m_1/m_0=G\tau
m˙solid=ρckvm3,Q=n˙solidΔHcryst\dot m_{solid}=\rho_c k_v m_3,\quad Q = -\dot n_{solid}\,\Delta H_{cryst}
BB
nucleation rate — how many new crystals appear [#/m³/s]
GG
linear growth rate [m/s]
kb,kgk_b, k_g
nucleation and growth rate constants
ΔC\Delta C
supersaturation [mol/m³] — the driving force for both nucleation and growth
C,CsatC, C_{sat}
actual and saturation concentrations [mol/m³]
b,gb, g
nucleation and growth order
m0..m3m_0..m_3
moments of the crystal size distribution (count, length, area, volume)
τ\tau
residence time [s]
Lˉ\bar L
mean crystal size [m] = Gτ — the product-quality number
ρc\rho_c
crystal density [kg/m³]
kvk_v
volume shape factor
m˙solid\dot m_{solid}
crystal mass production [kg/s]
ΔHcryst\Delta H_{cryst}
heat of crystallization [J/mol]
QQ
duty [W]

Parameters

solute [feed component]; MSMPR: csat [mol/L], feed_concentration [mol/L], residence_time [s], optional kb/b/kg/g kinetics, crystal_density [kg/m^3], kv, molar_mass [kg/mol] — OR short-cut yield [0..1]; optional temperature [K], heat_of_crystallization

Example flowsheets that use it

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