Initial Microstructure
The average precipitate morphologies in each sample were quantified with the use of ImageJ on the backscattered electron (BSE) micrographs shown in Fig. 4. Consistent with published literature on similar alloys [42, 55–58], all the microstructures were found to consist of nearly spheroidal low atomic number contrast precipitates in a higher atomic number contrast matrix. The precipitates appeared to be uniform within each specimen and, generally, increased in size and volume fraction with increased annealing time and decreased annealing temperature, respectively, as is shown in Fig. 5. These findings were consistent with the overall trend as calculated by TC-Prisma, with most predicted average diameters falling within the measured 95% confidence interval. To confirm the changes in microstructure, 30-point arrays of Vickers hardness indents were made and measured (Fig. 6). The median sample hardness increased with both increasing precipitate volume fraction and decreasing precipitate size as expected from first principles calculations [59].
To further characterize the chemistry of the precipitate phase, APT tips were extracted from bulk polished samples. Subsequent data reconstruction and analysis, summarized in Fig. 7 and Table 4, highlighted the significant chemical segregation between the matrix and intermetallic precipitates. The results are summarized in Table 4. The precipitate chemistry, determined by isoconcentration surfaces from selected regions of interest, was found to be similar to classic superalloy \({{\gamma }}^{{\prime }}\)-Ni3(Al,Ti) phase but with some additional substitution of cobalt [60, 61]. The matrix was found to contain high concentrations of nickel, cobalt, iron, and chromium. The configurational entropy, \(\varDelta {S}_{conf}\), of the total alloy and matrix phase as a function of mole fraction, \({X}_{i}\), was calculated using Eq. 1:
$$\varDelta {S}_{conf}=-R\sum _{i}{X}_{i}\text{ln}\left({X}_{i}\right)$$
1
Both the bulk alloy and matrix conform to the empirical high entropy alloy definition with an average \(\varDelta {S}_{conf}\) of \(1.54R\) and \(1.57R\), respectively [4]. Further analysis of the APT data showed little-to-no enrichment of the phase boundaries, indicating that the elements within the precipitates and matrix were uniformly distributed and that a variation of the inverse lever rule could be used to estimate each samples volume fraction more accurately. If the nominal and precipitate atomic concentrations are taken as \({X}_{\text{n}\text{o}\text{m}}\) and \({X}_{P}\), respectively, then (Eq. 2):
$${X}_{\text{n}\text{o}\text{m}}={\phi }_{P}{X}_{P}+\left(1-{\phi }_{P}\right){X}_{P}$$
2
where \({\phi }_{P}\) is the volume fraction of precipitate phase [60, 61]. Plotting \(\left({X}_{\text{n}\text{o}\text{m}}-{X}_{P}\right)\) vs. \(\left({X}_{P}-{X}_{\text{M}\text{a}\text{t}\text{r}\text{i}\text{x}}\right)\) for each atomic species should yield an approximately linear relationship, the slope of which is the volume fraction of precipitate phase (Fig. 8).
Table 4
Mean atomic concentrations for the original casting calculated by mass spectrum analysis of 1000L, 1000S, 800L, and 800S APT specimens
Ion
|
Average Phase Composition (at%)
|
Calculated Bulk Composition (at%)
|
PPT
|
Matrix
|
–
|
Ni
|
60.1 ± 2.4
|
32.8 ± 3.8
|
45.2 ± 7.4
|
Al
|
12.2 ± 0.7
|
5.8 ± 0.9
|
8.4 ± 1.6
|
Ti
|
11.9 ± 0.4
|
2.4 ± 0.9
|
6.7 ± 3.1
|
Cr
|
2.6 ± 0.8
|
19.9 ± 2.3
|
12.3 ± 5.2
|
Fe
|
3.1 ± 0.7
|
16.3 ± 1.6
|
10.4 ± 3.7
|
Co
|
9.9 ± 1.7
|
22.1 ± 1.8
|
16.9 ± 3.1
|
Comparison of the EDS and APT measurements as well as examination of the inter-precipitate spacing indicated that the sampling distribution of the APT may be over sampling the matrix. To obtain more representative volume fractions, \({X}_{nom}\) was taken as the atomic concentrations measured through bulk EDS for the specific sample group (Table 1), rather than what was measured from the total APT volumes. The resulting phase fractions correlated well the values determined from backscattered SEM images and those predicted by ThermoCalc. The bulk concentration measured by APT indicates that the samples were significantly leaner in nickel than originally synthesized which may indicate some sampling bias in the selected areas or poor separation of overlapped mass spectrum peaks.
Oxidation Behavior
Figure 9 shows a representative FIB trench milled from sample 800L. All the investigated samples formed a similar tri-layer oxide structure consisting of a faceted external layer consisting of
TiO2 + FexOy, an underlying continuous Cr2O3 layer, and a discontinuous internally oxidized Al2O3 subsurface region.
Figure 10 shows (a) nonisothermal-isothermal and (b) cyclic oxidation results for the samples annealed at 1000°C. No statistically significant change in overall mass gain was observed in the nonisothermal-isothermal TGA tests (Fig. 10a). However, there was a measurable difference in the parabolic oxidation rate between the two groups. The alloy that contained smaller precipitates (i.e., 1000S) exhibited a lower average \({K}_{p}\). Samples 1000L and 1000S behaved similarly in terms of total mass gain with each reaching steady oxidation, defined as the time until parabolic oxide growth is achieved, in approximately 10 hours. Under cyclic loading conditions, 1000S outperformed 1000L in terms of oxidation resistance and gained less specific mass overall (Fig. 10b).
Figure 11 shows BSE cross sections, EDS maps, and surface oxide morphologies for TGA specimens of 1000L and 1000S. Both groups displayed similar three-layered microstructures as described above, however, there was a difference in the Al2O3 morphology and subsequently a measurable change in the outer Cr2O3 layer thickness. Externally, no appreciable difference was observed in outer scale that formed after 100 hours of oxidation. Closer inspection of the metal-oxide interface and the internally grown Al2O3 showed that the internal oxide layer of sample 1000S was more continuous than the one in 1000L and had begun to grow a mixed titanium – aluminum oxide that was reminiscent of the beginnings of internal-to-external scale conversion (shown in the magnified inset). The Cr2O3 scale in sample 1000L was found to have a thickness of \(8.83\pm 0.97 {\mu }\text{m}\) while the scale in sample 1000S was found to be thinner at \(7.95\pm 1.23 {\mu }\text{m}\).
Figure 12 shows results from nonisothermal-isothermal and cyclic oxidation tests on samples annealed at 800°C. The smaller precipitate containing samples (i.e., 800S) gained less mass for the given oxidation time with lower parabolic oxidation rates in both nonisothermal-isothermal and cyclic tests. Examination in cross section, the nonisothermal-isothermal test specimens (Fig. 13) exhibited the same three-layered oxide structure as the 1000°C specimens, with sample 800S having a more contiguous Al2O3 layer than sample 800L. As observed in the other sample groups, a systematic variation in Cr2O3 layer thickness was observed with the larger precipitate containing material displaying a thicker Cr2O3 scale after 100 hours accounting for the observed increase in mass gain for the 800L samples.
To elucidate the microstructural evolution during earlier stages of oxidation, several nonisothermal-isothermal oxidation samples were removed from the TGA after approximately 30 minutes (see Fig. 14). The 800S samples formed a nearly continuous inner Al2O3 layer and a more cohesive outer TiO2 + FexOy scale as compared to the 800L samples after the same oxidation time. The 900°C annealed sample group (900L and 900S) displayed similar behavior to that of the 800L and 800S samples with smaller precipitate containing samples performing better. Several of the samples in this group were found to contain shrinkage cavities which subsequently caused some samples to crack during sample preparation. This led to an increase in total oxide mass gain over the prescribed time in nonisothermal-isothermal testing and produced inconsistent results in cyclic testing (Fig. 15).
While the inconsistent mass gains may have obscured the true oxidation reaction rates, the microstructural changes observed in the samples mirror those of previous specimens (Fig. 16). The samples containing smaller precipitates (i.e., 900S), consistently exhibited a more continuous Al2O3 sublayer than the larger precipitate containing samples (i.e., 900L). In addition, the 900L samples also exhibited more spallation and cracking of the oxide scale. As such, only small regions of those samples contained cohesive, measurable layers.
Oxidation Kinetics and Modelling
To understand the kinetic oxidation behavior of each sample group, a generalized reaction rate model, expressed in Eq. 3, was fit to the specific mass gains for each sample.
$${\left(\text{S}.\text{M}.\text{G}.\right)}^{n}={k}_{n}t$$
3
The degree of the reaction (i.e., the reaction kinetics exponent), \(n\), should vary between ~ 1 and ~ 3, representing the spread from linear through parabolic and logarithmic growth (Fig. 17) [20, 62]. The reaction products which formed and controlled further diffusion through each sample were similar for all sample groups. Therefore, \(n\) was also consistent across the sample groups. However, there was a consistent decrease in the reaction order with decreasing precipitate size for the 900°C and 800°C sample groups (i.e., 900L, 900S, 800L, and 800S).
To allow for comparison with other high temperature oxidation experiments, the parabolic reaction rate was calculated for the steady state portion of the mass gain curves (see Fig. 18). The parabolic oxidation rate constant, \({K}_{p}^{{\prime \prime }}\), was determined through calculation of the best fit curve of the mass gain data using Eq. 4 [63]:
$$\text{S}.\text{M}.\text{G}. =\sqrt{{K}_{p}^{{\prime \prime }}\bullet t}+{m}_{0}$$
4
To more accurately capture the transition from linear to parabolic growth and asses the variability in the parabolic rate constant throughout the experiment, a sliding window analysis approach was used to measure \({K}_{p}^{{\prime }{\prime }}\) as a function of the starting time (\(t\)) of the data range. Through an iterative process the range of data used to estimate the fitting parameters, \({K}_{p}^{{\prime \prime }}\) and \({m}_{0}\), was reduced by increasing the starting time considered. For example, iteration 1 considered the mass gain over the entire time range 1 minute to 100 hours, iteration 2 considered only 2 minutes to 100 hours, iteration 3 considered 3 minutes to 100 hours, and so on until only the final portion of the mass gain curve was considered. During the multi-oxide growth process, the initial oxidation front does not conform to the steady-state, diffusion-controlled region until some time, \({t}_{p}\), whereupon the parabolic growth model becomes valid. The determination of \({t}_{p}\) is still qualitative, but the data suggests that the time to reach steady state oxidation is reduced with decreasing precipitate size. The values for \({K}_{p}^{{\prime }{\prime }}\) that are most comparable to other literature are in the flat regions from approximately 40 hours onward and steadily decrease with both increased precipitate volume fraction and decreased precipitate size. In general the HESAs \({K}_{P}^{{\prime }{\prime }}\)rates are on the order of \({10}^{-11} {\text{g}}^{2}/\text{c}{\text{m}}^{4}\text{s}\) which falls in line with many other MPEAs and Ni-based superalloys [11, 62, 64, 65].