Solidification of colloidal suspensions

We present a mathematical model of the unidirectional solidification of a suspension of hard-sphere colloids. Similarity solutions are obtained for the volume fraction and temperature profiles ahead of a planar solidification front. The highly nonlinear functional dependence of the diffusion coefficient on the volume fraction gives rise to a range of behaviours. For small particles, Brownian diffusion dominates and the system behaviour is reminiscent of binary-alloy solidification. Constitutional supercooling occurs at the interface under certain conditions, leading potentially to an instability in the shape of the interface. For larger particles, Brownian diffusion is weak and the particles form a porous layer above the interface. In this case constitutional supercooling reaches a maximum near the surface of the layer, and the porous medium itself is potentially unstable. In stable systems there exists the possibility of secondary nucleation of ice.