This manuscript is devoted to the nonlinear dynamics of particulate assemblages in metastable liquids, caused by various dynamical laws of crystal growth and nucleation kinetics. First of all, we compare the quasi-steady-state and unsteady-state growth rates of spherical crystals in supercooled and supersaturated liquids. It is demonstrated that the unsteady-state rates transform to the steady-state ones in a limiting case of fine particles. We show that the real crystals evolve slowly in a more actual case of unsteady-state growth laws. Various growth rates of particles are tested against experimental data in metastable liquids. It is demonstrated that the unsteady-state rates describe the nonlinear behaviour of experimental curves with increasing the growth time or supersaturation. Taking this into account, the crystal-size distribution function and metastability degree are analytically found and compared with experimental data on crystallization in inorganic and organic solutions. It is significant that the distribution function is shifted to smaller sizes of particles if we are dealing with the unsteady-state growth rates. In addition, a complete analytical solution constructed in a parametric form is simplified in the case of small fluctuations in particle growth rates. In this case, a desupercooling/desupersaturation law is derived in an explicit form. Special attention is devoted to the biomedical applications for insulin and protein crystallization. This article is part of the theme issue ‘Patterns in soft and biological matters’.

The transition of a metastable liquid (supersaturated solution or supercooled melt) occurring from the intermediate stage (where the crystals nucleate and grow) to the concluding stage (where the larger particles evolve at the expense of the dissolution of smaller particles) is theoretically described, with allowance for various mass transfer mechanisms (reaction on the interface surface, volume diffusion, grain-boundary diffusion, diffusion along the dislocations) arising at the stage of Ostwald ripening (coalescence). The initial distribution function (its ‘tail’) for the concluding stage (forming as a result of the evolution of a particulate assemblage during the intermediate stage) is taken into account to determine the particle-size distribution function at the stage of Ostwald ripening. This modified distribution function essentially differs from the universal Lifshitz–Slyozov (LS) solutions for several mass transfer mechanisms. Namely, its maximum lies below and is shifted to the left in comparison with the LS asymptotic distribution function. In addition, the right branch of the particle-size distribution lies above and is shifted to the right of the LS blocking point. It is shown that the initial ‘tail’ of the particle-size distribution function completely determines its behaviour at the concluding stage of Ostwald ripening. The present theory agrees well with experimental data. This article is part of the theme issue ‘Patterns in soft and biological matters’.

The Ostwald ripening stage of a phase transformation process with allowance for synchronous operation of various mass transfer mechanisms (volume diffusion and diffusion along the block boundaries and dislocations) and the initial condition for the particle-radius distribution function is theoretically studied. The initial condition is taken from the analytical solution describing the intermediate stage of a phase transition process. The present theory focuses on relaxation dynamics from the beginning of the ripening process to its final asymptotic state, which is described by the previously constructed theories (Slezov VV. et al. 1978 J. Phys. Chem. Solids 39 , 705–709. ( doi:10.1016/0022-3697(78)90002-1 ) and Alexandrov & Alexandrova 2020 Phil. Trans. R. Soc. A 378 , 20190247. ( doi:10.1098/rsta.2019.0247 )). An evolutionary behaviour of particle growth rates dependent on various mass transfer mechanisms and time is analytically described. The boundaries of the transition layer, which surround the blocking point, are found. The fundamental and relaxation contributions to the particle-radius distribution function are derived for the simultaneous occurrence of various mass transfer mechanisms. The left branch of this function is shifted to smaller particle radii whereas its right branch extends to the right of the blocking point as compared with the asymptotic universal distribution function. The theory under consideration well agrees with experimental data. This article is part of the theme issue ‘Transport phenomena in complex systems (part 1)’.

Motivated by a wide range of applications in various fields of physics and materials science, we consider a generalized approach to the evolution of a polydisperse ensemble of spherical particles in metastable media. An integrodifferential system of governing equations, consisting of a kinetic equation for the particle-size distribution function (Fokker–Planck type equation) and a balance equation for the temperature (concentration) of a metastable medium, is formulated. The kinetic equation takes into account fluctuations in the growth/reduction rates of individual particles, the velocity of particles in a spatial direction, the withdrawal of particles of a given size from the metastable medium, and their source/sink term. The heat (mass) balance equation takes into account the growth/reduction of particles in a metastable system as well as heat (mass) exchange with the environment. A generalized system of equations describes various physical and chemical processes of phase transformations, such as the growth and dissolution of crystals, the evaporation of droplets, the boiling of liquids and the combustion of a polydisperse fuel. The ways of analytical solution of the formulated integrodifferential system of equations based on the saddle-point technique and the separation of variables method are considered. The theory can be applied when describing the evolution of an ensemble of particles at the initial and intermediate stages of phase transformation when the distances between the particles are large enough, and interactions between them can be neglected.

Evolution of crystal ensembles in supersaturated solutions is studied at the initial and intermediate stages of bulk crystallization. An integro-differential model includes fluctuations in crystal growth rates, initial crystal-size distribution and arbitrary nucleation and growth kinetics of crystals. Two methods based on variables separation and saddle-point technique for constructing a complete analytical solution to this model are considered. Exact parametric solutions based on these methods are derived. Desupersaturation dynamics is in good agreement with the experimental data for bovine and porcine insulin. The method based on variables separation has a strong physical limitation on exponentially decaying initial distribution and leads to the distribution function increasing with time. The method based on saddle-point technique leads to a dome-shaped crystal-size distribution function decreasing with time and has no strong physical limitations. The latter circumstance makes this method more reasonable for describing the kinetics of bulk crystallization in solutions and melts.