Size, model, and benchmark continuous Mixed-Suspension Mixed-Product-Removal (MSMPR) Crystallizers using the Randolph-Larson Population Balance Model (PBM). Computes crystal linear growth rates ((G)), secondary nucleation rates ((B_0)), dominant size ((L_D)), cumulative mass size distributions (CSD), magma suspension density ((M_T)), and fines dissolution loop improvements.
1. Growth & Nucleation Kinetics
2. Crystallizer Vessel Sizing
3. CSD Moments & Magma Density
Continuous MSMPR Crystallizer & Population Density Curve (CSD)
Real-time animated visualization displaying draft-tube agitator suspension, crystal growth kinetics, and logarithmic population density profile (ln n(L)) vs crystal size (L) illustrating dominant modal peak.
5 Fatal Traps & Industrial Pitfalls in Industrial Crystallizers
1. Secondary Contact Nucleation & Uncontrollable Fines Explosion
Increasing crystallizer impeller RPM to ensure solids suspension accelerates crystal-impeller and crystal-baffle collisions. Secondary collision nucleation birth rate scales as (B_0 propto N^h cdot M_T^j) (where (h approx 2 ext{ to }3)). A 25% increase in impeller speed can trigger a ten-fold explosion in microscopic crystal nuclei (<20 μm). These trillions of fine crystals consume all dissolved supersaturation before existing crystals can grow, causing product dominant size (L_D) to collapse, choking downstream centrifuge basket screens and blinding filter cloths.
2. Cold-Wall Encrustation & Thermal Transfer Choking
Chasing higher production rates by lowering the cooling jacket temperature introduces high local supersaturation at the internal stainless steel vessel wall. When the temperature difference between the bulk slurry and jacket coolant exceeds 3.5°C to 5.0°C, the boundary layer crosses the metastable limit, causing spontaneous primary nucleation and crystallization directly on the metal wall. A 3 mm thick crystalline crust has an exceptionally low thermal conductivity ((k approx 1.2 ext{ W/m}cdot ext{K})), cutting heat transfer by 80% and necessitating complete batch boil-out shutdowns.
3. High-Shear Draft Tube Cavitation & Crystal Cleavage Shattering
Draft-Tube Baffle (DTB) crystallizers utilize an axial flow marine propeller or hydrofoil impeller to circulate slurry up the draft tube. If the impeller tip speed exceeds 3.5 m/s or blade clearance is improper, localized hydrodynamic cavitation and high shear stress cleave large, fragile inorganic salt crystals along their crystalline lattice planes. The resulting shattered fragments distort the ideal exponential MSMPR size distribution into an irregular, bimodal distribution with poor flowability and severe cake agglomeration.
4. Magma Density Overload (>25 vol%) & Transfer Line Sanding
Operating with a slurry solids concentration exceeding 22% to 25% by volume sharply increases slurry apparent viscosity. If product discharge slurry pumps or gravity transfer piping drop below the critical Durand settling velocity ((v_{ ext{crit}} approx 1.5 ext{ to }2.2 ext{ m/s})), coarse crystals immediately drop out of suspension. The stationary bed quickly dams the pipe, causing instant hydraulic sanding and line blockage that requires manual dismounting and high-pressure water jetting.
5. Impurity Mother Liquor Inclusions at High Linear Growth Rates
Pushing linear growth rate (G) beyond critical limits ((>5 imes 10^{-8} ext{ m/s})) destabilizes planar crystal face advancement. Dendritic or step-bunching growth traps microscopic droplets of mother liquor inside internal crystal vacuoles ("inclusions"). Because the contaminated solvent is physically sealed inside the crystal interior, no amount of cake displacement washing on the centrifuge can remove residual mother liquor, resulting in severe chemical product purity specification failure.
Randolph-Larson Population Balance Equations & Moments
The steady-state 1D Population Balance Equation (PBE) for an unseeded, continuous mixed-suspension mixed-product-removal (MSMPR) crystallizer assuming size-independent growth (McCabe (Delta L) law) is:
1. Differential Population Balance Equation
d(G · n) / dL + n / tau = 0
Integrating with boundary condition at nuclei birth (L
ightarrow 0) ((n(0) = n_0 = B_0 / G)):
n(L) = n_0 · exp(-L / (G · tau))
where:
• n(L): Population density ((# / ( ext{m} cdot ext{m}^3)))
• G: Linear growth rate (dL/dt) (m/s)
• tau: Mean slurry residence time (V / Q) (seconds)
• B₀: Total nucleation rate ((# / ( ext{m}^3 cdot ext{s})))
2. Moments of the Distribution
The (j)-th moment of the crystal population density:
m_j = integral[0 to infinity] L^j · n(L) dL = j! · n_0 · (G · tau)^(j+1)
• Total crystal number: m_0 = n_0 · G · tau = B_0 · tau ((# / ext{m}^3))
• Total crystal length: m_1 = n_0 · (G · tau)² (( ext{m} / ext{m}^3))
• Total crystal surface area: A_c = k_a · m_2 = 2 · k_a · n_0 · (G · tau)³ (( ext{m}^2 / ext{m}^3))
• Total crystal volume: V_c = k_v · m_3 = 6 · k_v · n_0 · (G · tau)⁴ (( ext{m}^3 / ext{m}^3))
3. Slurry Magma Density & Dominant Size
Slurry suspension magma density (M_T) (kg/m³):
M_T = rho_c · k_v · m_3 = 6 · k_v · rho_c · n_0 · (G · tau)⁴ = 6 · k_v · rho_c · B_0 · G³ · tau⁴
Dominant crystal size (L_D) (peak of differential mass distribution):
L_D = 3 · G · tau
Median mass diameter (L_{50} = 3.67 cdot G cdot au).
4. Cumulative Mass Undersize Distribution W(L)
W(L) = 1 - exp(-L / (G·tau)) · [1 + L/(G·tau) + 1/2 · (L/(G·tau))² + 1/6 · (L/(G·tau))³]