Crystal Yield Optimization Calculation

Reference ID: MET-7F81 | Process Engineering Reference Sheets Calculation Guide

Introduction & Context

Crystal yield optimization is a critical unit operation in process engineering, particularly within the pharmaceutical, food, and specialty chemical industries, and it often incorporates techniques such as melt crystallization for purification to maximize the recovery of solid product from a supersaturated solution through controlled cooling. By employing multi‑stage crystallization, engineers can significantly increase the total recovery of solutes that would otherwise remain dissolved in the mother liquor after a single cooling step.

This calculation is essential for designing efficient separation trains, determining the economic viability of secondary recovery stages, and minimizing product loss in waste streams. It is typically applied in the production of high-value crystalline solids such as lactose, citric acid, and various active pharmaceutical ingredients (APIs).

Methodology & Formulas

The calculation follows a sequential mass balance approach. First, solubility data provided in mass-based ratios (e.g., g solute per 100 g solvent) must be converted to a mass fraction (w/w) to align with the total mass balance of the system. This methodology assumes solvent mass is conserved (no evaporation) and that each stage reaches equilibrium at the specified temperature.

The conversion from solubility S to mass fraction w_sat is defined as:

\[ w_{\text{sat}} = \frac{S}{100 + S} \]

For each stage, the theoretical mass of crystals m_{c,theoretical} is determined by the difference between the feed concentration w_feed and the saturation concentration w_sat at the operating temperature:

\[ m_{\text{c,theoretical}} = m_{\text{total}} \cdot \frac{w_{\text{feed}} - w_{\text{sat}}}{1 - w_{\text{sat}}} \]

The actual crystal recovery m_c,actual accounts for the mechanical separation efficiency η of the centrifuge or filtration unit:

\[ m_{\text{c,actual}} = m_{\text{c,theoretical}} \cdot \eta \]

The mother liquor exiting a stage becomes the feed for the subsequent stage. The solute mass balance for the mother liquor is calculated as:

\[ m_{\text{solute,ml}} = m_{\text{solute,feed}} - m_{\text{c,actual}} \]

Parameter	Condition / Constraint
Separation Efficiency (η)	0.0 ≤ η ≤ 1.0
Crystallization Driving Force	w_feed − w_sat > 0 (positive supersaturation required)
Solubility (w_sat)	w_sat > 0

The overall system performance is evaluated by the cumulative recovery across all stages:

\[ \text{Yield}_{\text{overall}} = \frac{\sum m_{\text{c,actual}}}{m_{\text{solute,initial}}} \]

To calculate the theoretical maximum yield, you must evaluate the solubility limits of the solute at the starting and ending temperatures. Follow these steps:

Determine the initial concentration of the solute in the feed solution.
Identify the final solubility of the solute at the target cooling temperature.
Calculate the difference between the initial solute mass and the final solute mass remaining in the mother liquor at equilibrium.
Verify that the mass balance closes and that solvent mass is conserved (unless evaporation data are explicitly included in the model).

Several process parameters influence the final output. The most critical variables include:

Cooling rate profile, which affects the supersaturation level.
Seeding density and the timing of seed addition.
Impurity concentration, which can inhibit crystal growth kinetics.
Agitation intensity, which influences secondary nucleation and crystal breakage.

If your yield is lower than expected, investigate the following potential causes:

Check for incomplete cooling or premature termination of the batch.
Verify the accuracy of the solubility curve data used in your model.
Analyze the mother liquor for residual solute to determine if the system reached equilibrium.
Inspect the filter cake for excessive solvent retention or crystal dissolution during the washing stage.

Worked Example: Two-Stage Lactose Crystal Recovery

Scenario: A dairy processing plant operates a two-stage continuous cooling crystallizer to recover lactose from whey permeate. The feed stream enters at 1000 kg/h containing 30 wt% lactose. The first stage cools to 20°C, and the second stage further cools the mother liquor to 5°C. The separation efficiency in the centrifuge for the first and second stages are 95% and 92%, respectively. The objective is to calculate the actual crystal mass from each stage, total actual crystals, overall yield, and final mother liquor loss.

Knowns:

Feed total mass flow: \(1000.0 \, \text{kg/h}\)
Feed solute (lactose) mass: \(300.0 \, \text{kg/h}\)
Feed water mass: \(700.0 \, \text{kg/h}\)
Feed mass fraction: \(0.3\)
Stage 1 temperature: 20°C
Solubility at 20°C (mass fraction): \(w_{\text{sat,1}} = 0.138\)
Stage 1 separation efficiency: \(\eta_1 = 0.95\)
Stage 2 temperature: 5°C
Solubility at 5°C (mass fraction): \(w_{\text{sat,2}} = 0.09\)
Stage 2 separation efficiency: \(\eta_2 = 0.92\)

Step-by-Step Calculation:

Stage 1 Theoretical Crystals:
\[ m_{\text{c,theoretical,1}} = m_{\text{total,1}} \cdot \frac{w_{\text{feed,1}} - w_{\text{sat,1}}}{1 - w_{\text{sat,1}}} \] Given \(m_{\text{total,1}} = 1000.0 \, \text{kg/h}\) and \(w_{\text{feed,1}} = 0.3\): \[ m_{\text{c,theoretical,1}} = 1000.0 \cdot \frac{0.3 - 0.138}{1 - 0.138} = 1000.0 \cdot \frac{0.162}{0.862} = 187.9 \, \text{kg/h} \]
Stage 1 Actual Crystals:
\[ m_{\text{c,actual,1}} = m_{\text{c,theoretical,1}} \cdot \eta_1 = 187.9 \cdot 0.95 = 178.5 \, \text{kg/h} \]
Mother Liquor After Stage 1:
Solute remaining: \(m_{\text{solute,ML,1}} = 300.0 - 178.5 = 121.5 \, \text{kg/h}\)
Water remains: \(m_{\text{water,ML,1}} = 700.0 \, \text{kg/h}\)
Total mother liquor: \(m_{\text{total,ML,1}} = 121.5 + 700.0 = 821.5 \, \text{kg/h}\)
Feed mass fraction to stage 2: \(w_{\text{feed,2}} = \dfrac{121.5}{821.5} = 0.1479\)
Stage 2 Theoretical Crystals:
\[ m_{\text{c,theoretical,2}} = m_{\text{total,2}} \cdot \frac{w_{\text{feed,2}} - w_{\text{sat,2}}}{1 - w_{\text{sat,2}}} \] Given \(m_{\text{total,2}} = 821.5 \, \text{kg/h}\), \(w_{\text{feed,2}} = 0.1479\), and \(w_{\text{sat,2}} = 0.09\): \[ m_{\text{c,theoretical,2}} = 821.5 \cdot \frac{0.1479 - 0.09}{1 - 0.09} = 821.5 \cdot \frac{0.0579}{0.91} = 52.3 \, \text{kg/h} \]
Stage 2 Actual Crystals:
\[ m_{\text{c,actual,2}} = m_{\text{c,theoretical,2}} \cdot \eta_2 = 52.3 \cdot 0.92 = 48.1 \, \text{kg/h} \]
Total Actual Crystals:
\[ m_{\text{c,total}} = 178.5 + 48.1 = 226.6 \, \text{kg/h} \]
Overall Yield:
\[ Y_{\text{overall}} = \frac{226.6}{300.0} = 0.755 \]
Final Mother Liquor Loss:
\[ m_{\text{loss,final}} = 300.0 - 226.6 = 73.4 \, \text{kg/h} \]

Final Answer: The two-stage crystallizer produces 226.6 kg/h of lactose crystals (Stage 1: 178.5 kg/h, Stage 2: 48.1 kg/h), achieving an overall yield of 0.755 (75.5%). The final mother liquor contains 73.4 kg/h of unrecovered lactose.

"Un projet n'est jamais trop grand s'il est bien conçu." — André Citroën

"La difficulté attire l'homme de caractère, car c'est en l'étreignant qu'il se réalise." — Charles de Gaulle

Crystal Yield Optimization Calculation

Introduction & Context

Methodology & Formulas

Worked Example: Two-Stage Lactose Crystal Recovery

🚀 Skip the Manual Math!