# A Simplified Scaling Law of Cell-Dendrite Transition in Directional Solidification

Authors: Yaochan Zhu; Hua Qiu; Zhijun Wang; Eckart Schnack
Journal: Advances in Condensed Matter Physics
 (2019)

---

## Abstract

To describe the cell-dendrite transition (CDT) during directional solidification, a new simplified scaling law is proposed and verified by quantitative phase field simulations. This scaling law bears clear physical foundation with consideration of the overall effects of primary spacing, pulling velocity, and thermal gradient on the onset of sidebranches. The analysis results show that the exponent parameters in this simplified scaling law vary within different systems, which mediates the discrepancy of exponent parameters in previous experiments. The scaling law also presents an explanation for the destabilizing mechanism of thermal gradient in sidebranching dynamics.


---

## Body

## 1. Introduction

As one of the typical problems of pattern formation, solidification microstructures are interesting for both scientific and technologic reasons [1–3]. Though remarkable progress in the study of microstructure evolution during solidification has been made, the mechanism of sidebranching dynamics is still unsatisfactorily understood [4–6]. During directional solidification, the onset of sidebranches means the occurrence of cell-dendrite transition (CDT). Over past few years, the identification of CDT has been widely performed experimentally [7–11], which showed that cell and dendrite coexist over a large range of control parameters. Given pulling velocity and thermal gradient, small primary spacing corresponds to cellular morphology in CDT region, while large primary spacing corresponds to dendrite. Accordingly, CDT significantly depends on the primary cellular/dendritic spacing, and a critical primary spacing λ
1
C
D corresponding to CDT should exist. In terms of the correlation between λ
1
C
D and control parameters, researchers have put much effort all along on finding the criteria of CDT by considering the spacing, pulling velocity V, and thermal gradient G. Based on abundant experiments [7–10], an empirical scaling law with the form of λ
1
C
D ∞ V
α
G
β has been summarized to describe the critical primary spacing λ
1
C
D in which α and β are exponent parameters.Although the form of scaling law has been proposed, some controversies still exist. On one hand, the values of the exponentα and β are inconsistent in different experiments. α
=
-
1
/
2 and β
=
-
1
/
8 were found in Gerogelin et al.’s experiments [7] and α
=
β
=
-
1
/
3 in Trivedi et al.’s experiments [8, 9]. On the other hand, this empirical scaling law was based on the data fitting without further physical foundation. Researchers have tried to expound the intrinsic physical foundation of the scaling law. In Trivedi et al.’s work [8], the critical primary spacing λ
1
C
D was given by the geometrical meaning of three characteristic lengths: solutal diffusion length l
D, thermal length l
T, and capillary length d
0. Gerogelin et al. [7] also presented a self-similar asymptotic regime about l
T
/
l
D, λ
1
C
D
/
d
0, and l
T
/
d
0 based on their experimental data. However, these analyses only focused on the assembly of characteristic lengths, not referring to the sidebranching dynamics, which characterizes the CDT. The exponent parameters selection and the physical foundation in the scaling law of CDT are still unclear and need further exploration.The crossover of CDT is usually defined by the occurring of sidebranches. Accordingly, sidebranching should be one of the typical characteristics in CDT, and the connections between the scaling law and sidebranching dynamics should exist. However, previous scaling law did not take sidebranching dynamics into account. Therefore, it should be more reasonable to characterize CDT by the onset of sidebranching instability. Furthermore, the empirical scaling law indicates that, beyond CDT, increasing the thermal gradientG will enhance the sidebranching dynamics. This destabilizing effect of thermal gradient G on the sidebranching dynamics has been observed in [7–10]. Noise amplification theory failed in describing the effect of thermal gradient G on sidebranching dynamics [7]. Therefore, to acquire a deep understanding of CDT and sidebranching, it is essential to propose more reasonable physical explanations on the scaling law and the destabilizing mechanism.Although sidebranching dynamics has received considerable attentions in free dendritic growth [1, 4–6], only some basic understanding about the sidebranching dynamics has been obtained in directional solidification [7, 10–15]. Gerogelin et al [7] presented a model to describe the sidebranching dynamics, in which the noise amplitude at the cellular tip was controlled by the feedback of sidebranches. However, it is difficult to determine the growth factor in their model. On the other hand, experimental results and phase field simulations indicated that the diffusion instability of dendritic trunk is the most possible reason for sidebranching dynamics [14–16] and revealed that the initial sidebranching spacing only depends on the pulling velocity, but the sidebranch amplitude is determined not only by pulling velocity V but also by the primary spacing and the thermal gradient G. With consideration of the overall effects of primary spacing, pulling velocity V, and thermal gradient G on the onset of sidebranches, the scaling law of CDT may gain more physical foundation.In this article, firstly, we briefly review the factors in determining the onset of sidebranches from previous experiments and quantitative phase field simulations. Then a new scaling law of CDT is proposed based on the sidebranching dynamics and the physical foundation of this new scaling law is presented. The effects of pulling velocity and thermal gradient on CDT will be analyzed according to this new scaling law. Finally, the quantitative phase field simulations will be used to testify this proposed scaling law.


## 2. Factors Determining Sidebranches during Directional Solidification

It has been widely accepted that the primary spacingλ
1, pulling velocity V, and thermal gradient G play an important role in determining the sidebranches dynamics. The effects of λ
1, V, and G on the dendritic growth have also been well studied. Here the related results and how the sidebranching dynamics is affected by the controlled parameters are recalled briefly.Primary spacingλ
1 affects the sidebranch dynamics significantly. In experiments, the microstructures near the CDT show dendritic array with large spacing and cellular array with small spacing [7–10]. Quantitative phase field simulations revealed the intrinsic reason [15, 17] that the smaller spacing suppresses the sidebranch growth due to the strong interdendritic solutal interaction. However, after the appearance of sidebranch, the location of first sidebranch and the initial sidebranch spacing are almost independent of primary spacing. Therefore, the primary spacing only influences the amplitude of the sidebranch but does not change the initial sidebranch spacing. Accordingly, the primary spacing only supplies spacing for sidebranch growth.Pulling velocityV is an important control parameter during directional solidification. According to experimental investigations and simulation, the initial sidebranch spacing λ
2 has a scaling law with the pulling velocity λ
2
∞
V
α [7, 14, 15, 18]. Dendrite trunk also depends on the pulling velocity greatly. The width of dendrite trunk decreases as pulling velocity increases. Within the same primary spacing, the decrease of pulling velocity enlarges the diffusion length, which enhances the interdendritic interaction and suppresses the sidebranch growth. As aforementioned, the interdendritic interaction does not change the initial sidebranch spacing. Therefore, the variation of initial sidebranch spacing is mainly attributed to the variation of pulling velocity. In previous investigations, the variation of initial sidebranch spacing λ
2 with pulling velocity satisfies λ
2
∞
V
-
0.59 [7, 14, 18].The role of thermal gradientG in sidebranch growth is a little bit complex. Previous investigation on the cell-dendrite transition indicated that positive thermal gradient promotes the generation of sidebranch [7], but the initial sidebranch spacing λ
2 was independent of thermal gradient G [14]. It is attributed to the remarkable interdendritic solutal interaction near CDT, where the thermal gradient G significantly affects the mush zone in directional solidification. However, the details on the thermal gradient effects are still absent and some confusion still exists. For example, the effects of thermal gradient on dendritic morphologies are relatively weaker for dendritic array growth, where the sidebranch amplitude is almost independent of thermal gradient [15].


## 3. Scaling Law of CDT from Sidebranching

From the above review, we can found that previous investigations have presented lots of details about sidebranching dynamics during directional solidification However, the intrinsic sidebranching mechanism during directional solidification is still not fully revealed. Compared with free growth, the interdendritic interaction will suppress the sidebranch growth, so the interdendrite interaction plays an important role in determining the sidebranches. It is still difficult to reveal CDT by directly deriving the sidebranching amplitude evolution from the basic diffusion equation and interface condition.Here we focus on the interdendritic solutal interaction to reveal the CDT and derive the scaling law of CDT in the following. Sketch of the dendrite with sidebranches and definition of parameters during the derivations are shown in Figure1, where ρ
t
i
p is the dendritic tip radius, z
t
i
p is the tip position along z-axis, λ
1 is primary spacing, λ
2 is the sidebranch spacing, and Λ
1 and Λ
2 are groove width. The parameter definitions are as follows: solutal diffusion length l
D
=
D
/
V, thermal length l
T
=
∆
T
0
/
G, and capillary length d
0
=
Γ
/
∆
T
0, where D is the liquid solutal diffusion coefficient, ∆
T
0
=
m
C
0
(
k
-
1
)
/
k, Γ is the Gibbs-Thomson coefficient [3], C
0 is the initial alloys composition, m is the liquid slope, and k is the equilibrium solute-partition coefficient. In the interdendrite region, we define the groove width 2
Λ and lateral diffusion length l
L
D
=
D
/
V
L, where V
L is the lateral interface velocity.Figure 1
Sketch of sidebranches and parameters definition in the derivation of scaling law, whereλ
1 is primary spacing, λ
2 is the sidebranch spacing, Λ
1 and Λ
2 are groove width, and z
t
i
p is the tip position along z-axis.Here, the CDT is characterized by the appearance of sidebranches. Sidebranches firstly appear near the dendritic tip and then are amplified in the groove by the lateral growth of dendritic trunk. Two assumptions are adopted in the derivation:(
1
) the occurrence of CDT corresponds with the visibility of sidebranches; (
2
) the intrinsic conditions for sidebranching are based on the fact that the lateral diffusion accelerates the lateral instability, while the groove width blocks the lateral instability. According to assumption (
1
), the ratio of the sidebranch amplitude A and the sidebranch wavelength λ
2 can be adopted to characterize the CDT. According to assumption (
2
), CDT is determined by the competition between the lateral diffusion length l
L
D and the groove width 2
Λ. Both l
L
D and 2
Λ have the length dimension, so the groove instability can be assumed to be proportional to the groove width 2
Λ and inversely proportional to the lateral diffusion length l
L
D. Therefore, the two dimensionless characteristics A
/
λ
2 and Λ
/
l
L
D can be the representative of dimensionless characterization and driving force of sidebranching dynamics, respectively. For smaller sidebranch amplitude, A
/
λ
2 is a function of Λ
/
l
L
D:(1)
A
λ
2
=
f
Λ
l
L
DNear CDT, the lateral diffusion length is very large and the groove width is small, which results in that theΛ
/
l
L
D approaches zero. With Taylor expansion of A
/
λ
2
≈
f
(
0
)
+
f
′
(
0
)
Λ
/
l
L
D
+
o
Λ
/
l
L
D
2, we can assume that the relationship between the driving force and the characterization of the sidebranch is linear near CDT, where both A
/
λ
2 and Λ
/
l
L
D are very small. Then the function in (1) can be approximated by a linear function: A
/
λ
2
~
Λ
/
l
L
D or A
~
λ
2
Λ
/
l
L
D. If we define certain finite amplitude A as the representation of sidebranches appearance, then this finite amplitude A leads to a constant λ
2
Λ
/
l
L
D at CDT for different pulling velocities and thermal gradients.According to the asymptotic analysis of dendrite growth in directional solidification by Spencer and Huppert [19], in the region of ρ
t
i
p
≪
z
t
i
p
-
z
≪
l
T as shown in Figure 1, we have(2)
Λ
=
λ
1
1
-
a
-
z
t
i
p
-
z
/
l
T
2where a is a modified parameter to represent the dendritic trunk width in the groove. Then, within a period of sidebranches, the lateral interface velocity V
L and the lateral diffusion length l
L
D can be represented as(3)
V
L
=
Λ
1
-
Λ
2
λ
2
V
(4)
l
L
D
∝
4
l
D
l
T
z
t
i
p
-
z
λ
1In experiments, the visible amplitude is most likely to appear atz
=
z
t
i
p
-
n
λ
2, where n is a constant [7, 8]. According to (1)–(4), a new simplified scaling law of CDT can be described as(5)
λ
1
C
D
~
λ
2
-
1
/
4
l
D
1
/
2
l
T
1
/
4
1
-
a
-
n
λ
2
1
/
2
l
T
-
1
/
2
-
1
/
2In formula (5), sidebranch spacing λ
2 and solutal diffusion length l
D
=
D
/
V are related to the pulling velocity Vwhile thermal length l
T
=
∆
T
0
/
G is related to the thermal gradient G. The first three terms on the right side of formula (5) are directly connected to the exponent parameters in the scaling law of λ
1
C
D
∞
V
α
G
β, while the contribution of the last term 1
-
α
-
n
λ
2
1
/
2
l
T
-
1
/
2
-
1
/
2 to the exponent parameters depends on λ
2
/
l
T. According to this new scaling law, the effects of pulling velocity V and thermal gradient G on the CDT can be well addressed.To validate this new scaling law, the data in [12] is adopted as a paradigm for further analysis. To describe the effect of thermal length l
T more conveniently, we define(6)
f
l
T
=
l
T
1
/
4
1
-
a
-
n
λ
2
1
/
2
l
T
-
1
/
2
-
1
/
2According to the data in [12], the relationship between sidebranch spacing λ
2 and pulling velocity V is λ
2 ≈ 350V
−0.6, the thermal length l
T is about 500~1400μm, and the parameters of a and n in (6) are 0.1 and 2, respectively.In the above new scaling law described by formula (5), the exponent parameter α for pulling velocity V is not constant but is related to the thermal gradient G. Figure 2 shows the effect of l
T on the exponent parameter α. When l
T increases from 500μm to 1400μm, α increases from -0.49 to -0.41 correspondingly. It shows that, within the range of l
T = 500~1400μm, the value of α consists with the scaling law proposed by Gerogelin et al. (α
=
-
1
/
2) [7], where α
=
-
0.46 in uniform cellular/dendritic array. When λ
2
/
l
T is small enough, we have λ
1
C
D
~
λ
2
-
1
/
4
l
D
1
/
2
l
T
1
/
4
~
V
-
0.35; then α will be -0.35, which agrees with the result of Trivedi et al.’s experiments [8, 9].Figure 2
The scaling law between the critical primary spacing and the pulling velocity with different thermal lengths in the cell-dendrite transition, described by formula (5), where λ
2 ≈ 350V
−0.6, a=0.1, and n=2.As to the exponent parameterβ for thermal gradient G in the scaling law, Figure 3 shows the variation of f
(
l
T
) with l
T. When l
T is small, the effect of the term 1
-
a
-
n
λ
2
1
/
2
l
T
-
1
/
2
-
1
/
2 cannot be overlooked. As shown in the inset, the power function fitting gives the exponent β as –0.122, which is consistent with the fitting results in Gerogelin et al.’s experiments (β
=
-
1
/
8) [7]. As the thermal gradient G decreases, l
T increases and the variation of term 1
-
a
-
n
λ
2
1
/
2
l
T
-
1
/
2
-
1
/
2 makes the exponent parameter β increase. As shown in Figure 3, large l
T ensures β
=
-
1
/
4 (λ
1
C
D
~
λ
2
-
1
/
4
l
D
1
/
2
l
T
1
/
4
~
G
-
1
/
4). Therefore, for a small thermal gradient corresponding to large l
T, the exponent parameter of the thermal gradient in the scaling law is very close to Trivedi et al.’s results [8, 9].Figure 3
The effect of the thermal length on the scaling law of formula (5). The inset presents the local region with power law fitting. The corresponding parameters are λ
2 ≈ 350V
−0.6 with V=50μm/s, a=0.1, and n=2.The form of this scaling law can return back to the empirical one and reconciles the difference of the exponent parameters in previous experiments. Analysis on the exponent parameters in the scaling law ofλ
1
C
D
∞
V
α
G
β indicates the variation of exponents α and β with different solidification systems. When λ
2
/
l
T
→
0, λ
1
C
D
~
λ
2
-
1
/
4
l
D
1
/
2
l
T
1
/
4
~
V
-
0.35
G
-
1
/
4, that is, α=-0.35 and β=-0.25, which is close to Trivedi et al.’s experimental results [8, 9]. Previous experiment also mentioned that the change of thermal gradient G influences the exponent parameters. Teng et al. [9] pointed out that the increase of thermal gradient G induces a larger systematic deviation from α
=
β
=
-
1
/
3. Similarly, in Gerogelin et al.’s experiments [7], the thermal gradient G was relatively large and α
=
-
0.5 and β
=
-
1
/
8. This agrees with the analysis presented here, which shows that the exponent parameters vary with thermal gradient G. For a large thermal gradient G, the exponent parameter may be selected as α
=
-
0.45 and β
=
-
0.122. Accordingly, the power law presented here could settle the argument in previous experiments about the discrepancy of exponent parameters.As to the destabilizing effect of thermal gradientG on the sidebranching dynamics, it also can be explained according to this new scaling law. The sidebranching instability is determined by the competition between lateral diffusion length l
L
D and the groove width 2
Λ. Both of l
L
D and Λ decrease with the increase of thermal gradient G. Small diffusion length implies the enhanced diffusion instability, while small groove width is helpful to stabilize the cell. However, from (2), (3), and (4), we can get Λ
/
l
L
D
~
1
-
a
/
l
T
z
t
i
p
-
z
-
1, so Λ
/
l
L
D decreases with the thermal length l
T; that is, Λ
/
l
L
D increases with the thermal gradient G. This indicates that the effect of thermal gradient G on the lateral diffusion length l
L
D predominates. The destabilizing effect of thermal gradient G comes from the decrease of lateral diffusion length l
L
D, which is similar to the destabilizing effect of increasing pulling velocity V in planar instability [1, 20, 21]. So the lateral diffusion length l
L
D bridges the thermal gradient G and the sidebranching instability, which can be used to explain the destabilizing nature of thermal gradient.The scaling law here is related to three parameters, the sidebranch spacing, the diffusion length, and the thermal length, which are related to the undercooling and thermal gradient, respectively. In experiments, all these parameters vary in different systems with different parameters. Accordingly, the validity of the scaling law can be well checked with designed experiments. In this research, the new scaling law will be validated by a benchmark from quantitative phase field simulation instead.Note that the noise effects on the sidebranching dynamics in dendrite growth have been a controversial issue for many years [22]. However, there is an agreement that the presence of stochastic noise will not affect the frequency of sidebranching, while the value of the amplitude is a function of noise intensity level. Just as predicted by [23], the transition of cellular to dendrite growth may be modified by the noise intensity level. In the derivation of the scaling law, the sidebranching amplitude is considered as a function of three length scales that are independent of the noise. Therefore, the noise may affect the transition points, but the scaling law between these lengths still exists. At an adequate noise level similar to that in the experiment, the scaling law is valid.


## 4. Validation of the Scaling Law by Phase Field Simulation

The development of phase field method makes it possible to quantitatively investigate microstructure evolution in solidification [18] and the quantitative phase field simulation has been widely used to investigate the sidebranching dynamics in crystal growth [6, 15, 16] and the primary spacing selection mechanism in directional solidification [17, 24–26]. To further validate the scaling law, the quantitative phase field method [18] is employed. For the directional solidification simulation by (7) and (8), assume one-sided diffusion and frozen temperature approximation, in which T
=
T
C
+
G
z and T
C is the temperature at the cooling end and G is the thermal gradient along z-axis. The dynamic evolution equations of phase field model in the moving frame with pulling velocity V are (7)
τ
0
1
-
1
-
k
z
+
m
c
∞
/
k
/
G
l
T
∂
t
ϕ
-
V
∂
z
ϕ
=
∇
W
θ
2
∇
ϕ
-
∂
x
W
θ
W
′
θ
∂
y
ϕ
+
∂
y
W
θ
W
′
θ
∂
x
ϕ
+
ϕ
-
ϕ
3
+
λ
1
-
ϕ
2
2
U
+
z
+
m
c
∞
/
k
/
G
l
T
(8)
∂
t
c
-
V
∂
z
c
=
∇
·
D
L
1
-
ϕ
1
+
k
-
1
-
k
ϕ
∇
c
+
D
L
1
-
ϕ
1
+
k
-
1
-
k
ϕ
+
W
0
2
∂
t
ϕ
-
V
∂
z
ϕ
∇
ϕ
c
1
-
k
1
+
k
-
1
-
k
ϕ
∇
ϕwith(9)
W
θ
=
W
0
1
+
γ
4
cos
⁡
4
θ
(10)
U
=
2
k
c
/
c
∞
/
1
+
k
-
1
-
k
ϕ
-
1
1
-
k
(11)
l
T
=
m
c
∞
1
-
k
k
Gwhere W
0 are the parameters of the interface thickness, τ
0 is the relaxation time for phase field model, and λ is the coupling constant, which are related to physical quantities by d
0
=
a
1
W
0
/
λ and τ
0
=
a
2
λ
W
0
2
/
D; γ
4 is the anisotropic intensity of the surface tension, θ is the angle between the normal vector of the interface and the preferred orientation, m is the liquidus slope,k is the partition coefficient, and c
∞ is the concentration in the far away field. Here d
0
=
Γ
/
(
m
c
∞
(
1
-
k
)
/
k
), a
1=0.8839, and a
2=0.6267, and Γ is the Gibbs-Thomson coefficient.The transparent alloy SCN-0.43wt%C152 is adopted, which has been widely used to investigate the evolution of dendritic pattern [27]. The chemical diffusion coefficient of the liquid phase is D ≈ 0.45×10−9
m
2/s, partition coefficient k=0.05, Gibbs-Thomson coefficient Γ = 6.48×10−8 K/m, and the slope of liquidus line m=-542 K/mol. The surface tension anisotropy intensity is assumed as γ
4=0.005.This quantitative phase field simulation on dendritic growth with the presented parameters has been widely performed in our previous investigations [15, 17]. The phase field simulation on the sidebranching indicates that λ
2=150V
−0.59 and the parameters a and n in (6) are about 0.4 and 5, respectively. Then the exponents in the new proposed scaling law for the simulation system are α
=
-
0.35 and β
=
-
0.21 according to formula (5), which are different from α
=
β
=
-
1
/
3 [8, 9] and α
=
-
1
/
2 and β
=
-
1
/
8 [7]. These three different groups of exponents will be quantitatively examined according to the phase field simulation results.Here, a benchmark is designed to directly compare the exponents in different scaling laws. The exponents are usually obtained by fitting plenty of simulation or experiment results with different control parameters. However, heavy workload and artificial judgment on the onset of sidebranching are required by this method. In a different way, we design a benchmark instead of finding the exponents in the simulation system as follows. The steady state of interface morphology in specific primary spacing with onset of sidebranching is firstly presented, where the thermal gradientG
0 and pulling velocity V
0 as well as the morphology will be the references. Then, the critical pulling velocity for the CDT can be extrapolated from different scaling laws along with the variation of thermal gradient G in the fixed primary spacing. The solid/liquid morphologies corresponding to the different critical velocities are obtained by phase field simulation. Finally, by comparing the simulation results with the reference morphology, the scaling law with different exponents is evaluated. Here G
0=20.2 K/cm and V
0=20μm/s are chosen as the referential control parameters. The critical primary spacing for CDT is λ
1=160μm and the referential morphology is shown in Figure 4(a). Withλ
1=160μm and G
0=5.2 K/cm, the different critical pulling velocities for CDT are V
I=77μm/s, V
I
I=45μm/s, and V
I
I
I=28μm/s, respectively, for three groups of exponent parameters (-1/3, -1/3), (-0.35, -0.21), and (-1/2, -1/8), as shown by the dot in Figure 4(b).Figure 4
The critical microstructure of cell-dendritic transition at the benchmark (a) and the criteria of cell-dendrite transition according to the different power laws based on the benchmark (b).
(a)
(b)Figure5 presents the cellular/dendrite morphologies for different pulling velocities with λ
1=160μm and G
0=5.2 K/cm. Compared with the referential morphology in Figure 4(a), the critical pulling velocity is around 35μm/s. The pulling velocity 77μm/s predicted by (-1/3, -1/3) is obviously larger than the critical pulling velocity for CDT, while the pulling velocity 28μm/s predicted by (-1/2, -1/8) is close to the critical pulling velocity but with cellular morphology. The pulling velocity 45μm/s from the new proposed scaling is also close to the critical pulling velocity of CDT.Figure 5
The cellular/dendritic morphologies with different pulling velocities whenG=5.2 K/cm.The critical pulling velocity of CDT can be further found within higher accuracy. In the simplified form of the new proposed scaling law (formula (5)), the contribution of the term λ
2
/
l
T on the critical pulling velocity in (6) is overlooked. Here, by submitting λ
2=150V
−0.59 into (6), the effect of the term λ
2
/
l
T can be revealed. Considering the overall effects of thermal gradient G and pulling velocity V, the variation of the left-hand side of formula (5) with pulling velocity for two different thermal gradients is presented in Figure 6. It shows that the exponent α is near -0.35 with small thermal gradient G. However, α deviates from -0.35 gradually as the thermal gradient G increases, which means that the influence of pulling velocity V on λ
2
/
l
T has significant impact on f
(
l
T
) with relatively small l
T. Considering the contribution of pulling velocity V on f
(
l
T
), Figure 6 gives the critical pulling velocity of CDT as 37μm/s when λ
1=160μm and G
0=5.2 K/cm in the simulation system. The interface morphology for 37μm/s is very close to the reference morphology in the benchmark.Figure 6
The variation of the value of right-hand side in (5) with pulling velocity for different thermal gradients.


## 5. Results and Discussion

To summarize, a simplified scaling law of CDT during directional solidification is derived with considering the sidebranching dynamics. The exponent parameters corresponding to the pulling velocity and thermal gradient in the new scaling law are discussed. The analysis shows that the exponent parameters in the scaling law vary with different solidification systems and reconcile the discrepancy in previous experimental results. The form of this scaling law can return back to the empirical one and reconciles the difference of the exponent parameters in previous experiments. The destabilizing mechanism of thermal gradient in the sidebranching dynamics can be revealed by lateral diffusion length. The new scaling law is also validated by a benchmark from quantitative phase field simulation. The appropriate experimental verification of the scaling law can be similar to that done by Teng et al [9]. With the apparatus, one can use different systems and parameters to check the cellular-to-dendrite transition.Furthermore, the proposed scaling law is more than the conciliation of the controversy in previous experiments. On one hand, compared with previous scaling law, the new scaling law is with more physical foundation related to the sidebranching dynamics. On the other hand, it indicates that the thermal gradient and pulling velocity are coupled together in describing the CDT within a large range parameter space. Only in local parameter space, the scaling law of CDT has a simple form consisted with previous experimental results where the thermal gradient and pulling velocity are decoupled. Therefore, the proposed scaling law is with more precision in predicting CDT in a large range parameter space compared with previous investigations.



---

*Source: 8767640-2019-06-02.xml*