🎓 Lesson 18
D5
Crystallizer Design: MSMPR Modeling and Residence Time Distribution
An MSMPR crystallizer is a tank where crystals grow steadily as fresh solution flows in and slurry flows out, keeping the crystal sizes evenly mixed over time.
🎯 Learning Objectives
- ✓ Calculate steady-state crystal size distribution (CSD) parameters—including mean crystal size and coefficient of variation—from MSMPR operating data
- ✓ Design an MSMPR crystallizer volume and residence time to achieve target crystal size and yield using population balance models
- ✓ Analyze residence time distribution (RTD) curves to diagnose mixing inefficiencies or bypassing in industrial crystallizers
- ✓ Apply Fick’s law and interfacial mass transfer principles to estimate growth rates from supersaturation profiles
- ✓ Explain how nucleation kinetics and residence time jointly determine crystal size and fines content in MSMPR output
📖 Why This Matters
In mining and mineral processing, crystallization is critical for recovering high-purity salts (e.g., potash, sodium sulfate) and metal compounds (e.g., nickel sulfate, lithium carbonate) from leach solutions. Poorly designed crystallizers produce off-spec crystals—too fine (hard to filter, high impurity carryover) or too large (slow dissolution in downstream processes). MSMPR modeling gives engineers predictive control over crystal size and shape before building expensive pilot plants—saving months of trial-and-error and millions in capital cost.
📘 Core Principles
MSMPR behavior rests on three pillars: (1) Perfect mixing ensures uniform supersaturation and temperature, eliminating concentration gradients; (2) Steady-state operation implies constant holdup, flow rates, and crystal size distribution; (3) Population balance modeling links nucleation (birth of new crystals) and growth (size increase) to RTD via the moment-based method. The zeroth moment (total crystal number) and first moment (total crystal length) govern key outputs: nucleation rate (B₀) and growth rate (G). RTD analysis—using tracer experiments—validates whether real equipment behaves like an ideal MSMPR; deviations signal dead zones or short-circuiting that distort CSD.
📐 Key Calculation
The mean crystal size (L̄) in an ideal MSMPR is directly proportional to residence time (τ) and growth rate (G), and inversely related to nucleation rate (B₀). This arises from solving the population balance under steady-state, constant-G assumptions.
Mean Crystal Size (MSMPR)
L̄ = G · τPredicts average crystal length under steady-state, constant growth rate assumption.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| L̄ | Mean crystal size | mm | Length of crystals at steady state |
| G | Linear growth rate | mm/h | Rate of crystal size increase per unit time |
| τ | Residence time | h | Average time slurry remains in crystallizer |
Typical Ranges:
Potash (KCl) MSMPR: 0.6 – 1.2 mm
Sodium sulfate decahydrate: 0.3 – 0.8 mm
💡 Worked Example
Problem: A potash MSMPR crystallizer operates at τ = 4.5 h. Lab tests show G = 0.18 mm/h and B₀ = 1.2 × 10⁶ #/(m³·h). Calculate L̄ and CV (coefficient of variation).
1.
Step 1: Apply L̄ = G·τ → L̄ = 0.18 mm/h × 4.5 h = 0.81 mm
2.
Step 2: For ideal MSMPR with constant G, CV = 0.5 (theoretical value from exponential CSD)
3.
Step 3: Verify consistency: B₀ = 4G / L̄² → RHS = 4×0.18/(0.81)² ≈ 1.098 × 10⁶ #/(m³·h), close to measured 1.2 × 10⁶ — acceptable within ±10% experimental error.
Answer:
The mean crystal size is 0.81 mm, with CV = 0.5 — indicating a broad, exponential distribution typical of MSMPR. This meets potash filtration specs (target L̄ > 0.7 mm, CV < 0.6).
🏗️ Real-World Application
At the Mosaic Co.’s Carlsbad potash facility, engineers redesigned their MSMPR crystallizers after field CSD measurements showed excessive fines (<0.3 mm) causing belt-filter blinding. RTD studies using NaCl tracer revealed E-curve skewness (θₚₑₐₖ/τ = 0.68), indicating 22% short-circuiting. By installing baffles and adjusting impeller tip speed from 2.1 to 2.8 m/s, mixing efficiency improved (θₚₑₐₖ/τ → 0.95), increasing L̄ from 0.52 mm to 0.79 mm and reducing fines by 65% — meeting ISO 14827 filtration standards without changing residence time.
📋 Case Connection
📋 Pharmaceutical API Purification via Crystallization
Polymorphic instability and residual solvent > ICH Q3C limits (e.g., acetone > 5000 ppm)