The Effects of Inhomogeneous Elasticity and Dislocation on Thermodynamics and the Kinetics of the Spinodal Decomposition of a Fe-Cr System: A Phase-Field Study

: The effects of inhomogeneous elasticity and dislocation on the microstructure evolution of α (cid:48) precipitate in a Fe-Cr system was investigated using a Computer Coupling of Phase Diagrams and Thermochemistry (CALPHAD)-type free energy incorporated phase-ﬁeld method. In order to simulate the precipitation behavior by phase-ﬁeld modeling in consideration of inhomogeneous elasticity, a Multiphysics Object-Oriented Simulation Environment (MOOSE) framework was used, which makes it easy to use powerful numerical means such as parallel computing and ﬁnite element method (FEM) solver. The effect of inhomogeneous elasticity due to the compositional inhomogeneity or the presence of dislocations affects the thermodynamic properties of the system was investigated, such as the lowest Cr concentration at which spinodal decomposition occurs. The effect of inhomogeneous elasticity on phase separation kinetics is also studied. Finally, we analyzed how inhomogeneous elasticity caused by compositional ﬂuctuation or dislocation affects microstructure characteristics such as ratio between maximum precipitate size with respect to the average on early stage and later stage, respectively.


Introduction
α precipitate hardening in a Fe-Cr system is an important research subject in terms of maintaining the mechanical integrity of ferritic steels [1][2][3][4][5][6][7]. Precipitation can occur through the spinodal decomposition and the nucleation and growth mechanism, and researchers reported that both mechanisms are involved in the precipitation of α [2,3]. The effects of inhomogeneous elasticity and dislocation on Fe-Cr spinodal decomposition have been investigated [6] using the phase-field method; however, the evaluation of the minimum Cr concentration at which spinodal decomposition occurs was not performed, and the quantitative microstructure analysis was also limited. In this study, we analyzed the effect of the dislocation on the kinetics of the spinodal decomposition and, in particular, the effect on the lowest nominal composition of Cr for the spinodal decomposition. Through this study, we improved our systematic understanding regarding α nucleation in the presence of the dislocation.
The role of the elasticity on the spinodal decomposition was analyzed using analytic approaches [8,9]. For cubic crystals [8], Cahn claimed that coherent elastic interaction increased the free energy, and therefore inhibited the spinodal decomposition. On the other hand, Eshelby predicted that spinodal decomposition would be accelerated when elastic interaction was considered based on the theory of continuum elastic mechanics [9]. At first glance, this appears to be a conflicting argument; however, Eshelby's predictions describe the rate at which spinodal decomposition occurs and Cahn's

Cahn-Hilliard Equation
The Cr concentration field will be evolved by solving the following Cahn-Hilliard equation [13,14], The molar free energy F(r, t) in Equation (1) is given by Equation (2), and M(r, t) is the mobility of the diffusion species. The mobility was assumed as a constant, and the value was set to 1.0 in this study. Moreover, κ is the gradient coefficient.

Free Energy of Fe-Cr System with Dislocation
The regular-solution type chemical free energy f (c) in Equation (2) is obtained from [15] where c is the composition of Cr, and G 0 Fe and G 0 Cr are the molar Gibbs free energies for the elements Fe and Cr, respectively. L FeCr is the interaction parameter of Fe and Cr. The detailed numbers are written in [11,15]. We set the temperature as T = 535 K. Moreover, we plot the phase diagram in Fe-Cr system in Figure 1 using the FactSage thermochemical software and the FSstel database, which is the FactSage steel alloy database [16].
We determined the elastic energy density f e (c, r) using Equation (2). We implemented Khachaturyan's strain interpolation scheme (KHS) [17]. The stiffness tensor and misfit strain were obtained by the interpolation and create a global stiffness tensor, which is a linear combination of end members Cr and Fe, and the global strain is KHS ij where tot ij (r) is the total strain, • ij,α is the eigenstrain of α phase, and • ij,d (r) is the eigenstrain induced by the dislocation. The elastic energy density is given by The mechanical equilibrium equation was solved to obtain the strain-stress state as below.

Computational Details
First, we normalized the free energy in Equation (3) by RT for computational efficiency.
Two options were considered for the normalized strain tensor (1) written in Table 1A [6], and (2) written in Table 1B. The stiffness tensor of Fe and Cr in Table 1A has a value at room temperature, and the values in Table 1B are dependent on the temperature. The stiffness tensor of Fe in Table 1B at 535K was calculated as quadratic approximation by extrapolation. As the stiffness of Cr has insufficient data to extrapolation, it was assumed to be the values at 500K. To examine how much the microstructure varies depending on which of the two different stiffness tensors in Table 1 is chosen, we evaluated microstructure of Fe-21Cr system with two different sets of stiffness tensor and the results of the average α radius and the number of α precipitation were plotted in Figures 2 and 3. As no significant difference was found in the microstructure characteristics in the two results, it was chosen to use a value at room temperature in Table 1A rather than the value calculated by extrapolation, even if there was a dependence on temperature in the stiffness. Table 1. The normalized stiffness tensor of Fe and Cr at 535K in Voigt notation [18][19][20].  Table 1A without dislocation, (c) with inhomogeneous elasticity having the stiffness in Table 1A and dislocation, (d) with inhomogeneous elasticity having the stiffness in Table 1B without dislocation, and (e) with inhomogeneous elasticity having the stiffness in Table 1B with dislocation.  The number of α precipitations with respect to the time for the case of (a) without elasticity, (b) with inhomogeneous elasticity having the stiffness in Table 1A without dislocation, (c) with inhomogeneous elasticity having the stiffness in Table 1A and dislocation, (d) with inhomogeneous elasticity having the stiffness in Table 1B without dislocation, and (e) with inhomogeneous elasticity having the stiffness in Table 1B with dislocation. We used κ = 4.25 in Equation (2). The eigenstrain of α phase • ij,α in Equation (4) is 0.006 when (4) is evaluated by the authors of [21].

Results and Discussions
The Effect of Elasticity and Dislocation on α Phase Precipitation Induced by the Spinodal Decomposition Mechanism We examined multiple compositions within the range of Fe-14Cr to Fe-21Cr to investigate the role of dislocation on α precipitate nucleation through the spinodal decomposition. To examine the role of the inhomogeneous elasticity and the dislocation, we performed simulations (1) without the f e term in Equation (3) (no elasticity case), (2) with the f e term without dislocation (with elasticity only arising from the concentration inhomogeneity considered), and (3) with the f e term with the dislocation.
Typically, the concentration of Cr in Fe-Cr steel is expressed as wt%; however, in this study, the value was set in terms of at% in the simulation. For example, in the Fe-13Cr system, we set the initial average composition of Cr was 13.86 at%. The initial Cr concentration spanned from 13.36 at% to 14.36 at%. The range of the fluctuation was ±0.5 at% for all the cases. First, we found that the presence of the inhomogeneous concentration did not alter the minimum Cr average concentration for the spinodal decomposition of Fe-16.7Cr; however, it altered the kinetics significantly. We observed the α precipitate from t = 14,300 when inhomogeneous elasticity is not incorporated; on the other hand, we observed α precipitate from t = 1732 when we considered inhomogeneous elasticity. This observation is consistent with the predictions made by Eshelby [9].
We also found that presence of the dislocation lowered the minimum Cr concentration for spinodal decomposition. As shown in Figure 4, we found that the α phase nucleated even in Fe-16.4Cr. As the Cr concentration increased in the vicinity of the dislocation when compared to the nominal composition, dislocation is a suitable location for the α nucleation as experimentally observed [22]. To confirm these findings, we performed four different sets of simulation with different random seeds.
We also analyzed the role of the inhomogeneous elasticity and dislocation on the microstructural evolution of Fe-21Cr in Figure 5. As shown in Figures 2 and 6, when elasticity was not considered, the microstructure evolution was the fastest. When elasticity only due to the compositional inhomogeneity was considered, the average precipitate size increased and the increased rate of the α phase fraction were the slowest among the three cases.
In summary, when elasticity was not considered, the system approached equilibrium the most rapidly. However, even if elasticity was considered, if the effects of structural inhomogeneity, such as dislocation, were considered simultaneously, the growth rate of the precipitate was faster in the case of the absence of the dislocation with consideration of the elasticity. At t = 70, as shown in Figure 7f, there was Cr concentration partitioning when a dislocation was present in the particle; therefore, a relatively large α precipitate existed in the vicinity of the dislocation.
In Figure 5, we observed overshootings in the R m /<R> plot at early stages for all three cases, where R m represents the maximum radius of the α phase and <R> is the average radius of the α precipitate. In the presence of elasticity only arising from compositional inhomogeneity, transient overshooting was lower than case of the absence of the elasticity. However, with the presence of dislocation, the largest overshooting existed at approximately t = 70, where R m /<R> exceeded 3.5. This indicates that the largest relative α precipitate existed with respect to the average value among all situations.
As time proceeded, <R> increased and, after a certain time (after t = 470 in Figure 5), R m /<R> demonstrated the maximum value with the presence of inhomogeneous elasticity without any dislocation. When a dislocation was present in the system, the morphology of the α precipitate deviated from the circular shape as shown in Figure 4; thereby, the driving force for the coarsening reduced [23]. Consequently, the original largest precipitates in the vicinity of the dislocation grew relatively slowly in the late stage, and the value of R m /<R> became smaller compared to the case of the presence of dislocation.
According to the analytic theory developed by Kawasaki and Enomoto [24], the size distribution of coarsened particles becomes wide when a repulsive elastic field is applied. As the eigenstrain is proportional to the Cr composition, we determined that repulsive elastic fields were generated between the Cr rich phases, and accordingly, the R m /<R> value was found to be larger in the presence of elasticity due to the concentration inhomogeneity without dislocation in the later stage.
In the following Figures 7 and 8, the α phase nucleation was more progressed in the Fe-21Cr system without elasticity (a) than in the cases with inhomogeneous elasticity at t = 70 (b,c). The Fe-Cr21 system with inhomogeneous elasticity and the dislocation showed more active α nucleation in the vicinity of the dislocation than at other sites. In this study, the effect of inhomogeneous elasticity upon precipitation of precipitates by spinodal decomposition was evaluated, but, in fact, precipitates can also be initiated and evolved by nucleation and growth mechanisms [25]. In fact, the authors of [25] state that the α precipitation mechanism changes on Fe-Cr at 28Cr. It is necessary to quantitatively analyze the effect of inhomogeneous elasticity when α phase is precipitated by the nucleation and growth mechanism. A more reliable explanation on the effect of elasticity on nucleation frequency or critical nuclei size is necessary to quantitatively predict phase separation kinetics induced by nucleation and growth mechanism.  (a) Fe-21Cr without elasticity (b) Fe-21Cr with inhomogeneous elasticity without dislocation (c) Fe-21Cr with inhomogeneous elasticity and dislocation Figure 5. R m /<R> with respect to time for the cases of (a) without elasticity, (b) with inhomogeneous elasticity without dislocation, and (c) with inhomogeneous elasticity and dislocation, where R m is the maximum radius of the α precipitate and <R> is the average radius of it. We performed four different sets of simulations to evaluate R m and <R> for enhanced statistical rigor of the analysis. The system size was 256 × 256 µm 2 . The composition of the system is Fe-21Cr. Cr segregation in the vicinity of dislocation (loop) has been consistently observed experimentally [22,26,27]; however, no studies have been found showing evidence that the segregation leads to α precipitation. It is possible to estimate that a locally high Cr concentration near the microstructure inhomogeneity, such as dislocation, will help in the formation of Cr-rich precipitate, and this study showed the anticipation is the case using the phase-field simulation. Studies on α behavior at the early stage of the phase separation in the Fe-Cr system have been conducted experimentally with relatively high Cr compositions, such as Fe-32 wt%Cr [5,28], Fe-36 wt%Cr [29] and Fe-25.28 wt%Cr, and Fe-30.42 wt%Cr and Fe-36.10 wt%Cr [30]. Their concentration ranges are clearly regions in which spinodal decomposition occurs, and the above studies [5,[28][29][30] are valuable studies that help to understand the early-stage microstructural behavior when the phase separation occurs by spinodal decomposition. Relatively recently, there has been a report that α is precipitated even in Fe-9 at%Cr during neutron irradiation. In this case, one cannot confidently claim that α precipitation is nucleated by spinodal decomposition, if no other effects are included other than the composition. From this experimental observation, we hypothesized that the dislocation loop generated by neutron irradiation caused Cr partitioning to lower the lowest average Cr concentration at which phase separation occurred, and phase-field modeling was used to verify it. Our results found that the dislocation lowered the lowest mean Cr concentration at which spinodal occurs, even though the effect is relatively marginal. Additionally, it was found that the incubation time required for α to precipitate was significantly reduced for Fe-16.4 wt%Cr by Cr partitioning by the dislocation. Moreover, we show that the model for Fe-Cr microstructure evolution with elastic effect proposed in [6] can be performed without advanced programming skill by utilizing MOOSE framework.  Microstructures of the Fe-21Cr system (a) without elasticity at t = 100, (b) with inhomogeneous elasticity without dislocation at t = 100, (c) with inhomogeneous elasticity and dislocation at t = 100, (d) without elasticity at t = 500, (e) with inhomogeneous elasticity without dislocation at t = 500, and (f) with inhomogeneous elasticity and dislocation at t = 500. The system size was 256 × 256 µm 2 . The microstructures were visualized by Cr composition.

Conclusions
Phase-field simulation using the MOOSE framework has been utilized to study the effect of inhomogeneous elasticity including dislocations on α precipitate nucleation through the spinodal decomposition in Fe-Cr alloy. Our study yields three main findings: We found that inhomogeneous elasticity induced by the compositional inhomogeneity of the α phase did not alter the minimum Cr composition for the spinodal decomposition; however, it altered the kinetics significantly. Moreover, we found that the presence of the dislocation lowered the minimum Cr composition for spinodal decomposition. At the early stage of the phase separation, the relative size of the α precipitate was maximized when a dislocation existed in the system. In the later stage, the relative size of the precipitate showed the maximum size when inhomogeneous elasticity was incorporated without any dislocation.