Continuous MSMPR Crystallizer & Population Balance Calculator
Design industrial Mixed Suspension Mixed Product Removal (MSMPR) crystallizers: crystal size distribution (CSD), population density n(L), dominant crystal size LD, nucleation rate B0, and production capacity.
1. Crystallizer Operating Conditions
2. Crystal & Fluid Physical Properties
Population Balance & CSD Metrics
Crystal Size Distribution w(L) & Population Density n(L)
5 Fatal Engineering Traps in Industrial Crystallizer Operation
1. Metastable Zone Width (MZW) Breach & Fines Shower Runaway
Pushing evaporation or cooling beyond the critical supersaturation threshold triggers catastrophic primary homogeneous nucleation. Instead of steady crystal growth onto existing seed magma, a milky cloud of sub-5 micron fines erupts instantaneously. This fine slurry clogs centrifuge baskets, produces wet un-washable filter cakes, and requires hours of steam dissolution to recover.
2. Impeller Tip Speed Shear & Attrition-Driven Secondary Nucleation
Agitator tip speeds exceeding 2.5–3.0 m/s collide violently with suspension crystals, fracturing edges and generating microscopic secondary nuclei. This artificially inflates nucleation rate (B0), forcing the dominant crystal size (LD) to plummet. Proper crystallizer design mandates low-speed axial marine impellers with draft tubes to circulate high magma volumes at minimal mechanical shear.
3. Non-Isokinetic Slurry Withdrawal & Population Balance Corruption
If discharge slurry velocity through the dip pipe does not match local vessel magma velocity, selective sampling occurs: large crystals are excluded while liquid and fines are preferentially pulled, or vice versa. This violates the foundational MSMPR assumption that discharge slurry composition mirrors bulk vessel slurry, corrupting measured residence times and particle size analyzers.
4. Heat Exchanger Encrustation & Thermal Shock Scaling
In cooling crystallizers, maintaining an excessive temperature delta between the bulk magma and the cooling jacket or external heat exchanger tubes causes local supersaturation to skyrocket at the tube wall. Solute precipitates directly onto metal heat transfer surfaces, forming a rock-hard crystalline scale layer that rapidly destroys thermal conductivity and forces costly chemical boil-outs.
5. Slurry Stagnation & Dead-Zone Gravity Settling Salting Out
Operating below the Zwietering minimum suspension speed (Njs) allows coarse crystals to settle into dish-bottom dead zones around discharge valves. Stagnant crystals consolidate into monolithic boulders that bind draft-tube impellers upon restart, burn out electric motors, and shear off drive shaft couplings.
Mathematical Derivations & Population Balance Equations (Randolph & Larson)
For an ideal continuous MSMPR crystallizer operating at steady state without crystal breakage or agglomeration, the one-dimensional population balance equation reduces to:
d[ n(L) ] / dL = - n(L) / (G · τ)
Integrating with boundary condition at zero size n(0) = n0 = B0 / G:
n(L) = n0 · exp( - L / (G · τ) )
The j-th moment of the population density distribution is defined as mj = ∫0∞ Lj n(L) dL = j! · n0 · (G · τ)j+1.
Total crystal mass concentration (magma density MT):
MT = ρc · kv · m3 = 6 · ρc · kv · n0 · (G · τ)4
The mass-weighted crystal size distribution frequency w(L) and dominant crystal size LD:
w(L) = [ ρc · kv · L3 · n(L) ] / MT = [ L3 / (6 · (G · τ)4) ] · exp( - L / (G · τ) )
dw/dL = 0 → LD = 3 · G · τ (Mode of mass distribution)