Next Article in Journal
The Variation of Schottky Barrier Height Induced by the Phase Separation of InAlAs Layers on InP HEMT Devices
Next Article in Special Issue
Compositional Glass: A State with Inherent Chemical Disorder, Exemplified by Ti-rich Ni3(Al,Ti)1 D024 Phase
Previous Article in Journal
Mechanical Behavior and Microstructure of Graphene Oxide Electrodeposited Carbon Fiber Reinforced Cement-Based Materials
Previous 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
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

A Stable Mode of Dendritic Growth in Cases of Conductive and Convective Heat and Mass Transfer

by
Liubov V. Toropova
1,2,*,
Peter K. Galenko
2,3 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, Ekaterinburg 620000, 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, Ekaterinburg 620000, Russia
*
Author to whom correspondence should be addressed.
Crystals 2022, 12(7), 965; https://doi.org/10.3390/cryst12070965
Submission received: 24 June 2022 / Revised: 6 July 2022 / Accepted: 8 July 2022 / Published: 11 July 2022
(This article belongs to the Special Issue Phase Transition in External Fields)

Abstract

:
In this paper, we develop a theory of stable dendritic growth in undercooled melts in the presence of conductive and convective heat and mass transfer boundary conditions at the solid/liquid interface of a dendrite. To simplify the matter and construct the analytical theory, conductive and convective mechanisms are considered separately. Namely, the laws for total undercooling and selection criterion defining the stable growth mode (dendrite tip velocity and diameter) are derived for conductive and convective boundary conditions. To describe the case of simultaneous occurrence of these heat and mass transfer mechanisms, we sew together conductive and convective laws using power stitching functions. The generalised selection theory is compared with experimental data for Al 24 Ge 76 and Ti 45 Al 55 undercooled melts.

1. Introduction

A crystal with a main trunk and lateral branches extending along the main crystallographic axes of the crystal lattice is called a dendrite [1,2,3,4]. Growing single dendritic crystals and the intertwined side branches of adjacent dendrites form one of the main types of a solid structure in solidifying materials. As this takes place, the features of such a structure (e.g., its dispersion or mean interdendritic distance) are determined by heat and mass transfer processes, hydrodynamic and convective melt currents throughout the solidifying system, and the kinetics of atom attachment to the interfacial crystal surfaces [5,6,7].
One of the actual problems allowing theoretical description of dendritic structure (e.g., transition between monocrystalline and polycrystalline patterns) is the problem of selection for a stable mode of dendritic crystal growth. By solving this problem, we may find the selection criterion (the equation relating the steady-state growth rate V of the dendrite, its tip diameter ρ and melt undercooling Δ θ ). This equation and the total undercooling balance (the second equation linking these quantities) represent a nonlinear system of two algebraic equations to determine the governing dependencies V ( Δ θ ) and ρ ( Δ θ ) [8,9,10,11,12]. This system of equations has been intensively studied in recent decades for conductive heat and mass fluxes at solid–liquid dendritic interfaces (see, among others, recent review papers [13,14]). However, experimental data [15,16] show that such a model does not work over the entire range of melt undercooling. In many cases comparing experimental data with theoretical modeling predictions, it has been shown that convection may influence the growth kinetics if the dendrite growth velocity is comparable by magnitude with the value of the flow velocity [17]. This occurs at relatively “small” crystal growth velocity (mm per second or cm per second) and, respectively, at low undercooling. So, for example, when undercooling is small enough, the model no longer describes the experimental data. This can be explained by the fact that at small undercoolings, the convective melt fluxes in the thin boundary layer near the dendrite surface become essential [18]. In this case, the selection criterion and the balance of total undercooling change its form and are entirely determined by the convective mechanism of heat and mass transfer near the dendrite surface. Therefore, to describe the whole range of melt undercooling, we need to use the theory of convective heat and mass transfer for small undercooling and conductive heat and mass transfer for moderate and large undercooling. How to combine these approaches using a unified selection theory and undercooling balance is described below, where we sew together the corresponding dependencies to better describe experimental data within a broad range of melt undercooling.

2. The Heat and Mass Transfer Model of Dendritic Growth

For simplicity, we first consider the steady-state growth of a dendritic crystal in an undercooled binary melt illustrated in Figure 1. The heat-transfer and diffusion equations in the solid and melt phases read as
θ s τ = D θ Δ θ s ( solid phase ) ,
θ l τ + ( v · ) θ l = D θ Δ θ l , σ τ + ( v · ) σ = D σ Δ σ ( liquid phase ) ,
where θ s and θ l are the temperatures in solid and liquid, σ is the solute concentration in melt, τ is the time, D θ is the temperature conductivity coefficient, D σ is the diffusivity coefficient, v is the velocity of oncoming melt, ∇ is the nabla operator, and Δ is the Laplacian.
We consider the viscous melt model in the Oseen hydrodynamic approximation [19,20,21], as the more general Navier–Stokes model has no analytical solutions.
u v z = 1 ρ m p + μ Δ v , · v = 0 , ( liquid phase ) ,
where u is the flow velocity far from the crystal at z , z is the spatial coordinate corresponding to the direction of dendrite growth, ρ m is the melt density, p is the hydrodynamic pressure, and μ is the kinematic viscosity. Note that the last equation in (3) describes the incompressible fluid model.
The system of hydrodynamic Equation (3) must be supplemented with appropriate boundary conditions for fluid adhesion to the dendrite surface. Since the hydrodynamic problem about the flow of a viscous melt around a dendrite does not depend on the thermal and concentration distributions, but is an independent one, below we use its solution constructed in Refs. [22,23].
Considering the case of plane-parallel melt flow far from the dendrite, we fix temperature θ and solute concentration σ far from the growing crystal, i.e.,
θ l θ , σ σ ( far from the crystal ) .
Below, we consider two possible regimes of melt flow around the growing crystal. The first of them occurs in slow laminar currents when thermal and mass fluxes are described by classical Fourier’s laws. Dealing with this case, we have
V · n = D θ κ p Q θ s θ l · n , ( 1 k e ) σ V · n = D σ σ · n ( at the interface ) ,
where V stands for the growth rate of solid/liquid interface, κ p is the specific heat, Q is the latent heat parameter, and k e is the equilibrium segregation coefficient.
The second regime occurs in cases of intense hydrodynamic currents near the solid/liquid surface of a dendrite. If this is really the case, the heat and mass fluxes become of convective-type and boundary conditions read as [24,25,26,27,28]
θ Q D θ V · n = θ s · n + α θ ρ m κ p u * λ s ( θ i θ ) , ( at the interface ) ( 1 k e ) σ i V · n = α σ u * ( σ i σ ) , ( at the interface ) ,
where θ Q = Q / κ p is the adiabatic temperature, α θ and α σ are the convective heat and mass coefficients, and λ s is the temperature conductivity in the solid. The friction velocity u * , defining the intensity of flow near the dendrite, can be expressed in terms of shear stress s s as u * = ( s s / ρ m ) 1 / 2 [29]. In addition, the ratio of convective heat and mass coefficients is expressed as α θ / α σ = ( D θ / D σ ) n ¯ , 2 / 3 < n ¯ < 4 / 5 [24,30,31,32,33].
The last boundary condition defines the interfacial temperature in the form of
θ s = θ l = θ * m e σ θ Q d K β ˜ V · n , ( at the interface ) ,
where θ * is the phase transition temperature for a single-component melt ( σ = 0 ), m e is the equilibrium liquidus line slope, and d and β ˜ stand for the capillary length and kinetic coefficient, which due to the anisotropic growth, are the functions of polar angle θ [34] shown in Figure 1
d ( θ ) = d 0 1 χ d cos [ n ( θ θ d ) ] , β ˜ ( θ ) = β 0 θ Q 1 χ β cos [ n ( θ θ β ) ] .
Here, d 0 and β 0 are the constant factors. Note that expressions (8) work in the case of n-fold dendritic symmetry. The constants χ d and χ β are assumed to be much smaller than unity. It is also significant to note that the crystal curvature K equates to 2 / ρ and 4 / ρ in the 2D and 3D geometries, respectively. Here, ρ stands for the dendrite tip diameter. Following papers [14,18] we assume that the misalignment angles θ d and θ β are negligibly small.
It is convenient to solve the model (1)–(8) using special curvilinear coordinate systems (Figure 2). So, we use the parabolic cylinder reference frame in the 2D case and paraboloid of revolution reference frame in the 3D case. These curvilinear coordinates are introduced as follows
z = ρ ζ ξ 2 , x = ρ ξ ζ in 2 D , z = ρ ζ ξ 2 , x = ρ ξ ζ cos ϕ , y = ρ ξ ζ sin ϕ in 3 D .
Here, ρ represents the dendrite tip radius, and x, y, and z are the Cartesian coordinates. The coordinate ζ = 1 is introduced to be constant at the crystal surface and the surfaces ζ = c o n s t represent isothermal (isoconcentration) surfaces surrounding the growing dendrite.
Rewriting the model (1)–(8) in the coordinates (9) and considering that θ l = θ l ( ζ ) and σ = σ ( ζ ) (the temperature θ s in the solid phase is fixed), we obtain the steady-state solutions in the form of
θ l ( ζ ) = θ i + ( θ θ i ) I θ ( ζ ) I θ ( ) , σ ( ζ ) = σ i + ( σ σ i ) I σ ( ζ ) I σ ( ) .
Here, θ i and σ i are the temperature and solute concentration at the dendrite surface ζ = 1 , and I θ ( ζ ) and I σ ( ζ ) are the temperature and concentration integrals. These functions depend on the heat and mass transfer mechanism (conductive or convective) near the crystal and are given in Section 3.1 and Section 3.2.

3. The Total Undercooling Balance

The undercooling is a measure of the deviation of the temperature at any point of a system from the equilibrium temperature of phase co-existence. The full and experimentally measurable undercooling is a sum of a number of processes: heat and mass transport in the bulk and at the interface as well as the solid–liquid interface stability due to interface energy and attachment of particles (atoms, molecules, or clusters) to the interface. All these contributions give their own influence on the growth kinetics that is expressed by the following undercooling balance.
The melt undercooling Δ θ = θ * θ m e σ represents the driving force of dendritic growth and is given by the following terms
Δ θ = Δ θ θ + Δ θ σ + Δ θ ρ + Δ θ k ,
where Δ θ θ is the thermal undercooling, Δ θ σ is the concentration undercooling, Δ θ ρ is the undercooling connected with the curved phase transition interface, and Δ θ k is the kinetic undercooling. The first two of them are substantially dependent on the heat and mass transfer mechanism at the dendrite interface. The last two contributions are only dependent on dendrite tip diameter ρ and tip velocity V.
So, Δ θ ρ is given by
Δ θ ρ ( ρ ) = 2 d 0 θ Q ρ ( 2 D ) and Δ θ ρ ( ρ ) = 4 d 0 θ Q ρ ( 3 D ) .
The kinetic undercooling Δ θ k reads as
Δ θ k ( V ) = V ν k j ,
where ν k represents the kinetic coefficient and j is the constant derived from experimental data or numerical calculations.

3.1. Conductive Heat and Mass Transfer

Dealing with conductive heat and mass transfer near the dendritic surface, we come to the following thermal and concentration undercoolings [13,18]
Δ θ θ ( ρ , V ) = θ i θ = θ Q Iv θ * ( ρ , V ) , Δ θ σ ( ρ , V ) = m e ( σ i σ ) = m e σ ( 1 k e ) Iv σ * ( ρ , V ) 1 ( 1 k e ) Iv σ * ( ρ , V ) ,
where the modified (thermal and concentration) Ivantsov functions take the form
Iv θ * ( ρ , V ) = P θ ( ρ , V ) exp [ P 0 ( ρ , V ) ] I θ ( ) ,
Iv σ * ( ρ , V ) = P θ ( ρ , V ) D θ D σ exp P 0 ( ρ , V ) D θ D σ I σ ( ) .
Here, the thermal ( P θ ) and flow ( P u ) Péclet numbers are defined as
P θ ( ρ , V ) = ρ V 2 D θ , P u ( ρ ) = ρ u 2 D θ and P 0 ( ρ , V ) = P θ ( ρ , V ) + P u ( ρ ) .
The interfacial temperature θ i and solute concentration σ i , as well as the integrals I θ and I σ , are given by [13,18]
θ i = θ + θ Q P θ exp ( P 0 ) I θ ( ) ,
σ i = σ 1 ( 1 k e ) exp ( P 0 D θ / D σ ) P θ I σ ( ) D θ / D σ ,
I θ ( ζ ) = 1 ζ exp ( k 1 ) P u 1 ζ 1 G ( ζ 2 ) d ζ 2 ζ 2 P 0 ζ 1 d ζ 1 ζ 1 ( k 1 ) / 2 ,
I σ ( ζ ) = 1 ζ exp ( k 1 ) P u D θ D σ 1 ζ 1 G ( ζ 2 ) d ζ 2 ζ 2 P 0 D θ D σ ζ 1 d ζ 1 ζ 1 ( k 1 ) / 2 ,
where k = 2 and k = 3 in cases of the 2D and 3D dendritic geometry.
The function G ( ζ ) entering these expressions is determined from the hydrodynamic solution of the Oseen viscous flow model and is given by [22,23] in the 2D and 3D cases
G ( ζ ) = 2 π erfc / 2 exp 2 exp ζ 2 + erfc ζ / 2 erfc / 2 ζ in 2 D ,
G ( ζ ) = exp ( / 2 ) exp ( ζ / 2 ) ζ E 1 ( / 2 ) + ζ E 1 ( ζ / 2 ) 2 E 1 ( / 2 ) in 3 D ,
where = ρ u / μ stands for the Reynolds number, and E 1 ( w ) = w λ 1 exp ( λ ) d λ .
As a result, the contributions in the total undercooling balance (11) in the case of conductive heat and mass transfer mechanism are defined by expressions (12)–(14).

3.2. Convective Heat and Mass Transfer

Dealing with intense fluid motions and mixing near the dendrite surface, we come to the following temperature and concentration undercoolings for convective heat and mass transfer mechanism [18]
Δ θ θ ( V ) = θ i θ = θ Q V λ s α θ ρ m κ p u * D θ , Δ θ σ ( V ) = m e ( σ i σ ) = ( 1 k e ) V m e σ α σ u * ( 1 k e ) V .
The integrals I θ ( ζ ) and I σ ( ζ ) entering in (10) in the case of convective transfer mechanism are given by
I θ ( ζ ) = 1 ζ exp ( P θ ζ 1 ) ζ 1 ( k 1 ) / 2 d ζ 1 , I σ ( ζ ) = 1 ζ exp ( P σ ζ 1 ) ζ 1 ( k 1 ) / 2 d ζ 1 , P σ = ρ V 2 D σ ,
and, as before, k = 2 and k = 3 describe the 2D and 3D growth geometries. Thus, the contributions in the total undercooling balance (11) in the case of convective heat and mass transfer mechanism are defined by expressions (12), (13) and (15).

3.3. Sewing Together Undercooling Balances

In this section, we sew together the undercooling balances for conductive and convective heat and mass transfer to describe mixed transfer mechanisms near the dendritic surface. To do this, we introduce the corresponding superscripts “cond” and “conv”, i.e., we write the melt undercoolings as Δ θ c o n d for conductive mechanism and Δ θ c o n v for convective mechanism. To sew together these functions of Δ θ let us write down the following law
Δ θ = Δ θ s e w e d = Δ θ c o n v B c o n v ( Δ θ ) + Δ θ c o n d B c o n d ( Δ θ ) B c o n v ( Δ θ ) + B c o n d ( Δ θ ) ,
where B c o n d ( Δ θ ) and B c o n v ( Δ θ ) are the sewing functions. These functions should satisfy the following limiting cases: B c o n d ( Δ θ ) 0 for Δ θ 0 and B c o n v ( Δ θ ) 0 for Δ θ Δ θ c , where Δ θ c is a characteristic undercooling (e.g., initial undercooling). In other words, Δ θ s e w e d Δ θ c o n v for small and Δ θ s e w e d Δ θ c o n d for large melt undercoolings. The sewing functions can be chosen in various ways. So, for example, when describing dendritic tip shapes, the authors of Refs. [35,36] used the exponential sewing functions. On the other hand, when describing dynamics of intracellular clusters of nanoparticles, the authors of Ref. [37] used power functions to stitch together different asymptotic solutions.

4. Selection Criterion

Equation (16) expresses ρ , V and Δ θ . To obtain ρ ( Δ θ ) and V ( Δ θ ) independently, we need one more equation between these variables. Such an equation, known as the selection criterion, follows from the solvability theory and stability analysis [13,18,34,38,39,40].

4.1. Solvability Condition

The solvability condition defines the family of solutions for temperature and solute concentration distributions lying in the vicinity of parabolic/paraboloidal Ivantsov solutions. This condition takes the form [41]
A [ Ξ 0 ( ) ] m ( ) d = 0 , m ( ) = exp i 0 κ m ( 1 ) d 1 .
Here, A represents a curvature operator, κ m is the wavenumber found from stability theory (see Section 4.2), Ξ 0 ( ) is the solution of corresponding boundary-value problem, and i is the unit imaginary number. Equation (17) will be used below to obtain the second equation connecting ρ , V and Δ θ .

4.2. Conductive Heat and Mass Transfer

Now we consider the conductive-type boundary conditions at the growing dendritic surface. Following the theory [10], we expand the flow velocity components u ξ and u ζ in series at the point ζ = 1 for the 2D crystal. The result reads as
u ξ = ξ ξ + 1 1 / 2 V + b ( ) u ( ζ 1 ) , u ζ = V ( ξ + 1 ) 1 / 2 ,
where
b ( ) = 2 π 1 / 2 exp ( / 2 ) erfc ( / 2 ) in 2 D and b ( ) = exp ( / 2 ) E 1 ( / 2 ) in 3 D .
Again, following the theory [10], we use the local Cartesian reference frame ( x l , y l ) fixed at the dendrite surface and shifting together with the crystal. Keeping this in mind, we rewrite Formulas (18) as (here the axes x l and y l point out the tangent and normal directions)
u ˜ = y l b ( ) u 2 ρ sin ( 2 θ ) V sin θ , v ˜ = V cos θ ,
where u ˜ and v ˜ represent the tangent and normal velocities near dendritic surface.
The boundary conditions (5) determine the following derivatives at the solid/liquid surface y l = 0
d θ ˜ l d y l = θ Q V cos θ D θ , d σ ˜ d y l = ( k e 1 ) V σ i cos θ D σ , y l = 0 .
Here, and below the “tilde” above, the symbols indicate the steady-state solutions. Keeping this in mind, we have the temperature and solute concentration near the solid/liquid surface
θ ˜ l = θ i y l θ Q V cos θ D θ , σ ˜ = σ i + y l ( k e 1 ) V σ i cos θ D σ .
To find the wavenumber κ m in (17), we perturb the steady-state solutions of (i) fluid flow near the growing dendrite as u and v , (ii) temperature field in the solid and liquid phases as θ s and θ l , (iii) solute concentration in liquid as σ , (iv) hydrodynamic pressure as p , and (v) solid/liquid interface as Σ . All perturbations are assumed to be small compared to the corresponding steady-state functions. Equations governing a behaviour of perturbations follow from (1)–(3) and take the form
θ s τ + u ˜ θ s x l + v ˜ θ s y l + v d θ ˜ s d y l = D θ Δ θ s , θ l τ + u ˜ θ l x l + v ˜ θ l y l + v d θ ˜ l d y l = D θ Δ θ l , σ τ + u ˜ σ x l + v ˜ σ y l + v d σ ˜ d y l = D σ Δ σ , p = μ ρ m Δ v , · v = 0 .
Here, v means the vector of hydrodynamic perturbations.
Expanding expressions (5) and (7) in series near dendrite surface y l = 0 , we arrive at
θ s + θ l θ Q V cos θ D θ Σ = 0 , θ s θ Q d ( θ ) 2 Σ y l 2 + β ˜ ( θ ) Σ τ m e σ + m e σ i ( 1 k e ) V cos θ D σ Σ = 0 , u = 0 , v + Σ τ = 0 , θ Q D θ Σ τ θ s y l + θ l y l + θ Q V 2 cos 2 θ D θ 2 Σ = 0 .
Let us especially underline that all perturbations are proportional to
exp ϖ τ + i κ x l ε κ y l
within the framework of linear stability theory (see, among others, [10,13,18]). Here, ϖ and κ stand for the frequency and wavenumber of perturbations, and | ε | = 1 defines the convergence of perturbations for y l .
The method for constructing the analytical expressions for perturbations is described in detail in several papers devoted to the morphological stability theory (see, among others, [10,13,18]). Therefore, to save room, we mention only the main idea for finding such solutions. Considering all perturbations to be proportional to the exponential function (24), we find equations for their amplitudes from the boundary conditions (23). Then, to obtain a nontrivial solution for a system of linear equations for such amplitudes, we equate the determinant of this linear system to zero and come to the following cubic equation for the wavenumber κ = κ m
κ m 3 κ m V 2 D θ d ( θ ) + m e ( 1 k e ) σ i V θ Q D σ d ( θ ) exp ( i θ ) i κ m b ( ) u sin ( 2 θ ) 16 ρ D θ i κ m b ( ) u sin ( 2 θ ) 8 ρ D σ + i κ m 2 V sin θ 2 D θ V 2 cos θ exp ( i θ ) 4 D θ 2 d ( θ ) i κ m 2 β ˜ ( θ ) V sin θ θ Q d ( θ ) = 0 .
In deriving Equation (25), it was taken into account that the moving coordinate system associated with the dendrite tip is travelling in space at velocity i κ V sin θ . As this takes place, the perturbation evolves with the velocity ϖ ( κ ) . Therefore, the real velocity of a perturbation represents the sum of these velocities and looks like ϖ ( κ ) i V κ sin θ . Again, let us mention that the growth rate V of the dendritic tip is assumed to be constant throughout this paper. Here, we refer the interested reader to the theory [42] of how fast the dendrite reaches its steady-state velocity V. Furthermore, we have considered the neutral stability surface, where κ = κ m ( ϖ = 0 ), put ε = 1 , and changed i by i accordingly to the theory of Ref. [10]. An important point is that expression (25) corresponds to previously developed theories [8,9,10,13,41,43,44,45,46,47] in various limiting cases.
Next, combining (17) and (25), we get the selection criterion for conductive heat and mass transfer boundary conditions
σ * c o n d 2 d 0 D θ ρ 2 V 1 = σ 0 n χ d 7 / n B n 7 / n 1 + β ν ¯ n s n 1 1 + b 1 n χ d 2 / n B n 2 / n P θ 1 + Δ 0 β 0 D θ / d 0 2 + 2 m e ( 1 k e ) σ i D θ 1 + b 2 n χ d 2 / n B n 2 / n P σ 1 + Δ 0 β 0 D σ / d 0 C D 2 D σ θ Q ,
where
ν ¯ n = χ d 3 / n B n 3 / n b ( ) u d 0 4 V ρ R + b ( ) u d 0 D θ 2 V ρ R D σ , s n = 7 + n 2 ( 3 + n ) , R = 1 + 2 m e ( 1 k e ) σ i D θ D σ θ Q ,
d 0 C D = d 0 θ Q 2 m e ( 1 k e ) σ i , B n = 1 2 3 n / 4 j = 0 n n j i n j cos π ( n j ) 2 ,
b 1 n = 8 σ 0 n 7 1 / 2 3 56 3 / 8 B n 3 / 2 n χ d 3 / 2 n 3 / 8 , b 2 n = 2 b 1 n , Δ 0 = 20 7 24 1 / 2 56 3 3 / 8 .
Here, σ 0 n and β are the selection constants. The first of them defines the stable mode of dendritic growth in the absence of hydrodynamic flow, whereas the second one is responsible for hydrodynamic flows. The selection criterion (26) determines the second equation connecting ρ , V and Δ θ for conductive heat and mass transfer boundary conditions. Now, coupling (11) and (26), one can find the functions ρ ( Δ θ ) and V ( Δ θ ) . Note that criterion (26) has limiting transitions to previously developed theories [8,9,10,13,41,43,44,45,46,47]. The theory developed can be used for the three-dimensional dendritic growth as well. In this case, we have to replace the solution of the two-dimensional heat and mass transfer and hydrodynamic equations with the corresponding three-dimensional solutions. Therefore, the system (11) and (26) works in the 3D case too.

4.3. Convective Heat and Mass Transfer

Now, we consider the convective-type boundary conditions at the growing dendritic surface. In this case, the theory considered in Section 4.2 should be modified. At first, now expressions (19) look like u ˜ = V sin θ and v ˜ = V cos θ . The boundary conditions (23) become
θ l + ( f 1 + m e f 2 ) Σ + m e σ + d θ Q 2 Σ y l 2 β ˜ Σ τ = 0 , θ s m e f 2 Σ m e σ d θ Q 2 Σ y l 2 + β ˜ Σ τ = 0 , θ Q D θ Σ τ θ s y l + 2 h f 1 Σ + 2 h θ l , 1 k e α σ u * V cos θ ( σ + f 2 Σ ) + σ i Σ τ + σ + f 2 Σ = 0 ,
where h = α θ ρ m κ p u * / ( 2 λ s ) .
Equating to zero the determinant of coefficients at perturbations following from Equation (27), we obtain the marginal mode of the wavenumber
κ m 2 + κ m 2 h i β ˜ V sin θ θ Q d i Y 2 sin θ Y 1 d i 2 h β ˜ V sin θ θ Q d i V sin θ D θ d i 2 h Y 2 sin θ Y 1 d = 0 ,
where
Y 1 = 1 + ( 1 k e ) V cos θ α σ u * , Y 2 = ( 1 k e ) m e σ i V α σ θ Q u * .
Now, combining (17) and (28), we get the selection criterion for convective heat and mass transfer boundary conditions
σ * c o n v 2 d 0 D θ ρ 2 V 1 = 2 σ 0 n c χ d 2 / n B n 2 / n D θ β 1 ρ + σ 0 n c χ d 5 / n B n 5 / n 1 + b 1 β 1 D θ 1 + r n ν ˜ 1 n ( n + 5 ) / 2 ( n 1 ) 2 1 + l n χ d 3 / n B n 3 / n b 1 ρ + 3 χ d 1 / 4 B n 2 / n P θ D θ β 1 / ( 2 1 / 4 d 0 ) 2 ,
where
b 1 = α θ u * ρ m κ p 2 λ s , β 1 = β 0 + m e ( 1 k e ) σ i α σ u * θ Q , l n 2 = 2 9 / 2 25 27 σ 0 n c , ν ˜ 1 n = b 1 2 χ d 1 / n B n 1 / n ρ d 0 2 1 / 4 P θ ( 1 + b 1 β 1 D θ ) .
Here, σ 0 n c and r n represent selection constants dependent on the crystalline symmetry (subscript n).
The selection criterion (29) defines the second equation connecting ρ , V and Δ θ for convective heat and mass-transfer boundary conditions. Now, coupling (11) and (29), one can find the functions ρ ( Δ θ ) and V ( Δ θ ) . As before, the theory developed can be used for the three-dimensional dendritic growth as well.

4.4. Sewing Together Selection Criteria

Let us now sew together selection criteria (26) and (29) for conductive and convective heat and mass transfer near a dendrite surface
σ s e w e d * = σ * c o n v B c o n v ( Δ θ ) + σ * c o n d B c o n d ( Δ θ ) B c o n v ( Δ θ ) + B c o n d ( Δ θ ) .
As before, σ s e w e d * σ * c o n v for smaller and σ s e w e d * σ * c o n d for larger melt undercoolings.This means that criterion (30) contains limiting transitions to the conductive and convective mechanisms of heat and mass transfer. If both mechanisms take place simultaneously, we have two nonlinear Equations (16) and (30) to find ρ ( Δ θ ) and V ( Δ θ ) . An important point is that these equations should describe both undercooling limits: low Δ θ (convective mechanism prevails) and moderate Δ θ (conductive mechanism prevails). Note that a behaviour of ρ ( Δ θ ) and V ( Δ θ ) will substantially depend on the sewing functions B c o n d ( Δ θ ) and B c o n v ( Δ θ ) , which can be chosen in different ways [35,36,37,48,49]. The essential conditions are only those B c o n d ( Δ θ ) 0 for Δ θ 0 and B c o n v ( Δ θ ) 0 for Δ θ Δ θ c . In the next section, we analyse the sewed solutions (16) and (30) for different undercooled melts.

5. Behaviour of Sewed Functions

First of all, let us choose the stitching functions entering in the sewed undercooling balance (16) and selection criterion (30). Following previous studies [35,36,37], we chose power functions as follows
B c o n v ( Δ θ ) = ϵ c o n v θ Q Δ θ j c o n v , B c o n d ( Δ θ ) = ϵ c o n d Δ θ θ * j c o n d .
Here, ϵ c o n v and ϵ c o n d are the stitching constants, and j c o n v 1 and j c o n d 1 are the parameters determining the smoothness of stitching functions. These values should also provide the correct behaviour of B c o n v ( Δ θ ) and B c o n d ( Δ θ ) for small and large melt undercoolings. Namely, B c o n v ( Δ θ ) 0 at large Δ θ ( Δ θ θ Q ) and B c o n d ( Δ θ ) 0 at small Δ θ ( Δ θ θ * ). Taking these limiting cases into account, we see that
Δ θ Δ θ c o n v and σ * σ * c o n v , small undercooling ,
Δ θ Δ θ c o n d and σ * σ * c o n d , large undercooling .
Such behaviour provides a correct smooth transition of the sewed solution compared with purely convective and conductive branches shown in Figure 3 for the Al 24 Ge 76 melt. The solidification experiments with Al 24 Ge 76 alloy (see for details [15]) were carried out using a horizontal sample alignment, minimising gravity-induced fluid flow. The driving force for solidification is the total undercooling Δ θ at the dendrite tip. It is determined as the sum of solutal and curvature (Gibbs–Thomson effect) undercooling; kinetic and thermal undercooling were neglected (i.e., Δ θ θ = 0 and Δ θ k = 0 ), because fast thermal diffusion was assumed.
Comparing the theory with the experimental data, it is easy to see that each of these branches moves away from the experimental points into the region of its invalidity: the convective branch (solid line in Figure 3a) does not work for larger Δ θ , while the conductive one (dashed line in Figure 3a) does not describe experiments for smaller Δ θ . Both branches sewed together (solid line in Figure 3b) well satisfy experimental points in the whole diapason of melt undercooling.
To check the theory functionality for another system, we compare the sewed solution for the Ti 45 Al 55 droplets processed in electromagnetic levitators [13] in Figure 4. The sewed model is able to quantitatively reproduce the experimental results obtained by Hartmann et al. (see [16]). Note that that the intensity of convection and flow velocity are not measured in experiments of [16] directly. In this sense, the theoretical model uses the fluid velocity u * as the parameter, which relates to the process of viscous flow in levitated droplets realistically. Thus, convective boundary conditions well describe two experimental points lying near 40 K and conductive boundary conditions covering the intermediate region, as well as the whole undercooling diapason. As is easily seen, the sewed “tip velocity—melt undercooling” curve describes both limiting regimes with smaller and larger Δ θ , as well as the intermediate regime between them. In other words, the whole diapason of experimental measurements is covered.
The reason for this behaviour, from a practical viewpoint, is that vigorous currents of liquid could occur in the vicinity of growing dendrites. This in turn is associated with a laminar-to-turbulent transition when undercooled droplets solidify in electromagnetic levitators [50]. This transition may be driven as a result of vigorous fluid swirls, resulting in a convective heat and mass transport near the surfaces of crystals. Therefore, taking all of the above into account, we conclude that convective boundary conditions might outline experimental data on melt solidification for small undercooling by analogy with ice dendrites grown in sea water [24,25,26,27].

6. Conclusions

In summary, we developed a generalised theory of stable dendritic growth in undercooled binary melts in the presence of convective flows. First of all, we considered two kinds of solidification conditions: (i) forced melt convection when conductive-type heat and mass transfer boundary conditions take place, and (ii) intense melt convection (e.g., caused by turbulent flows near the dendrite surface or various peculiarities of hydrodynamics currents) when convective-type heat and mass transfer boundary conditions occur. The first of these cases describes moderate and high melt undercoolings Δ θ , whereas the second one corresponds to smaller values of the driving force Δ θ . Let us especially underline that small driving forces Δ θ happen at P θ 10 3 , moderate Δ θ occur at 10 3 P θ 10 2 , and high Δ θ correspond to P θ 10 2 [35]. Sewing together these cases, we described a wider range of melt undercoolings (both of these cases) until local-equilibrium crystallisation conditions are fulfilled. To do this, we derived undercooling balances and stability criteria for cases (i) and (ii) independently. Then these expressions were sewed together using the generalised Formulas (16) and (30), which represent the main outcome of our theory. Analysing these laws, we saw that the present solution contained two branches, convective and conductive. The first of them describes low melt undercoolings, while the second one is valid for moderate and high driving forces. The stitching functions B c o n v ( Δ θ ) and B c o n d ( Δ θ ) ensure a smooth transition between these branches. These functions can only be chosen by comparing theory and experimental data. Our simulations based on Formulas (16) and (30) demonstrated that the theory and experiments are in good agreement. Let us especially highlight in this conclusion that the stable growth mode of dendrites in local nonequilibrium crystallisation conditions P θ 1 can also be selected using the hyperbolic diffusion equation, modified mass balance condition, and nonequilibrium expressions for the segregation coefficient and liquidus line slope. Such a theory was developed in our previous studies [51,52].
Note that it is convenient to use the power laws “dendrite tip velocity—undercooling” and “dendrite tip diameter—undercooling” to compare the theory and experimental data. Such laws can be derived analytically in some limiting cases (see, among others, Ref. [46]). When studying different asymptotics to derive similar power laws for different undercoolings, Péclet and Reynolds numbers represent one of the important tasks for future investigations.
The theory under consideration could be further developed to describe the anomalous (U-shape) behaviour of the velocity–undercooling relationship found in experiments aboard the International Space Station [53]. To describe these experiments where dendritic growth and crystal nucleation occur simultaneously ahead of the recalescence front, we also need to use nucleation theory at initial and intermediate stages (see, among others, recently published papers [54,55,56,57]). This will be the subject of further investigations in the near future.

Author Contributions

Conceptualisation, D.V.A. and L.V.T.; methodology, D.V.A., P.K.G. and L.V.T.; software, L.V.T. and P.K.G.; validation, L.V.T. and P.K.G.; formal analysis, D.V.A. and P.K.G.; investigation, D.V.A.; resources, P.K.G.; writing—original draft preparation, D.V.A. and L.V.T.; writing—review and editing, D.V.A., L.V.T. and P.K.G.; visualisation, L.V.T.; supervision, P.K.G. and D.V.A.; project administration, P.K.G.; funding acquisition, P.K.G. and D.V.A. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Russian Science Foundation (project No. 21-19-00279).

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. Chernov, A.A. Modern Crystallography III: Crystal Growth; Springer: Berlin, Germany, 1984. [Google Scholar]
  2. Huang, S.-C.; Glicksman, M.E. Overview 12: Fundamentals of dendritic solidification—I. Steady-state tip growth. Acta Metall. 1981, 29, 701–715. [Google Scholar] [CrossRef]
  3. Trivedi, R.; Kurz, W. Dendritic growth. Int. Mater. Rev. 1994, 39, 49–74. [Google Scholar] [CrossRef]
  4. Galenko, P.K.; Zhuravlev, V.A. Physics of Dendrites; World Scientific: Singapore, 1994. [Google Scholar]
  5. Kurz, W.; Fisher, D.J. Fundamentals of Solidification; Trans Tech Publ.: Aedermannsdorf, Switzerland, 1989. [Google Scholar]
  6. Huppert, H.E. The fluid mechanics of solidification. J. Fluid Mech. 1990, 212, 209–240. [Google Scholar] [CrossRef] [Green Version]
  7. Herlach, D.; Galenko, P.; Holland-Moritz, D. Metastable Solids from Undercooled Melts; Elsevier: Amsterdam, The Netherlands, 2007. [Google Scholar]
  8. Pelcé, P. Dynamics of Curved Fronts; Academic Press: Boston, MA, USA, 1988. [Google Scholar]
  9. Ben Amar, M.; Pelcé, P. Impurity effect on dendritic growth. Phys. Rev. A 1989, 39, 4263–4269. [Google Scholar] [CrossRef] [PubMed]
  10. Bouissou, P.; Pelcé, P. Effect of a forced flow on dendritic growth. Phys. Rev. A 1989, 40, 6637–6680. [Google Scholar] [CrossRef]
  11. Brener, E.A.; Mel’nikov, V.I. Pattern selection in two-dimensional dendritic growth. Adv. Phys. 1991, 40, 53–97. [Google Scholar] [CrossRef]
  12. Alexandrov, D.V.; Galenko, P.K. Selected mode of dendritic growth with n-fold symmetry in the presence of a forced flow. EPL 2017, 119, 16001. [Google Scholar] [CrossRef]
  13. Alexandrov, D.V.; Galenko, P.K. Dendrite growth under forced convection: Analysis methods and experimental tests. Physics-Uspekhi 2014, 57, 771–786. [Google Scholar] [CrossRef]
  14. Alexandrov, D.V.; Galenko, P.K. A review on the theory of stable dendritic growth. Phil. Trans. R. Soc. A 2021, 379, 20200325. [Google Scholar] [CrossRef]
  15. Becker, M.; Klein, S.; Kargl, F. Free dendritic tip growth velocities measured in Al-Ge. Phys. Rev. Mat. 2018, 2, 073405. [Google Scholar] [CrossRef] [Green Version]
  16. Hartmann, H.; Galenko, P.K.; Holland-Moritz, D.; Kolbe, M.; Herlach, D.M.; Shuleshova, O. Non-equilibrium solidification in undercooled Ti45Al55 melts. J. Appl. Phys. 2008, 103, 073509. [Google Scholar] [CrossRef]
  17. Galenko, P.K.; Reuther, K.; Kazak, O.V.; Alexandrov, D.V.; Rettenmayr, M. Effect of convective transport on dendritic crystal growth from pure and alloy melts. Appl. Phys. Lett. 2017, 111, 031602. [Google Scholar] [CrossRef]
  18. Alexandrov, D.V.; Galenko, P.K.; Toropova, L.V. Thermo-solutal and kinetic modes of stable dendritic growth with different symmetries of crystalline anisotropy in the presence of convection. Phil. Trans. R. Soc. A 2018, 376, 20170215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Lamb, H. Hydrodynamics; Dover Publications: New York, NY, USA, 1945. [Google Scholar]
  20. Kochin, N.E.; Kibel, I.A.; Roze, N.V. Theoretical Hydromechanics; Interscience: New York, NY, USA, 1964. [Google Scholar]
  21. Buyevich, Y.A.; Alexandrov, D.V.; Zakharov, S.V. Hydrodynamics. Examples and Problems; Begell House: New York, NY, USA, 2001. [Google Scholar]
  22. Alexandrov, D.V.; Galenko, P.K. Dendritic growth in an inclined viscous flow. Part 1. Hydrodynamic solutions. AIP Conf. Proc. 2017, 1906, 200003. [Google Scholar]
  23. Alexandrov, D.V.; Galenko, P.K. Dendritic growth in an inclined viscous flow. Part 2. Numerical examples. Aip Conf. Proc. 2017, 1906, 200004. [Google Scholar]
  24. McPhee, M.G.; Maykut, G.A.; Morison, J.H. Dynamics and thermodynamics of the ice/upper ocean system in the marginal ice zone of the Greenland sea. J. Geophys. Res. 1987, 92, 7017–7031. [Google Scholar] [CrossRef]
  25. Notz, D.; McPhee, M.G.; Worster, M.G.; Maykut, G.A.; Schlünzen, K.H.; Eicken, H. Impact of underwater-ice evolution on Arctic summer sea ice. J. Geophys. Res. 2003, 108, 3223. [Google Scholar] [CrossRef]
  26. Alexandrov, D.V.; Nizovtseva, I.G. Nonlinear dynamics of the false bottom during seawater freezing. Dokl. Earth Sci. 2008, 419, 359–362. [Google Scholar] [CrossRef]
  27. Alexandrov, D.V.; Nizovtseva, I.G. To the theory of underwater ice evolution, or nonlinear dynamics of ‘false bottoms’. Int. J. Heat Mass Trans. 2008, 51, 5204–5208. [Google Scholar] [CrossRef]
  28. Alexandrov, D.V.; Malygin, A.P. Convective instability of directional crystallization in a forced flow: The role of brine channels in a mushy layer on nonlinear dynamics of binary systems. Int. J. Heat Mass Trans. 2011, 54, 1144–1149. [Google Scholar] [CrossRef]
  29. Tritton, D.J. Physical Fluid Dynamics; Clarendon Press: Oxford, UK, 1988. [Google Scholar]
  30. Feltham, D.L.; Worster, M.G.; Wettlaufer, J.S. The influence of ocean flow on newly forming sea ice. J. Geophys. Res. 2002, 107, 3009. [Google Scholar] [CrossRef] [Green Version]
  31. Owen, P.R.; Thomson, W.R. Heat transfer across rough surfaces. J. Fluid Mech. 1963, 15, 321–334. [Google Scholar] [CrossRef]
  32. Yaglom, A.M.; Kader, B.A. Heat and mass transfer between a rough wall and turbulent flow at high Reynolds and Peclet numbers. J. Fluid Mech. 1974, 62, 601–623. [Google Scholar] [CrossRef]
  33. 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. (Met.) 2019, 2019, 787–794. [Google Scholar] [CrossRef]
  34. Barbieri, A.; Langer, J.S. Predictions of dendritic growth rates in the linearized solvability theory. Phys. Rev. A 1989, 39, 5314–5325. [Google Scholar] [CrossRef]
  35. Alexandrov, D.V.; Galenko, P.K. The shape of dendritic tips. Phil. Trans. R. Soc. A 2020, 378, 20190243. [Google Scholar] [CrossRef] [Green Version]
  36. Alexandrov, D.V.; Toropova, L.V.; Titova, E.A.; Kao, A.; Demange, G.; Galenko, P.K.; Rettenmayr, M. The shape of dendritic tips: A test of theory with computations and experiments. Phil. Trans. R. Soc. A 2021, 379, 20200326. [Google Scholar] [CrossRef]
  37. Alexandrov, D.V.; Korabel, N.; Currell, F.; Fedotov, S. Dynamics of intracellular clusters of nanoparticles. Cancer Nanotechnol. 2022, 13, 15. [Google Scholar] [CrossRef]
  38. Ben Amar, M. Theory of needle-crystal. Phys. D 1988, 31, 409–423. [Google Scholar] [CrossRef]
  39. Brener, E.A.; Mel’nikov, V.I. Two-dimensional dendritic growth at arbitrary Peclet number. J. Phys. France 1990, 51, 157–166. [Google Scholar] [CrossRef]
  40. Brener, E.A. Effects of surface energy and kinetics on the growth of needle-like dendrites. J. Cryst. Growth 1990, 99, 165–170. [Google Scholar] [CrossRef]
  41. Pelcé, P.; Bensimon, D. Theory of dendrite dynamics. Nucl. Phys. B 1987, 2, 259–270. [Google Scholar] [CrossRef]
  42. Titova, E.A.; Galenko, P.K.; Alexandrov, D.V. Method of evaluation for the non-stationary period of primary dendritic crystallization. J. Phys. Chem. Solids 2019, 134, 176–181. [Google Scholar] [CrossRef]
  43. Langer, J.S.; Hong, D.C. Solvability conditions for dendritic growth in the boundary-layer model with capillary anisotropy. Phys. Rev. A 1986, 34, 1462–1471. [Google Scholar] [CrossRef]
  44. Alexandrov, D.V.; Galenko, P.K. Selection criterion of stable dendritic growth at arbitrary Péclet numbers with convection. Phys. Rev. E 2013, 87, 062403. [Google Scholar] [CrossRef] [Green Version]
  45. Alexandrov, D.V.; Galenko, P.K. Thermo-solutal and kinetic regimes of an anisotropic dendrite growing under forced convective flow. Phys. Chem. Chem. Phys. 2015, 17, 19149–19161. [Google Scholar] [CrossRef]
  46. Alexandrov, D.V.; Galenko, P.K. Dendritic growth with the six-fold symmetry: Theoretical predictions and experimental verification. J. Phys. Chem. Solids 2017, 108, 98–103. [Google Scholar] [CrossRef]
  47. Alexandrov, D.V.; Galenko, P.K. Thermo-solutal growth of an anisotropic dendrite with six-fold symmetry. J. Phys. Condens. Matter 2018, 30, 105702. [Google Scholar] [CrossRef]
  48. Toropova, L.V. Shape functions for dendrite tips of SCN and Si. Eur. Phys. J. Spec. Top. 2022, 231, 1129–1133. [Google Scholar] [CrossRef]
  49. 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. Condens. Matter Phys. 2022, 34, 094002. [Google Scholar] [CrossRef]
  50. Hyers, R.W.; Trapaga, G.; Abedian, B. Laminar-turbulent transition in an electromagnetically levitated droplet. Metall. Mater. Trans. B 2003, 34B, 29–36. [Google Scholar] [CrossRef]
  51. Alexandrov, D.V.; Galenko, P.K. Selection criterion of stable mode of dendritic growth with n-fold symmetry at arbitrary Péclet numbers with a forced convection. In Proceedings of the IUTAM Symposium on Recent Advances in Moving Boundary Problems in Mechanics, Christchurch, New Zealand, 12–15 February 2019; pp. 203–215. [Google Scholar]
  52. Alexandrov, D.V.; Galenko, P.K. Selected mode for rapidly growing needle-like dendrite controlled by heat and mass transport. Acta Mater. 2017, 137, 64–70. [Google Scholar] [CrossRef]
  53. Reinartz, M.; Kolbe, M.; Herlach, D.M.; Rettenmayr, M.; Toropova, L.V.; Alexandrov, D.V.; Galenko, P.K. Study on anomalous rapid solidification of Al-35 at%Ni in microgravity. JOM 2022, 74, 2420–2427. [Google Scholar] [CrossRef]
  54. 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. (Met.) 2017, 2017, 127–135. [Google Scholar] [CrossRef]
  55. 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]
  56. 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]
  57. Toropova, L.V.; Alexandrov, D.V. Dynamical law of the phase interface motion in the presence of crystals nucleation. Sci. Rep. 2022, 12, 10997. [Google Scholar] [CrossRef]
Figure 1. A three-dimensional sketch of dendritic crystal. (a) Laminar flow with conductive boundary conditions, (b) local convection in front of dendrite tip with convective boundary conditions, where u * represents friction velocity appearing as a result of intensive heat and mass transfer.
Figure 1. A three-dimensional sketch of dendritic crystal. (a) Laminar flow with conductive boundary conditions, (b) local convection in front of dendrite tip with convective boundary conditions, where u * represents friction velocity appearing as a result of intensive heat and mass transfer.
Crystals 12 00965 g001
Figure 2. A sketch of dendritic crystal streamlined by fluid flow in 2D (panel a, parabolic cylinder coordinates ζ and ξ ) and 3D (panel b, paraboloid of revolution coordinates ζ , ξ and ϕ ). The phase transition boundary corresponds to ζ = 1 in panels a and b.
Figure 2. A sketch of dendritic crystal streamlined by fluid flow in 2D (panel a, parabolic cylinder coordinates ζ and ξ ) and 3D (panel b, paraboloid of revolution coordinates ζ , ξ and ϕ ). The phase transition boundary corresponds to ζ = 1 in panels a and b.
Crystals 12 00965 g002
Figure 3. Tip velocity V as a function of the melt undercooling Δ θ for Al 24 Ge 76 [15] undercooled melt. (a) Convective (solid line) and conductive (dashed line) boundary conditions. Material and calculation parameters are given in Table 1. (b) Sewed model. Stitching parameters: ϵ c o n v = 4.5 × 10 10 , ϵ c o n d = 0.48 , j c o n v = j c o n d = 2 .
Figure 3. Tip velocity V as a function of the melt undercooling Δ θ for Al 24 Ge 76 [15] undercooled melt. (a) Convective (solid line) and conductive (dashed line) boundary conditions. Material and calculation parameters are given in Table 1. (b) Sewed model. Stitching parameters: ϵ c o n v = 4.5 × 10 10 , ϵ c o n d = 0.48 , j c o n v = j c o n d = 2 .
Crystals 12 00965 g003
Figure 4. Tip velocity V as a function of the melt undercooling Δ θ for Ti 45 Al 55 [16] undercooled melt. Sewed model, material and calculation parameters are given in Table 1. Stitching parameters: ϵ c o n v = 3 × 10 9 , ϵ c o n d = 0.03 , j c o n v = j c o n d = 3 .
Figure 4. Tip velocity V as a function of the melt undercooling Δ θ for Ti 45 Al 55 [16] undercooled melt. Sewed model, material and calculation parameters are given in Table 1. Stitching parameters: ϵ c o n v = 3 × 10 9 , ϵ c o n d = 0.03 , j c o n v = j c o n d = 3 .
Crystals 12 00965 g004
Table 1. Material and calculation parameters for the Al 24 Ge 76 [15] and Ti 45 Al 55 [16] melts.
Table 1. Material and calculation parameters for the Al 24 Ge 76 [15] and Ti 45 Al 55 [16] melts.
ParameterAl 24 Ge 76 Ti 45 Al 55 Units
Liquidus slope, m e 10.410.72 K / at . %
Hypercooling, θ Q 353272K
Liquidus temperature, θ * 7321748K
Solute diffusion coefficient, D σ 7 × 10 9 8 × 10 9 m 2 s 1
Initial composition, σ 2455 at . %
Capillary constant, d 0 4.3 × 10 10 7.8 × 10 10 m
Thermal diffusivity, D θ 4 × 10 5 7.5 × 10 6 m 2 s 1
Liquid density, ρ m 1.26 × 10 3 2.46 × 10 3 kg m 3
Heat capacity, κ p 5501237 J kg 1 K 1
Thermal conductivity in the solid, λ s 29.2229.22 W m 1 K 1
Friction velocity of flow, u * 3.3 × 10 4 2 m s 1
Segregation coefficient, k e 0.110.86-
Surface energy stiffness, χ d 0.0260.030-
Solvability constant, σ 0 n / σ 0 n c 0.09/0.090.02/0.15-
Convective coefficient of heat, α θ 00.25-
Convective coefficient of mass, α σ 2.881-
Order of crystalline symmetry, n44-
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.; Galenko, P.K.; Alexandrov, D.V. A Stable Mode of Dendritic Growth in Cases of Conductive and Convective Heat and Mass Transfer. Crystals 2022, 12, 965. https://doi.org/10.3390/cryst12070965

AMA Style

Toropova LV, Galenko PK, Alexandrov DV. A Stable Mode of Dendritic Growth in Cases of Conductive and Convective Heat and Mass Transfer. Crystals. 2022; 12(7):965. https://doi.org/10.3390/cryst12070965

Chicago/Turabian Style

Toropova, Liubov V., Peter K. Galenko, and Dmitri V. Alexandrov. 2022. "A Stable Mode of Dendritic Growth in Cases of Conductive and Convective Heat and Mass Transfer" Crystals 12, no. 7: 965. https://doi.org/10.3390/cryst12070965

APA Style

Toropova, L. V., Galenko, P. K., & Alexandrov, D. V. (2022). A Stable Mode of Dendritic Growth in Cases of Conductive and Convective Heat and Mass Transfer. Crystals, 12(7), 965. https://doi.org/10.3390/cryst12070965

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