3.1. CALPHAD Calculations
The plot of the thermodynamically stable phases in terms of oxygen activity (
aO) for all investigated steel grades at 650 °C is shown in
Figure 1,
Figure 2,
Figure 3 and
Figure 4. The increase in oxygen activity is related to the penetration of oxygen into the steel matrix and the subsequent oxidation process. The highest oxygen activity corresponds to the outermost surface of the oxide scale, while the inner oxide layers are characterised by lower oxygen activities.
Figure 1 illustrates the results for the steel 16Mo3 at 650 °C. Within the range of oxygen activity (
aO) up to about 10
−13, the stable oxide (Fe, Cr)
2O
3 (95.9 wt.%) predominates, followed by (Fe, Mn)
2SiO
4 (2.8 wt.%), and the remaining 1.3 wt.% is (Fe, Cr, Mn, Mo)
3O
4. In the range
aO ≈ 10
−13–10
−23, (Fe, Cr, Mn, Mo)
3O
4 is stable and predominant (96.7 wt.%), followed by (Fe, Mn)
2SiO
4 (3.3 wt.%), and in the range
aO ≈ 10
−22–10
−23, the (Fe, Mn, Ni)O content increases to 90.7 wt.%, with the remaining components being (Fe, Cr, Mn, Mo)
3O
4 and (Fe, Mn)
2SiO
4 (although their content varies with
aO due to changes in wüstite content). Also noteworthy is the occurrence of internal oxidation at 650 °C, with (Fe, Cr, Mn, Mo)
3O
4 and (Fe, Mn)
2SiO
4 remaining stable below the steel surface. Internal oxidation is assessed by examining the content of the α-ferrite phase, which is 97.5 wt.% at equilibrium at a temperature of 650 °C. Therefore, we can determine whether internal oxidation takes place or not. If the α-ferrite content increases (black line), it should immediately increase to 97.5 wt.%; if not, internal oxidation takes place, as shown in
Figure 1. Thus, we can determine the occurrence of internal oxidation. When the α-ferrite content increases (indicated by the black line), it should immediately reach 97.5 wt.%. If this threshold is not reached, it indicates that internal oxidation is taking place, as shown in
Figure 1.
Figure 2 shows the results for 13Cr steel at 650 °C. In the range of oxygen activity (
aO) up to about 10
−13, the stable oxide (Fe, Cr, Mn)
2O
3 (95.6 wt.%) predominates, followed by Fe
2SiO
4 (4.4 wt.%). In the
aO range of about 10
−13–10
−22, (Fe, Cr, Mn, Mo)
3O
4 becomes stable and predominant (96.7 wt.%), followed by Fe
2SiO
4 (3.3 wt.%). In the
aO range of about 10
−22–10
−23, the (Fe, Mn, Ni)O content increases up to 83 wt.% while the remainder consists of (Fe, Cr, Mn, Mo)
3O
4 and Fe
2SiO
4. It is important to note that the content of these two components changes with the
aO content in line with the variations in wüstite content. Furthermore, the results indicate the presence of internal oxidation at 650 °C, where (Fe, Cr, Mn, Mo)
3O
4 and Fe
2SiO
4 remain stable under the steel surface. The determination of internal oxidation is based on the content of the α-ferrite phase, which is calculated to be 97 wt.% in equilibrium at 650 °C. Any deviation from this value indicates the occurrence of internal oxidation. If the α-ferrite content increases (black line), it should immediately increase to 97 wt.%. Thus, we can determine the occurrence of internal oxidation. When the α-ferrite content increases (indicated by the black line), it should immediately reach 97 wt.%. If this threshold is not reached, it indicates that internal oxidation is taking place, as shown in
Figure 2.
Figure 3 shows the results for T24 steel at 650 °C. Within the range of oxygen activity (
aO) up to approximately 10
−13, the stable oxide (Fe, Cr, V)
2O
3 (94.3 wt.%) predominates, followed by (Fe, Mn)
2SiO
4 (5.7 wt.%). In the
aO range of about 10
−13–10
−22, (Fe, Mo, Mn, Cr)
3O
4 becomes stable and predominant (96.7 wt.%), followed by (Fe, Mn)
2SiO
4 (3.3 wt.%). In the
aO range of about 10
−22–10
−23, (Fe, Cr, V)O content increases up to 76.4 wt.% while the remainder consists of (Fe, Mo, Mn, Cr)
3O
4 and (Fe, Mn)
2SiO
4. It is important to note that the content of these two components changes with the
aO content in line with the variations in wüstite content. Furthermore, the results indicate the presence of internal oxidation at 650 °C, where (Fe, Mo, Mn, Cr)
3O
4, (Fe, Mn)
2SiO
4 and (Fe, Cr, V)
2O
3 remain stable under the steel surface. The determination of internal oxidation is based on the content of the α-ferrite phase, which is calculated to be 98.5 wt.% in equilibrium at 650 °C. Any deviation from this value indicates the occurrence of internal oxidation. When the α-ferrite content increases (indicated by the black line), it should immediately reach 98.5 wt.%. If this threshold is not reached, it indicates that internal oxidation is taking place, as shown in
Figure 3.
Figure 4 shows the results for the steel P91 at 650 °C. In the range of oxygen activity (
aO) up to about 10
−14, the stable oxide (Cr, V)
2O
3 (94.2 wt.%) predominates, followed by (Fe, Cr)
2SiO
4 (5.8 wt.%). In the
aO range of about 10
−14–10
−22, (Cr, Mn, V, Fe)
3O
4 becomes stable and predominant (96.1 wt.%), followed by (Fe, Cr)
2SiO
4 (3.9 wt.%). In the
aO range of about 10
−23–10
−24, (Fe, Cr, V, Mn)O content increases up to 45.8 wt.% while the remainder consists of (Cr, Mn, V, Fe)
3O
4 and (Fe, Cr)
2SiO
4. It is important to note that the content of these two components changes with the
aO content in line with the variations in wüstite content. Furthermore, the results indicate the presence of internal oxidation at 650 °C, where (Cr, Mn, V, Fe)
3O
4, (Fe, Cr)
2SiO
4 and (Cr, V)
2O
3 remain stable under the steel surface. The determination of internal oxidation is based on the content of the α-ferrite phase, which is calculated to be 97.8 wt.% in equilibrium at 650 °C. Any deviation from this value indicates the occurrence of internal oxidation. When the α-ferrite content increases (indicated by the black line), it should immediately reach 97.8 wt.%. If this threshold is not reached, it indicates that internal oxidation is taking place, as shown in
Figure 4.
Based on the CALPHAD results, we can gain some initial insight into the equilibrium composition of the oxide layers that would form on the steels studied at an oxidation temperature of 650 °C. Of course, these data refer to the equilibrium state, so they cannot be directly related to the composition of the oxide layer formed after 48 h, as equilibrium had not yet been reached. We also get information about the alloying elements in the individual oxide sublayers, which is useful for further metallographic analyses. Another very useful piece of information is the wüstite content in each equilibrium oxide layer formed on the steels studied. The results show that 16Mo3 steel has up to 90.7 wt.% wüstite in the equilibrium oxide layer formed, 13Cr steel has up to 83 wt.% wüstite in the oxide layer formed, followed by T24 steel, which has up to 76.4 wt.% wüstite, and P91 which has up to 76.4 wt.% wüstite in the oxide layer formed. Based on the wüstite content, we can assume that the oxidation kinetics are fastest for 16Mo3 (the steel with the highest wüstite content), followed by 13Cr, T24 and P91 (the steel with the lowest wüstite content), and that the oxidation kinetics are slowest for the P91 steel. This is due to the well-known fact that diffusion through wüstite is the fastest [
20,
21].
3.2. Thermogravimetric Analysis
Thermogravimetric analysis (TGA) served as the primary tool for monitoring the weight variations in steel samples during exposure to elevated temperatures. TGA is the preferred technique for evaluating oxidation rates as it allows for continuous or intermittent data collection [
22]. To effectively model the TGA results, we used an iterative rectangular distance regression algorithm. Three mathematical functions were used to explain the TGA curves: a two-phase exponential growth function (Equation (1)), a second-degree polynomial (Equation (2)—the so-called parabolic law) and a third-degree polynomial (Equation (3)—the so-called cubic law). These functions were used to accurately describe the observed TGA results.
In all cases, t is the time in s, Δm is the change in weight in mg, A is the specific surface area of the sample in cm2, while the other coefficients depend on the temperature and the chemical composition of the steel.
In each case, t stands for the time in s, Δm for the change in weight mg and A for the specific surface area of the sample in cm2, while the other coefficients depend on both the temperature and the chemical composition of the steel.
However, since the analysis of oxidation kinetics typically revolves around the change in weight over time, a parabolic law can be succinctly formulated by using the Pilling–Bedworth equation [
20,
23]:
Here,
W stands for the change in weight per unit area due to the oxidation of iron.
kp (=
) stands for the parabolic constant in g
2 cm
−4 s
−1, while
W0 corresponds to the initial weight at the beginning of the parabolic oxidation of parabolic oxidation (
t = 0). It is noteworthy that in the original Pilling–Bedworth equation [
20,
23],
W0 is set to zero. The most appropriate approach for deriving the equation to calculate the rate constant for steel oxidation under continuous heating or cooling was presented by Kofstad [
24]. According to Kofstad’s interpretation, the oxidation behaviour, regardless of whether it follows a linear, parabolic or cubic law, can be expressed as follows:
In this context,
W is defined as the change in weight per unit area with respect to time
t. The variable
n takes on constant values of 1, 2 or 3, corresponding to the linear, parabolic and cubic laws, respectively. In addition,
kn represents a time-independent rate constant and is expressed as follows:
Here, B stands for the constant, T means the absolute temperature, R stands for the gas constant and Q denotes the activation energy.
The rate constants for all the samples studied were calculated after deriving Equation (5) by a simple linear regression approach. The basic equation for calculating the rate constant (
kn) is expressed as follows:
In this context, Δm stands for the weight change in mg, A for the specific surface area in cm2, t for the oxidation time in s and kn for rate constant. The value n can be 1, 2 or 3, which corresponds to a linear, parabolic or cubic law. For the exponential law, the value of n is in the range 1 < n < 3. The designation of the index n is based on the type of the equation, which leads to the exponential law being associated with the index e and the cubic law with the index c.
The following graph (
Figure 5) shows the TGA results of the steels examined. The results show an increase in mass over time. Since a constant airflow was introduced, the increase in mass is due to the chemical bonding of oxygen, i.e., the formation of an oxide layer/oxidation of the samples. After 4 h, the oxidation kinetics of steel 16Mo3 starts to follow a parabolic law and the TGA results can be expressed with a parabolic equation, but overall, the best fit is still the exponential equation. For steel 13Cr, the same trend is observed, but it seems to be faster as the oxidation kinetics starts to follow a parabolic law after 2 h. On the other hand, the oxidation kinetics of steel T24 follows the cubic law. The oxidation kinetics of steel P91 can be described by the cubic law, as the TGA results agree best with the cubic equation. In a study on the high-temperature oxidation behaviour of P91 steel [
18], it was suggested that it follows a parabolic law (at 600 °C and 700 °C). We cannot confirm this, but on the other hand, its oxidation time was 1000 h. Furthermore, the authors suggest a more severe oxidation, with a high increase in the oxide layer. Our experimental data suggest otherwise, as we did not find any evidence of substantial P91 oxidation at 650 °C. We assume that this is the same trend as for the 16Mo3 and 13Cr steels in this study.
The equation types with the corresponding coefficients that best fit the TGA results are listed in
Table 2.
The calculated results of the rate constants are shown in
Table 3. The results of the rate constants are a clear indication that the oxidation kinetics is fastest for 13Cr steel, followed by 16Mo3, T24 and P91. This means that 13Cr steel oxidises the most and P91 the least.
Based on the TGA results, we can confirm that the CALPHAD assumption about the correlation between the amount of wüstite in the oxide layer formed and the oxidation kinetics is true. This is because 16Mo3 has the fastest oxidation kinetics, followed by 13Cr, T24 and P91. The same trend was observed in the CALPHAD results with respect to the wüstite content, with 13Cr having the highest amount of wüstite in the oxide layer, followed by 13Cr, T24 and P91.
3.3. Microscopy
The oxide layers formed were further analysed using SEM (
Figure 6). It can be seen that the thickest oxide layer was formed on steel 16Mo3, followed by 13 Cr, T24 and P91. The results are in agreement with the analysis of TG (
Figure 5) and consequently with the calculated rate constants (
Table 3).
In addition, the thickness of the oxide layers formed on the steels examined was measured. The following table (
Table 4) shows the average of the measured thickness. The results show that the thickest oxide layer was formed on 13Cr steel, followed by 16Mo3, T24 and P91 (with no oxide layer observed on SEM). For 13Cr steel, there are deviations from the TG analysis results, but if you look at
Figure 6, it is obvious that the oxide layer on 13Cr steel is cracked and chipped, which is why these results differ. In addition, as already mentioned, no oxide layer was found on the surface of P91 steel after high-temperature oxidation, which is also in agreement with the results of the TG analysis (
Figure 5).
The oxide layers formed were also examined using EBSD analysis. The EBSD analysis also gives an insight into the crystal grain size, and the samples with a more substantial oxide layer (16Mo3, 13Cr and T24) show that the inner layer is usually composed of very small grains, while the outer layers are composed of larger crystal grains. The outer layer has grown for a longer time and the grains have undergone coalescence. The smaller grains in the inner layer present a difficulty for proper EBSD indexation.
Figure 7 shows the oxide layer formed after high-temperature oxidation of 16Mo3 steel. It is obvious that the outer and middle oxide layers consist mainly of hematite. On the other hand, the inner oxide sublayer is mainly magnetite. IPF-Z mapping also showed that the grain size in the outer and inner oxide layers is smaller than in the middle oxide sublayer. However, EDS analysis showed no differences within the oxide layer. Fe, O and Cr are essentially homogeneously distributed throughout the oxide layer formed. In the oxide layer (mainly in the middle and inner oxide sublayers) there are some yellow-coloured grains, which are wüstite. Since the wüstite crystal grains are barely visible, we have enlarged the image of the EBSD phase analysis (
Figure 8), on which the wüstite crystal grains are clearly visible. In the meantime, the matrix/steel has a ferritic (bcc—body-centred cubic) crystal structure, which is true for all samples examined (
Figure 7,
Figure 8,
Figure 9,
Figure 10,
Figure 11,
Figure 12 and
Figure 13).
Figure 9 shows the oxide layer formed after high-temperature oxidation of 13Cr steel. It is obvious that the outer and inner oxide sublayer consists mainly of hematite. On the other hand, the inner oxide sublayer is mainly magnetite. IPF-Z mapping also shows that the grain size in the outer and inner oxide sublayer is smaller than in the middle oxide sublayer. The EDS analysis of the surface distribution of the elements shows that there is an increased Cr content in the inner oxide sublayer. However, as far as O and Fe are concerned, they are essentially homogeneously distributed throughout the oxide layer formed. In the oxide layer (mainly in the middle and inner oxide sublayers), there are some yellow-coloured grains which are not ferrite but wüstite, as both have the same body-centred cubic crystal structure (bcc). In this case, there are fewer wüstite crystal grains in the oxide layer than in the 16Mo3 steel. Since the wüstite crystal grains are hardly visible, we have enlarged the image of the EBSD phase analysis (
Figure 10), on which the wüstite crystal grains are clearly visible. This is also consistent with the results of the TG analysis (
Figure 5), as diffusion is fastest in the wüstite, which means that the oxidation rate should be higher for 16Mo3 steel than for 13Cr and the calculated rate constant (
Table 3) should also be lower than that calculated for 16Mo3 steel.
The last oxide layer analysed was one that had formed on T24 steel after high-temperature oxidation (
Figure 11). In this case, although an oxidation layer was found, the entire sublayer consisted mainly of hematite. Some magnetite crystal grains were found in the inner oxide sublayer and some wüstite grains were also found in the middle oxide sublayer. There are only a few wüstite grains in the oxide layer, which are hardly visible, compared to the steels 16Mo3 and 13Cr, where the amount of wüstite crystal grains was higher. Since the wüstite crystal grains are hardly visible, we have enlarged the image of the EBSD phase analysis (
Figure 12), on which the few wüstite crystal grains are clearly visible. As far as the size of the grains is concerned, the same trend as in the other analysed oxide layers is shown, i.e., the outer and inner oxide sublayers consist of smaller crystal grains than the inner oxide sublayer. EDS analysis of the surface distribution of the elements shows that the Cr content is increased in the inner oxide sublayer, while the Fe content decreases in the inner oxide sublayer. O, on the other hand, is essentially homogeneously distributed over the entire oxide layer formed. The results again agree well with the analytical results from TG (
Figure 5), where T24 steel shows the second lowest weight increase during high-temperature oxidation and the second lowest calculated rate constant (
Table 3).
Figure 13 shows the surface of P91 steel after high-temperature oxidation. In this case, there was no oxide layer that could be analysed with EBSD. The result was to be expected as the TG analysis also showed that there was only a minimal increase in weight during oxidation. EDS analysis was also only carried out to show that there is no oxide layer that can be analysed (there is no increased amount of oxygen on the steel surface); it only shows that there are some Cr-based carbides in the matrix (areas with increased Cr). This is also consistent with the results of the TG analysis (
Figure 5), as P91 steel has the lowest weight increase during high-temperature oxidation and the lowest calculated rate constant (
Table 3).
The Cr content is the decisive factor in high-temperature oxidation. Even an average content of 2.6 wt.% lowers the oxidation rate considerably (as can be seen in
Table 3); the Cr content also changes the oxidation curve from exponential to cubic. This means that steels with low Cr content such as 16Mo3 and 13Cr are susceptible to severe oxidation, even in the presence of short-term, transient overheating. The strong weight gain due to oxidation is already visible after one hour (see
Figure 5). The thicker oxide layers in the 16Mo3 and 13Cr steels show a clear distinction between hematite and magnetite, as shown by the EBSD analysis in
Figure 7,
Figure 8,
Figure 9 and
Figure 10. The formation of FeO increases the oxidation kinetics. Therefore, the less stable samples such as 16Mo3 should have the highest amount of FeO, followed by 13 Cr and T24. This is consistent with CALPHAD calculations and EBSD analysis, as well as TGA results. Furthermore, such oxidation changes the surface of the steel and can affect structural integrity. At higher Cr contents, such as in steel T24, magnetite should form, but the crystal structure is difficult to assess with EBSD due to the smaller grain size, so only the Fe-rich hematite can be clearly distinguished. The surface is certainly affected by the short-term overheating, but it is not as severe and would not compromise the integrity of the components as severely as in the case of 16Mo3 and 13Cr. The high Cr contents in P91 prevent a much thicker oxide layer formation so that the oxide layer cannot be observed (
Figure 13). This is also consistent with the TG curves (
Figure 5) and the calculated constant rates (
Table 3). The oxide layers of 16Mo3 and 13Cr are brittle and crumble during sample preparation, which leads to unreliable measurements of the oxide layer thickness, so we focused on the thermogravimetric measurements.