Next Article in Journal
Annealing Engineering in the Growth of Perovskite Grains
Next Article in Special Issue
Evolution of a Polydisperse Ensemble of Spherical Particles in a Metastable Medium with Allowance for Heat and Mass Exchange with the Environment
Previous Article in Journal
Correction: Crenn et al. Influence of Anodized Titanium Surfaces on the Behavior of Gingival Cells in Contact with: A Systematic Review of In Vitro Studies. Crystals 2021, 11, 1566
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Nucleation and Growth of an Ensemble of Crystals during the Intermediate Stage of a Phase Transition in Metastable Liquids

by
Liubov V. Toropova
1,2,*,
Eugenya V. Makoveeva
3,
Sergei I. Osipov
4,
Alexey P. Malygin
4,
Yang Yang
5 and
Dmitri V. Alexandrov
3
1
Laboratory of Mathematical Modeling of Physical and Chemical Processes in Multiphase Media, Department of Theoretical and Mathematical Physics, Ural Federal University, Lenin Ave., 51, 620000 Ekaterinburg, Russia
2
Otto-Schott-Institut für Materialforschung, Friedrich-Schiller-Universität-Jena, 07743 Jena, Germany
3
Laboratory of Multi-Scale Mathematical Modeling, Department of Theoretical and Mathematical Physics, Ural Federal University, Lenin Ave., 51, 620000 Ekaterinburg, Russia
4
Institute of Natural Sciences and Mathematics, Ural Federal University, Lenin Ave., 51, 620000 Ekaterinburg, Russia
5
Key Laboratory of Polar Materials and Devices and Physics Department, East China Normal University, Shanghai 200241, China
*
Author to whom correspondence should be addressed.
Crystals 2022, 12(7), 895; https://doi.org/10.3390/cryst12070895
Submission received: 6 June 2022 / Revised: 21 June 2022 / Accepted: 22 June 2022 / Published: 24 June 2022
(This article belongs to the Special Issue Phase Transition in External Fields)

Abstract

:
In this paper, an analytical method of solving the integro-differential system of kinetic and balance equations describing the evolution of an ensemble of crystals during the intermediate phase of the bulk crystallization process is described. The theory is developed for kinetic equations of the first- and second order corresponding to the absence and presence of fluctuations in particle growth rates. The crystal-size distribution function as well as the dynamics of metastability reduction in a supercooled melt (supersaturated solution) are analytically found using the saddle-point and the Laplace transform methods. The theory enables us to obtain the crystal-size distribution function that establishes in a supercooled (supersaturated) liquid at the beginning of the final stage of a phase transformation process when Ostwald ripening, coagulation and fragmentation of crystals are able to occur.

1. Introduction

The phase transformations from a metastable to a thermodynamically stable state are accompanied by the growth of new phase particles on microscopic inclusions or impurities. In the early stages of a phase transformation, the evolution of nuclei can be considered independently of each other, whereas in the later stages of the process it is necessary to take into account the interaction between growing crystallites, which affects the degree of liquid metastability (its supercooling or supersaturation) [1,2,3,4,5]. The dynamic process of evolution of a metastable zone filled with the crystals at all stages of the phase transition is important from the theoretical and applied points of view. In particular, this applies to crystallization from supercooled melts or supersaturated solutions, where the phase transition is maintained by fluctuations in the growth rates of particles, and the contribution from impurity (representing nucleation sites) is rather small [6,7]. Note that crystal nucleation and growth processes control the size distributions obtained in laboratory and industrial crystallizers [8,9,10,11,12,13]. As this occurs, the growth rate of solid phase particles also plays a key role in the evolution of the metastable system and, thus, in the physico-chemical properties of the produced materials [14,15].
Generally speaking, the process of bulk phase transformation from a supercooled melt or supersaturated solution may be described by four main stages [16,17]. The supercooled (supersaturated) state is established in the preliminary stage. Then, in the initial stage, the formation of new phase nuclei occurs (in a heterogeneous or homogeneous manner). Next, nucleation, crystal growth, and partial supercooling (or desupersaturation) occur in the intermediate stage. A final stage describes the relaxation processes of Ostwald ripening, coagulation, and particle fragmentation [18]. The study of the initial stages of the phase transition is generally based on kinetic nucleation theory and on the assumption that the growth of individual particles is independent [19,20,21]. The final stage of a phase transformation is analyzed according to the Lifshitz-Slyozov theory, which is based on the assumptions that no particles of a new phase appear and that large crystals grow at the expense of small ones [22,23,24,25,26,27,28].
The present paper is focused on the study of the intermediate stage when the processes of nucleation of the new phase elements and the growth of already existing crystals are of same importance. The general theory of this stage is incomplete, however the intermediate stage often takes most of the time of the whole phase transition process. This is caused, in particular, since the corresponding mathematical model of the process is integro-differential, and some of the boundary conditions are put on the moving growth boundaries. Note that complete information on the dynamics of metastability removal and evolution of the new phase may be derived by analyzing just this stage. The absence of universal methods for solving the integrodifferential equations is the reason why the dynamics of desupercooling (or desupersaturation) is completely ignored in several evolutionary models [29,30]. Note that the necessity of simultaneous accounting for the processes of nucleation and growth of a new phase is described, for example, in Refs. [31,32]. Nevertheless, many models are based on stationary approximations (see, e.g., [33,34]), which may occur only at the very initial stages or in the case of particular system parameters.
The earliest attempt to overcome the above described difficulties of formulating and solving a mathematical model was undertaken in Ref. [35] using the saddle-point method of calculating a Laplace-type integral. A similar approach was then used to solve problems about the evolution of aggregates in magnetic fluids and metastable colloids [36,37]. The approach and solutions [35] were then applied to describe protein crystallization [38] and to verify von Weimarn theory describing the average size of crystals [39]. It is important to note that the analytical solutions were obtained only for the main contribution to the Laplace integral [35]. This contribution leads to the zero-order approximation of the Laplace integral in the neighborhood of the saddle point. In this case, the dimensionless Gibbs number p included in the nucleation frequency I = I * exp p / υ 2 formally tends to infinity, where υ and I * denote the dimensionless supercooling/supersaturation and pre-exponential factor. Since p varies from ∼ 10 1 to ∼ 10 3 for real supercooled melts or supersaturated solutions, the theory of Ref. [35] must be developed to describe such systems. It is also important to note that even in the case of sufficiently large p 10 2 10 3 the zero-order solution of Ref. [35] differs significantly from the complete solution containing the following terms. Here we develop the nucleation and growth theory at the intermediate stage with allowance for additional terms to the aforementioned zero-order solution and possible fluctuations in crystal growth rates.

2. The Model of Nucleation and Crystal Growth

The molecular-kinetic theory is based on the assumption that micronuclei of a new phase may form as a result of density fluctuations in metastable melts or solutions [19,20,40,41,42]. If the characteristic size of such nuclei exceeds some critical value r * , they are able to grow further and it will take work to form a new spherical particle of radius r in a one-component system
Y ( r ) = 4 π r 2 γ i 4 3 π r 3 ρ s μ l μ s ,
where ρ s is the solid phase density, γ i is the surface tension, μ s and μ l are the chemical potentials of the solid and liquid phases, respectively, and μ l > μ s for a metastable melt. However, the critical radius r * of a particle can be achieved when Y ( r ) reaches a maximum. In this case, particals with radii r = r * are unstable, and the work Y * of their formation can be written as
Y * = 16 π γ i 3 3 ρ s 2 ( μ l μ s ) 2 .
The chemical potential difference may be calculated using the Gibbs-Helmholtz equation and has a form [43]
μ l μ s = L Δ T T p , Δ T = T p T l ,
where L represents the latent heat of the phase transition, and T p and T l are the phase transition and melt temperatures, respectively. Nucleation of a new phase occurs if the energy barrier Y * that prevents nucleation is overcome. In this case, the nucleation frequency can be represented as an exponential function of the energy barrier height [16,42,44]. Thus, we obtain
I = I * exp Y * k B T l = I * exp 16 π γ i 3 T p 2 3 ρ s 2 L 2 Δ T 2 k B T l ,
where k B stands for the Boltzmann constant. The pre-exponential factor I * weakly depends on supercooling/supersaturation (metastability degree) [20,40,41,45] and can be a function of r [33]. However, in the present paper I * is assumed to be a constant, because the nucleation frequency is included only in the boundary condition at r = r * . Considering Δ T T p , expression (4) can be rewritten as
I = I * exp 16 π γ i 3 T p 3 ρ s 2 L 2 Δ T 2 k B .
The model is applicable for single-component melts and should be modified for binary systems. In that case, the phase transition temperature T p must be replaced by the function T p ( C l ) , as it depends on the impurity concentration C l . Equation (5) can be rewritten as
I = I * exp p ( Δ T / Δ T 0 ) 2 , p = 16 π γ i 3 T p 3 ρ s 2 L 2 Δ T 0 2 k B ,
where Δ T 0 is the characteristic supercooling. Expression (6) defines the nucleation frequency as a function of dimensionless supercooling υ = Δ T Δ T 0 , and in the case of supersaturated solutions has a form [46]
I = I * exp p ln 2 ( C l / C p ) , p = 16 π γ i 3 M s 2 3 ρ s 2 R g 2 T s 3 k B .
Here C p is the supersaturation concentration, M s is the molecular weight, R g is the universal gas constant, and T s is the solution temperature.
Further, we will use a quasi-stationary approximation to solve a Stefan-type problem with a moving phase transition boundary, which occurs as a result of nucleus growth. The temperature/concentration field T around the growing spherical particle and a growth rate d r d τ take a form
2 T = 0 , ρ > r ( τ ) , d r d τ = λ l ρ s L T ρ = β * ( T p T ) , ρ = r ( τ ) , ρ / r ( τ ) 1 , T T l ,
where ρ represents the spherical coordinate, λ l is the thermal conductivity coefficient of the liquid phase, β * is the kinematic parameter, and τ stands for time. The solution of Equation (8) defines the particle growth rate as
d r d τ = β * Δ T 1 + β * q r , q = ρ s L λ l .
Expression (9) represents the growth rate as a linear function of supercooling that has two regimes. When the crystal size is small, r β * q 1 , growth regime is called kinetic, since it is entirely determined by surface processes. The other regime occurs when the crystals significantly exceed the value of β * q 1 , and the growth rate is controlled by heat dissipation velocity. In the case of liquid solutions, Equation (9) can be rewritten as [47]
d r d τ = β * Δ C 1 + β * q C r , q C = C p ( k 0 1 ) D ,
where k 0 is the partition coefficient, D and Δ C represent the diffusion coefficient and the liquid supersaturation, respectively.
Next, consider a system of spherical solid crystals in a macroscopically homogeneous one-component supercooled (supersaturated) system. The evolution of the system can be described by the Fokker-Planck type kinetic equation. Neglecting random fluctuations in the crystal growth rate and considering a temperature distribution in the system as homogeneous, we get
ϕ τ + r d r d τ ϕ = 0 , r > r * ,
ρ m C m d Δ T d τ = 4 π ρ s L r * r 2 d r d τ ϕ d r ,
where ϕ is a number size density, C m and ρ m are heat capacity and constant density of a melt, respectively. We also suppose that the liquid instantaneously cools down below the crystallization temperature by Δ T 0 at the initial moment. The boundary conditions for Equations (11) and (12) can be written as
τ = 0 , Δ T = Δ T 0 , ϕ = 0 ,
r = r * , d r d τ ϕ = I ( Δ T ) .
Expressions (4)–(7) and (9)–(14) represent a complete system of equations describing the intermediate stage of bulk phase transitions.

3. Kinetics of Supercooling (Supersaturation) Removal

Let us introduce the following dimensionless parameters and variables
ζ = r l 0 , t = τ τ 0 , υ = Δ T Δ T 0 , Φ = l 0 4 ϕ ,
τ 0 = β * 3 Δ T 0 3 I 0 1 / 4 , l 0 = β * Δ T 0 1 / 4 I 0 1 / 4 , α * = β * ρ s L l 0 λ l , B 1 = 4 π ρ s L ρ m C m Δ T 0 ,
where I 0 = I ( Δ T 0 ) . Rewriting the statement of the problem ( r * 0 ) in dimensionless variables (16), we arrive at
Φ t + υ ζ Φ 1 + α * ζ = 0 , ζ > 0 ,
d υ d t = B 1 υ 0 ζ 2 Φ ( t , ζ ) 1 + α * ζ d ζ ,
Φ = 0 , t = 0 , υ = 1 ,
Φ = 1 υ exp p g ( υ ) , ζ = 0 .
Note that the dimensionless model (17)–(20) defines an intermediate stage of phase transitions both for supercooled melts and supersaturated solutions. The function g ( υ ) determines the kinetic law of the process in Equation (20) that can be used for both homogeneous and heterogeneous nucleation [48].
In the case of supercooled melt (lower index “sm”) and Weber-Volmer-Frenkel-Zel’dovich kinetics (WVFZ), we obtain the function g ( υ ) from Equation (6) as
g ( υ ) = g s m = 1 υ 2 .
Using the empirical law I = I * ( Δ T ) p for the nucleation frequency, (see [49]), we have
g ( υ ) = g s m = ln υ .
In the case a supersaturated solution (lower index “ss”), υ represents dimensionless supersaturation, and dimensionless parameters can be defined as
υ = Δ C Δ C 0 , α * = β * C p l 0 D , B 1 = 4 π C p Δ C 0 , τ 0 = β * 3 Δ C 0 3 I 0 1 / 4 , l 0 = β * Δ C 0 1 / 4 I 0 1 / 4 ,
where Δ C 0 is the initial supersaturation.
Further, for a supersaturated solution with WVFZ kinetics, we obtain from expression (7)
g ( υ ) = g s s = ln 2 1 + 1 υ p ln 2 1 + υ υ p ,
where υ p = C p Δ C 0 .
In the case of a supersaturated solution with Meirs kinetics I = I * ( Δ C ) p [44], we have
g ( υ ) = g s s = ln υ .
The kinetics of metastability removal is described by the nonlinear integro-differential Equations (17) and (18). The analytical solution of the model is obtained from the saddle-point method for calculating a Laplace-type integral [50,51]. The solution to Equation (17) with boundary and initial conditions (19) and (20) is given by
Φ ( t , ζ ) = 1 + α * ζ φ ξ ( t ) η ( ζ ) H e a v ξ ( t ) η ( ζ ) ,
where
φ ( ξ ( t ) ) = 1 υ ( t ) exp p g ( υ ( t ) ) , ξ ( t ) = 0 t υ ( t ) d t , η ( ζ ) = 0 ζ ( 1 + α * ζ ) d ζ = ζ + α * ζ 2 2 ,
and H e a v ( ξ ) denotes the Heaviside function.
Integration of the dimensionless growth rate from Equations (9) and (10) using variables (16) and (23) we arrive at
d ζ d t = υ ( t ) 1 + α * ζ ,
and with the initial condition ζ = 0 at t = ν , we obtain
α * ζ = 1 + 2 α * ξ ( t ) ξ ( ν ) 1 .
Expression (27) determines the crystal radius ζ m ( t ) in a dimensionless form (at ν = 0 ) appearing as a function of time t.
Now we replace the variable ζ at any constant time t by the new variable ν as
ξ ( ν ) = ξ ( t ) η ( ζ ) .
Note that the limits ζ = ζ m ( t ) and ζ = 0 of the integration over ζ correspond to the limits ν = t and ν = 0 of the integration over ν , respectively. Hence, by substituting Equations (26) and (27) into expression (7), we arrive at
d υ d t = B 1 υ 0 t h ( ν , t ) exp p g ( ν ) d ν ,
where
h ( ν , t ) = 1 + 2 α * ξ ( t ) ξ ( ν ) 1 2 α * 2 1 + 2 α * ξ ( t ) ξ ( ν ) ,
and g ( ν ) = g ( υ ( ν ) ) is defined by Equations (21), (22), (24) or (25).
The integro-differential Equation (28) can be simplified if the dimensionless Gibbs number p is higher than unity and the nucleation is sufficiently intense. In this case, the Laplace-type integral (28) can be calculated using the saddle-point method.
Equations (21), (22), (24) and (25) introduce d g d ν < 0 for all kinetic mechanisms. This represents that the function g ( ν ) reaches the maximum at the boundary ν = 0 . Using the Equation (28) and calculating the derivatives of the function υ over ν , we arrive at first three derivatives that are equal to zero at ν = 0 , and g ( n ) | ν = 0 0 at n 4 .
The first four derivatives are given in Table 1, where
χ = 6 B 1 1 / 4 , κ = 1 1 + υ p ln 3 1 + υ p 1 .
Thus, the saddle point method for an integral of Laplace type (28) gives the solution to the problem as follows [50]
0 t h ( ν , t ) exp p g ( ν ) d ν k = 0 p ( k + 1 ) / 4 a k ( ξ ( t ) ) ,
a k ( ξ ( t ) ) = ( 1 ) k + 1 4 k k ! Γ k + 1 4 H k ( 0 ) d k d ν k h ( ν , t ) H ( ν ) ν = 0 ,
where Γ represents the Euler gamma function, and
H ( ν ) = g ( ν ) 3 / 4 g ( ν ) .
The first four coefficients in expression (31) have the form ( z = z ( ξ ( t ) ) )
a 0 ( z ) = Γ 1 4 h ( 0 , t ) H ( 0 ) , a 1 ( z ) = 4 Γ 1 2 H ( 0 ) H ( 0 ) h ν ν = 0 + h ( 0 , t ) H ( 0 ) ,
a 2 ( z ) = 8 Γ 3 4 H 2 ( 0 ) 2 H ( 0 ) h ν ν = 0 + H ( 0 ) 2 h ν 2 ν = 0 + h ( 0 , t ) H ( 0 ) ,
a 3 ( z ) = 32 3 H 3 ( 0 ) 3 H ( 0 ) h ν ν = 0 + h ( 0 , t ) H ( 0 ) + 3 H ( 0 ) 2 h ν 2 ν = 0 + H ( 0 ) 3 h ν 3 ν = 0 ,
where H and its derivatives at the point ν = 0 are given in Table 1, and
h 0 , t = z 1 2 α * 2 z , h ν ν = 0 = z 1 2 α * z 3 / 2 2 z 1 α * z ,
2 h ν 2 ν = 0 = 3 z 5 / 2 1 z 3 / 2 , 3 h ν 3 ν = 0 = 15 α * z 7 / 2 3 α * z 5 / 2 .
The functions ϕ 0 ( z ) , ϕ 1 ( z ) , ϕ 2 ( z ) , and ϕ 3 ( z ) in Table 1, have the form
ϕ 0 = Γ 1 4 B 1 χ z 1 3 12 α * 3 p 1 / 4 , ϕ 1 = 6 π B 1 z 1 2 z z + 5 20 α * 2 p z ,
ϕ 2 = Γ 3 4 B 1 χ 3 200 α * p 3 / 4 z 3 / 2 25 + 5 z 66 z 3 / 2 + 48 z 2 18 z 5 / 2 + 6 z 3 ,
ϕ 3 = 75 + 20 z 99 z 2 + 77 z 5 / 2 111 z 3 + 57 z 7 / 2 19 z 4 100 p z 5 / 2 .
Note that
z = z ( ξ ( t ) ) = 1 + 2 α * ξ ( t ) , ξ ( t ) = 0 t υ ( t ) d t .
Substituting now
d υ d t = d υ d ξ d ξ d t = υ d υ d ξ
to the left side of Equation (28) and taking expression (30) into account, we obtain the differential equation for dimensionless metastability degree
d υ d ξ = B 1 k = 0 p ( k + 1 ) / 4 a k ( ξ ) .
By integrating it, we obtain the solution for υ as
υ ( ξ ) = B 1 p 1 / 4 0 ξ k = 0 p k / 4 a k ( ξ ) d ξ + 1 .
Further, we find an expression for t ( ξ ) to determine υ as a function of time t. Substituting υ ( ξ ) from (32) into the relation υ = d ξ / d t , we obtain
t ( ξ ) = 0 ξ d ξ υ ( ξ ) .
The Equations (32) and (33) determine the dimensionless metastability degree υ as a function of time t in parametric form. Figure 1 shows this function for the two mechanisms of kinetic crystal growth from supercooled melt, where υ = υ 0 is the main contribution to the sum in expression (32) and υ = υ i ( i = 1 , 2 , 3 ) represents the i-approximation of the complete solution for υ
υ = υ i = B 1 p 1 / 4 0 ξ k = 0 i p k / 4 a k ( ξ ) d ξ + 1 .
Figure 1 illustrates dimensionless supercooling υ to the dimensionless time t. If the melt was initially supercooled (see Figure 1a,b), the rate of supercooling removal is high (this can be explained mathematically by the inflection point in the Figure 1b for small t). Otherwise, υ 2 and υ 3 in Figure 1b,c show that the initially fast rate of metastability removal decreases with time as nucleation slows and coarsening of existing crystallites continues.
The given solutions allow us to determine an equation for the density function of crystal size distribution. Considering the maximum size of crystals appearing at the zero moment of time
ζ m ( ξ ( t ) ) = 1 + 2 α * ξ ( t ) 1 α * ,
we arrive at
Φ ( ξ ( t ) , ζ ) = 1 + α * ζ E 1 1 ( t , ζ ) exp p ( 1 E 1 2 ) E 1 2 ,
where
E 1 ( ξ ( t ) , ζ ) = 1 B 1 0 ξ ( t ) ζ α * ζ 2 / 2 k = 0 p ( k + 1 ) / 4 a k ( ξ ) d ξ .
Here ζ is limited by ζ m , ξ ( t ) is defined by (33) as an inverse function.
Figure 2 represents the evolution of the dimensionless density of the distribution function Φ . It can be seen that weak inflection points appear as time increases. This means that the high nucleation rate is typical for the initial stage of the process when supercooling/supersaturation is close to the value υ = 1 . The nucleation process is followed by coarsening of the crystallites, which reduces the metastability degree of the system, slows down the intensity of nucleation, and leads to inflection points on the density graphs of the distribution function. It is important to note that such behavior was previously observed in experiments (see, for example, [55,56,57,58]).
The total number of crystals in the metastable system is an important parameter for the intermediate stage of phase transitions
N ( ξ ( t ) ) = 1 l 0 3 0 ζ m ( ξ ( t ) ) Φ ( ξ ( t ) , ζ ) d ζ ,
where ζ m ( ξ ) and ξ ( t ) are defined by expressions (33) and (34). The other important parameter of evolving particulate ensembles is the average size of crystals
L ¯ ( ξ ( t ) ) = l 0 0 ζ m ( ξ ( t ) ) ζ Φ d ζ 0 ζ m ( ξ ( t ) ) Φ d ζ 1 .
Equations (36) and (37) are complete solutions written in parametric form.
The dimensionless number of crystals n i ( t ) and average crystal size u i ( t ) , corresponding to the i-approximation of the solution (35), are shown in Figure 3 (functions u i and n i are defined below). It is easy to see that approximations n 2 and u 2 correspond to the third-order solutions n 3 and u 3 . It is important to note the stability of the characteristics with increasing time at the intermediate stage of the phase transformation process, while the main contributions n 0 and u 0 previously found in [35], give only a rough approximation to the solution behavior.
The number N i of crystals and the average crystal size L ¯ i , corresponding to the i-approximation of the solution can be written as
N i ( ξ ( t ) ) = 1 l 0 3 0 ζ m ( ξ ( t ) ) Φ i ( ξ ( t ) , ζ ) d ζ , L ¯ i ( ξ ( t ) ) = l 0 0 ζ m ( ξ ( t ) ) ζ Φ i d ζ 0 ζ m ( ξ ( t ) ) Φ i d ζ 1 ,
where
Φ i ( ξ ( t ) , ζ ) = 1 + α * ζ E 1 i 1 ( t , ζ ) exp p ( 1 E 1 i 2 ) E 1 i 2 ,
E 1 i ( ξ ( t ) , ζ ) = 1 B 1 0 ξ ( t ) ζ α * ζ 2 / 2 k = 0 i p ( k + 1 ) / 4 a k ( ξ ) d ξ .
The relative number of particles n i and the average crystal radius u i , shown in Figure 3, can be written in the form
n i ( ξ ( t ) ) = N i ( ξ ( t ) ) N 3 ( ξ * ) , u i ( ξ ( t ) ) = L ¯ i ( ξ ( t ) ) L ¯ 3 ( ξ * ) ,
where, the parameter ξ * corresponds to the time t * when supercooling/supersaturation of the system is completely removed. This dependence is defined by Equation (33), where t * = t ( ξ * ) .

4. The Model with Fluctuations in Crystal Growth Rates

Experiments demonstrate that the crystal growth rate has random fluctuations [59,60,61,62] occuring a random crystal radius r ( τ ) . Thus, the kinetic equation for the density ϕ can be written as
ϕ τ + r d r d τ ϕ = r D ϕ r , τ > 0 , r > r * ,
where D is a function determining fluctuation rate. Then, equating the crystal flux of minimum size r * to the nucleation frequency I, the boundary conditions are found as
d r d τ ϕ D ϕ r = I Δ T , r = r * , τ > 0 ,
ϕ = 0 , r ,
where
d r d τ = β * Δ T ,
D = d 1 d r d τ = β * d 1 Δ T .
Parameter d 1 is the proportionality constant. Note that this theoretical assumption is also confirmed experimentally [63]. Thus, the equations and conditions (12), (13), (38)–(42) represent the complete statement of the problem.
Let us now rewrite this problem statement in dimensionless variables (16) ( Φ = Φ ( t , ζ ) , r * 0 , υ = υ ( t ) ) as
Φ t + υ Φ ζ = u 0 υ 2 Φ ζ 2 , t > 0 , ζ > 0 ,
Φ = 0 , ζ ,
υ = 1 B 0 ζ 3 Φ ( t , ζ ) d ζ , t > 0 ,
J ( υ ) 1 υ exp p φ ( υ ) = Φ u 0 Φ ζ , ζ = 0 .
Here B = B 1 / 3 , u 0 = d 1 / l 0 , Equation (45) is time-integrated with (38). Initial conditions (19) in the case of WVFZ kinetics can be written as φ ( υ ) = 1 υ 2 and in the case of Meirs kinetic as φ ( υ ) = ln ( υ ) .
Entering the modified time
ξ ( t ) = 0 t υ ( t 1 ) d t 1 ,
we rewrite Equation (43) in the following form ( υ = υ ( ξ ) , Φ = Φ ( ξ , ζ ) )
Φ ξ + Φ ζ = u 0 2 Φ ζ 2 , ξ > 0 , ζ > 0 .
The equation may be solved using the Laplace integral method (the initial condition gives ξ = 0 and Φ = 0 ).
Considering the boundary conditions (44) and (46), we write the final solution of Equation (47) as
Φ ( ξ , ζ ) = exp ζ 2 u 0 0 ξ J [ υ ( ξ η ) ] u 0 exp η 4 u 0 1 π η exp ζ 2 4 u 0 η
1 2 u 0 exp ζ 2 u 0 + η 4 u 0 erfc ζ 2 u 0 η + η 2 u 0 d η .
Substituting the Equation (48) into (45) and evaluating the integral over ζ , we obtain
υ ( ξ ) = 1 B 0 ξ J [ υ ( ξ η ) ] h ( η ) d η ,
where
h ( η ) = 2 π exp η 4 u 0 1 2 u 0 1 / 2 η 5 / 2 + 7 2 u 0 3 / 2 η 3 / 2 3 u 0 5 / 2 η 1 / 2
+ π 9 2 u 0 η 2 + 3 u 0 3 + η 3 2 erfc η 2 u 0 6 u 0 3 + 9 2 u 0 η 2 + η 3 2 ,
J [ υ ( ξ η ) ] = 1 υ ( ξ η ) exp p φ [ υ ( ξ η ) ] .
Next, we replace the variable by the inverse function η = η ( υ ) to solve the Equation (49) and arrive at
υ = 1 B 1 υ J ( υ 1 ) h [ η ( υ 1 ) ] d η d υ 1 d υ 1 .
Differentiating the Equation (50) over the variable υ gives an equation with separate variables, the integration of which determines the implicit expression for υ ( ξ ) as
1 υ d υ 1 J ( υ 1 ) = 1 υ υ 1 exp p φ ( υ 1 ) d υ 1 = B 0 ξ h ( η ) d η H ( ξ ) .
Substituting the expressions for φ , we find υ ( ξ ) for WVFZ kinetics in a form
υ 2 exp p 1 1 υ 2 1 p exp p E p ( υ ) = 2 H ( ξ ) ,
and for Meirs kinetics
υ ( ξ ) = ( 2 p ) H ( ξ ) + 1 1 / ( 2 p ) , p 2 ; υ ( ξ ) = exp H ( ξ ) , p = 2 .
Here E p ( ζ ) denotes the integral function
E p ( υ ) = p p / υ 2 exp ( v ) d v v .
The analytical solutions (51)–(53) satisfy the integral Equation (49).
Given the density of the distribution function (48), find the number N of crystals as a function of the parameter ξ
N = 0 ϕ ( τ , r ) d r = 1 l 0 3 0 Φ ( ξ , ζ ) d ζ = 1 l 0 3 0 ξ J υ ( ξ η ) d η .
Now rewriting ξ ( t ) , we obtain the solution with parameter ξ as
t ( ξ ) = 0 ξ d ξ 1 υ ( ξ 1 ) .
Equations (48), (51)–(55) are exact analytical solutions of the integro-differential model describing nucleation and growth of crystals at the intermediate stage of phase transitions. The solution is expressed in a general form that is valid for various kinetic mechanisms.
Figure 4 and Figure 5 show an exact analytical solutions according to expressions (48), (53), and (55) for Meirs nucleation kinetics. It can be seen that the removal of supercooling slows down with time when the nucleation process is not intense enough (when ξ is greater than or near unity in Figure 4), while the density of the distribution function becomes larger with increasing time (when ξ increases in Figure 5).
Figure 6 demonstrates the transition from a bell-shaped distribution to a distribution with a break point as the parameter u 0 decreases. When u 0 is small enough, the distribution function has the behavior described, for example, in [64], where the early stages of the Ostwald ripening process were observed. It shows that the diffusion term in Equation (38) is important in the initial stages of the intermediate phase transition. Figure 7 compares the theory with experimental data on crystallization kinetics of lysozyme and canavalin. The metastability reduction is in good agreement with experimental data for a broad range of time changes, as can easily be seen. Significantly, the analytical method defining the intermediate stage of phase transitions can be applied to determine the initial crystal distribution at the final stage of phase transformation and the Ostwald ripening stage (when the supercooling of the system tends to zero).

5. Conclusions

This paper deals with the intermediate stage of a phase transformation in one-component metastable systems (melts or solutions). The main results are as follows:
(i)
A complete analytic solution of the integro-differential model that describes the intermediate phase transition in one-component melts and solutions without taking into account fluctuations in the crystal growth rate has been found. Within the framework of this model, an exact analytical solution of the kinetic equation is obtained, and the density of the crystal size distribution function is found. An integro-differential equation for the metastability degree has been derived. A complete analytical solution to this equation is constructed based on the saddle-point method for calculating the Laplace-type integral.
(ii)
Our studies show that the found second approximation of the asymptotic solution is sufficient for quantitative analysis of the solution to the problem. In other words, our study demonstrates the sign-variable convergence of solutions to the second approximation, which indicates the fundamental importance of taking into account the first and second corrections in the asymptotic solution of the problem. The convergence of the found solution is ensured by a large value of the dimensionless Gibbs number p and/or by the smallness of the parameter α * , which is typical for a wide class of real systems. An explicit form of the first four coefficients of the found asymptotic solutions for the Weber-Volmer-Frenkel-Zel’dovich and Meiers nucleation kinetics is derived.
(iii)
Within the framework of the obtained solution, it is shown that the metastability degree and density of the distribution function have inflection points (observed experimentally) responsible for different stages of the phase transition process. At the initial times, when supercooling/supersaturation is large, the appearance of nuclei is the dominant process; then this process is followed by the growth of crystallites, which becomes dominant at the final stages of the intermediate phase transition when the metastability of the system is sufficiently small and nucleation proceeds unintensively.
(iv)
An exact analytical solution has been found for an integro-differential model describing the intermediate stage of the phase transition in one-component melts and solutions taking fluctuations in the crystal growth rate into account. Namely, the density of the crystal size distribution function has been determined and an implicit expression for metastability degree of the system in the case of different nucleation kinetics has been driven. The parametric solution to the problem constructed in this work is detailed for the frequently encountered Weber-Volmer-Frenkel-Zel’dovich and Meirs nucleation kinetics.
(v)
It is shown that the exact analytical solution to the integro-differential model with allowance for the “diffusion” term in the Fokker-Planck equation passes into a complete solution of the “non-diffusion” model (when the “diffusion” term in the kinetic equation is small enough). The analytical solutions, describing the evolution of the phase transformation at the intermediate stage, can be used as initial conditions at the final stage of the phase transition (e.g., at the stage of Ostwald ripening or coagulation of crystals).
In general, the paper shows methods for solving the kinetic and balance equations for the intermediate stage of phase transformation with and without considering fluctuations in crystal growth rates (respectively in the case of the kinetic equation of second and first orders). The solutions found can be used as initial conditions for the final phase of the phase transformation, which is described by different models with consideration of Ostwald ripening of the particles, their possible coagulation, and fragmentation [66,67,68,69,70]. The theory developed in this paper may be generalized to the case of a simultaneous phase transformation in the volume of the metastable phase and in a given spatial direction (directional phase transition), which is caused by the presence of temperature gradients. This situation occurs when crystals grow in the moving region of the phase transformation—the two-phase layer. A theory combining volumetric and directional phase transformation can be developed in the spirit of Refs. [71,72,73,74,75,76,77,78,79,80].

Author Contributions

Conceptualization, L.V.T., E.V.M. and D.V.A.; methodology, D.V.A.; software, L.V.T., E.V.M. and S.I.O.; validation, E.V.M. and A.P.M.; formal analysis, A.P.M. and E.V.M.; investigation, E.V.M.; resources, L.V.T.; writing—original draft preparation, L.V.T., Y.Y. and D.V.A.; writing—review and editing, Y.Y. and D.V.A.; visualization, L.V.T.; supervision, E.V.M. and D.V.A.; project administration, L.V.T.; funding acquisition, L.V.T., E.V.M. and S.I.O. All authors have read and agreed to the published version of the manuscript.

Funding

Authors gratefully acknowledge financial support from the Russian Science Foundation (project no. 21-79-10012).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analysed during this study are included in this published article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Madras, G.M.; McCoy, B.J. Transition from nucleation and growth to Ostwald ripening. Chem. Eng. Sci. 2002, 57, 3809–3818. [Google Scholar] [CrossRef]
  2. Zhang, T.H.; Liu, X.Y. Nucleation: What happens at the initial stage. Angew. Chem. Int. Ed. 2009, 48, 1308–1312. [Google Scholar] [CrossRef] [PubMed]
  3. Uwaha, M.; Koyama, K. Transition from nucleation to ripening in the classical nucleation model. J. Cryst. Growth 2010, 312, 1046–1054. [Google Scholar] [CrossRef]
  4. Makoveeva, E.V.; Alexandrov, D.V. On the theory of phase transformation process in a binary supercooled melt. Eur. Phys. J. Spec. Top. 2020, 229, 375–382. [Google Scholar] [CrossRef]
  5. Alexandrov, D.V.; Alexandrova, I.V. On the theory of the unsteady-state growth of spherical crystals in metastable liquids. Phil. Trans. R. Soc. A 2019, 377, 20180209. [Google Scholar] [CrossRef] [Green Version]
  6. Chernov, A.A. Modern Crystallography III; Springer: Berlin/Heidelberg, Germany, 1984. [Google Scholar]
  7. Alexandrova, I.V.; Alexandrov, D.V. Dynamics of particulate assemblages in metastable liquids: A test of theory with nucleation and growth kinetics. Phil. Trans. R. Soc. A 2020, 378, 20190245. [Google Scholar] [CrossRef] [Green Version]
  8. Chalmers, B. Physical Metallurgy; Wiley: New York, NY, USA, 1959. [Google Scholar]
  9. Todes, O.M.; Seballo, V.A.; Gol’tsiker, L.D. Bath Crystallization from Solutions; Khimia: Leningrad, Russia, 1984. [Google Scholar]
  10. Poon, J.M.-H.; Immanuel, C.D.; Doyle, F.J., III; Litser, J.D. A threedimensional population balance model of granulation with a mechanistic representation of the nucleation and aggregation phenomena. Chem. Eng. Sci. 2008, 63, 1315–1329. [Google Scholar] [CrossRef]
  11. Poon, J.M.-H.; Ramachandran, R.; Sanders, C.F.W.; Glaser, T.; Immanuel, C.D.; Doyle, F.J., III; Litser, J.D.; Stepanek, F.; Wang, F.Y.; Cameron, I.T. Experimental validation studies on a multi-scale and multi-dimensional population balance model of batch granulation. Chem. Eng. Sci. 2009, 64, 775–786. [Google Scholar] [CrossRef]
  12. Alexandrov, D.V.; Dubovoi, G.Y.; Malygin, A.P.; Nizovtseva, I.G.; Toropova, L.V. Solidification of ternary systems with a nonlinear phase diagram. Russ. Metall. 2017, 2017, 127–135. [Google Scholar] [CrossRef]
  13. Makoveeva, E.V.; Alexandrov, D.V.; Ivanov, A.A. Mathematical modeling of crystallization process from a supercooled binary melt. Math. Meth. Appl. Sci. 2021, 44, 12244–12251. [Google Scholar] [CrossRef]
  14. Thompson, C.V.; Spaepen, F. Homogeneous crystal nucleation in binary metallic melts. Acta Metall. 1983, 31, 2021–2027. [Google Scholar] [CrossRef]
  15. Kelton, K.F.; Greer, A.L. Nucleation in Condensed Matter: Applications in Materials and Biology; Elsevier: Amsterdam, The Netherlands, 2010. [Google Scholar]
  16. Zettlemoyer, A.C. Nucleation; Dekker: New York, NY, USA, 1969. [Google Scholar]
  17. Alexandrov, D.V.; Nizovtseva, I.G.; Alexandrova, I.V. On the theory of nucleation and nonstationary evolution of a polydisperse ensemble of crystals. Int. J. Heat Mass Trans. 2019, 128, 46–53. [Google Scholar] [CrossRef]
  18. Alexandrov, D.V. Kinetics of particle coarsening with allowance for Ostwald ripening and coagulation. J. Phys. Condens. Matter 2016, 28, 035102. [Google Scholar] [CrossRef] [PubMed]
  19. Volmer, M. Kinetik der Phasenbildung; Steinkopf: Dresden, Germany, 1939. [Google Scholar]
  20. Frenkel, J. Kinetic Theory of Liquids; Dover: New York, NY, USA, 1955. [Google Scholar]
  21. Velmurugan, J.; Noel, J.-M.; Nogala, W.; Mirkin, M.V. Nucleation and growth of metal on nanoelectrodes. Chem. Sci. 2012, 3, 3307–3314. [Google Scholar] [CrossRef]
  22. Lifshitz, I.M.; Slyozov, V.V. The kinetics of precipitation from supersaturated solid solutions. J. Phys. Chem. Solids 1961, 19, 35–50. [Google Scholar] [CrossRef]
  23. Tokuyama, M.; Kawasaki, M.; Enomoto, V. Kinetic equations for Ostwald ripening. Physics A 1986, 134, 323–338. [Google Scholar] [CrossRef]
  24. Akaiwa, N.; Meiron, D.I. Numerical simulation of two-dimensional latestage coarsening for nucleation and growth. Phys. Rev. E 1995, 51, 5408–5421. [Google Scholar] [CrossRef] [Green Version]
  25. Farjoun, Y.; Neu, J.C. Aggregation according to classical kinetics: From nucleation to coarsening. Phys. Rev. E 2011, 83, 051607. [Google Scholar] [CrossRef] [Green Version]
  26. Alexandrova, I.V.; Alexandrov, D.V.; Makoveeva, E.V. Ostwald ripening in the presence of simultaneous occurrence of various mass transfer mechanisms: An extension of the Lifshitz-Slyozov theory. Phil. Trans. R. Soc. A 2021, 379, 20200308. [Google Scholar] [CrossRef]
  27. Alexandrov, D.V. On the theory of Ostwald ripening: Formation of the universal distribution. J. Phys. A Math. Theor. 2015, 48, 035103. [Google Scholar] [CrossRef]
  28. Alexandrov, D.V. Relaxation dynamics of the phase transformation process at its ripening stage. J. Phys. A Math. Theor. 2015, 48, 245101. [Google Scholar] [CrossRef]
  29. Avrami, M. Granulation, phase change, and microstructure kinetics of phase change. III. J. Chem. Phys. 1991, 9, 177–184. [Google Scholar] [CrossRef]
  30. Aastuen, D.J.W.; Clark, N.A.; Swindal, J.C.; Muzny, C.D. Determination of the colloidal crystal nucleation rate density. Phase Trans. 1990, 21, 139–155. [Google Scholar] [CrossRef]
  31. Binder, K.; Stauffer, D. Statistical theory of nucleation, condensation and coagulation. Adv. Phys. 1976, 25, 343–396. [Google Scholar] [CrossRef]
  32. Langer, J.S.; Schwartz, A.J. Kinetics of nucleation in near-critical Fluids. Phys. Rev. A 1980, 21, 948–958. [Google Scholar] [CrossRef]
  33. Shneidman, V.A. Transient nucleation with a monotonically changing barrier. Phys. Rev. E 2010, 82, 031603. [Google Scholar] [CrossRef] [Green Version]
  34. Shneidman, V.A. Time-dependent distributions in self-quenching nucleation. Phys. Rev. E 2011, 84, 031602. [Google Scholar] [CrossRef] [Green Version]
  35. Buyevich, Y.A.; Mansurov, V.V. Kinetics of the intermediate stage of phase transition in batch crystallizatione. J. Cryst. Growth 1990, 104, 861–867. [Google Scholar] [CrossRef]
  36. Buyevich, Y.A.; Ivanov, A.O. Kinetics of phase separation in colloids II. Non-linear evolution of a metastable colloid. Physics A 1993, 193, 221–240. [Google Scholar] [CrossRef]
  37. Ivanov, A.O.; Zubarev, A.Y. Non-linear evolution of a system of elongated droplike aggregates in a metastable magnetic Fluid. Physics A 1998, 251, 348–367. [Google Scholar] [CrossRef]
  38. Barlow, D.A. Theory of the intermediate stage of crystal growth with applications to protein crystallization. J. Cryst. Growth 2009, 311, 2480–2483. [Google Scholar] [CrossRef]
  39. Barlow, D.A.; Baird, J.K.; Su, C.-H. Theory of the von Weimarn rules governing the average size of crystals precipitated from a supersaturated solution. J. Cryst. Growth 2004, 264, 417–423. [Google Scholar] [CrossRef]
  40. Volmer, M.; Weber, A. Nukleation in übersättigten Strukturen. Z. Phys. Chem. 1926, 119, 277–301. [Google Scholar] [CrossRef]
  41. Zel’dovich, J.B. On the theory of formation of new phases: Cavitation. J. Exp. Theor. Phys. 1942, 12, 525–538. [Google Scholar]
  42. Lifshitz, E.M.; Pitaevskii, L.P. Physical Kinetic; Pergamon Press: Oxford, UK, 1981. [Google Scholar]
  43. Landau, L.D.; Lifshitz, E.M. Statistical Physics; Pergamon Press: Oxford, UK, 1980. [Google Scholar]
  44. Mullin, J.W. Crystallization; Butterworths: London, UK, 1972. [Google Scholar]
  45. Gherras, N.; Fevotte, G. Comparison between approaches for the experimental determination of metastable zone width: A case study of the batch cooling crystallization of ammonium oxalate in water. J. Cryst. Growth 2012, 342, 88–98. [Google Scholar] [CrossRef] [Green Version]
  46. Kidyarov, B.I. Kinetics of Crystal Formation from the Liquid Phase; Nauka Publishing House: Novosibirsk, Russia, 1979. [Google Scholar]
  47. Alexandrov, D.V.; Alexandrova, I.V.; Ivanov, A.A.; Malygin, A.P.; Starodumov, I.O.; Toropova, L.V. On the Theory of the Nonstationary Spherical Crystal Growth in Supercooled Melts and Supersaturated Solutions. Russ. Metall. 2019, 2019, 787–794. [Google Scholar] [CrossRef]
  48. Nyvlt, J.; Sohnel, O.; Matuchova, M.; Broul, M. The Kinetics of Industrial Crystallization; Elsevier: Amsterdam, The Netherlands, 1985. [Google Scholar]
  49. Randolph, A.; Larson, M. Theory of Particulate Processes; Academic Press: New York, NY, USA, 1988. [Google Scholar]
  50. Fedoruk, M.V. Saddle-Point Method; Nauka: Moscow, Russia, 1977. [Google Scholar]
  51. Alexandrov, D.V. Nonlinear dynamics of polydisperse assemblages of particles evolving in metastable media. Eur. Phys. J. Spec. Top. 2020, 229, 383–404. [Google Scholar] [CrossRef]
  52. Skripov, V.P. Methastable Liquids; Wiley: New York, NY, USA, 1974. [Google Scholar]
  53. Porter, D.A.; Easterling, K.E. Phase Transformations in Metalls and Alloys; Van Nostrand Reinhold: New York, NY, USA, 1981. [Google Scholar]
  54. Hamza, M.A.; Berge, B.; Mikosch, W.; Ruhl, E. Homogeneous nucleation of supersaturated KCl-solutions from single levitated microdroplets. Phys. Chem. Chem. Phys. 2004, 6, 3484–3489. [Google Scholar] [CrossRef]
  55. Janse, A.H. Nucleation and Crystal Growth in Batch Crystallizers; Delft University of Technology: Delft, The Netherlands, 1977. [Google Scholar]
  56. Pot, A. Industrial Sucrose Crystallization; Delft University of Technology: Delft, The Netherlands, 1980. [Google Scholar]
  57. Leubner, I.H. Balanced nucleation and growth model for controlled crystal size distribution. J. Disp. Sci. Technol. 2002, 23, 577–590. [Google Scholar] [CrossRef]
  58. Hanhoun, M.; Montastruc, L.; Azzaro-Pantel, C.; Biscans, B.; Freche, M.; Pibouleau, L. Simultaneous determination of nucleation and crystal growth kinetics of struvite using a thermodynamic modeling approach. Chem. Eng. J. 2013, 215–216, 903–912. [Google Scholar] [CrossRef] [Green Version]
  59. Zumstein, R.C.; Rousseau, R.W. Growth rate dispersion by initial growth rate distributions and growth rate fluctuations. AIChE J. 1987, 33, 121–129. [Google Scholar] [CrossRef]
  60. Buyevich, Y.A.; Alexandrov, D.V.; Mansurov, V.V. Macrokinetics of Crystallization; Begell House: New York, NY, USA, 2001. [Google Scholar]
  61. Toropova, L.V.; Aseev, D.L.; Osipov, S.I.; Ivanov, A.A. Mathematical modeling of bulk and directional crystallization with the moving phase transition layer. Math. Meth. Appl. Sci. 2021. [Google Scholar] [CrossRef]
  62. Toropova, L.V.; Alexandrov, D.V.; Rettenmayr, M.; Liu, D. Microstructure and morphology of Si crystals grown in pure Si and Al-Si melts. J. Phys. Condens. Matter 2022, 34, 094002. [Google Scholar] [CrossRef] [PubMed]
  63. Melikhov, I.V.; Beloussova, T.Y.; Uludev, N.A.; Blyudev, N.T. Fluctuain the rate of growth of monocrystals. Crystallography 1974, 19, 1263–1268. [Google Scholar]
  64. Shneidman, V. Early stages of Ostwald ripening. Phys. Rev. E 2013, 88, 010401. [Google Scholar] [CrossRef] [Green Version]
  65. Caraballo, K.G.; Baird, J.K.; Ng, J.D. Comparison of the crystallization kinetics of canavalin and lysozyme. Cryst. Growth Des. 2006, 6, 874–888. [Google Scholar] [CrossRef]
  66. Friedlander, S.K. Smoke, Dust and Haze: Fundamentals of Aerosol Behaviour; Wiley: New York, NY, USA, 1977. [Google Scholar]
  67. Herlach, D.M. Phase Transformations in Multicomponent Melts; Wiley-VCH: Weinheim, Germany, 2008. [Google Scholar]
  68. Slezov, V.V. Kinetics of First-Order Phase Transitions; Wiley: Weinheim, Germany, 2009. [Google Scholar]
  69. Alexandrov, D.V. On the theory of fragmentation process with initial particle volume. Commun. Theor. Phys. 2017, 68, 269–271. [Google Scholar] [CrossRef]
  70. Makoveeva, E.V.; Alexandrov, D.V. Effects of external heat/mass sources and withdrawal rates of crystals from a metastable liquid on the evolution of particulate assemblages. Eur. Phys. J. Spec. Top. 2019, 228, 25–34. [Google Scholar] [CrossRef]
  71. Huppert, H.E.; Worster, M.G. Dynamic solidification of a binary melt. Nature 1985, 314, 703–707. [Google Scholar] [CrossRef]
  72. Kerr, R.C.; Woods, A.W.; Worster, M.G.; Huppert, H.E. Solidification of an alloy cooled from above Part 1. Equilibrium growth. J. Fluid Mech. 1990, 216, 323–342. [Google Scholar] [CrossRef] [Green Version]
  73. Peppin, S.S.L.; Aussillous, P.; Huppert, H.E.; Worster, M.G. Steady-state mushy layers: Experiments and theory. J. Fluid Mech. 2007, 570, 69–77. [Google Scholar] [CrossRef] [Green Version]
  74. Peppin, S.S.L.; Huppert, H.E.; Worster, M.G. Steady-state solidification of aqueous ammonium chloride. J. Fluid Mech. 2008, 599, 465–476. [Google Scholar] [CrossRef] [Green Version]
  75. Toropova, L.V.; Galenko, P.K.; Alexandrov, D.V.; Rettenmayr, M.; Kao, A.; Demange, G. Non-axisymmetric growth of dendrite with arbitrary symmetry in two and three dimensions: Sharp interface model vs. phase-field model. Eur. Phys. J. Spec. Top. 2020, 229, 2899–2909. [Google Scholar] [CrossRef]
  76. Toropova, L.V. Shape functions for dendrite tips of SCN and Si. Eur. Phys. J. Spec. Top. 2022, 231, 1129–1133. [Google Scholar] [CrossRef]
  77. Alexandrov, D.V.; Malygin, A.P. Analytical description of seawater crystallization in ice fissures and their influence on heat exchange between the ocean and the atmosphere. Dokl. Earth Sci. 2006, 411, 1407–1411. [Google Scholar] [CrossRef]
  78. Toropova, L.V.; Alexandrov, D.V. Dynamical law of the phase interface motion in the presence of crystals nucleation. Sci. Rep. 2022. [Google Scholar] [CrossRef]
  79. Alexandrov, D.V.; Bashkirtseva, I.A.; Ryashko, L.B. Nonlinear dynamics of mushy layers induced by external stochastic fluctuations. Phil. Trans. R. Soc. A 2018, 376, 20170216. [Google Scholar] [CrossRef] [Green Version]
  80. Nizovtseva, I.G.; Alexandrov, D.V. The effect of density changes on crystallization with a mushy layer. Phil. Trans. R. Soc. A 2020, 378, 20190248. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Dimensionless supercooling υ = Δ T Δ T 0 as a function of dimensionless time t = τ τ 0 : (a) WVFZ kinetics at Δ T 0 = 373 K, p = 66.9 ; (b) WVFZ kinetics at Δ T 0 = 573 K, p = 7.4 ; (c) Meirs kinetics at Δ T 0 = 573 K, p = 4 . The physical parameters are typical for supercooled melts [52,53,54]: γ i = 0.3 J m 2 , L V = 7 × 10 9 J m 3 , T p = 1273 K, ρ s = 7.8 × 10 3 kg m 3 , ρ m = ρ l = 7 × 10 3 kg m 3 , β * = 0.5 × 10 7 m/s 1 K 1 , λ l = 63 J m 1 K 1 s 1 , I 0 = 10 11 m 3 s 1 , C m = 840 Jkg 1 K 1 .
Figure 1. Dimensionless supercooling υ = Δ T Δ T 0 as a function of dimensionless time t = τ τ 0 : (a) WVFZ kinetics at Δ T 0 = 373 K, p = 66.9 ; (b) WVFZ kinetics at Δ T 0 = 573 K, p = 7.4 ; (c) Meirs kinetics at Δ T 0 = 573 K, p = 4 . The physical parameters are typical for supercooled melts [52,53,54]: γ i = 0.3 J m 2 , L V = 7 × 10 9 J m 3 , T p = 1273 K, ρ s = 7.8 × 10 3 kg m 3 , ρ m = ρ l = 7 × 10 3 kg m 3 , β * = 0.5 × 10 7 m/s 1 K 1 , λ l = 63 J m 1 K 1 s 1 , I 0 = 10 11 m 3 s 1 , C m = 840 Jkg 1 K 1 .
Crystals 12 00895 g001
Figure 2. The dimensionless distribution function Φ ( t , ζ ) as a function of ζ = r l 0 at different times t = τ τ 0 for supercooled melts with WVFZ kinetics: (a) p = 669 , Δ T 0 = 373 K, (1) t = 0.05 , ξ = 0.05 , (2) t = 0.254 , ξ = 0.25 , (3) t = 0.503 , ξ = 0.45 , and (4) t = 1.513 , ξ = 0.591 ; (b) p = 7.4 , Δ T 0 = 573 K, (1) t = 0.051 , ξ = 0.05 , (2) t = 0.263 , ξ = 0.25 , (3) t = 0.496 , ξ = 0.45 , and (4) t = 1.95 , ξ = 0.71 . The maximal size ζ m of crystals is shown by vertical lines at different moments of time. The curves denoted by simbole (4) corresponding to the moment when supercooling is removed ( υ = 0 ).
Figure 2. The dimensionless distribution function Φ ( t , ζ ) as a function of ζ = r l 0 at different times t = τ τ 0 for supercooled melts with WVFZ kinetics: (a) p = 669 , Δ T 0 = 373 K, (1) t = 0.05 , ξ = 0.05 , (2) t = 0.254 , ξ = 0.25 , (3) t = 0.503 , ξ = 0.45 , and (4) t = 1.513 , ξ = 0.591 ; (b) p = 7.4 , Δ T 0 = 573 K, (1) t = 0.051 , ξ = 0.05 , (2) t = 0.263 , ξ = 0.25 , (3) t = 0.496 , ξ = 0.45 , and (4) t = 1.95 , ξ = 0.71 . The maximal size ζ m of crystals is shown by vertical lines at different moments of time. The curves denoted by simbole (4) corresponding to the moment when supercooling is removed ( υ = 0 ).
Crystals 12 00895 g002
Figure 3. The dimensionless number of crystals n i = N i ( ξ ) N 3 ( ξ * ) (a) and average crystal size u i = L ¯ i ( ξ ) L ¯ 3 ( ξ * ) (b) as a function of dimensionless time t = τ τ 0 for supercooled melts with WVFZ kinetics, p = 7.4 , Δ T 0 = 573 K. The vertical lines show the time of the process when supercooling is totally removed.
Figure 3. The dimensionless number of crystals n i = N i ( ξ ) N 3 ( ξ * ) (a) and average crystal size u i = L ¯ i ( ξ ) L ¯ 3 ( ξ * ) (b) as a function of dimensionless time t = τ τ 0 for supercooled melts with WVFZ kinetics, p = 7.4 , Δ T 0 = 573 K. The vertical lines show the time of the process when supercooling is totally removed.
Crystals 12 00895 g003
Figure 4. The dimensionless supercooling υ = Δ T Δ T 0 and time t = τ τ 0 depending on the modified time ξ . The physical parameters are typical for metallic alloys: u 0 = 10 2 , B = 14.92 , Δ T 0 = 300 °C, p = 2.2 .
Figure 4. The dimensionless supercooling υ = Δ T Δ T 0 and time t = τ τ 0 depending on the modified time ξ . The physical parameters are typical for metallic alloys: u 0 = 10 2 , B = 14.92 , Δ T 0 = 300 °C, p = 2.2 .
Crystals 12 00895 g004
Figure 5. The dimensionless density of the distribution function Φ as a function of the crystal radius at different time moments.
Figure 5. The dimensionless density of the distribution function Φ as a function of the crystal radius at different time moments.
Crystals 12 00895 g005
Figure 6. The dimensionless density of the distribution function Φ as a function of the crystal radius at different values of the parameter u 0 ( ξ = 1 ).
Figure 6. The dimensionless density of the distribution function Φ as a function of the crystal radius at different values of the parameter u 0 ( ξ = 1 ).
Crystals 12 00895 g006
Figure 7. Comparison of theory (solid line) with experiment (dots) [65]. Calculated parameters in Equation (53): u 0 = 10 1 , p = 4.5 .
Figure 7. Comparison of theory (solid line) with experiment (dots) [65]. Calculated parameters in Equation (53): u 0 = 10 1 , p = 4.5 .
Crystals 12 00895 g007
Table 1. Functions and parameters of the analytic solution for various kinetics. Here z = z ( ξ ( t ) ) = 1 + 2 α * ξ ( t ) , and ϕ 0 , ϕ 1 , ϕ 2 and ϕ 3 are given below.
Table 1. Functions and parameters of the analytic solution for various kinetics. Here z = z ( ξ ( t ) ) = 1 + 2 α * ξ ( t ) , and ϕ 0 , ϕ 1 , ϕ 2 and ϕ 3 are given below.
Supercooled MeltsSupersaturated Solutions
WVFZ KineticsMeirs KineticsWVFZ KineticsMeirs Kinetics
I = I * exp p Δ T / Δ T 0 2 I = I * Δ T p I = I * exp p ln 2 C l / C p I = I * Δ C p
Values(6) and (16)(16)(7) and (23)(23)
g ( υ ) (21)(22)(24)(25)
g ( IV ) | ν = 0 4 B 1 2 B 1 4 κ B 1 2 B 1
g ( V ) | ν = 0 24 B 1 α * 12 B 1 α * 24 κ B 1 α * 12 B 1 α *
g ( VI ) | ν = 0 180 B 1 α * 2 90 B 1 α * 2 180 κ B 1 α * 2 90 B 1 α * 2
g ( VII ) | ν = 0 1680 B 1 α * 3 840 B 1 α * 3 1680 κ B 1 α * 3 840 B 1 α * 3
H | ν = 0 χ / 4 2 1 / 4 χ / 4 κ 1 / 4 χ / 4 2 1 / 4 χ / 4
H | ν = 0 3 α * χ / 20 3 α * 2 1 / 4 χ / 20 3 κ 1 / 4 α * χ / 20 3 α * 2 1 / 4 χ / 20
H | ν = 0 9 α * 2 χ / 50 9 α * 2 2 1 / 4 χ / 50 9 κ 1 / 4 α * 2 χ / 50 9 α * 2 2 1 / 4 χ / 50
H | ν = 0 57 α * 3 χ / 100 57 α * 3 2 1 / 4 χ / 100 57 κ 1 / 4 α * 3 χ / 100 57 α * 3 2 1 / 4 χ / 100
υ 0 ( z ) 1 ϕ 0 ( z ) 1 2 1 / 4 ϕ 0 ( z ) 1 κ 1 / 4 ϕ 0 ( z ) 1 2 1 / 4 ϕ 0 ( z )
υ 1 ( z ) υ 0 ( z ) + ϕ 1 ( z ) υ 0 ( z ) + 2 ϕ 1 ( z ) υ 0 ( z ) + κ 1 / 2 ϕ 1 ( z ) υ 0 ( z ) + 2 ϕ 1 ( z )
υ 2 ( z ) υ 1 ( z ) + ϕ 2 ( z ) υ 1 ( z ) + 2 3 / 4 ϕ 2 ( z ) υ 1 ( z ) + κ 3 / 4 ϕ 2 ( z ) υ 1 ( z ) + 2 3 / 4 ϕ 2 ( z )
υ 3 ( z ) υ 2 ( z ) + ϕ 3 ( z ) υ 2 ( z ) + 2 ϕ 3 ( z ) υ 2 ( z ) + κ 1 ϕ 3 ( z ) υ 2 ( z ) + 2 ϕ 3 ( z )
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Toropova, L.V.; Makoveeva, E.V.; Osipov, S.I.; Malygin, A.P.; Yang, Y.; Alexandrov, D.V. Nucleation and Growth of an Ensemble of Crystals during the Intermediate Stage of a Phase Transition in Metastable Liquids. Crystals 2022, 12, 895. https://doi.org/10.3390/cryst12070895

AMA Style

Toropova LV, Makoveeva EV, Osipov SI, Malygin AP, Yang Y, Alexandrov DV. Nucleation and Growth of an Ensemble of Crystals during the Intermediate Stage of a Phase Transition in Metastable Liquids. Crystals. 2022; 12(7):895. https://doi.org/10.3390/cryst12070895

Chicago/Turabian Style

Toropova, Liubov V., Eugenya V. Makoveeva, Sergei I. Osipov, Alexey P. Malygin, Yang Yang, and Dmitri V. Alexandrov. 2022. "Nucleation and Growth of an Ensemble of Crystals during the Intermediate Stage of a Phase Transition in Metastable Liquids" Crystals 12, no. 7: 895. https://doi.org/10.3390/cryst12070895

APA Style

Toropova, L. V., Makoveeva, E. V., Osipov, S. I., Malygin, A. P., Yang, Y., & Alexandrov, D. V. (2022). Nucleation and Growth of an Ensemble of Crystals during the Intermediate Stage of a Phase Transition in Metastable Liquids. Crystals, 12(7), 895. https://doi.org/10.3390/cryst12070895

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop