Characterizing Short-Time Aging Precipitation Behavior of a Novel Nickel–Iron-Based Alloy via Electrical Performance

In this paper, the precipitation behavior and its effect on resistivity in a new type of nickel–iron-based alloy during short-term aging were investigated. During the aging process, the γ′ phase increases in average size and decreases in number, with its area fraction fluctuating over time. This fluctuation is caused by the mismatch in the redissolution and growth rates of the γ′ phase. As the area fraction of the γ′ phase increases, the content of solute atoms in the matrix that scatter electrons decreases, lowering the resistivity of the alloy. Additionally, the continuous precipitation of M23C6 at grain boundaries during aging causes the resistivity to gradually increase. This paper explains the fluctuation in the total amount of γ′ phase during short-term aging and proposes a new method for characterizing the precipitation behavior of the γ′ phase in the novel alloy using the relative trend of resistivity changes.


Introduction
Nickel-based and iron-nickel-based superalloys are widely utilized in advanced ultrasupercritical (A-USC) steam power generation technology due to their outstanding mechanical properties at high temperatures, including tensile strength, yield strength, and creep resistance [1][2][3][4].These properties are intrinsically linked to factors such as the size, volume fraction, and distribution of the precipitated phases [5][6][7][8][9].As one of the common strengthening phases in nickel-based alloys, the precipitation of γ ′ phases, including Ni 3 Al and Ni 3 Ti, significantly influences the mechanical properties of nickel-based alloys.For instance, with the progression of aging, the size of the γ ′ phases gradually increases, and in certain alloys the γ ′ phases may transition from a spherical to a cubic shape.These morphological changes can significantly influence the material's mechanical properties, including its tensile and fatigue behaviors [7,[10][11][12][13][14][15][16].Consequently, ideal properties can be achieved through the optimization of the γ ′ phase shape and distribution [17].The alloy studied in this research is a novel γ ′ precipitation-strengthened nickel-iron-based alloy.To ensure its long-term service reliability, it is necessary to investigate the microstructural evolution of this iron-nickel-based superalloys under different heat treatment conditions.
To characterize the morphology and size evolution of the precipitates in nickel-ironbased alloys during the aging process, current research methods include direct characterization techniques such as scanning electron microscopy (SEM) [18,19] and transmission electron microscopy (TEM) [20,21], as well as indirect characterization techniques like tensile tests [22,23], impact tests [24,25], creep tests [26,27], and differential scanning calorimetry (DSC) [28].Additionally, some researchers have used neutron diffraction to characterize the elastic lattice strains in the γ matrix, γ ′ phases, and γ ′′ phases during creep tests, stress relaxation tests, and stress-free aging experiments in Inconel 718 alloys [6].However, among these precipitation characterization techniques, SEM, TEM, and DTA have higher requirements for preparation.Besides, due to the tiny sample size used in TEM, it is challenging to obtain the overall characteristics of macroscopic precipitates.Meanwhile, it is difficult to conduct large-scale mechanical and thermodynamic experiments in actual production due to their long experimental cycles.Therefore, there is an urgent need for a characterization technique that has relatively low sample preparation requirements and a short experimental cycle, while remaining sensitive to microstructural changes.
The microstructure of the alloy, including precipitates and the matrix, changes during the aging process, which in turn affects the electrical properties.Hence, the resistivity test holds significant value for characterizing the precipitation evolution during the aging process.Tian et al. [29] found that the scattering effect of Fe solute atoms in the Cu matrix is the primary cause of changes in the electrical properties of the Cu-Fe composite during aging.The precipitation of Fe-rich phases leads to a decrease in the resistivity of the composite, while the re-dissolution of Fe atoms into the matrix increases the resistivity.Short-rangeordered (SRO) structures and solute atom-enriched regions, i.e., Guinier Preston zones (GP zones), also significantly impact the electrical properties [30][31][32].Lee, SH et al. [33] characterized the precipitation behavior of solute atoms in Al-Zn-Mg-Cu alloys during natural aging by monitoring changes in resistivity, observing a rapid rise in resistivity due to the formation of GP zones as natural aging-time increased.Zhang et al. [34] investigated the influence of Ag and Ge elements on the aging behavior of AA6061 alloys.It was found that Ge atoms more readily form secondary phases with Mg atoms, reducing the concentration of Mg atoms in the solid solution and thus the alloy's resistivity.In the overaging stage, the resistivity increased due to the aggregation of Ag atoms and their inhibitory effect on the ordering of precipitates.These studies indicate that the resistivity tests can be used to characterize the precipitation behavior of alloys during heat treatment.However, the resistivity tests are primarily applied to conductors and metal films with few element types, such as Cu and Al, and there is limited research on the characterization of the aging behavior of complex nickel-iron-based alloys.
This paper aims to investigate the relationship between the resistivity changes and the precipitate evolution in the new type of nickel-iron-based alloy during short-term aging.Through resistivity measurements and SEM observations, the variations in resistivity, size, quantity, and distribution of γ ′ phases at different aging times were analyzed.The precipitation of M 23 C 6 at grain boundaries was also investigated.Combining solute atom diffusion theory, the reasons for the fluctuating total precipitation amount of γ ′ phases were analyzed.Based on the insight into the precipitation behavior of nickel-iron-based alloy, a new method for determining the variation trend of γ ′ phase density during short-term aging was proposed.

Materials and Heat Treatment
The composition of the novel nickel-iron-based alloy was analyzed using the electron probe microanalyzer analyzer (EPMA, JXA8230, JEOL, Tokyo, Japan) with a surface scan performed at a voltage of 20 kV, as detailed in Table 1.Due to the necessity of eliminating the dendritic and precipitated phases inherent in the material before aging [9,35], a complete solution treatment was required prior to the aging process.Three 50 mm × 30 mm × 25 mm samples were solution-treated at 1080 • C. To avoid the oxidation and hardened layer, a uniform 2 mm-thickness layer was removed from all exposed sides.Then, samples were cut into 15 mm × 10 mm × 8 mm specimens.After resistivity measurement, the solutiontreated specimens were placed in a muffle furnace for aging at 750 • C, 775 • C, 800 • C, 825 • C, and 850 • C for 0.5 h to 20 h.The heat treatment process is illustrated in Figure 1.layer, a uniform 2 mm-thickness layer was removed from all exposed sides.Then, samples were cut into 15 mm × 10 mm × 8 mm specimens.After resistivity measurement, the solution-treated specimens were placed in a muffle furnace for aging at 750 °C, 775 °C, 800 °C 825 °C, and 850 °C for 0.5 h to 20 h.The heat treatment process is illustrated in Figure 1.

Electrical Resistivity Measurement
To characterize the changes in resistivity brought by aging treatment, the resistivity of each specimen was measured before and during the aging process.Before each resistivity measurement, the surfaces of the specimen were polished to become scratch-free to diminish the impact of the oxidation layer and surface roughness on the resistivity.A direct current resistance measuring instrument (QJ36, Shanghai, China) was used to measure and calculate the resistivity with the four-point method.To minimize the influence of probe-contact force on the measurement results, probes were soldered to the specimen surface using tin solder.To eliminate measurement errors caused by contact potential and thermoelectric potential, the positive and negative currents were measured 5 times, respectively, and the average value was taken as the resistivity of the sample.The relative change in resistivity was calculated using the following formula: where  and  represent the electrical resistivity of the same specimen in the solution-treated state and the aged state, respectively.A solution-treated specimen was selected as the control, which only underwent repeated measurements accompanying the testing specimens without the aging treatment.

Microstructure Characterization
(1) Scanning Electron Microscopy (SEM) Observation Specimens were cut into 15 mm × 10 mm × 3 mm dimensions for scanning electron microscopy (SEM) observation using an SU8220 field-emission SEM (Tokyo, Japan) at 5 kV, with their surfaces polished to become scratch-free.The specimens were electrolytic etched for 30 s under a current of 1.2 A in the etching solution mixed with 65 mL H PO + 15 mL H SO + 12 mLH O + 3 g H O + 5 g CrO .One grain was randomly selected from five different regions on the surface.The attributes of the γ′ phase were statistically analyzed in interior regions of each grain distant from grain boundaries and carbides.The aggregate of γ′ particles counted in each specimen was not less than 1.3 × 10 3 .The tota

Electrical Resistivity Measurement
To characterize the changes in resistivity brought by aging treatment, the resistivity of each specimen was measured before and during the aging process.Before each resistivity measurement, the surfaces of the specimen were polished to become scratch-free to diminish the impact of the oxidation layer and surface roughness on the resistivity.A direct current resistance measuring instrument (QJ36, Shanghai, China) was used to measure and calculate the resistivity with the four-point method.To minimize the influence of probe-contact force on the measurement results, probes were soldered to the specimen surface using tin solder.To eliminate measurement errors caused by contact potential and thermoelectric potential, the positive and negative currents were measured 5 times, respectively, and the average value was taken as the resistivity of the sample.The relative change in resistivity was calculated using the following formula: where ρ 0 and ρ Ageing represent the electrical resistivity of the same specimen in the solutiontreated state and the aged state, respectively.A solution-treated specimen was selected as the control, which only underwent repeated measurements accompanying the testing specimens without the aging treatment.

Microstructure Characterization
(1) Scanning Electron Microscopy (SEM) Observation Specimens were cut into 15 mm × 10 mm × 3 mm dimensions for scanning electron microscopy (SEM) observation using an SU8220 field-emission SEM (Tokyo, Japan) at 5 kV, with their surfaces polished to become scratch-free.The specimens were electrolytic etched for 30 s under a current of 1.2 A in the etching solution mixed with 65 mL H 3 PO 4 + 15 mL H 2 SO 4 + 12 mL H 8 O 3 + 3 g H 2 O + 5 g CrO 3 .One grain was randomly selected from five different regions on the surface.The attributes of the γ ′ phase were statistically analyzed in interior regions of each grain distant from grain boundaries and carbides.The aggregate of γ ′ particles counted in each specimen was not less than 1.3 × 10 3 .The total area of the γ ′ phase among each region was measured using ImageJ version 1.54f [36], and the total area fraction was calculated by: where A i represents the total area of the γ ′ phase in region i, and S i represents the total area of region i.
(2) Electron Backscatter Diffraction (EBSD) Tests EBSD tests were used to assess the effect of aging treatment on the microstructure.Specimens aged at 850 • C for 1 h and 10 h were mechanically polished, followed by ionbeam milling with an ion-beam mill (Leica EM TIC 3X, Wetzlar, Germany).During milling, the argon ion beam was accelerated at 5.5 kV with a current of 1.5 mA for a duration of 2.5 h.EBSD tests were conducted using a scanning electron microscope (TESCAN S9000X, Brno, Czech Republic) equipped with an EBSD detector.The electron beam was accelerated at 15 kV.A 250 µm × 250 µm region was scanned for phase distribution with a step size of 0.5 µm, and a 10 µm × 10 µm region was scanned in detail with a step size of 0.025 µm.
(3) X-ray Diffraction (XRD) Analysis XRD analysis was performed on the specimens aged at 850 • C using an X-ray diffractometer (Rigaku D/max-2500 diffractometer, Tokyo, Japan) with a scanning rate of 1 • /min from 2θ = 40 • to 125 • .A Cu target with its Kα1 wavelength equal to 1.542 Å was employed.A Cu target was employed with its Kα1 wavelength equal to 1.542 Å.

Resistivity Changes Induced by Aging Treatment
The measured resistivity of the novel nickel-iron-based alloy in the solution-treated state was approximately 55.36 ± 4.59 µΩ•cm.The resistivities of both specimens aged at 800 • C and 850 • C exhibited a fluctuating pattern in the ∆ρ/ρ 0 over time with a general upward trend as time increased, as illustrated in Figure 2. In comparison, the ∆ρ/ρ 0 of the control specimen varied within a tiny scale from −0.53% to 0.62%, indicating a fine accuracy and reliability of the resistivity measurement.
Materials 2024, 17, x FOR PEER REVIEW 4 of 21 area of the γ′ phase among each region was measured using ImageJ version 1.54f [36], and the total area fraction was calculated by: where  represents the total area of the γ′ phase in region , and  represents the total area of region .
(2) Electron Backscatter Diffraction (EBSD) Tests EBSD tests were used to assess the effect of aging treatment on the microstructure.Specimens aged at 850 °C for 1 h and 10 h were mechanically polished, followed by ionbeam milling with an ion-beam mill (Leica EM TIC 3X, Wetzlar, Germany).During milling, the argon ion beam was accelerated at 5.5 kV with a current of 1.5 mA for a duration of 2.5 h.EBSD tests were conducted using a scanning electron microscope (TESCAN S9000X, Brno, Czech Republic) equipped with an EBSD detector.The electron beam was accelerated at 15 kV.A 250 µm × 250 µm region was scanned for phase distribution with a step size of 0.5 µm, and a 10 µm × 10 µm region was scanned in detail with a step size of 0.025 µm.
(3) X-ray Diffraction (XRD) Analysis XRD analysis was performed on the specimens aged at 850 °C using an X-ray diffractometer (Rigaku D/max-2500 diffractometer, Tokyo, Japan) with a scanning rate of 1°/min from 2 = 40° to 125°.A Cu target with its Kα1 wavelength equal to 1.542 Å was employed.A Cu target was employed with its Kα1 wavelength equal to 1.542 Å.

Resistivity Changes Induced by Aging Treatment
The measured resistivity of the novel nickel-iron-based alloy in the solution-treated state was approximately 55.36 ± 4.59 μΩ • cm.The resistivities of both specimens aged at 800 °C and 850 °C exhibited a fluctuating pattern in the ∆/ over time with a general upward trend as time increased, as illustrated in Figure 2. In comparison, the ∆/ of the control specimen varied within a tiny scale from −0.53% to 0.62%, indicating a fine accuracy and reliability of the resistivity measurement.

Evolution of 𝛾′ Precipitate Morphology
Figure 3 and Figure 4 show the morphology of the γ′ phase in the solution-treated alloy after aging at 850 °C and 800 °C over time, respectively, and the average diameter of the γ′ phase is given in Figure 5a.The results indicate that, after aging for the same

Evolution of γ ′ Precipitate Morphology
Figures 3 and 4 show the morphology of the γ ′ phase in the solution-treated alloy after aging at 850 • C and 800 • C over time, respectively, and the average diameter of the γ ′ phase is given in Figure 5a.The results indicate that, after aging for the same duration, the particle diameters at 850 • C were larger than those at 800 • C.During aging from 1 h to 20 h at 850 • C, the average particle radius of the γ ′ phase increased from 23.17 nm to 51.40 nm.At 800 • C, the aging process from 0.5 h to 20 h saw the average particle radius of the γ ′ phase grow from 14.77 nm to 34.34 nm.The γ ′ phase maintained a spherical morphology throughout the aging process at both 850 • C and 800 • C.
Materials 2024, 17, x FOR PEER REVIEW 5 of 21 duration, the particle diameters at 850 °C were larger than those at 800 °C.During aging from 1 h to 20 h at 850 °C, the average particle radius of the γ′ phase increased from 23.17 nm to 51.40 nm.At 800 °C, the aging process from 0.5 h to 20 h saw the average particle radius of the γ′ phase grow from 14.77 nm to 34.34 nm.The γ′ phase maintained a spherical morphology throughout the aging process at both 850 °C and 800 °C.duration, the particle diameters at 850 °C were larger than those at 800 °C.During aging from 1 h to 20 h at 850 °C, the average particle radius of the γ′ phase increased from 23.17 nm to 51.40 nm.At 800 °C, the aging process from 0.5 h to 20 h saw the average particle radius of the γ′ phase grow from 14.77 nm to 34.34 nm.The γ′ phase maintained a spherical morphology throughout the aging process at both 850 °C and 800 °C.As shown in Figure 5b, the number of γ′ phases per unit area showed an overall decreasing trend.At the same aging duration, the number of γ′ phases per unit area at 800 °C was significantly higher than at 850 °C.This is because at higher temperatures, the diffusion rates of solute atoms and vacancies are greater, resulting in a faster dissolution rate of smaller-sized γ′ phases.Therefore, the higher the temperature, the lower the particle density at the same aging time.Considering the increase in the particle size and the overall decrease in the particle number over time, it can be inferred that that both precipitation and Ostwald ripening phenomena occurred during the short-term aging process.In other words, smaller particles dissolved, releasing solute atoms that then re-precipitate onto larger particles, which leads to an increase in the average particle size.As shown in Figure 5b, the number of γ ′ phases per unit area showed an overall decreasing trend.At the same aging duration, the number of γ ′ phases per unit area at 800 • C was significantly higher than at 850 • C.This is because at higher temperatures, the diffusion rates of solute atoms and vacancies are greater, resulting in a faster dissolution rate of smaller-sized γ ′ phases.Therefore, the higher the temperature, the lower the particle density at the same aging time.As shown in Figure 5b, the number of γ′ phases per unit area showed an overall decreasing trend.At the same aging duration, the number of γ′ phases per unit area at 800 °C was significantly higher than at 850 °C.This is because at higher temperatures, the diffusion rates of solute atoms and vacancies are greater, resulting in a faster dissolution rate of smaller-sized γ′ phases.Therefore, the higher the temperature, the lower the particle density at the same aging time.Considering the increase in the particle size and the overall decrease in the particle number over time, it can be inferred that that both precipitation and Ostwald ripening phenomena occurred during the short-term aging process.In other words, smaller particles dissolved, releasing solute atoms that then re-precipitate onto larger particles, which leads to an increase in the average particle size.Considering the increase in the particle size and the overall decrease in the particle number over time, it can be inferred that that both precipitation and Ostwald ripening phenomena occurred during the short-term aging process.In other words, smaller particles dissolved, releasing solute atoms that then re-precipitate onto larger particles, which leads to an increase in the average particle size.

Evolution of Carbide with Aging Time
M C is also one of the main strengthening phases of the novel nickel-iron-based alloy.As shown in Figure 7, the precipitation of M C at the grain boundaries gradually increases in number and size with the aging time at 850 °C.The morphology of these carbides also transformed from a particulate to a continuous cellular morphology along the grain boundaries.Additionally, the decomposition of the carbide (MC + γ → M C + γ′) was observed at 850 °C for 4 h and 10 h, which is consistent with phenomena observed by other researchers [8,9,37].Consequently, finer γ′ particles were formed around the MC carbides.As shown in Figure 8, when aged at 850 °C for 10 h, the size of γ′ particles near MC was 15 nm, whereas the average size of γ′ particles distant from grain boundaries and MC was 41.85 nm.

Evolution of Carbide with Aging Time
M 23 C 6 is also one of the main strengthening phases of the novel nickel-iron-based alloy.As shown in Figure 7, the precipitation of M 23 C 6 at the grain boundaries gradually increases in number and size with the aging time at 850 • C. The morphology of these carbides also transformed from a particulate to a continuous cellular morphology along the grain boundaries.Additionally, the decomposition of the carbide ( MC + γ → M 23 C 6 + γ ′ ) was observed at 850 • C for 4 h and 10 h, which is consistent with phenomena observed by other researchers [8,9,37].Consequently, finer γ ′ particles were formed around the MC carbides.As shown in Figure 8, when aged at 850 • C for 10 h, the size of γ ′ particles near MC was 15 nm, whereas the average size of γ ′ particles distant from grain boundaries and MC was 41.85 nm.

Evolution of Carbide with Aging Time
M C is also one of the main strengthening phases of the novel nickel-iron-based alloy.As shown in Figure 7, the precipitation of M C at the grain boundaries gradually increases in number and size with the aging time at 850 °C.The morphology of these carbides also transformed from a particulate to a continuous cellular morphology along the grain boundaries.Additionally, the decomposition of the carbide (MC + γ → M C + γ′) was observed at 850 °C for 4 h and 10 h, which is consistent with phenomena observed by other researchers [8,9,37].Consequently, finer γ′ particles were formed around the MC carbides.As shown in Figure 8, when aged at 850 °C for 10 h, the size of γ′ particles near MC was 15 nm, whereas the average size of γ′ particles distant from grain boundaries and MC was 41.85 nm.

Causes of the Resistivity Fluctuation
Several factors can influence electrical properties of the alloy, including the grain boundary behavior, dislocations, and precipitates [38].In Huang's resea HT700P alloy, no significant change in grain morphology and size was found d aging process [9].Unlike cubic γ′ phases, spherical γ′ phases exhibit a sma misfit with the matrix [39], thus the scattering effect of phase boundaries on ele be disregarded.G. Kelekanjeri's study on Waspaloy alloy suggested that, in th tion stage, the influence of dislocations, vacancies, and other factors on resis subtle compared with the solute concentration in the matrix [40,41].The solute c tion in the matrix is closely related to precipitation behavior.Therefore, the fl and upward trend of the ∆/ −  curve during the aging process of the nov iron-based alloy can be explained from two aspects, namely the precipitation and grain boundary behavior.
Roebuck's study [42] on CMSX4 single-crystal, nickel-based alloy shows t sistivity of the precipitated γ′ phase is higher than that of the γ matrix, resu positive correlation between the resistivity of this alloy and the change in the tot of precipitation [43].However, the resistivity of the alloy in the solution-treated measured to be 55.36 ± 4.59 μΩ • cm, and studies have shown that the resistivit phase is about 44.2 ± 0.5 μΩ • cm [44], which is lower than that of the γ phase.ally, change in the volume fraction of the precipitates should also be considered cipitation of the γ′ phase can decrease the solute concentration in the γ phase tionship between the resistivity of the γ phase and the solute concentration ca trated by the following equation [45][46][47]: where  is the resistivity of the alloy matrix without solute atoms,  is the caused by the solid solution atoms,  is the atomic fraction of the solute ele ∆ / represents the change in resistivity caused by a sole solute atom.This that the resistivity of the γ phase is positively correlated with the solute conc As a result, the precipitation of γ′ phases can decrease the overall resistivity of Taking the fluctuation of resistivity and the fluctuation of γ′ phase area 850 °C as an example, the corresponding relationship between resistivity and ar fluctuations can be analyzed based on Figure 9a.
(1) During the stages of 2-4 h, 5-10 h, and 15-20 h, the area fraction of the decreased, leading to the redissolution of solute atoms into the matrix.Hence, concentration in the matrix increased.On one hand, the reduction in the total a

Causes of the Resistivity Fluctuation
Several factors can influence electrical properties of the alloy, including the grain size, grain boundary behavior, dislocations, and precipitates [38].In Huang's research on the HT700P alloy, no significant change in grain morphology and size was found during the aging process [9].Unlike cubic γ ′ phases, spherical γ ′ phases exhibit a smaller lattice misfit with the matrix [39], thus the scattering effect of phase boundaries on electrons can be disregarded.G. Kelekanjeri's study on Waspaloy alloy suggested that, in the maturation stage, the influence of dislocations, vacancies, and other factors on resistivity was subtle compared with the solute concentration in the matrix [40,41].The solute concentration in the matrix is closely related to precipitation behavior.Therefore, the fluctuation and upward trend of the ∆ρ/ρ 0 − t curve during the aging process of the novel nickel-ironbased alloy can be explained from two aspects, namely the precipitation behavior and grain boundary behavior.
Roebuck's study [42] on CMSX4 single-crystal, nickel-based alloy shows that the resistivity of the precipitated γ ′ phase is higher than that of the γ matrix, resulting in a positive correlation between the resistivity of this alloy and the change in the total amount of precipitation [43].However, the resistivity of the alloy in the solution-treated state was measured to be 55.36 ± 4.59 µΩ•cm, and studies have shown that the resistivity of the γ ′ phase is about 44.2 ± 0.5 µΩ•cm [44], which is lower than that of the γ phase.Additionally, change in the volume fraction of the precipitates should also be considered.The precipitation of the γ ′ phase can decrease the solute concentration in the γ phase.The relationship between the resistivity of the γ phase and the solute concentration can be illustrated by the following equation [45][46][47]: where ρ 0 is the resistivity of the alloy matrix without solute atoms, ρ imp is the resistivity caused by the solid solution atoms, x imp is the atomic fraction of the solute element, and ∆ρ imp /a represents the change in resistivity caused by a sole solute atom.This indicates that the resistivity of the γ phase is positively correlated with the solute concentration.As a result, the precipitation of γ ′ phases can decrease the overall resistivity of the alloy.
Taking the fluctuation of resistivity and the fluctuation of γ ′ phase area fraction at 850 • C as an example, the corresponding relationship between resistivity and area fraction fluctuations can be analyzed based on Figure 9a.
(2) During the stages of 1-2 h, 4-5 h, and 10-15 h, the area fraction of the γ′ phase increased.The increased amount of γ′ phases led to an overall reduction in resistivity.In addition, the decrease in solute concentration within the matrix and the reduction in lattice distortion helped restore the periodicity of the lattice in the matrix.According to Bloch's theorem, the movement of electrons will experience reduced energy consumption over a larger lattice range, thereby lowering the resistivity.As for the upward trend in resistivity, research has shown that the pinning of dislocations by M C carbides leads to the enrichment of carbon atoms, thereby increasing the resistivity [49].The kernel average misorientation (KAM) value from EBSD tests is deemed to be positively correlated with the dislocation density [50].To be specific, the geometrically necessary dislocation (GND) density of the sample can be calculated using the following formula [51]: where  represents the Burgers vector.∆ is the scanning step size. is the KAM angle.As shown in Figure 10a,b, after being aged at 850 °C for 1 h, no M C phase was observed in the phase map, and the KAM value near the grain boundaries distributed uniformly.However, as shown in Figure 10c,d, after 10 h of aging, Cr C carbides precipitated at grain boundaries.Higher KAM values could be observed near the carbides, indicating a greater dislocation density and higher degree of lattice distortion around the carbides, compared to other regions.(1) During the stages of 2-4 h, 5-10 h, and 15-20 h, the area fraction of the γ ′ phase decreased, leading to the redissolution of solute atoms into the matrix.Hence, the solute concentration in the matrix increased.On one hand, the reduction in the total amount of γ ′ phases, which have lower resistivity, resulted in an overall increase in resistivity.On the other hand, the outer electrons of Al and Ti atoms tended to overlap with the 3d orbitals of transition elements including Ni, Cr, and Fe, and the redistribution of the electrons could form short-range-ordered structures similar to intermetallic compounds [48], increasing the resistivity of the alloy [30].Additionally, the increased solute concentration in the matrix led to a greater lattice distortion, disrupting the periodicity of electron scattering, thus hindering electron scattering and increasing the resistivity.
(2) During the stages of 1-2 h, 4-5 h, and 10-15 h, the area fraction of the γ ′ phase increased.The increased amount of γ ′ phases led to an overall reduction in resistivity.In addition, the decrease in solute concentration within the matrix and the reduction in lattice distortion helped restore the periodicity of the lattice in the matrix.According to Bloch's theorem, the movement of electrons will experience reduced energy consumption over a larger lattice range, thereby lowering the resistivity.
As for the upward trend in resistivity, research has shown that the pinning of dislocations by M 23 C 6 carbides leads to the enrichment of carbon atoms, thereby increasing the resistivity [49].The kernel average misorientation (KAM) value from EBSD tests is deemed to be positively correlated with the dislocation density [50].To be specific, the geometrically necessary dislocation (GND) density of the sample can be calculated using the following formula [51]: where b represents the Burgers vector.∆d is the scanning step size.θ KAM is the KAM angle.As shown in Figure 10a,b, after being aged at 850 • C for 1 h, no M 23 C 6 phase was observed in the phase map, and the KAM value near the grain boundaries distributed uniformly.However, as shown in Figure 10c,d, after 10 h of aging, Cr 23 C 6 carbides precipitated at grain boundaries.Higher KAM values could be observed near the carbides, indicating a greater dislocation density and higher degree of lattice distortion around the carbides, compared to other regions.The precipitation of the γ′ phase can cause lattice contraction, facilitating the dislocation multiplication.Strengthening phases such as M C carbides can pin these dislocations.As shown in Figure 7, as aging time was prolonged, there was an increasing amount of M C precipitates at the grain boundaries, leading to an increased number of pinned dislocations at these sites.This resulted in severe lattice distortion around the M C phases, attracting more solute atoms and causing an increase in resistivity at the grain boundaries.

The Short-Term Aging Coarsening Kinetics of the 𝛾′ Phase
The frequency distribution histograms for γ′ phases of different sizes during aging at 850 °C are presented in Figure 11.The calculation formula for the skewness coefficient (SK) is as follows: where,  denotes the mean of the sample,  represents the median of the sample, and  signifies standard deviation of the sample.When SK = 0, the particle diameter follows a normal distribution; SK > 0 indicates a higher proportion of larger-sized γ′ phases, and vice versa.With the aging time prolongation, SK changed from positive to negative, The precipitation of the γ ′ phase can cause lattice contraction, facilitating the dislocation multiplication.Strengthening phases such as M 23 C 6 carbides can pin these dislocations.As shown in Figure 7, as aging time was prolonged, there was an increasing amount of M 23 C 6 precipitates at the grain boundaries, leading to an increased number of pinned dislocations at these sites.This resulted in severe lattice distortion around the M 23 C 6 phases, attracting more solute atoms and causing an increase in resistivity at the grain boundaries.

The Short-Term Aging Coarsening Kinetics of the γ ′ Phase
The frequency distribution histograms for γ ′ phases of different sizes during aging at 850 • C are presented in Figure 11.The calculation formula for the skewness coefficient (SK) is as follows: where, µ D denotes the mean of the sample, M represents the median of the sample, and σ signifies standard deviation of the sample.When SK = 0, the particle diameter follows a normal distribution; SK > 0 indicates a higher proportion of larger-sized γ ′ phases, and vice versa.With the aging time prolongation, SK changed from positive to negative, indicating that small size γ ′ particles gradually dissolved, and the proportion of larger-sized γ ′ particles increased.This indicates that ripening occurs during short-term aging.indicating that small size γ′ particles gradually dissolved, and the proportion of largersized γ′ particles increased.This indicates that ripening occurs during short-term aging.The classical formula for particle size evolution based on the LSW theory is given by ̅ −  = , where ̅ represents the average radius of the γ′ phase,  is the critical nucleus radius, and  is the coarsening rate constant.Although the classical LSW theory assumes that the precipitated phase is in equilibrium [52][53][54][55], the results indicate that the cube of the γ′ phase radius presents a relatively high fitting goodness, as shown in Figure 12.This suggests that the particle size of the γ′ phase during short-term aging can be predicted by the LSW theory.Meanwhile, the growth rates of the γ′ phase during short-term aging at 850 °C and 800 °C were 986.25 nm /h and 262.30nm /h, respectively, showing that a higher temperature results in a greater growth rate.The classical formula for particle size evolution based on the LSW theory is given by r 3 − r 3 0 = kt, where r represents the average radius of the γ ′ phase, r 0 is the critical nucleus radius, and k is the coarsening rate constant.Although the classical LSW theory assumes that the precipitated phase is in equilibrium [52][53][54][55], the results indicate that the cube of the γ ′ phase radius presents a relatively high fitting goodness, as shown in Figure 12.This suggests that the particle size of the γ ′ phase during short-term aging can be predicted by the LSW theory.Meanwhile, the growth rates of the γ ′ phase during short-term aging at 850 • C and 800 • C were 986.25 nm 3 /h and 262.30nm 3 /h, respectively, showing that a higher temperature results in a greater growth rate.
The activation energy for the growth of the γ ′ phase can be calculated using the equation ln(kT) = C − Q/RT, where k is the growth rate, T is the absolute temperature, C is a constant, and R = 8.314 J/(mol•K) is the ideal gas constant.The calculated coarsening activation energy for the γ ′ phase was 277.59 kJ/mol.In the Ni matrix, the diffusion activation energy of Al was approximately 260 kJ/mol [56]~269.2kJ/mol [57], and that of Ti in the Ni matrix was 257 kJ/mol [58,59].Meanwhile, the diffusion activation energies of Ni and Fe in the Ni matrix were 284 kJ/mol [60,61] and 281 kJ/mol [60], respectively.Therefore, the solid-state diffusion capabilities of Al and Ti atoms are greater than those of Fe and Ni atoms.In this new type of nickel-iron-based superalloy, the coarsening activation energy of the γ ′ phase is slightly higher than that of Al and Ti, but lower than that of Ni and Fe.Therefore, under the given temperature and time conditions, the bulk diffusion of Al and Ti elements is more feasible, indicating that the growth of the γ ′ phase was primarily achieved through the bulk diffusion of Al and Ti elements, and the coarsening activation energy was dominated by the mole fractions of Al and Ti atoms in the matrix.
The classical formula for particle size evolution based on the LSW theory is given by ̅ −  = , where ̅ represents the average radius of the γ′ phase,  is the critical nucleus radius, and  is the coarsening rate constant.Although the classical LSW theory assumes that the precipitated phase is in equilibrium [52][53][54][55], the results indicate that the cube of the γ′ phase radius presents a relatively high fitting goodness, as shown in Figure 12.This suggests that the particle size of the γ′ phase during short-term aging can be predicted by the LSW theory.Meanwhile, the growth rates of the γ′ phase during short-term aging at 850 °C and 800 °C were 986.25 nm /h and 262.30nm /h, respectively, showing that a higher temperature results in a greater growth rate.

Analysis of the Fluctuation in the Total Amount of γ ′ Phase Precipitation
The transformation process of γ ′ phase is driven by the minimization of Gibbs free energy.In this process, the driving force is the difference in chemical free energy between the matrix and γ ′ phases, and the resistance is the interfacial energy and elastic strain energy.The fluctuations in the area fraction of the γ ′ phases are inherently the result of the competition between the driving and resisting forces of the phase transformation.The competition results in the rates of precipitation, redissolution, and growth during the aging process not being stable and being in competition with each other.
In the first stage of aging, the chemical free energy of the γ ′ phase was lower than that of the γ matrix.To minimize the Gibbs free energy, the γ ′ phase began to precipitate and grow, reducing the solute concentration in the matrix and increasing the total amount of the γ ′ phase.Ostwald ripening occurred in the second stage of aging.To minimize the total interfacial energy, the solute concentration at the surfaces of different particles varied due to differences in chemical potential, leading to strong interactions between γ ′ phases of different sizes.The interactions between particles involved the dissolution of smaller-sized γ ′ phases into the matrix, as well as the growth of larger-sized γ ′ phases via absorbing solute atoms from the matrix.However, the rates of two such processes may differ, which mainly accounts for the fluctuation in the area fraction during the aging process.
To explain the changes in solute concentration within the matrix, two types of equilibrium concentrations were first defined: ➀ At a specific moment, t 0 , during a very short aging period, i.e., when ∆t| t=t 0 → 0 , the solute concentration at a specific site in the matrix remained almost constant during the exchange of solute atoms between the γ ′ particles and the matrix.This concentration is defined as C e , representing the equilibrium solute concentration in the matrix.➁ When the aging time t → ∞ , the precipitation process reached equilibrium, and the equilibrium concentration of solute atoms in the matrix at this point is defined as C ∞ .The following assumptions are made specific to the short-term aging: ➀ the distribution of γ ′ phases in the matrix was random; ➁ the equilibrium concentration of solute atoms, C e , was equal among all particles in the matrix.
(1) Stage I: the precipitation and growth of the γ ′ phase.
The solute concentration distribution within the matrix in the solid solution state during the initial stage of aging is depicted in Figure 13.In this stage, γ ′ phases precipitated and grew rapidly.At this point, the equilibrium concentration C e within the γ matrix was higher than the solute concentration in the matrix near the surface of the γ ′ phase.Solute atoms diffused from regions of high concentration to the surface of the γ ′ phase with lower concentration, causing a decrease in the solute concentration in the matrix.The growth of the γ ′ phase led to an increase in the area fraction.
Solute atoms diffused from regions of high concentration to the surface of the γ′ phase with lower concentration, causing a decrease in the solute concentration in the matrix.The growth of the γ′ phase led to an increase in the area fraction.(2) Stage II: the redissolution of smaller-sized γ′ particles and growth of larger-sized γ′ particles.
As shown in Figure 14a, when the equilibrium concentration  in the matrix decreased to  ,  <  <  , solute atoms at the surface of smaller-sized γ′ particles dissolved back into the matrix, while larger-sized γ′ particles continued to grow.To compare the rates of redissolution and growth, the following definitions are made: the radius of the smaller-sized γ′ particles is  , the solute atom concentration in the matrix near the surface of the smaller-sized particles is  , the radius of the larger-sized γ′ particles is  , and the solute atom concentration in the matrix near the surface of the larger-sized γ′ particles is  .Since the driving force for Ostwald ripening is the difference in chemical potential between the smaller-sized and larger-sized γ′ particles, the chemical potential at the equilibrium concentration,  , is defined as  , the chemical potential of solute atoms near the surface of the larger-sized γ′ particles in the matrix is defined as  , and the chemical potential of solute atoms near the surface of the smaller-sized γ′ particles in the matrix is defined as  .According to Raoult's law: where  is the chemical potential of component  at a concentration of 100%,  = 8.314 J/(mol • K) is the ideal gas constant,  is the absolute temperature, and  is the concentration of component .
We assume that the chemical driving force for the dissolution of small-sized γ′ and the growth of large-sized γ′ is identical, i.e.,  −  =  −  .In this case, the equilibrium concentration between the two particles is  =  •  .
According to the spherical symmetry steady-state diffusion equation: where is the flow rate of solute atoms through the spherical shell of the γ′ particle with a radius of  due to diffusion,  is the volume of the new phase formed by a mole of solute atoms, and  is the diffusion coefficient of solute atoms in the matrix near the surface of the γ′ particle with a radius of  .Since = 2 • , Equation ( 8) can be transformed into: (2) Stage II: the redissolution of smaller-sized γ ′ particles and growth of larger-sized γ ′ particles.
As shown in Figure 14a, when the equilibrium concentration C e in the matrix decreased to C 2 , C 3 < C e < C 1 , solute atoms at the surface of smaller-sized γ ′ particles dissolved back into the matrix, while larger-sized γ ′ particles continued to grow.To compare the rates of redissolution and growth, the following definitions are made: the radius of the smaller-sized γ ′ particles is r S , the solute atom concentration in the matrix near the surface of the smaller-sized particles is C S , the radius of the larger-sized γ ′ particles is r B , and the solute atom concentration in the matrix near the surface of the larger-sized γ ′ particles is C B .Since the driving force for Ostwald ripening is the difference in chemical potential between the smaller-sized and larger-sized γ ′ particles, the chemical potential at the equilibrium concentration, C e , is defined as µ e , the chemical potential of solute atoms near the surface of the larger-sized γ ′ particles in the matrix is defined as µ B , and the chemical potential of solute atoms near the surface of the smaller-sized γ ′ particles in the matrix is defined as µ S .According to Raoult's law: where µ 0 i is the chemical potential of component i at a concentration of 100%, R = 8.314 J/(mol•K) is the ideal gas constant, T is the absolute temperature, and C i is the concentration of component i.
Considering that the concentration at the particle surface is the smaller the radius of the γ′ particle, the higher the solute concentration is in the surrounding matrix, i.e.,  <  <  .Combining Equation ( 10), it can be found that > 0, indicating that the larger-sized γ′ particles continue to grow, and < 0, indicating that the smaller-sized γ′ particles will dissolve back into the matrix.The total area change of the large and small particles during the competitive growth process is given by: The diffusion coefficient,  , of solute atoms in the matrix is concentration-dependent.Within the temperature range of 800 °C to 900 °C, when the Al content in the Fe-based matrix was less than 20%, the interdiffusion coefficient  increased with the concentration of Al solute atoms [63].Similarly, the diffusion coefficient of Al in Ni-based matrices also increased with the concentration of Al [64].This implies that  >  and   −   < 0 .That is, for two adjacent γ′ particles of different sizes, when competitive growth occurs, it can be found from Equation ( 11) that < 0, meaning that when competitive growth occurs, the dissolution rate of the smaller-sized γ′ exceeds the growth rate of the larger-sized γ′ under the same chemical driving force.As a result, the dissolution of the smaller-sized γ′ dominated at this stage, leading to a decrease in the area fraction.Such a process maintained until the moment when the smaller-sized γ′ phase disappeared, as shown in Figure 14b.The dissolution rate fell to 0, while residual solute atoms kept contributing to the growth of the adjacent γ′ phase.In other words, Figure 14b represents the inflection point where the area fraction changed from decreasing to increasing.(3) Stage III: a new loop of growth and redissolution process of the γ′ phase.
As shown in Figure 14b, with the complete dissolution of the smaller-sized γ′ particles, the interaction between adjacent γ′ particles transformed from γ′ ↔ γ′ and γ′ ↔ γ′ to γ′ ↔ γ′ .At this point, the solute concentration in the matrix reached a maximum.Subsequently, the process shown in Figure 15a began, where solute atoms diffused from high-concentration areas to low-concentration areas.The growth of γ′ phases led to an increase in the area fraction and a decrease in the solute concentration in the matrix.
After a period of aging, the equilibrium concentration,  , in the matrix further decreased to  <  <  , when a new loop of competitive growth occurred.The redissolution rate was still greater than the growth rate, leading to a continuous fluctuation in the area fraction.We assume that the chemical driving force for the dissolution of small-sized γ ′ and the growth of large-sized γ ′ is identical, i.e., µ S − µ e = µ e − µ B .In this case, the equilibrium concentration between the two particles is C e = √ C S •C B .According to the spherical symmetry steady-state diffusion equation: where dm i dt is the flow rate of solute atoms through the spherical shell of the γ ′ particle with a radius of r i due to diffusion, Ω is the volume of the new phase formed by a mole of solute atoms, and D i is the diffusion coefficient of solute atoms in the matrix near the surface of the γ ′ particle with a radius of r i .Since dA dt = 2πr• dr dt , Equation ( 8) can be transformed into: Considering that the concentration at the particle surface is RTr i [62], the smaller the radius of the γ ′ particle, the higher the solute concentration is in the surrounding matrix, i.e., C B < C e < C S .Combining Equation (10), it can be found that dA B dt > 0, indicating that the larger-sized γ ′ particles continue to grow, and dA S dt < 0, indicating that the smallersized γ ′ particles will dissolve back into the matrix.The total area change of the large and small particles during the competitive growth process is given by: The diffusion coefficient, D i , of solute atoms in the matrix is concentration-dependent.Within the temperature range of 800 • C to 900 • C, when the Al content in the Fe-based matrix was less than 20%, the interdiffusion coefficient D increased with the concentration of Al solute atoms [63].Similarly, the diffusion coefficient of Al in Ni-based matrices also increased with the concentration of Al [64].This implies that D S > D B and That is, for two adjacent γ ′ particles of different sizes, when competitive growth occurs, it can be found from Equation ( 11) that dA dt < 0, meaning that when competitive growth occurs, the dissolution rate of the smaller-sized γ ′ exceeds the growth rate of the larger-sized γ ′ under the same chemical driving force.As a result, the dissolution of the smaller-sized γ ′ dominated at this stage, leading to a decrease in the area fraction.Such a process maintained until the moment when the smaller-sized γ ′ phase disappeared, as shown in Figure 14b.The dissolution rate fell to 0, while residual solute atoms kept contributing to the growth of the adjacent γ ′ phase.In other words, Figure 14b represents the inflection point where the area fraction changed from decreasing to increasing.
(3) Stage III: a new loop of growth and redissolution process of the γ ′ phase.
As shown in Figure 14b, with the complete dissolution of the smaller-sized γ ′ particles, the interaction between adjacent γ ′ particles transformed from γ At this point, the solute concentration in the matrix reached a maximum.Subsequently, the process shown in Figure 15a began, where solute atoms diffused from high-concentration areas to low-concentration areas.The growth of γ ′ phases led to an increase in the area fraction and a decrease in the solute concentration in the matrix.In summary, during aging for 1-2 h, 4-5 h, and 10-15 h at 850 °C, and during aging for 0.5-1 h, 1.5-2 h, and 7-15 h at 800 °C, the area fraction showed an upward trend, indicating that the growth of the larger-sized γ′ phase dominated in these periods, and the solute concentration in the matrix decreased.In contrast, during aging for 2-4 h, 5-10 h, and 15-20 h at 850 °C, and during aging for 1-1.5 h, 2-7 h, and 15-20 h at 800 °C, the area After a period of aging, the equilibrium concentration, C e , in the matrix further decreased to C 3 < C e < C 2 , when a new loop of competitive growth occurred.The redissolution rate was still greater than the growth rate, leading to a continuous fluctuation in the area fraction.
In summary, during aging for 1-2 h, 4-5 h, and 10-15 h at 850 • C, and during aging for 0.5-1 h, 1.5-2 h, and 7-15 h at 800 • C, the area fraction showed an upward trend, indicating that the growth of the larger-sized γ ′ phase dominated in these periods, and the solute concentration in the matrix decreased.In contrast, during aging for 2-4 h, 5-10 h, and 15-20 h at 850 • C, and during aging for 1-1.5 h, 2-7 h, and 15-20 h at 800 • C, the area fraction exhibited a downward trend, suggesting that the redissolution of the smaller-sized γ ′ phase was the predominant process, and the solute concentration in the matrix increased.
It is known that the lattice constant of the γ ′ phase was smaller compared to the γ matrix.Additionally, the decrease in solute atom concentration in the γ matrix led to a reduction in the degree of lattice distortion in the γ matrix.Therefore, when the total amount of the γ ′ phase increased, the lattice constant of the alloy decreased, corresponding to a shift in the diffraction peak to a higher degree, and vice versa.The change in the average lattice constant can be taken as a reflection of the fluctuation in the total amount of the γ ′ phase and, as shown in Figure 16, the lattice constant increased during aging for 2-4 h and 5-10 h at 850 • C, corresponding to a decrease in the total amount of the γ ′ phase.In contrast, the lattice constant decreased during aging for 0-2 h, 4-5 h, and 10-15 h, indicating an increase in the total amount of the γ ′ phase.The fluctuation in the lattice constant obtained by XRD was consistent with the SEM results and the fluctuation pattern of resistivity.In summary, during aging for 1-2 h, 4-5 h, and 10-15 h at 850 °C, and during aging for 0.5-1 h, 1.5-2 h, and 7-15 h at 800 °C, the area fraction showed an upward trend, indicating that the growth of the larger-sized γ′ phase dominated in these periods, and the solute concentration in the matrix decreased.In contrast, during aging for 2-4 h, 5-10 h, and 15-20 h at 850 °C, and during aging for 1-1.5 h, 2-7 h, and 15-20 h at 800 °C, the area fraction exhibited a downward trend, suggesting that the redissolution of the smallersized γ′ phase was the predominant process, and the solute concentration in the matrix increased.
It is known that the lattice constant of the γ′ phase was smaller compared to the γ matrix.Additionally, the decrease in solute atom concentration in the γ matrix led to a reduction in the degree of lattice distortion in the γ matrix.Therefore, when the total amount of the γ′ phase increased, the lattice constant of the alloy decreased, corresponding to a shift in the diffraction peak to a higher degree, and vice versa.The change in the average lattice constant can be taken as a reflection of the fluctuation in the total amount of the γ′ phase and, as shown in Figure 16, the lattice constant increased during aging for 2-4 h and 5-10 h at 850 °C, corresponding to a decrease in the total amount of the γ′ phase.In contrast, the lattice constant decreased during aging for 0-2 h, 4-5 h, and 10-15 h, indicating an increase in the total amount of the γ′ phase.The fluctuation in the lattice constant obtained by XRD was consistent with the SEM results and the fluctuation pattern of resistivity.

Characterization of Resistivity and Precipitation during Aging at Other Temperatures
According to the prior analysis, an increase in ∆/ indicates that the redissolution of small-sized γ′ particles was the dominant process, leading to a decrease in the total

Characterization of Resistivity and Precipitation during Aging at Other Temperatures
According to the prior analysis, an increase in ∆ρ/ρ 0 indicates that the redissolution of small-sized γ ′ particles was the dominant process, leading to a decrease in the total amount of the γ ′ phase.Conversely, a decrease in ∆ρ/ρ 0 suggests that the growth of the γ ′ phase is the dominant process, resulting in an increase in the total amount of the γ ′ phase.Therefore, the inflection points in the ∆ρ/ρ 0 − t curve can be regarded as the transition points between alternate dominance of the redissolution and growth processes.Figure 17a displays the ∆ρ/ρ 0 − t curves of the novel nickel-iron-based alloy aged at 750 • C, 775 • C, 800 • C, 825 • C, and 850 • C for various durations, which all exhibited fluctuation.Four curves can be obtained by connecting the 1st, 2nd, 3rd, and 4th extreme values in order of temperature, as shown in Figure 17b.
"minimum total content of γ′ phase curve", corresponding to the second and fourth curves from the left in Figure 17b.These curves can be used to determine whether the total content of the γ′ phase increases or decreases during a specific aging period at the given temperature.For example, at 825 °C the total content of the γ′ phase decreased during aging for 0.5-1 h, increased during aging for 1-2 h, and then decreased again during aging for 2-7 h.As observed in Figure 17b, there is a general trend of increasing interval between the extrema curves of the total content of the γ′ phase.In line with the prior analysis, the fluctuation in the total volume of the γ′ phase is attributed to the alternating dominance of the redissolution process of small-sized γ′ phases and the growth process of largesized γ′ phases, both of which were driven by the bulk diffusion of solute elements.Consequently, the interval between the extrema curves is associated with the time taken for solute atom concentrations to reach equilibrium.According to Fick's Second Law, ( , ) =  ( , ) , it is evident that the concentration distribution of solute elements in the matrix was influenced by both time and co-ordinates.
Assuming that the gaps between particles in unit area are arranged in a honeycomb pattern, as depicted in Figure 18a, and the distance between the centers of two particles is ℎ, the area of a regular hexagon can be calculated by  = √ ℎ .Further, it can be obtained based on the number of particles per unit area N that: The curve formed by the first and third minimum points of the ∆ρ/ρ 0 − t curve in Figure 17a is defined as the "maximum total content of γ ′ phase curve", corresponding to the first and third curves from the left in Figure 17b.Likewise, the curve formed by the second and fourth maximum points of the ∆ρ/ρ 0 − t curve in Figure 17a is defined as the "minimum total content of γ ′ phase curve", corresponding to the second and fourth curves from the left in Figure 17b.These curves can be used to determine whether the total content of the γ ′ phase increases or decreases during a specific aging period at the given temperature.For example, at 825 • C the total content of the γ ′ phase decreased during aging for 0.5-1 h, increased during aging for 1-2 h, and then decreased again during aging for 2-7 h.
As observed in Figure 17b, there is a general trend of increasing interval between the extrema curves of the total content of the γ ′ phase.In line with the prior analysis, the fluctuation in the total volume of the γ ′ phase is attributed to the alternating dominance of the redissolution process of small-sized γ ′ phases and the growth process of large-sized γ ′ phases, both of which were driven by the bulk diffusion of solute elements.Consequently, the interval between the extrema curves is associated with the time taken for solute atom concentrations to reach equilibrium.According to Fick's Second Law, ∂C(x,t) ∂t = D ∂ 2 C(x,t) ∂x 2 , it is evident that the concentration distribution of solute elements in the matrix was influenced by both time and co-ordinates.
Assuming that the gaps between particles in unit area are arranged in a honeycomb pattern, as depicted in Figure 18a, and the distance between the centers of two particles is h, the area of a regular hexagon can be calculated by S RH = √ 3 2 h 2 .Further, it can be obtained based on the number of particles per unit area N that: Hence, the average distance d between two γ ′ particles at a certain moment is represented as: where r is the average diameter of the γ ′ particles shown in Figure 5a.
Hence, the average distance  between two γ′ particles at a certain moment is rep-resented as: where ̅ is the average diameter of the γ′ particles shown in Figure 5a.As shown in Figure 18b, the increase in the average distance between γ′ particles indicates a prolonged time required for solute atoms to reach equilibrium through diffusion, which corresponds to increased intervals between the extrema curves of the total content of the γ′ phase.Previous research has shown that the redissolution of the γ′ phase can lead to the precipitation of new γ′ in oversaturated regions of elements [65].The extensive occurrence of such a process during aging for 4 h and 10 h at 850 °C, and for 2 h, 4 h, and 10 h at 800 °C brought about a decrease in particle spacing.

Conclusions
(1) The resistivity exhibits fluctuation over aging time and is negatively correlated with the area fraction of the γ′ phase.This indicates that the solute concentration in the matrix is the primary factor for the resistivity change.The enrichment of solute atoms led to lattice distortion, facilitating the formation of short-range-ordered structures, which increase electron scattering and thus elevate the resistivity.Additionally, the precipitation and growth of M C at grain boundaries also contributed to the gradual increase in resistivity over aging time.(2) The γ′ phase exhibited pronounced ripening during the aging process.The LSW ripening theory successfully predicted the size evolution of the γ′ phase, and its coarsening activation energy was determined to be 277.59kJ/mol, indicating that the coarsening of γ was primarily due to bulk diffusion of Al and Ti elements.Differences in the chemical potential of solute atoms in the matrix near the surfaces of differently sized γ′ particles led to interactions between adjacent γ′ particles.This resulted in the redissolution of small-sized γ′ particles and the growth of large-sized γ′ particles.Since the growth rate was not consistent with the dissolution rate, fluctuations appeared in the area fraction of the γ′ phase.As shown in Figure 18b, the increase in the average distance between γ ′ particles indicates a prolonged time required for solute atoms to reach equilibrium through diffusion, which corresponds to increased intervals between the extrema curves of the total content of the γ ′ phase.Previous research has shown that the redissolution of the γ ′ phase can lead to the precipitation of new γ ′ in oversaturated regions of elements [65].The extensive occurrence of such a process during aging for 4 h and 10 h at 850 • C, and for 2 h, 4 h, and 10 h at 800 • C brought about a decrease in particle spacing.

Conclusions
(1) The resistivity exhibits fluctuation over aging time and is negatively correlated with the area fraction of the γ ′ phase.This indicates that the solute concentration in the matrix is the primary factor for the resistivity change.The enrichment of solute atoms led to lattice distortion, facilitating the formation of short-range-ordered structures, which increase electron scattering and thus elevate the resistivity.Additionally, the precipitation and growth of M 23 C 6 at grain boundaries also contributed to the gradual increase in resistivity over aging time.
(2) The γ ′ phase exhibited pronounced ripening during the aging process.The LSW ripening theory successfully predicted the size evolution of the γ ′ phase, and its coarsening activation energy was determined to be 277.59kJ/mol, indicating that the coarsening of γ was primarily due to bulk diffusion of Al and Ti elements.Differences in the chemical potential of solute atoms in the matrix near the surfaces of differently sized γ ′ particles led to interactions between adjacent γ ′ particles.This resulted in the redissolution of small-sized γ ′ particles and the growth of large-sized γ ′ particles.Since the growth rate was not consistent with the dissolution rate, fluctuations appeared in the area fraction of the γ ′ phase.(3) By connecting the inflection points of the resistivity change curves at different temperatures, the extrema curves of the total content of the γ ′ phase were obtained.These curves can be used to determine the total content changes of the γ ′ phase during aging at a specific temperature and time.This provides a new method for characterizing the fluctuation of the total content of the γ ′ phase.

Figure 1 .
Figure 1.Temperature-time curve for solution and aging heat treatment.

Figure 5 .
Figure 5. (a) Evolution of the average diameter of the γ′ phase with aging time at 800 °C and 850 °C; (b) evolution of particle density with aging time at 800 °C and 850 °C.

Figure 6
Figure 6 depicts the relationship between the area fraction of the γ′ phase and aging time at 850 °C and 800 °C.The area fraction of the γ′ phase fluctuates with aging time, indicating that the precipitated phase does not reach equilibrium during short-term aging.Considering the increase in the particle size and the overall decrease in the particle number over time, it can be inferred that that both precipitation and Ostwald ripening phenomena occurred during the short-term aging process.In other words, smaller particles dissolved, releasing solute atoms that then re-precipitate onto larger particles, which leads to an increase in the average particle size.

Figure 5 .
Figure 5. (a) Evolution of the average diameter of the γ′ phase with aging time at 800 °C and 850 °C; (b) evolution of particle density with aging time at 800 °C and 850 °C.

Figure 6
Figure 6 depicts the relationship between the area fraction of the γ′ phase and aging time at 850 °C and 800 °C.The area fraction of the γ′ phase fluctuates with aging time, indicating that the precipitated phase does not reach equilibrium during short-term aging.Considering the increase in the particle size and the overall decrease in the particle number over time, it can be inferred that that both precipitation and Ostwald ripening phenomena occurred during the short-term aging process.In other words, smaller particles dissolved, releasing solute atoms that then re-precipitate onto larger particles, which leads to an increase in the average particle size.

Figure 5 .
Figure 5. (a) Evolution of the average diameter of the γ ′ phase with aging time at 800 • C and 850 • C; (b) evolution of particle density with aging time at 800 • C and 850 • C.

Figure 6
Figure 6 depicts the relationship between the area fraction of the γ ′ phase and aging time at 850 • C and 800 • C. The area fraction of the γ ′ phase fluctuates with aging time, indicating that the precipitated phase does not reach equilibrium during short-term aging.Considering the increase in the particle size and the overall decrease in the particle number over time, it can be inferred that that both precipitation and Ostwald ripening phenomena occurred during the short-term aging process.In other words, smaller particles dissolved, releasing solute atoms that then re-precipitate onto larger particles, which leads to an increase in the average particle size.

Figure 6 .
Figure 6.The variation of the γ′ phase area fraction with time during short-term aging at (a) 850 °C and (b) 800 °C.

Figure 7 .
Figure 7. Evolution of M C at the grain boundaries during aging at 850 °C for different durations: (a) 4 h; (b) 10 h.

Figure 6 .
Figure 6.The variation of the γ ′ phase area fraction with time during short-term aging at (a) 850 • C and (b) 800 • C.

21 Figure 6 .
Figure 6.The variation of the γ′ phase area fraction with time during short-term aging at (a) 850 °C and (b) 800 °C.

Figure 7 .
Figure 7. Evolution of M C at the grain boundaries during aging at 850 °C for different durations: (a) 4 h; (b) 10 h.

Figure 8 .
Figure 8. (a,b) The γ ′ particles around the MC carbide when carbide decomposition occurs at 850 • C for 10 h of aging; (c) a magnified image of finer γ ′ particles.

Figure 9 .
Figure 9.The relationship between the relative change rate of resistivity and the area fraction of the precipitated γ′ phase during aging at (a) 850 °C, (b) 800 °C.

Figure 9 .
Figure 9.The relationship between the relative change rate of resistivity and the area fraction of the precipitated γ ′ phase during aging at (a) 850 • C, (b) 800 • C.

Figure 10 .
Figure 10.Phase and KAM value distribution maps after aging at 850 °C.(a) Phase map and (b) KAM value distribution map of the same location after 1 h of aging.(c) Phase map and (d) KAM value distribution map near the grain boundary after 10 h of aging.

Figure 10 .
Figure 10.Phase and KAM value distribution maps after aging at 850 • C. (a) Phase map and (b) KAM value distribution map of the same location after 1 h of aging.(c) Phase map and (d) KAM value distribution map near the grain boundary after 10 h of aging.

Figure 12 .
Figure 12.Fitting results of the relationship between γ ′ phase particle size and aging time according to the LSW ripening theory at 850 • C and 800 • C.

Figure 13 .
Figure 13.(a) The solid solution matrix with a solute concentration of  ; (b) the precipitation and growth stages of the γ′ phase.

Figure 13 .
Figure 13.(a) The solid solution matrix with a solute concentration of C 0 ; (b) the precipitation and growth stages of the γ ′ phase.

Figure 14 .
Figure 14.Solute concentration distribution at (a) the stage when the redissolution process of smaller-sized γ′ particles was dominant; (b) the moment when the smaller-sized γ′ phase disappeared.

Figure 14 .
Figure 14.Solute concentration distribution at (a) the stage when the redissolution process of smallersized γ ′ particles was dominant; (b) the moment when the smaller-sized γ ′ phase disappeared.

Materials 2024 , 21 Figure 15 .
Figure 15.The new round of fluctuation in the total amount of γ′ phase: (a) the growth process of γ′ phase dominates in the new loop; (b) the redissolution process of γ′ phase dominates in the new loop.

Figure 15 .
Figure 15.The new round of fluctuation in the total amount of γ ′ phase: (a) the growth process of γ ′ phase dominates in the new loop; (b) the redissolution process of γ ′ phase dominates in the new loop.

Figure 15 .
Figure 15.The new round of fluctuation in the total amount of γ′ phase: (a) the growth process of γ′ phase dominates in the new loop; (b) the redissolution process of γ′ phase dominates in the new loop.

Figure 16 .
Figure 16.(a) XRD diffraction patterns corresponding to different aging times at 850 °C; (b) changes in the average lattice constant calculated from 2θ.

Figure 16 .
Figure 16.(a) XRD diffraction patterns corresponding to different aging times at 850 • C; (b) changes in the average lattice constant calculated from 2θ.

Figure 18 .
Figure 18.(a) Schematic diagram of close-packed particles; (b) variation of surface spacing between two adjacent γ′ phases with aging time at 800 °C and 850 °C.

Figure 18 .
Figure 18.(a) Schematic diagram of close-packed particles; (b) variation of surface spacing between two adjacent γ ′ phases with aging time at 800 • C and 850 • C.
Figure 1.Temperature-time curve for solution and aging heat treatment.