Crystallization, Nucleation and Crystal Growth Kinetics

Engineering context

A batch cooling profile produces a high apparent crystal yield but highly variable filtration time. Sudden nucleation may create fines even when the final equilibrium solubility predicts strong recovery. A controlled seed charge and slower cooling ramp can move crystal growth onto existing surfaces, producing a narrower size distribution and improving downstream separation. The trade-off is batch time and possible mother-liquor losses. Thermal history and impurity effects are as important as final solution temperature.

Unit learning outcomes

1. Develop and check supersaturation, solubility and metastable zones using S=C/C_sat(T); identify the principal assumptions and explain an experimental validation path.
2. Develop and check primary nucleation, secondary nucleation and growth using G=k_g(DeltaC)^g; identify the principal assumptions and explain an experimental validation path.
3. Develop and check batch cooling and crystal mass yield using m_crystal=m_solvent(C_initial-C_final); identify the principal assumptions and explain an experimental validation path.
4. Develop and check population balance and crystal-size distribution using n(L)=n_0 exp[-L/(G tau)]; identify the principal assumptions and explain an experimental validation path.

Concept map

1. Supersaturation, solubility and metastable zones

Supersaturation is the driving force for nucleation and growth. The equilibrium solubility curve depends on solvent composition and impurities; metastable zone width depends on cooling rate, mixing, vessel surfaces and history. A supersaturation ratio above one indicates a thermodynamic tendency to precipitate, but it does not prove an immediate nucleation event. Rapid cooling can create many small crystals, while controlled seeding can favor growth on existing crystals.

Central relation: S=C/C_sat(T). Operating boundary: Linearized solubility curves are illustrations only; real systems require measured equilibrium data.

2. Primary nucleation, secondary nucleation and growth

Primary nucleation can occur without crystals already present, while secondary nucleation is promoted by particles, collisions or attrition. In a mixed suspension, the number of nuclei affects the eventual particle-size distribution and downstream filtration performance. Growth kinetics often depend on supersaturation, temperature and impurity adsorption on specific crystal faces. A growth-only law cannot represent agglomeration or breakage.

Central relation: G=k_g(DeltaC)^g. Operating boundary: The example is for a fictitious compound with illustrative kinetics, not a general law for all crystals.

3. Batch cooling and crystal mass yield

A cooling crystallizer may not reach equilibrium before filtration. The mother liquor also retained within a filter cake lowers apparent recovery and contaminates product. A valid material balance should account for evaporated solvent, dissolved impurities and possible polymorphic transformations. A large theoretical yield can be economically poor if submicron crystals require energy-intensive separation.

Central relation: m_crystal=m_solvent(C_initial-C_final). Operating boundary: The calculation treats the crystallizing species as the only solute and assumes no solvent change.

4. Population balance and crystal-size distribution

In an ideal MSMPR crystallizer with size-independent growth, the steady-state number population decreases exponentially with size. The mean characteristic size scale is G tau, but that is not directly equal to the volume-weighted mean diameter. Crystal-size moments have different physical meanings: number, length, surface area and volume are weighted by powers of L. Measurements of chord length require care because they are not necessarily direct crystal-size distributions.

Central relation: n(L)=n_0 exp[-L/(G tau)]. Operating boundary: This ideal MSMPR population lacks aggregation, breakage, classification and shape dependence.

What to test experimentally

Establish the basis. A useful calculation begins with a consistent basis: mass rate, molar rate, solution volume, dry adsorbent mass or solvent mass. For equipment sizing a rate must not be confused with a cumulative inventory. Record every conversion as a separate line. Does the result represent an instantaneous rate, a time-integrated amount, a concentration, or a dimensionless ratio?

Identify equilibrium and rate assumptions. Equilibrium constrains an attainable composition or loading; kinetics determines how quickly a process approaches that constraint. Where an empirical coefficient is used, report its calibration range, units and experimental provenance. Which measured observation would show that an assumed equilibrium relation is not applicable?

Connect the model to a unit operation. A local model should be embedded within balances, operating constraints and product targets. Report useful recovery, purity and waste flows on comparable bases rather than optimizing a single local response. Which outlet becomes the limiting waste stream when product recovery increases?

Test uncertainty and reversal. If two alternatives are close, measurement error may reverse their apparent ranking. Estimate the sensitivity of performance to the largest uncertain inputs, and check that the conclusion survives a credible parameter range. Could an operating choice preferred at the reference feed become inferior when feed quality changes?

Laboratories supporting this unit

• [Lab 15: Cooling crystallization supersaturation](../labs/lab-15.html)
• [Lab 16: Crystal yield from batch cooling](../labs/lab-16.html)

Further work

• [Unit 5 assessment](../assessments/unit-05.html)
• [Unit 5 engineering project](../projects/unit-05.html)

Unit 5, Lesson 1 — Supersaturation, solubility and metastable zones

Purpose and operating question

Crystallization begins with a supersaturated solution, not simply a low temperature. Why this matters: A batch cooling profile produces a high apparent crystal yield but highly variable filtration time. Sudden nucleation may create fines even when the final equilibrium solubility predicts strong recovery. A controlled seed charge and slower cooling ramp can move crystal growth onto existing surfaces, producing a narrower size distribution and improving downstream separation. The trade-off is batch time and possible mother-liquor losses. Thermal history and impurity effects are as important as final solution temperature.

Physical interpretation

Supersaturation is the driving force for nucleation and growth. The equilibrium solubility curve depends on solvent composition and impurities; metastable zone width depends on cooling rate, mixing, vessel surfaces and history. A supersaturation ratio above one indicates a thermodynamic tendency to precipitate, but it does not prove an immediate nucleation event. Rapid cooling can create many small crystals, while controlled seeding can favor growth on existing crystals.

The principal mathematical relationship used in this lesson is

S=C/C_sat(T)

The equation is a compact model rather than an instruction to ignore transport resistance, property dependence or other equipment conditions. Before using it, translate every mathematical term into a measurable quantity. Identify whether the driving force is a local value at an interface, an average over equipment, or a measured inlet-to-outlet difference. Units and operating conditions must accompany any quoted property coefficient.

Derivation strategy and limitations

Start from the model: S=C/C_sat(T). Define the control volume or interface. Separate a measurable inlet/outlet quantity from an unobservable local state; do not fit a transport parameter to a mixed-stream observation without stating the model.

The physical assumptions matter more than a successful algebraic manipulation. If the property varies with temperature or composition, calculate locally or integrate over the changing state rather than using a single value without justification. For time-dependent problems, distinguish a true transient solution from a sequence of pseudo-steady snapshots. For multiphase systems, check that a concentration or composition basis has not silently changed from mole to mass fraction.

Linearized solubility curves are illustrations only; real systems require measured equilibrium data.

Worked investigation

A hypothetical solute has solubility 180 g/kg solvent at 50 C and 105 g/kg at 25 C. A batch containing 155 g/kg is cooled from 50 to 25 C. Calculate the initial and final S and explain when seeded crystallization becomes feasible. How does seed surface area influence induction time and crystal size?

Solution approach: First list the information provided, and convert all values to consistent units. Substitute only the values that are supported by the statement; identify any additional properties or conditions that would need measurement for a full numerical result. Use the central relation S=C/C_sat(T) as the mathematical starting point, then evaluate a limiting or reference case. The final conclusion should distinguish an explicitly calculated response from a hypothetical plant-level extrapolation. Do not invent the unspecified flow, area, equilibrium or mass-transfer parameters to make a single definitive answer appear possible.

Engineering decision and diagnostic testing

Measure a full material balance over a representative interval and identify sensors whose uncertainty can invert the claimed trend. Measure a full material balance over a representative interval and identify sensors whose uncertainty can invert the claimed trend.

A meaningful experiment uses paired samples or simultaneous inlet/outlet measurements at a steady operating condition. Record the sampling interval and the instrument calibration or reference standard. If the observed response changes as throughput rises, test whether the change originates in the material property, the boundary layer, equilibrium constraints or residence-time distribution. When comparing equipment, express the outcome as recovery, product purity, energy or material consumption on the same basis.

Failure mechanisms and how to investigate them

• Inconsistent basis: Recalculate using one mass or molar basis and check whether every outflow is physically nonnegative. In a composition balance, sum the species to unity after solving.
• A parameter outside its fitted range: Compare the new operating condition with the measurements that established the coefficient. An apparently precise result may be unreliable under a different feed chemistry.
• An unmodeled resistance: Determine whether external transfer, internal diffusion, hydraulic pressure losses or phase equilibrium changes could constrain the response.
• An ambiguous measurement: Confirm the sample location and timing. An apparent difference in purity may be caused by mixing product samples across a transition.
• A neglected waste or utility duty: Draw the full separation train and account for regenerant, wash solvent, backwash, concentrate or mechanical power as applicable.

Sensitivity study

Retain the reference case from the worked investigation, then change the principal driving force by +10% while holding the model's other independent parameters fixed. Repeat for the strongest resistance or equilibrium coefficient. If the response is strongly nonlinear, a percentage change in the input does not imply the same percentage change in product. Where the result is dimensionless, verify that it remains within physical bounds. Where the result is a cumulative amount, confirm the time and inventory basis.

Fully worked numerical demonstration

At 50 °C the specified solubility is 180 g/kg solvent. A 155 g/kg solution has S=155/180=0.861 and is undersaturated at this temperature. At 25 °C the solubility is 105 g/kg and S=155/105=1.476, so crystallization is thermodynamically favored. Yet the existence of supersaturation does not establish nucleation time: metastable zone width, mixing, seeds and impurities govern kinetics. An experimental cooling schedule can use a seed charge when supersaturation enters a controlled region to avoid uncontrolled fine-crystal formation.

Questions for independent work

1. Explain the physical meaning and units of every symbol in S=C/C_sat(T) and identify one property that would need experimental calibration.
2. Solve the worked investigation with a justified perturbation to its largest uncertain input, then report the direction of the result and a dimensional check.
3. Describe a realistic counterexample where the relation would fail even though all arithmetic was correctly implemented.
4. Develop a plant-level material balance that includes both recovered product and rejected or retained material.
5. Compare the model's predicted effect of increasing the driving force with an actual equipment constraint and state which measurement would reveal the discrepancy.

Further study

Read the related unit chapter and complete the unit assessment. Use the relevant simulation in the [laboratory directory](../index.html#labs) to investigate a parameter change. Keep the model limitations visible when drawing an engineering conclusion.

Unit 5, Lesson 2 — Primary nucleation, secondary nucleation and growth

Purpose and operating question

Nucleation determines number density while growth changes individual crystal mass. Why this matters: A batch cooling profile produces a high apparent crystal yield but highly variable filtration time. Sudden nucleation may create fines even when the final equilibrium solubility predicts strong recovery. A controlled seed charge and slower cooling ramp can move crystal growth onto existing surfaces, producing a narrower size distribution and improving downstream separation. The trade-off is batch time and possible mother-liquor losses. Thermal history and impurity effects are as important as final solution temperature.

Physical interpretation

Primary nucleation can occur without crystals already present, while secondary nucleation is promoted by particles, collisions or attrition. In a mixed suspension, the number of nuclei affects the eventual particle-size distribution and downstream filtration performance. Growth kinetics often depend on supersaturation, temperature and impurity adsorption on specific crystal faces. A growth-only law cannot represent agglomeration or breakage.

The principal mathematical relationship used in this lesson is

G=k_g(DeltaC)^g

The equation is a compact model rather than an instruction to ignore transport resistance, property dependence or other equipment conditions. Before using it, translate every mathematical term into a measurable quantity. Identify whether the driving force is a local value at an interface, an average over equipment, or a measured inlet-to-outlet difference. Units and operating conditions must accompany any quoted property coefficient.

Derivation strategy and limitations

Start from the model: G=k_g(DeltaC)^g. Choose a driving-force definition that matches the constitutive equation. A concentration gradient, pressure difference and electrochemical potential difference are not interchangeable.

The physical assumptions matter more than a successful algebraic manipulation. If the property varies with temperature or composition, calculate locally or integrate over the changing state rather than using a single value without justification. For time-dependent problems, distinguish a true transient solution from a sequence of pseudo-steady snapshots. For multiphase systems, check that a concentration or composition basis has not silently changed from mole to mass fraction.

The example is for a fictitious compound with illustrative kinetics, not a general law for all crystals.

Worked investigation

Assume growth velocity is 0.06 mm/h at a supersaturation driving force of 10 g/kg, with first-order dependence. Estimate growth velocity at 6 and 15 g/kg. For a constant 0.06 mm/h and six-hour batch, calculate ideal change in crystal length. Discuss why a real crystal population does not all grow by the same amount.

Solution approach: First list the information provided, and convert all values to consistent units. Substitute only the values that are supported by the statement; identify any additional properties or conditions that would need measurement for a full numerical result. Use the central relation G=k_g(DeltaC)^g as the mathematical starting point, then evaluate a limiting or reference case. The final conclusion should distinguish an explicitly calculated response from a hypothetical plant-level extrapolation. Do not invent the unspecified flow, area, equilibrium or mass-transfer parameters to make a single definitive answer appear possible.

Engineering decision and diagnostic testing

Hold all parameters but the candidate driving force constant; then document whether the observed output changes linearly, monotonically or nonmonotonically. Hold all parameters but the candidate driving force constant; then document whether the observed output changes linearly, monotonically or nonmonotonically.

A meaningful experiment uses paired samples or simultaneous inlet/outlet measurements at a steady operating condition. Record the sampling interval and the instrument calibration or reference standard. If the observed response changes as throughput rises, test whether the change originates in the material property, the boundary layer, equilibrium constraints or residence-time distribution. When comparing equipment, express the outcome as recovery, product purity, energy or material consumption on the same basis.

Failure mechanisms and how to investigate them

• Inconsistent basis: Recalculate using one mass or molar basis and check whether every outflow is physically nonnegative. In a composition balance, sum the species to unity after solving.
• A parameter outside its fitted range: Compare the new operating condition with the measurements that established the coefficient. An apparently precise result may be unreliable under a different feed chemistry.
• An unmodeled resistance: Determine whether external transfer, internal diffusion, hydraulic pressure losses or phase equilibrium changes could constrain the response.
• An ambiguous measurement: Confirm the sample location and timing. An apparent difference in purity may be caused by mixing product samples across a transition.
• A neglected waste or utility duty: Draw the full separation train and account for regenerant, wash solvent, backwash, concentrate or mechanical power as applicable.

Sensitivity study

Retain the reference case from the worked investigation, then change the principal driving force by +10% while holding the model's other independent parameters fixed. Repeat for the strongest resistance or equilibrium coefficient. If the response is strongly nonlinear, a percentage change in the input does not imply the same percentage change in product. Where the result is dimensionless, verify that it remains within physical bounds. Where the result is a cumulative amount, confirm the time and inventory basis.

Fully worked numerical demonstration

At first-order crystal-growth dependence, scaling the supersaturation driving force from 10 to 6 g/kg lowers velocity from 0.06 to 0.036 mm/h; increasing it to 15 g/kg raises velocity to 0.09 mm/h. If velocity stays at the reference 0.06 mm/h for six hours, crystal length grows by 0.36 mm. The result refers to length growth along a chosen crystal dimension, not necessarily a radius or diameter. A growing population competes for solute, so actual supersaturation may decrease and the integrated length increase will be smaller than the constant-rate estimate.

Questions for independent work

1. Explain the physical meaning and units of every symbol in G=k_g(DeltaC)^g and identify one property that would need experimental calibration.
2. Solve the worked investigation with a justified perturbation to its largest uncertain input, then report the direction of the result and a dimensional check.
3. Describe a realistic counterexample where the relation would fail even though all arithmetic was correctly implemented.
4. Develop a plant-level material balance that includes both recovered product and rejected or retained material.
5. Compare the model's predicted effect of increasing the driving force with an actual equipment constraint and state which measurement would reveal the discrepancy.

Further study

Read the related unit chapter and complete the unit assessment. Use the relevant simulation in the [laboratory directory](../index.html#labs) to investigate a parameter change. Keep the model limitations visible when drawing an engineering conclusion.

Unit 5, Lesson 3 — Batch cooling and crystal mass yield

Purpose and operating question

Solute balances convert solubility change into possible isolated crystal mass. Why this matters: A batch cooling profile produces a high apparent crystal yield but highly variable filtration time. Sudden nucleation may create fines even when the final equilibrium solubility predicts strong recovery. A controlled seed charge and slower cooling ramp can move crystal growth onto existing surfaces, producing a narrower size distribution and improving downstream separation. The trade-off is batch time and possible mother-liquor losses. Thermal history and impurity effects are as important as final solution temperature.

Physical interpretation

A cooling crystallizer may not reach equilibrium before filtration. The mother liquor also retained within a filter cake lowers apparent recovery and contaminates product. A valid material balance should account for evaporated solvent, dissolved impurities and possible polymorphic transformations. A large theoretical yield can be economically poor if submicron crystals require energy-intensive separation.

The principal mathematical relationship used in this lesson is

m_crystal=m_solvent(C_initial-C_final)

The equation is a compact model rather than an instruction to ignore transport resistance, property dependence or other equipment conditions. Before using it, translate every mathematical term into a measurable quantity. Identify whether the driving force is a local value at an interface, an average over equipment, or a measured inlet-to-outlet difference. Units and operating conditions must accompany any quoted property coefficient.

Derivation strategy and limitations

Start from the model: m_crystal=m_solvent(C_initial-C_final). Estimate uncertainty and parameter sensitivity before investing in a more elaborate model. Distinguish uncertainty in measured properties from uncertainty in model structure.

The physical assumptions matter more than a successful algebraic manipulation. If the property varies with temperature or composition, calculate locally or integrate over the changing state rather than using a single value without justification. For time-dependent problems, distinguish a true transient solution from a sequence of pseudo-steady snapshots. For multiphase systems, check that a concentration or composition basis has not silently changed from mole to mass fraction.

The calculation treats the crystallizing species as the only solute and assumes no solvent change.

Worked investigation

Charge a crystallizer with 800 kg solvent and 190 g solute per kg solvent at a warm condition. Cool to a condition where equilibrium solubility is 118 g/kg. Estimate equilibrium mass of dry crystal available from the solubility difference. If only 85% of this mass is recovered and the product filter cake retains 20 kg solution, quantify separate crystal recovery and mother-liquor contamination terms.

Solution approach: First list the information provided, and convert all values to consistent units. Substitute only the values that are supported by the statement; identify any additional properties or conditions that would need measurement for a full numerical result. Use the central relation m_crystal=m_solvent(C_initial-C_final) as the mathematical starting point, then evaluate a limiting or reference case. The final conclusion should distinguish an explicitly calculated response from a hypothetical plant-level extrapolation. Do not invent the unspecified flow, area, equilibrium or mass-transfer parameters to make a single definitive answer appear possible.

Engineering decision and diagnostic testing

At the reference operating point, perturb two high-impact inputs separately by about ten percent and record the response; interpret direction and magnitude. At the reference operating point, perturb two high-impact inputs separately by about ten percent and record the response; interpret direction and magnitude.

A meaningful experiment uses paired samples or simultaneous inlet/outlet measurements at a steady operating condition. Record the sampling interval and the instrument calibration or reference standard. If the observed response changes as throughput rises, test whether the change originates in the material property, the boundary layer, equilibrium constraints or residence-time distribution. When comparing equipment, express the outcome as recovery, product purity, energy or material consumption on the same basis.

Failure mechanisms and how to investigate them

• Inconsistent basis: Recalculate using one mass or molar basis and check whether every outflow is physically nonnegative. In a composition balance, sum the species to unity after solving.
• A parameter outside its fitted range: Compare the new operating condition with the measurements that established the coefficient. An apparently precise result may be unreliable under a different feed chemistry.
• An unmodeled resistance: Determine whether external transfer, internal diffusion, hydraulic pressure losses or phase equilibrium changes could constrain the response.
• An ambiguous measurement: Confirm the sample location and timing. An apparent difference in purity may be caused by mixing product samples across a transition.
• A neglected waste or utility duty: Draw the full separation train and account for regenerant, wash solvent, backwash, concentrate or mechanical power as applicable.

Sensitivity study

Retain the reference case from the worked investigation, then change the principal driving force by +10% while holding the model's other independent parameters fixed. Repeat for the strongest resistance or equilibrium coefficient. If the response is strongly nonlinear, a percentage change in the input does not imply the same percentage change in product. Where the result is dimensionless, verify that it remains within physical bounds. Where the result is a cumulative amount, confirm the time and inventory basis.

Fully worked numerical demonstration

At the initial concentration, solute inventory is 800×0.190=152 kg. At the final equilibrium solubility, the liquid holds 800×0.118=94.4 kg. The theoretical dry-crystal mass is 57.6 kg. At 85% solids recovery, the recovered dry crystal is 48.96 kg. A filter cake retaining 20 kg of mother liquor carries approximately 20×0.118/(1+0.118)=2.11 kg dissolved solute if the given solubility is on a solute-per-solvent basis. This entrained dissolved solute must be distinguished from dry recovered crystal and may cause contamination on drying.

Questions for independent work

1. Explain the physical meaning and units of every symbol in m_crystal=m_solvent(C_initial-C_final) and identify one property that would need experimental calibration.
2. Solve the worked investigation with a justified perturbation to its largest uncertain input, then report the direction of the result and a dimensional check.
3. Describe a realistic counterexample where the relation would fail even though all arithmetic was correctly implemented.
4. Develop a plant-level material balance that includes both recovered product and rejected or retained material.
5. Compare the model's predicted effect of increasing the driving force with an actual equipment constraint and state which measurement would reveal the discrepancy.

Further study

Read the related unit chapter and complete the unit assessment. Use the relevant simulation in the [laboratory directory](../index.html#labs) to investigate a parameter change. Keep the model limitations visible when drawing an engineering conclusion.

Unit 5, Lesson 4 — Population balance and crystal-size distribution

Purpose and operating question

A crystal-size distribution links nucleation, growth, residence time and separation behavior. Why this matters: A batch cooling profile produces a high apparent crystal yield but highly variable filtration time. Sudden nucleation may create fines even when the final equilibrium solubility predicts strong recovery. A controlled seed charge and slower cooling ramp can move crystal growth onto existing surfaces, producing a narrower size distribution and improving downstream separation. The trade-off is batch time and possible mother-liquor losses. Thermal history and impurity effects are as important as final solution temperature.

Physical interpretation

In an ideal MSMPR crystallizer with size-independent growth, the steady-state number population decreases exponentially with size. The mean characteristic size scale is G tau, but that is not directly equal to the volume-weighted mean diameter. Crystal-size moments have different physical meanings: number, length, surface area and volume are weighted by powers of L. Measurements of chord length require care because they are not necessarily direct crystal-size distributions.

The principal mathematical relationship used in this lesson is

n(L)=n_0 exp[-L/(G tau)]

The equation is a compact model rather than an instruction to ignore transport resistance, property dependence or other equipment conditions. Before using it, translate every mathematical term into a measurable quantity. Identify whether the driving force is a local value at an interface, an average over equipment, or a measured inlet-to-outlet difference. Units and operating conditions must accompany any quoted property coefficient.

Derivation strategy and limitations

Start from the model: n(L)=n_0 exp[-L/(G tau)]. Associate every ideal parameter with a realistic location in a piece of equipment: a membrane surface, packed bed, crystallizer slurry or filter cake.

The physical assumptions matter more than a successful algebraic manipulation. If the property varies with temperature or composition, calculate locally or integrate over the changing state rather than using a single value without justification. For time-dependent problems, distinguish a true transient solution from a sequence of pseudo-steady snapshots. For multiphase systems, check that a concentration or composition basis has not silently changed from mole to mass fraction.

This ideal MSMPR population lacks aggregation, breakage, classification and shape dependence.

Worked investigation

Take G=0.08 mm/h and mean residence time tau=3 h. Compute the exponential size scale G tau and the ratio n(0.5 mm)/n(0.1 mm). Repeat at tau=5 h. Explain how agglomeration and growth-rate dispersion would change the simple exponential population shape.

Solution approach: First list the information provided, and convert all values to consistent units. Substitute only the values that are supported by the statement; identify any additional properties or conditions that would need measurement for a full numerical result. Use the central relation n(L)=n_0 exp[-L/(G tau)] as the mathematical starting point, then evaluate a limiting or reference case. The final conclusion should distinguish an explicitly calculated response from a hypothetical plant-level extrapolation. Do not invent the unspecified flow, area, equilibrium or mass-transfer parameters to make a single definitive answer appear possible.

Engineering decision and diagnostic testing

Write an equipment-level feed/product/residue table. Do not claim whole-plant recovery from a local equilibrium or point-transport relation. Write an equipment-level feed/product/residue table. Do not claim whole-plant recovery from a local equilibrium or point-transport relation.

A meaningful experiment uses paired samples or simultaneous inlet/outlet measurements at a steady operating condition. Record the sampling interval and the instrument calibration or reference standard. If the observed response changes as throughput rises, test whether the change originates in the material property, the boundary layer, equilibrium constraints or residence-time distribution. When comparing equipment, express the outcome as recovery, product purity, energy or material consumption on the same basis.

Failure mechanisms and how to investigate them

• Inconsistent basis: Recalculate using one mass or molar basis and check whether every outflow is physically nonnegative. In a composition balance, sum the species to unity after solving.
• A parameter outside its fitted range: Compare the new operating condition with the measurements that established the coefficient. An apparently precise result may be unreliable under a different feed chemistry.
• An unmodeled resistance: Determine whether external transfer, internal diffusion, hydraulic pressure losses or phase equilibrium changes could constrain the response.
• An ambiguous measurement: Confirm the sample location and timing. An apparent difference in purity may be caused by mixing product samples across a transition.
• A neglected waste or utility duty: Draw the full separation train and account for regenerant, wash solvent, backwash, concentrate or mechanical power as applicable.

Sensitivity study

Retain the reference case from the worked investigation, then change the principal driving force by +10% while holding the model's other independent parameters fixed. Repeat for the strongest resistance or equilibrium coefficient. If the response is strongly nonlinear, a percentage change in the input does not imply the same percentage change in product. Where the result is dimensionless, verify that it remains within physical bounds. Where the result is a cumulative amount, confirm the time and inventory basis.

Fully worked numerical demonstration

For an ideal MSMPR crystallizer with constant G=0.08 mm/h and residence time tau=3 h, the characteristic size scale is G tau=0.24 mm. The population ratio n(0.5 mm)/n(0.1 mm) is exp[−(0.5−0.1)/0.24]≈0.189. At tau=5 h the characteristic scale is 0.40 mm and the same ratio is exp(−1)≈0.368. Longer residence can shift the number-based size distribution toward larger particles in this special model, but it does not capture agglomeration, breakage, shape-dependent growth or size-dependent removal.

Questions for independent work

1. Explain the physical meaning and units of every symbol in n(L)=n_0 exp[-L/(G tau)] and identify one property that would need experimental calibration.
2. Solve the worked investigation with a justified perturbation to its largest uncertain input, then report the direction of the result and a dimensional check.
3. Describe a realistic counterexample where the relation would fail even though all arithmetic was correctly implemented.
4. Develop a plant-level material balance that includes both recovered product and rejected or retained material.
5. Compare the model's predicted effect of increasing the driving force with an actual equipment constraint and state which measurement would reveal the discrepancy.

Further study

Read the related unit chapter and complete the unit assessment. Use the relevant simulation in the [laboratory directory](../index.html#labs) to investigate a parameter change. Keep the model limitations visible when drawing an engineering conclusion.

Unit 5: engineering project — Crystallization and controlled particle formation

Design brief

A batch cooling profile produces a high apparent crystal yield but highly variable filtration time. Sudden nucleation may create fines even when the final equilibrium solubility predicts strong recovery. A controlled seed charge and slower cooling ramp can move crystal growth onto existing surfaces, producing a narrower size distribution and improving downstream separation. The trade-off is batch time and possible mother-liquor losses. Thermal history and impurity effects are as important as final solution temperature.

Deliverables and evidence

1. A stream table with all independent feeds, products, residues and recycle flows; define a common mass or molar basis.
2. A labeled process diagram naming the physical location of the driving force and the location of the main measurement.
3. A one-page model statement with every coefficient, unit, validity range and inferred rather than measured property.
4. A nominal calculation and at least four parameter-variation cases, including a counterexample outside the simplifying assumptions.
5. A safety and environmental note describing rejected material, regenerant, solvent or concentrate management.
6. A short conclusion that states what the simulation can support, what it cannot, and the next experimental validation needed.

Technical focus

• Supersaturation, solubility and metastable zones — employ S=C/C_sat(T) only within its assumptions. Linearized solubility curves are illustrations only; real systems require measured equilibrium data.
• Primary nucleation, secondary nucleation and growth — employ G=k_g(DeltaC)^g only within its assumptions. The example is for a fictitious compound with illustrative kinetics, not a general law for all crystals.
• Batch cooling and crystal mass yield — employ m_crystal=m_solvent(C_initial-C_final) only within its assumptions. The calculation treats the crystallizing species as the only solute and assumes no solvent change.
• Population balance and crystal-size distribution — employ n(L)=n_0 exp[-L/(G tau)] only within its assumptions. This ideal MSMPR population lacks aggregation, breakage, classification and shape dependence.

Assessment rubric

• Material and energy accounting: 25%.
• Correct governing equations and units: 25%.
• Experimental test strategy: 20%.
• Sensitivity, uncertainty and limits: 20%.
• Clear engineering communication: 10%.

Safety boundary: No classroom model in this volume substitutes for validated pressure-equipment, chemical-handling or environmental-compliance design.

Laboratory 09 — Linear-driving-force uptake

Unit: 3. Model boundaries: Constant fluid equilibrium concentration and zero initial loading.

Learning objective

Identify the physical driving force, calculate the principal response, and compare the effect of two independent inputs. Explain the difference between an equation that fits measurements and an equipment design that accounts for full mass and energy balances.

Adjustable variables

Equilibrium loading (mg/g), Kinetic constant (min⁻¹), Contact time (min).

Experiment instructions

1. Record the reference-case output, its units and the assumed basis; confirm that the plotted point matches the numerical readout.
2. Change Contact time (min) to its lower and upper tested values. Export the CSV response, inspect monotonicity and identify any warnings or nonphysical outcomes.
3. Change one of the other parameters while holding the x-axis variable fixed. Explain the response in terms of the governing model and its dominant resistance or equilibrium contribution.
4. Estimate a practical equipment limitation omitted by this model, and describe a measurement that could reveal it. Do not claim equipment safety, product specification or regulatory compliance from this simulator alone.

Interpretation and model limits

Constant fluid equilibrium concentration and zero initial loading. This is a transparent educational approximation with explicit parameter inputs. It is not a validated design program. The UI may omit temperature changes, feed composition dynamics, axial dispersion, fouling, energy duties or multicomponent coupling depending on the selected phenomenon. Use the experimental log to identify whether those omitted effects would change the conclusion.

Three questions

1. Which manipulated parameter most strongly changes the displayed response over the chosen operating range, and why?
2. Where can the response be nonphysical or the model cease to represent the real process? State a specific boundary.
3. How would you validate one point of the predicted curve using laboratory or vendor measurements?

Laboratory 10 — Fixed-bed breakthrough front

Unit: 3. Model boundaries: Logistic experimental front; dimensionless effluent fraction, no competing contaminants.

Learning objective

Identify the physical driving force, calculate the principal response, and compare the effect of two independent inputs. Explain the difference between an equation that fits measurements and an equipment design that accounts for full mass and energy balances.

Adjustable variables

50% breakthrough time (min), Front steepness (min⁻¹), Elapsed time (min).

Experiment instructions

1. Record the reference-case output, its units and the assumed basis; confirm that the plotted point matches the numerical readout.
2. Change Elapsed time (min) to its lower and upper tested values. Export the CSV response, inspect monotonicity and identify any warnings or nonphysical outcomes.
3. Change one of the other parameters while holding the x-axis variable fixed. Explain the response in terms of the governing model and its dominant resistance or equilibrium contribution.
4. Estimate a practical equipment limitation omitted by this model, and describe a measurement that could reveal it. Do not claim equipment safety, product specification or regulatory compliance from this simulator alone.

Interpretation and model limits

Logistic experimental front; dimensionless effluent fraction, no competing contaminants. This is a transparent educational approximation with explicit parameter inputs. It is not a validated design program. The UI may omit temperature changes, feed composition dynamics, axial dispersion, fouling, energy duties or multicomponent coupling depending on the selected phenomenon. Use the experimental log to identify whether those omitted effects would change the conclusion.

Three questions

1. Which manipulated parameter most strongly changes the displayed response over the chosen operating range, and why?
2. Where can the response be nonphysical or the model cease to represent the real process? State a specific boundary.
3. How would you validate one point of the predicted curve using laboratory or vendor measurements?

Unit 5: assessment and answer guidance

These questions require a declared basis, explicit assumptions and units; no safety-critical design decision is authorized by classroom examples.

Topic 1: Supersaturation, solubility and metastable zones

Q1.1 — Interpretation: Starting from S=C/C_sat(T), identify the independent driving force and the dependent response. Give their dimensions and explain which parameter must be measured.

Q1.2 — Calculation: A hypothetical solute has solubility 180 g/kg solvent at 50 C and 105 g/kg at 25 C. A batch containing 155 g/kg is cooled from 50 to 25 C. Calculate the initial and final S and explain when seeded crystallization becomes feasible. How does seed surface area influence induction time and crystal size?

Q1.3 — Model validation: Linearized solubility curves are illustrations only; real systems require measured equilibrium data. Explain an experimental plan to identify whether this limitation is important at the proposed operating point.

Q1.4 — Decision: Compare a higher-throughput operating condition with a lower-throughput but more selective or stable condition. Identify the missing measurements needed to rank them for a plant.

Q1.5 — Fault diagnosis: Suppose an observed result deviates by 25% from the relation. Distinguish data-unit mistakes, instrument bias, a boundary condition change and genuinely invalid model assumptions.

Answer-check framework: A correct solution cites S=C/C_sat(T), uses consistent concentration or mass units, explicitly honors the limit (“Linearized solubility curves are illustrations only; real systems require measured equilibrium data.”), and reports a physically sensible direction of change. Numerical calculations must show substitutions and a final unit. Alternative answers can be valid when a different, clearly documented operating assumption is made.

Topic 2: Primary nucleation, secondary nucleation and growth

Q2.1 — Interpretation: Starting from G=k_g(DeltaC)^g, identify the independent driving force and the dependent response. Give their dimensions and explain which parameter must be measured.

Q2.2 — Calculation: Assume growth velocity is 0.06 mm/h at a supersaturation driving force of 10 g/kg, with first-order dependence. Estimate growth velocity at 6 and 15 g/kg. For a constant 0.06 mm/h and six-hour batch, calculate ideal change in crystal length. Discuss why a real crystal population does not all grow by the same amount.

Q2.3 — Model validation: The example is for a fictitious compound with illustrative kinetics, not a general law for all crystals. Explain an experimental plan to identify whether this limitation is important at the proposed operating point.

Q2.4 — Decision: Compare a higher-throughput operating condition with a lower-throughput but more selective or stable condition. Identify the missing measurements needed to rank them for a plant.

Q2.5 — Fault diagnosis: Suppose an observed result deviates by 25% from the relation. Distinguish data-unit mistakes, instrument bias, a boundary condition change and genuinely invalid model assumptions.

Answer-check framework: A correct solution cites G=k_g(DeltaC)^g, uses consistent concentration or mass units, explicitly honors the limit (“The example is for a fictitious compound with illustrative kinetics, not a general law for all crystals.”), and reports a physically sensible direction of change. Numerical calculations must show substitutions and a final unit. Alternative answers can be valid when a different, clearly documented operating assumption is made.

Topic 3: Batch cooling and crystal mass yield

Q3.1 — Interpretation: Starting from m_crystal=m_solvent(C_initial-C_final), identify the independent driving force and the dependent response. Give their dimensions and explain which parameter must be measured.

Q3.2 — Calculation: Charge a crystallizer with 800 kg solvent and 190 g solute per kg solvent at a warm condition. Cool to a condition where equilibrium solubility is 118 g/kg. Estimate equilibrium mass of dry crystal available from the solubility difference. If only 85% of this mass is recovered and the product filter cake retains 20 kg solution, quantify separate crystal recovery and mother-liquor contamination terms.

Q3.3 — Model validation: The calculation treats the crystallizing species as the only solute and assumes no solvent change. Explain an experimental plan to identify whether this limitation is important at the proposed operating point.

Q3.4 — Decision: Compare a higher-throughput operating condition with a lower-throughput but more selective or stable condition. Identify the missing measurements needed to rank them for a plant.

Q3.5 — Fault diagnosis: Suppose an observed result deviates by 25% from the relation. Distinguish data-unit mistakes, instrument bias, a boundary condition change and genuinely invalid model assumptions.

Answer-check framework: A correct solution cites m_crystal=m_solvent(C_initial-C_final), uses consistent concentration or mass units, explicitly honors the limit (“The calculation treats the crystallizing species as the only solute and assumes no solvent change.”), and reports a physically sensible direction of change. Numerical calculations must show substitutions and a final unit. Alternative answers can be valid when a different, clearly documented operating assumption is made.

Topic 4: Population balance and crystal-size distribution

Q4.1 — Interpretation: Starting from n(L)=n_0 exp[-L/(G tau)], identify the independent driving force and the dependent response. Give their dimensions and explain which parameter must be measured.

Q4.2 — Calculation: Take G=0.08 mm/h and mean residence time tau=3 h. Compute the exponential size scale G tau and the ratio n(0.5 mm)/n(0.1 mm). Repeat at tau=5 h. Explain how agglomeration and growth-rate dispersion would change the simple exponential population shape.

Q4.3 — Model validation: This ideal MSMPR population lacks aggregation, breakage, classification and shape dependence. Explain an experimental plan to identify whether this limitation is important at the proposed operating point.

Q4.4 — Decision: Compare a higher-throughput operating condition with a lower-throughput but more selective or stable condition. Identify the missing measurements needed to rank them for a plant.

Q4.5 — Fault diagnosis: Suppose an observed result deviates by 25% from the relation. Distinguish data-unit mistakes, instrument bias, a boundary condition change and genuinely invalid model assumptions.

Answer-check framework: A correct solution cites n(L)=n_0 exp[-L/(G tau)], uses consistent concentration or mass units, explicitly honors the limit (“This ideal MSMPR population lacks aggregation, breakage, classification and shape dependence.”), and reports a physically sensible direction of change. Numerical calculations must show substitutions and a final unit. Alternative answers can be valid when a different, clearly documented operating assumption is made.