Introduction & Context
Recrystallization, specifically through the mechanism of Ostwald ripening, is a critical phenomenon in process engineering involving crystal suspensions. It describes the process where larger crystals grow at the expense of smaller ones due to differences in surface energy and solubility. In industrial applications—such as the storage of ice cream, pharmaceutical suspensions, or chemical slurries—temperature fluctuations act as a catalyst for this process. As temperature cycles, the solubility of the solute changes, driving the dissolution of smaller particles and the subsequent deposition of mass onto larger particles. For a deeper quantitative insight, see our Ostwald ripening rate estimation methodology, which helps predict how quickly these size changes occur and supports maintaining product quality, texture, and stability over the shelf life of a product.
Methodology & Formulas
The evolution of the crystal size distribution is modeled using the Lifshitz-Slyozov-Wagner (LSW) theory. This approach assumes that the process is diffusion-limited and that the system maintains a self-similar size distribution over time. The calculation tracks the third moment of the crystal size distribution, which is proportional to the cube of the mean diameter.
The rate of change for the mean cube diameter is defined by the following expression:
\[ \frac{d\langle L^3 \rangle}{dt} = \frac{64 \cdot \sigma \cdot D \cdot C_{\mathrm{eq}} \cdot V_{\mathrm{m}}^2}{9 \cdot n \cdot R \cdot T} \]
To determine the final mean crystal size after a specific duration, the equation is integrated over time:
\[ \langle L^3 \rangle_{\mathrm{final}} = \langle L^3 \rangle_{\mathrm{initial}} + \left( \frac{64 \cdot \sigma \cdot D \cdot C_{\mathrm{eq}} \cdot V_{\mathrm{m}}^2}{9 \cdot n \cdot R \cdot T} \right) \cdot \Delta t \]
The final mean diameter is then derived by taking the cube root of the resulting third moment:
\[ L_{\mathrm{final}} = \sqrt[3]{\langle L^3 \rangle_{\mathrm{final}}} \]
| Parameter |
Condition/Threshold |
Implication |
| Solid Volume Fraction (φ) |
φ < 0.1 |
LSW theory is strictly valid; higher values require empirical correction. |
| Crystal Size (L) |
L > 1 nm |
Continuum thermodynamics apply; below this, molecular effects dominate. |
| Temperature (T) |
T > 0 K |
Absolute temperature must be positive for physical validity. |
Recrystallization occurs when thermal cycling causes the dissolution and subsequent growth of crystals. When the temperature fluctuates, the solubility of the solute changes, leading to the following phenomena:
- Small crystals dissolve during warmer phases due to higher solubility.
- Solute molecules migrate toward larger crystals during cooler phases.
- The system minimizes surface energy, resulting in the growth of larger particles at the expense of smaller ones, often referred to as Ostwald ripening.
Worked Example: Ice Crystal Growth Due to Temperature Fluctuation in Ice Cream
Scenario: A batch of ice cream (solid volume fraction 0.05) is stored in a freezer with daily temperature fluctuations. The average storage temperature is –18 °C (255 K). The initial mean ice crystal diameter is 30.0 μm. Below are the system parameters.
- Initial mean crystal diameter: L0 = 3.000×10–5 m
- Average storage temperature: Tavg = 255.000 K
- Storage duration: t = 2.592×106 s (30 days)
- Interfacial tension: σ = 0.025 J/m²
- Molar volume of solid: Vm = 1.800×10–5 m³/mol
- Diffusivity of solute in liquid: D = 1.000×10–12 m²/s
- Equilibrium concentration (flat surface): Ceq = 1.000×103 mol/m³
- Number of ions per molecule: n = 1.0
- Gas constant: R = 8.314 J/(mol·K)
- Convert initial crystal diameter to meters.
L0 = 30.0 μm = 3.000×10–5 m.
- Compute the LSW ripening rate constant K.
\[
K = \frac{64 \sigma D C_{\mathrm{eq}} V_{\mathrm{m}}^2}{9 n R T_{\mathrm{avg}}}
\]
Substituting the known values yields K = 2.717×10–23 m³/s (from numerical evaluation).
- Calculate the change in the cube of the mean diameter over the storage period.
\[
\Delta\langle L^3 \rangle = K \cdot t = (2.717\times10^{-23} \text{ m}^3/\text{s}) \times (2.592\times10^6 \text{ s}) = 7.042\times10^{-17} \text{ m}^3.
\]
- Determine the initial cube of the mean diameter.
\[
L_0^3 = (3.000\times10^{-5} \text{ m})^3 = 2.700\times10^{-14} \text{ m}^3.
\]
- Add the increment to obtain the final cube of the mean diameter.
\[
\langle L^3 \rangle_{\mathrm{final}} = L_0^3 + \Delta\langle L^3 \rangle = 2.700\times10^{-14} + 7.042\times10^{-17} = 2.707\times10^{-14} \text{ m}^3.
\]
- Take the cube root to find the final mean diameter.
\[
L_{\mathrm{final}} = \left( \langle L^3 \rangle_{\mathrm{final}} \right)^{1/3} = \left( 2.707\times10^{-14} \text{ m}^3 \right)^{1/3} = 3.003\times10^{-5} \text{ m}.
\]
- Convert to micrometers.
Lfinal = 3.003×10–5 m × 106 μm/m = 30.03 μm.
Final Answer: The mean ice crystal diameter after 30 days of storage under temperature fluctuations is Lfinal = 30.03 μm (to four significant figures). This represents a small increase from the initial 30.00 μm, consistent with the low solid volume fraction assumed in this example; real ice cream (solid fraction 0.3–0.5) would exhibit larger ripening due to enhanced diffusion-limited growth.
