Several phenomena can lead to grain size evolution in a metallic material during hot forming processes and heat treatments. When a material accumulates plastic deformation, dislocation networks develop within the grains. Subsequently, several thermal activated processes tend to reduce the energy of the system. These processes include recovery, that annihilates and rearranges the dislocations1, grain boundary migration due to capillary effects (the so-called grain growth) and stored energy gradient (also called strain-induced boundary migration), and recrystallization (when the stored energy promotes the formation of new grains which nucleate mainly at grain boundaries).
Different criteria must be locally achieved to promote the evolution of a cluster of atoms into a nucleus. The first is the formation of a high-angle mobile grain boundary during the nucleation process. The second is the presence of a strong gradient of stored energy across the interface, which provides sufficient driving pressure to counteract the capillary effects applied to the nucleus.
DIGIMU® relies on a Finite-Elements – Level-Set (FE-LS) formalism2 for modelling microstructural evolutions at the polycrystal level. Grain boundary kinetics are piloted by a thermo-dependent grain boundary mobility Mb, the capillary pressure depending on the grain boundary energy
where
Discontinuous dynamic recrystallization (DRX) is modelled using a critical nucleation energy
The nuclei are created with the critical nucleus radius r*of the « Bailey-Hirsch » criterion5, which is required to compensate for capillary effects exerted by neighbouring grains. This condition is approximated by this criterion, such that:
where
Furthermore, the nucleation rate is calculated according to a variant of the « Peczak & Luton »6 proportional nucleation model. At a given time
Although this model effectively describes DRX, the fact that for a given temperature and strain rate, all nuclei appear of the same radius have shown some limits for post-dynamic recrystallization (PDRX) evolutions. If nuclei sizes are too similar in size, capillary effects are compensated between neighbouring nuclei and grain growth is slowed down. Experimental micrographs show clearly a size dispersion around the theoretical value, due to local intragranular energy distributions (note that most of these micrographs are in 2D, meaning the observed sizes correspond to lateral sections of the nuclei). Taking that into account in the simulation maintains a certain competition between the nuclei and generates some grain growth due to capillary effects7.
This article demonstrates how integrating a nuclei size dispersion in DIGIMU® will maintain this competition between them, and lead to more realistic grain size evolutions and distributions during PDRX.
In a study conducted in collaboration with AUBERT & DUVAL, a comparison of experimentally measured grain size distribution, in an Inconel 718 alloy following DRX/PDRX, with that obtained by DIGIMU® for a constant nucleus radius, showed a slight divergence (see Figure 3). To remedy this discrepancy, we considered a log-normal distribution of the nuclei size around the calculated Bailey-Hirsch value and studied the impact of such a distribution on the resulting microstructure at the end of DRX/PDRX. Four different distributions, shown in Figure 1, were tested in our calculations to identify the one closest to the experimental results.
Figure 1: Representation of the tested lognormal nucleus size distributions. The value 1 in the x-axis represents the mean nucleus size (calculated with the Bailey-Hirsch criterion).
By following the microstructure evolution for each case, at the end of DRX and during PDRX, we compared the average equivalent diameter (Deq) and the recrystallized fraction (FRX) of the grains. We observed that the broader the nucleus size distributions, the larger the grain size at the end of DRX and the lower the recrystallized fraction. The results are summarized in the Figure 2, at the beginning and after 300 s of PDRX. A broadening of Deq values is observed compared with the constant radius case, particularly towards larger values around 20 µm. Among the different distributions, the 4th distribution appears to be the closest to the experimental result. However, a discrepancy in the recrystallized fraction, at the beginning of the PDRX, emerges between the simulation with constant nucleus size (100%) and the 4th distribution (73%). This is an undesired effect, since the FRX values obtained with the constant nucleus size were calibrated in previous studies performed on the same alloy, resulting in satisfactory outcomes.
Figure 2: Simulation results obtained at the end of DRX, which marks the start of PDRX (tPDRX=0), and after a duration of 300 s. The colored images show the microstructures during PDRX for the four nucleus size distributions, and the histograms represent the corresponding mean equivalent grain diameter distributions.
In order to make up for this discrepancy, we can adjust the parameters
To validate our improvements, the same material (Inconel 718 - 10/11 ASTM) was studied again in collaboration with AUBERT & DUVAL. The material was preheated (1000°C, 15 min), deformed (975°C, 5 min, 0.02 s-1, ε = 1.3) and then cooled to 860°C at a rate of 140°C/min. The evolution of the microstructure was then examined at the end of DRX and during PDRX. The histograms in Figure 3 show the results obtained for the different simulations, in comparison with the experimental results.
Figure 3: Comparison between the grain size distribution (in µm and ASTM) in an Inconel 718 alloy measured experimentally (in orange) at the end of DRX/PDRX, and the one obtained by DIGIMU® for: a constant nucleus radius (in green) and the 4th distribution with parameter adjustment (in blue).
First, we note that the simulation’s result without a nucleus size distribution reveals grain sizes centred around the mean value, in comparison with the experimental distribution which is much flatter. In fact, the simulation without a nucleus size distribution gives rise to a peak at 12 ASTM which is not present in the experimental result. This can be explained by an evolution of grain sizes in the experiment towards higher values (10-9 ASTM).
Nevertheless, for the lower classes (15-14 ASTM), we observe results relatively close to reality, although the peak at 15 ASTM is absent in the experimental distribution. This can be explained by the experimental characterization of the microstructure, which ignores smaller grains.
On the other hand, when we deploy DIGIMU®'s new ability to insert nuclei following a size distribution, the final evolution of the microstructure shows a distribution closer to experimental observations. In fact, the resulting grain size distribution after DRX/PDRX is flatter than the one that does not consider a nucleus size distribution. This result can be explained by the difference in the size of the nuclei at their creation, which gives rise to different curvatures, thus generating a driving force that favours competition between nuclei and influences their subsequent evolution.
Note that the absence of 9 ASTM grains can be attributed to a slower evolution in the simulation. This could result from various factors, such as slightly reduced mobility, minor temperature fluctuations during the experiment (a few degrees below 980°C), or a combination of both.
To highlight the benefits of incorporating a nucleus size distribution, we present an animation below that compares the evolution of a microstructure with and without nucleus size distribution. In the upper simulation, which uses the 4th distribution, the histogram reveals a broader spread of grain diameters and a more flattened profile. This highlights a heterogeneous grain size distribution that can be seen in the final microstructure. In contrast, the lower simulation, which does not account for nucleus size distribution, displays a homogeneous grain size with a histogram centred around the calculated Bailey-Hirsch value. This comparison clearly demonstrates how the adoption of a nucleus size distribution allows to obtain realistic microstructural simulations.
Figure 4: Animation comparing the evolution of a microstructure with and without the distribution of nuclei sizes.
In the latest updates of DIGIMU®, exciting functionalities were introduced to refine microstructural simulations and align them closer with experimental results. Among the highlights is the incorporation of nucleus size distribution during nucleation, enhancing the accuracy of simulated final grain size distributions as demonstrated in the present article for an Inconel 718. Additional features include heterogeneous hardening and updates to material files. These tools are ready to help you unlock better predictions and optimize your material properties.
Take your simulations to the next level with DIGIMU® 5.1!
References