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.

1. Introduction

As one of the typical problems of pattern formation, solidification microstructures are interesting for both scientific and technologic reasons [13]. Though remarkable progress in the study of microstructure evolution during solidification has been made, the mechanism of sidebranching dynamics is still unsatisfactorily understood [46]. 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 [711], 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 corresponding to CDT should exist. In terms of the correlation between and control parameters, researchers have put much effort all along on finding the criteria of CDT by considering the spacing, pulling velocity , and thermal gradient . Based on abundant experiments [710], an empirical scaling law with the form of has been summarized to describe the critical primary spacing 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. and were found in Gerogelin et al.’s experiments [7] and 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 was given by the geometrical meaning of three characteristic lengths: solutal diffusion length , thermal length , and capillary length . Gerogelin et al. [7] also presented a self-similar asymptotic regime about , , and 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 gradient will enhance the sidebranching dynamics. This destabilizing effect of thermal gradient on the sidebranching dynamics has been observed in [710]. Noise amplification theory failed in describing the effect of thermal gradient 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, 46], only some basic understanding about the sidebranching dynamics has been obtained in directional solidification [7, 1015]. 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 [1416] and revealed that the initial sidebranching spacing only depends on the pulling velocity, but the sidebranch amplitude is determined not only by pulling velocity but also by the primary spacing and the thermal gradient . With consideration of the overall effects of primary spacing, pulling velocity , and thermal gradient 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 , pulling velocity , and thermal gradient play an important role in determining the sidebranches dynamics. The effects of , , and 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 affects the sidebranch dynamics significantly. In experiments, the microstructures near the CDT show dendritic array with large spacing and cellular array with small spacing [710]. 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 velocity is an important control parameter during directional solidification. According to experimental investigations and simulation, the initial sidebranch spacing has a scaling law with the pulling velocity [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 with pulling velocity satisfies [7, 14, 18].

The role of thermal gradient 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 was independent of thermal gradient [14]. It is attributed to the remarkable interdendritic solutal interaction near CDT, where the thermal gradient 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 Figure 1, where is the dendritic tip radius, is the tip position along -axis, is primary spacing, is the sidebranch spacing, and and are groove width. The parameter definitions are as follows: solutal diffusion length , thermal length , and capillary length , where is the liquid solutal diffusion coefficient, , Γ is the Gibbs-Thomson coefficient [3], is the initial alloys composition, is the liquid slope, and is the equilibrium solute-partition coefficient. In the interdendrite region, we define the groove width and lateral diffusion length , where is the lateral interface velocity.

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: the occurrence of CDT corresponds with the visibility of sidebranches; 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 , the ratio of the sidebranch amplitude and the sidebranch wavelength can be adopted to characterize the CDT. According to assumption , CDT is determined by the competition between the lateral diffusion length and the groove width . Both and have the length dimension, so the groove instability can be assumed to be proportional to the groove width and inversely proportional to the lateral diffusion length . Therefore, the two dimensionless characteristics and can be the representative of dimensionless characterization and driving force of sidebranching dynamics, respectively. For smaller sidebranch amplitude, is a function of :

Near CDT, the lateral diffusion length is very large and the groove width is small, which results in that the approaches zero. With Taylor expansion of , we can assume that the relationship between the driving force and the characterization of the sidebranch is linear near CDT, where both and are very small. Then the function in (1) can be approximated by a linear function: or . If we define certain finite amplitude as the representation of sidebranches appearance, then this finite amplitude leads to a constant 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 as shown in Figure 1, we havewhere is a modified parameter to represent the dendritic trunk width in the groove. Then, within a period of sidebranches, the lateral interface velocity and the lateral diffusion length can be represented as

In experiments, the visible amplitude is most likely to appear at , where is a constant [7, 8]. According to (1)–(4), a new simplified scaling law of CDT can be described as

In formula (5), sidebranch spacing and solutal diffusion length are related to the pulling velocity while thermal length is related to the thermal gradient . The first three terms on the right side of formula (5) are directly connected to the exponent parameters in the scaling law of , while the contribution of the last term to the exponent parameters depends on . According to this new scaling law, the effects of pulling velocity and thermal gradient 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 more conveniently, we define

According to the data in [12], the relationship between sidebranch spacing and pulling velocity is ≈ 350V−0.6, the thermal length is about 500~1400 μm, and the parameters of and in (6) are 0.1 and 2, respectively.

In the above new scaling law described by formula (5), the exponent parameter for pulling velocity is not constant but is related to the thermal gradient . Figure 2 shows the effect of on the exponent parameter . When increases from 500 μm to 1400 μm, increases from -0.49 to -0.41 correspondingly. It shows that, within the range of = 500~1400 μm, the value of consists with the scaling law proposed by Gerogelin et al. () [7], where in uniform cellular/dendritic array. When is small enough, we have ; then will be -0.35, which agrees with the result of Trivedi et al.’s experiments [8, 9].

As to the exponent parameter for thermal gradient in the scaling law, Figure 3 shows the variation of with . When is small, the effect of the term 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 () [7]. As the thermal gradient decreases, increases and the variation of term makes the exponent parameter increase. As shown in Figure 3, large ensures (). Therefore, for a small thermal gradient corresponding to large , the exponent parameter of the thermal gradient in the scaling law is very close to Trivedi et al.’s results [8, 9].

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 indicates the variation of exponents and with different solidification systems. When , , 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 influences the exponent parameters. Teng et al. [9] pointed out that the increase of thermal gradient induces a larger systematic deviation from . Similarly, in Gerogelin et al.’s experiments [7], the thermal gradient was relatively large and and . This agrees with the analysis presented here, which shows that the exponent parameters vary with thermal gradient . For a large thermal gradient , the exponent parameter may be selected as and . 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 gradient 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 and the groove width . Both of and decrease with the increase of thermal gradient . 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 , so decreases with the thermal length ; that is, increases with the thermal gradient . This indicates that the effect of thermal gradient on the lateral diffusion length predominates. The destabilizing effect of thermal gradient comes from the decrease of lateral diffusion length , which is similar to the destabilizing effect of increasing pulling velocity in planar instability [1, 20, 21]. So the lateral diffusion length bridges the thermal gradient 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, 2426]. 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 and is the temperature at the cooling end and is the thermal gradient along -axis. The dynamic evolution equations of phase field model in the moving frame with pulling velocity are withwhere are the parameters of the interface thickness, is the relaxation time for phase field model, and is the coupling constant, which are related to physical quantities by and ; is the anisotropic intensity of the surface tension, is the angle between the normal vector of the interface and the preferred orientation, is the liquidus slope, k is the partition coefficient, and is the concentration in the far away field. Here , =0.8839, and =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 ≈ 0.45×10−9 m2/s, partition coefficient =0.05, Gibbs-Thomson coefficient = 6.48×10−8 K/m, and the slope of liquidus line =-542 K/mol. The surface tension anisotropy intensity is assumed as =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 =150V−0.59 and the parameters and in (6) are about 0.4 and 5, respectively. Then the exponents in the new proposed scaling law for the simulation system are and according to formula (5), which are different from [8, 9] and and [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 gradient and pulling velocity 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 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 =20.2 K/cm and =20 μm/s are chosen as the referential control parameters. The critical primary spacing for CDT is =160 μm and the referential morphology is shown in Figure 4(a). With λ1=160 μm and =5.2 K/cm, the different critical pulling velocities for CDT are =77 μm/s, =45 μm/s, and =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 5 presents the cellular/dendrite morphologies for different pulling velocities with =160 μm and =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.

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 on the critical pulling velocity in (6) is overlooked. Here, by submitting =150V−0.59 into (6), the effect of the term can be revealed. Considering the overall effects of thermal gradient and pulling velocity , 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 . However, deviates from -0.35 gradually as the thermal gradient increases, which means that the influence of pulling velocity on has significant impact on with relatively small . Considering the contribution of pulling velocity on , Figure 6 gives the critical pulling velocity of CDT as 37μm/s when =160 μm and =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.

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.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The work was supported by the fund of the State Key Laboratory of Solidification Processing in NWPU (Grant no. SKLSP201725).