Numerical Simulation of Microstructure Evolution in Solidification Process of Ferritic Stainless Steel with Cellular Automaton

In order to study the basic principles of vibration-excited liquid metal nucleation technology, a coupled model to connect the temperature field calculated by ANSYS Fluent and the dendritic growth simulated by cellular automaton (CA) algorithm was proposed. A two-dimensional CA model for dendrite growth controlled by solute diffusion and local curvature effects with random zigzag capture rule was developed. The proposed model was applied to simulate the temporal evolution of solidification microstructures under different degrees of surface undercooling and vibration frequency of the crystal nucleus generator conditions. The simulation results showed that the predicted columnar dendrites regions were more developed, the ratio of interior equiaxed dendrite reduced and the size of dendrites increased with the increase of the surface undercooling degrees on the crystal nucleus generator. It was caused by a large temperature gradient formed in the melt. The columnar-to-equiaxed transition (CET) was promoted, and the refined grains and homogenized microstructure were also achieved at the high vibration frequency of the crystal nucleus generator. The influences of the different process parameters on the temperature gradient and cooling rates in the mushy zone were investigated in detail. A lower cooling intensity and a uniform temperature gradient distribution could promote nucleation and refine grains. The present research has guiding significance for the process parameter selection in the actual experimental.


Introduction
Ferritic stainless steel is easy to form developed columnar dendrites during solidification [1]; in the rolling process, it is easy to form a strong texture and this is the main reason for wrinkles, cracks, and other defects [2], which restrict its yield and wide application. Increasing the ratio of the equiaxed dendrite is an effective method to improve its forming defects. A new vibration-excited liquid metal nuclear technique has been proposed to improve the solidification structure of ferritic stainless steel, and it has been proved that this technology is beneficial to increase the equiaxed grains [3][4][5]. At present, it is necessary to reveal the quantitative relationship between the process variables and the solidification microstructure, but it is difficult to visually observe the dendrite growth process in the molten metal by conventional experiment. Numerical simulations were developed to understand the dendritic growth and microstructure formation in solidification in the last twenty years [6][7][8][9]. The cellular automaton (CA) method has been rapidly developed as an attractive and efficient technique to predict the evolution of solidification structures [10][11][12] and study the mechanical properties of composite in the field of materials science [13,14].
Zhu and Hong [15,16] first developed the modified cellular automaton model (MCA) and studied dendritic growth at the microscopic scale. The model was able to reproduce

Nucleation Model
The continuous nucleation model proposed by Thevoz [43] was used to describe the heterogeneous nucleation process in the melt. The relationship between the increase in nucleus density dn and the undercooling d(∆T) in the Thevoz model follows the Gaussian distribution, (1) where ∆T is the nucleation undercooling, ∆T n and ∆T σ are the average value and standard deviation of the nucleation undercooling, n max is the maximum density of nuclei, T liq is liquidus temperature, m l is liquidus slope, T * is the local temperature, and C l is the local actual liquid concentration. The nucleation undercooling of the bulk liquid that is considered the phenomenon of solute enrichment will alter the actual liquidus temperature; thus, the change in the solute concentration affects the degree of undercooling and feeds back to the nucleation and growth processes. During one time step, δt, the undercooling of a given CA cell increases by an amount δ(∆T). Accordingly, the density of new grains, δn, which are nucleated within the volume of the melt is given by where V CA is the volume of the CA cell, ∆x and ∆x 2 (µm and µm 2 ) corresponding to the mold wall and in the bulk liquid, respectively, and P n is nucleation probability in the cell. During the time step, nucleation will be activated in a liquid cell that needs to meet two conditions-first, the total undercooling in the cell is greater than the critical nucleation undercooling, and second, the nucleation probability is greater than a random number that ranges from 0 to 1. After nucleation, the cell becomes active and grows with a preferential direction corresponding to the randomly crystallographic orientation.

Solute Field Calculation
During the solidification of alloys, the dendritic growth is mainly controlled by solute transport. Without considering the natural and forced convection in the liquid, the governing equation for diffusion in the domain can be described as follows: where C l and C s are the solute concentration in liquid and solid phases, respectively, D l is the solute diffusion coefficient in liquid phases, and k e is partition coefficient. The third term on the right-hand side is included only at the interface cell and denotes the amount of solute rejected at the interface. Additionally, the solute partition between liquid and solid at solid-liquid interface (SL interface) is calculated by Equation (6). A central finite difference method was used to calculate the liquid concentration of cells as described by Equation (5). The method [21] of setting the equivalent composition C e was used to calculate the local liquid composition in the SL interface at each time interval.

Interface Growth Kinetics
The interface growth kinetics used in the present study was proposed by Zhu and Stefanescu [21]. The fundamental concept of the dendritic growth kinetics was considered to be driven by the difference between the local equilibrium composition and the local actual liquid composition. According to the thermodynamic concept of local equilibrium between the liquid and solid phases, the interface equilibrium composition C * l and the increase in the solid fraction ∆F s of an interface cell can be calculated by where T liq is the equilibrium liquidus temperature at the initial composition, T * is the local temperature in the dendrite tip, C * l is the interface equilibrium composition, C 0 is the initial equilibrium concentration and m l is liquidus slope, Γ is the Gibbs-Thomson coefficient, K is the curvature of the SL interface, and f (ϕ, θ) is a function describing the anisotropy of the interface energy, where θ is crystallographic orientation, ϕ is the angle of the normal to the interface with the coordinate axis, ε is the interface anisotropy intensity parameter, which is 0.04 in present work. The term ΓK f (ϕ, θ) in Equation (7) indicates the effect of interfacial energy on the process of dendrite growth; only by accurately calculating the curvature of the interface, the effect of the interface energy can be characterized.
For each interface cell with a non-zero solid fraction, the compositions of liquid and solid are calculated and stored simultaneously. The solid and liquid composition is calculated as follows: where C new l and C new s are the liquid and solid composition in the interface cell at the current moment. Similarly, C old l and C old s represent the liquid and solid composition at the previous moment, respectively. F s is the solid fraction inside the cell. When the interface cell is completely solidified (i.e., F s = 1), the excess solute components are distributed to eight neighbor liquid cells according to the concentration differences, which aim at maintaining the solute conservation conditions of the system.

Calculation of the SL Interface Curvature and the Capture Rules
According to Wei's research [39,44], the accurate result of the interface curvature played a critical role in the CA algorithm to correct the anisotropy of the grid. The lack of accuracy in curvature calculation significantly influenced the accuracy of the CA model. The conventional curvature calculation method includes the counting cell method and the variation of the unit vector normal (VUVN) to the solid-liquid interface. The calculation method of curvature and interface growth angle is as follows: However, the Equations (12) and (13) needs to calculate the derivatives of solid fractions, which is difficult to be accurately calculated in a single interface cell. Thus, the capture rules in this article adopt the random zigzag capture rule and calculated the SL interface curvature by the bilinear interpolation VUVN method, which aims to avoid the mesh anisotropy introduced by Von Neumann and Moore's capture rule; the calculation details can be found in the reference [37]. Figure 1 presents the basic principles of the vibration-excited liquid metal nucleation technique. The main process [4,42] can be summarized as follows: a cylindrical bar with a cooling structure and high-frequency vibration was inserted into the superheated liquid metal, and the distribution of temperature and flow field in the melt correspondingly changed significantly. The function of the crystal nucleus generator could be divided into two important parts. It could maintain a large degree of undercooling between the nucleus generator surface and liquid metal, which could change the temperature gradient distribution and cooling intensity in the melt, as is shown in Figure 1a. Moreover, the violent vibration effect triggered a convection stirring process in the melt, which enhanced the convection capacity in the liquid melt, as is shown in Figure 1b. As a result, the effect of fluid flow could create uniform temperature gradient distribution and break off the primary dendrites arm, and the fragments of dendrites were drawn into the melt, which increased the potential heterogeneous nucleation site. However, the Equations (12) and (13) needs to calculate the derivatives of solid fractions, which is difficult to be accurately calculated in a single interface cell. Thus, the capture rules in this article adopt the random zigzag capture rule and calculated the SL interface curvature by the bilinear interpolation VUVN method, which aims to avoid the mesh anisotropy introduced by Von Neumann and Moore's capture rule; the calculation details can be found in the reference [37]. Figure 1 presents the basic principles of the vibration-excited liquid metal nucleation technique. The main process [4,42] can be summarized as follows: a cylindrical bar with a cooling structure and high-frequency vibration was inserted into the superheated liquid metal, and the distribution of temperature and flow field in the melt correspondingly changed significantly. The function of the crystal nucleus generator could be divided into two important parts. It could maintain a large degree of undercooling between the nucleus generator surface and liquid metal, which could change the temperature gradient distribution and cooling intensity in the melt, as is shown in Figure 1a. Moreover, the violent vibration effect triggered a convection stirring process in the melt, which enhanced the convection capacity in the liquid melt, as is shown in Figure 1b. As a result, the effect of fluid flow could create uniform temperature gradient distribution and break off the primary dendrites arm, and the fragments of dendrites were drawn into the melt, which increased the potential heterogeneous nucleation site.  The present dendrite growth algorithm does not consider the remelting effect, the breaking process of dendrite fragments, and the influence of the melt convection on the quantity and position of crystal nuclei because these required duplicated model construction and huge computational efficiency. Future research can learn from the methods in computational fluid dynamics [45,46] to optimize and expand the model. The contribution of dendrite fragmentation due to the vibration-induced convection was incorporated into the nucleation model by increasing the probability of nucleation event in liquid. Based on the previous experimental research, the nucleation probability P n in 25% increments as a simulation parameter, when the vibration frequency increases. The thermal properties of the ferritic stainless steel Cr17 used in the present simulation are listed in Table 1 [47,48]. Table 1. The thermal properties of the ferritic stainless steel Cr17 used in the simulation.

Definition and Symbols Values
Liquidus temperature, Figure 2 shows the process of a single equiaxed dendrite growing at the undercooled melt with different undercooling. A square calculation domain was divided into 1000 × 1000 cells with a cell size of 0.5 µm. The whole calculation area was maintained at a constant undercooling degree and the cell state was initialized to liquid phase with the solute concentration of C 0 . A round nucleus with the composition k·C 0 and the radius of 10 grid size was placed at the center of the domain. Zero-flux solute boundary conditions were applied; the effect of convection on solute transport was ignored in the model.

Equiaxed Dendritic Growth in Undercooled Melt
Comparing the two sets of pictures, it can be derived that the dendrite exhibits the branchless needle morphology at lower undercooling; on the contrary, the dendrites growth is fast with developed branches at larger undercooling. Figure 3 shows the characteristics of the solute distribution and the growth rate of the dendrite tip obtained by simulation at the conditions of ∆T = 5 K and mesh size ∆x = 0.5 µm. Figure 3a reveals that the excess solute is released by the dendrite arm and accumulates at the surrounding liquid phase at the initial solidification stage. There is an exponential drop in solute concentration at the SL interface, which represents the solute enrichment region and corresponds to the red envelope, shown in Figure 2. As the growth continues, the liquid composition at the front of the interface gradually stabilizes, as shown by the curves at 0.063 s and 0.072 s in Figure 3b. This indicates that the liquid composition at the front of the interface has reached the equilibrium liquid composition, and dendrites enter steady-state growth stage.

The Evolution of Microstructure
A computational domain with 1.6 × 4.0 mm 2 was meshed into CA cells with a size of ∆x = 2 μm and was initially full of the solute concentration of C 0 at the liquidus temperature. The transient thermal conditions of the microscopic CA grid are interpolated from the result calculated by ANSYS Fluent 17.0, and the dendritic evolution under the vibration-excited liquid metal nucleation technology is predicted with the 2D parallel CA-FD model, as shown in Figure 4.
Nucleation occurs first on the chilling surface of the nucleus generator. Columnar dendrites grow along the direction of the highest temperature gradient until the equiaxed grains are formed in the undercooled melt when the degree of local liquid undercooling is greater than the critical nucleation undercooling ∆ . The nucleation and growth processe of equiaxed grains are determined by the local thermal and solute fields according to Equation (2). Finally, the growth of columnar dendrites is hindered by equiaxed crystals, forming a tendency for columnar dendrites and equiaxed dendrites to compete with each other. Figure 5 depicts the grain structures obtained by simulation result and experiment, respectively. It can be seen that the grain morphology simulated by the CA model is in good agreement with the experimental result. The average grain size of equiaxed grain obtained by simulation and experiment are 103 μm and 123 μm, respectively. Compared with the simulation results, the average grain size distribution obtained in the experiment is not uniform.    Figure 4.
Nucleation occurs first on the chilling surface of the nucleus generator. Columnar dendrites grow along the direction of the highest temperature gradient until the equiaxed grains are formed in the undercooled melt when the degree of local liquid undercooling is greater than the critical nucleation undercooling ∆T N . The nucleation and growth processe of equiaxed grains are determined by the local thermal and solute fields according to Equation (2). Finally, the growth of columnar dendrites is hindered by equiaxed crystals, forming a tendency for columnar dendrites and equiaxed dendrites to compete with each other.
In the following part, the location of columnar-to-equiaxed transition (CET), the numbers of newly formed dendrites, and the average grain size are analyzed quantitatively. The effects of the undercooling degree and the vibration frequency on the thermal distribution conditions and microstructure evolution are discussed in detail.

Influence of the Surface Undercooling Degree
In this set of simulations, the vibration frequency of the crystal nucleus generator was set to a constant value of 1000 Hz, and the structure characteristic formed in the melt under different degrees of undercooling on the surface of the crystal nucleus generator was explored. Figure 6 shows the dendritic structures under the different undercooling conditions with 200-700 K, which means ∆T = T melt -T surf , where T melt and T surf is the temperature of the melt and the temperature on the crystal nucleus generator surface at the initial stage, respectively. When the undercooling increases from 300 to 700 K, the columnar dendrites regions are more developed, but the area of the equiaxed dendrites regions decreases significantly, the number of equiaxed dendrites decreases, and the average grain size increases. The undercooling on the surface of the nucleus generator directly influenced the thermal gradient and cooling intensity in the melt. In order to analyze the influence mechanism of the nucleus generator surface undercooling, the cooling rate and temperature gradient distribution curves of the two points A and B (as illustrated in Figure 4a) were extracted. The temperature field characteristics can be clearly observed in Figure 7.
At the initial stage, the temperature gradient between A and B is relatively large and gradually decreases to a steady value. When the degree of undercooling increases, the temperature gradient distribution in the melt is significantly improved; the value of the stable temperature gradient is 1.73, 2.45, 2.89 K/mm, which can be observed from Figure  7a. The temperature gradient at the front of the solid-liquid interface directly affects the solidification manner of the dendrites. Figure 8 shows the solute field during the dendrite growth process when the undercooling is 300 K and describes the solidification process in the microscopic area. When the temperature gradient at melt is uniform, a wider mush zone is formed, creating a condition for the interaction between dendrite growth and solute pile-up, and hence the wider mushy zone increases the nucleation probability of equiaxed dendrites. In the following part, the location of columnar-to-equiaxed transition (CET), the numbers of newly formed dendrites, and the average grain size are analyzed quantitatively. The effects of the undercooling degree and the vibration frequency on the thermal distribution conditions and microstructure evolution are discussed in detail.

Influence of the Surface Undercooling Degree
In this set of simulations, the vibration frequency of the crystal nucleus generator was set to a constant value of 1000 Hz, and the structure characteristic formed in the melt under different degrees of undercooling on the surface of the crystal nucleus generator was explored. Figure 6 shows the dendritic structures under the different undercooling conditions with 200-700 K, which means ∆T = T melt − T sur f , where T melt and T sur f is the temperature of the melt and the temperature on the crystal nucleus generator surface at the initial stage, respectively. When the undercooling increases from 300 to 700 K, the columnar dendrites regions are more developed, but the area of the equiaxed dendrites regions decreases significantly, the number of equiaxed dendrites decreases, and the average grain size increases. The undercooling on the surface of the nucleus generator directly influenced the thermal gradient and cooling intensity in the melt. In order to analyze the influence mechanism of the nucleus generator surface undercooling, the cooling rate and temperature gradient distribution curves of the two points A and B (as illustrated in Figure 4a) were extracted. The temperature field characteristics can be clearly observed in Figure 7. Crystals 2021, 11, x FOR PEER REVIEW 12 of 20  At the initial stage, the temperature gradient between A and B is relatively large and gradually decreases to a steady value. When the degree of undercooling increases, the temperature gradient distribution in the melt is significantly improved; the value of the stable temperature gradient is 1.73, 2.45, 2.89 K/mm, which can be observed from  Figure 9 is a solute field diagram during dendrite growth when the undercooling degree is 500 K. The temperature gradient in the melt is relatively large at this condition, which is conducive to the directional growth of columnar dendrites. In addition, a narrow mush zone is formed and the nucleation of equiaxed dendrites appears in the region B between the columnar dendrites firstly, which can be clearly seen in Figure 8a. From the solidification process described in Figure 8, it can be seen that the equiaxed dendrites nucleation area is local and small, which exhibited as the layer-by-layer advancement of equiaxed dendrites. The regions A and C with lower solute enrichment evolve into the regions of equiaxed dendrites nucleation and growth, as shown in regions A' and C' (as illustrated in Figure 9a Figure 8 shows the solute field during the dendrite growth process when the undercooling is 300 K and describes the solidification process in the microscopic area. When the temperature gradient at melt is uniform, a wider mush zone is formed, creating a condition for the interaction between dendrite growth and solute pile-up, and hence the wider mushy zone increases the nucleation probability of equiaxed dendrites. The temperature gradient distribution in the melt is uniform, and the process showing volumetric solidification trend, the growth of columnar dendrites is blocked by the nucleation of equiaxed crystals, which corresponds to the situation in Figure 6a-c. According to Figure 8, it can be estimated that the enriched solute near the columnar dendrite's tip can affects the nucleation undercooling distribution. The area with low solute enrichment is more conducive to nucleation and growth. For example, in Figure 8a, the solute concentration in regions A, B, C, and D are relatively low; therefore, the nucleation and growth of equiaxed dendrites preferentially appear in A', B', C', and D' in Figure 8b. Similarly, the regions E, F, and G are also priority regions for equiaxed dendrite nucleation, as shown in Figure 8c. Figure 9 is a solute field diagram during dendrite growth when the undercooling degree is 500 K. The temperature gradient in the melt is relatively large at this condition, which is conducive to the directional growth of columnar dendrites. In addition, a narrow mush zone is formed and the nucleation of equiaxed dendrites appears in the region B between the columnar dendrites firstly, which can be clearly seen in Figure 8a. From the solidification process described in Figure 8, it can be seen that the equiaxed dendrites nucleation area is local and small, which exhibited as the layer-by-layer advancement of equiaxed dendrites. The regions A and C with lower solute enrichment evolve into the regions of equiaxed dendrites nucleation and growth, as shown in regions A' and C' (as illustrated in Figure 9a,b). In the same way, the regions D, E, and F gradually evolve into the regions D', E', F' in Figure 9b,c. The result indicates that with the increase of the undercooling on the surface of the crystal nucleus generator, the undercooled zone width and the nucleation probability of equiaxed dendrites decrease.  The cooling rates in the mushy zone at different locations are plotted in Figure 7b,c, respectively. The cooling rate of point A gradually decreases from a high value to a low value as the solidification proceeds, which corresponds to the large discrepancy between the surface temperature of the crystal nucleus generator and the temperature in the melt at the initial stage, but the cooling rate of point A decreases visibly with the increase of the heat transfer as the solidification progresses. The cooling rate of point B rises from a lower value to a higher and steady value. The steady values are 0.56, 0.63, and 0.72 K/s, approximately corresponding to the minimum rate at point A at the same time, which illustrates that the solidification process gradually transitions from point A to point B. The large undercooling could provide a high cooling rate, which can promote the nucleation and growth of equiaxed grains; however, it is found that the large undercooling could create a strong temperature gradient distribution at the front of columnar grains, it was a favorable condition for columnar dendrites growth. Crystals 2021, 11, x FOR PEER REVIEW 15 of 20 Figure 9. The solute field during the dendrite growth process when the undercooling of the crystal nucleus generator is 500 K. The solidification time from (a)-(c) is 6.8, 7.0, and 7.2 s, respectively.
The cooling rates in the mushy zone at different locations are plotted in Figure 7b,c, respectively. The cooling rate of point A gradually decreases from a high value to a low value as the solidification proceeds, which corresponds to the large discrepancy between the surface temperature of the crystal nucleus generator and the temperature in the melt at the initial stage, but the cooling rate of point A decreases visibly with the increase of the heat transfer as the solidification progresses. The cooling rate of point B rises from a lower value to a higher and steady value. The steady values are 0.56, 0.63, and 0.72 K/s, approximately corresponding to the minimum rate at point A at the same time, which illustrates that the solidification process gradually transitions from point A to point B. The large undercooling could provide a high cooling rate, which can promote the nucleation and growth of equiaxed grains; however, it is found that the large undercooling could create a strong temperature gradient distribution at the front of columnar grains, it was a favorable condition for columnar dendrites growth. Figure 7d shows the average grain size distribution under different degrees of undercooling on the surface of the bar, which is measured from the simulation results. Another noticeable phenomenon can be observed, i.e., the average grain size is formed when the undercooling degree of the vibrating bar surface at 300 K and 400 K is finer than that at 200 K, which reflects that the influence of undercooling factors on microstructure evolution is not a simple linear relationship. The reason for the above phenomenon is probably that when the degree of undercooling on the surface of the crystal nucleus generator  Figure 7d shows the average grain size distribution under different degrees of undercooling on the surface of the bar, which is measured from the simulation results. Another noticeable phenomenon can be observed, i.e., the average grain size is formed when the undercooling degree of the vibrating bar surface at 300 K and 400 K is finer than that at 200 K, which reflects that the influence of undercooling factors on microstructure evolution is not a simple linear relationship. The reason for the above phenomenon is probably that when the degree of undercooling on the surface of the crystal nucleus generator is relatively lower, the crystal nucleus generator does not significantly affect the distribution of the temperature field and cooling intensity in the melt, and it is not conducive to creating conditions for heterogeneous nucleation of central equiaxed dendrites; therefore, the grain size of newly formed dendrites is relatively large.
The simulation results show that different temperature field distribution conditions will significantly affect the final microstructure. It can be concluded that under the vibration frequency of 1000 Hz, the undercooling on the surface of the crystal nucleus generator was maintained at 300~400 K; it was beneficial to the formation and refinement of the central equiaxed crystal structure by affecting the temperature field distribution conditions in the melt.

Influence of the External Vibration Effect
In this group of simulations, the undercooling of the crystal nucleus generator was a certain value of 400 K, and the structure formation under different vibration frequencies of the crystal nucleus generator was explored. Figure 10 shows the dendritic structures at the vibration frequency varies from 500 to 1800 Hz. With the increase in the vibration frequency, the area of the equiaxed dendrites regions is obviously extended, the development of columnar crystal regions is hindered, the numbers of equiaxed dendrites increase, and the average grain size decreases.  Figure 11b describes the average grain size under different vibration frequency conditions. As the vibration frequency increases, the average grain size decreases obviously. The reasons can be summarized in the following two aspects. As the vibration frequency As analyzed above, the possibility of dendrite fragmentation, separation, and settlement will increase with the convection effects caused by vibration are intensified, and thus the nucleation probability P n in the melt are increased to correspond to different conditions. Moreover, the enhancement of the vibration intensity can improve the convection stirring ability of the melt. It was found in previous studies that the vibration produced a circulation flow around the crystal nucleus generator, which drives the melt near the wall to collide with the mold wall and form a turbulent flow effect. As solidification progresses, the turbulent flow effect in the melt intensifies significantly.
The distribution of temperature gradient between A and B with the solidification time was investigated, as shown in Figure 11a. It is clearly shown that as solidification progresses, the temperature gradient between the two points drops significantly from a maximum value to a steady-state value. With the increase in the vibration frequency, the temperature gradient curve drops faster and reaches a relatively lower steady-state value. The large value of temperature gradient at the initial stage is due to the chilling crystal generator that is inserted into the melt; as the vibration starts, the strong fluid flow could improve the efficiency of the heat dissipation in the melt, which resulted in uniform temperature gradient distribution. Increasing the vibration frequency could significantly reduce the temperature gradient-the vibration frequency is increased from 500 to 1800 Hz, and the corresponding temperature gradient stability values in the melt after 6 s are 3.12 K/mm, 2.55 K/mm, 2.30 K/mm, 1.80 K/mm, 1.50 K/mm, respectively. increased, the thermal conductivity of the melt increased gradually, which resulted in the uniform distribution of temperature gradient, the superheat in front of the columnar crystal was reduced, and a larger part of the region entered an undercooled state, creating the necessary conditions for the nucleation of equiaxed crystals. Furthermore, the number of heterogeneous nucleation sites increased due to vibration force that breaks the primary dendrites and crushes them into fragments. According to the above analysis result, it can be found that selecting the appropriate process parameters exerts an important influence on the grain structure evolution, which can effectively increase the proportion of the central equiaxed crystal area, shorten the area of the columnar crystal area, and achieve refined grains and homogenized microstructure.

Conclusions
In this study, a macro-micro coupled model for simulating the process of nucleation and competitive dendrite growth was developed to study the morphological evolution of columnar-to-equiaxed transition in the vibration-excited liquid metal nucleation technique. The microstructure evolution of ferritic stainless steel Cr17 was predicted and the influence of the different process parameters on the CET location, number of newly formed dendrites, and average grain size were investigated. The main conclusions are summarized as follows: 1. With the increase of surface undercooling on the crystal nucleus generator, the columnar dendrites regions are more developed, and the equiaxed dendrites regions decrease visibly, the number of new-born equiaxed dendrites decreases, and the average grain size increases. The quantitative analysis of the temperature field distribution in the melt shows that the temperature gradient at the columnar crystal front played a critical role in the nucleation and growth of the center equiaxed grain; 2. With the vibration frequency increases, the area of the equiaxed dendrites regions is obviously extended, and the refined grains and homogenized microstructure are achieved. With the vibration intensified, the possibility of dendrite fragmentation is increased by increasing the probability of nucleation in the melt; the convection stirring ability of the melt is improved, resulting in uniform temperature gradient dis- Figure 11b describes the average grain size under different vibration frequency conditions. As the vibration frequency increases, the average grain size decreases obviously. The reasons can be summarized in the following two aspects. As the vibration frequency increased, the thermal conductivity of the melt increased gradually, which resulted in the uniform distribution of temperature gradient, the superheat in front of the columnar crystal was reduced, and a larger part of the region entered an undercooled state, creating the necessary conditions for the nucleation of equiaxed crystals. Furthermore, the number of heterogeneous nucleation sites increased due to vibration force that breaks the primary dendrites and crushes them into fragments.
According to the above analysis result, it can be found that selecting the appropriate process parameters exerts an important influence on the grain structure evolution, which can effectively increase the proportion of the central equiaxed crystal area, shorten the area of the columnar crystal area, and achieve refined grains and homogenized microstructure.

Conclusions
In this study, a macro-micro coupled model for simulating the process of nucleation and competitive dendrite growth was developed to study the morphological evolution of columnar-to-equiaxed transition in the vibration-excited liquid metal nucleation technique. The microstructure evolution of ferritic stainless steel Cr17 was predicted and the influence of the different process parameters on the CET location, number of newly formed dendrites, and average grain size were investigated. The main conclusions are summarized as follows: 1.
With the increase of surface undercooling on the crystal nucleus generator, the columnar dendrites regions are more developed, and the equiaxed dendrites regions decrease visibly, the number of new-born equiaxed dendrites decreases, and the average grain size increases. The quantitative analysis of the temperature field distribution in the melt shows that the temperature gradient at the columnar crystal front played a critical role in the nucleation and growth of the center equiaxed grain; 2.
With the vibration frequency increases, the area of the equiaxed dendrites regions is obviously extended, and the refined grains and homogenized microstructure are achieved. With the vibration intensified, the possibility of dendrite fragmentation is increased by increasing the probability of nucleation in the melt; the convection stirring ability of the melt is improved, resulting in uniform temperature gradient distribution and making a large area in the melt to reach the critical nucleation undercooling.