Low Alloy Steel Tempering Kinetics

When a quenched low alloy steel is tempered, two reactions run through it in sequence and partly overlap. A dilatometer records both of them as changes in the length of the specimen, and the curve can be separated into the two reactions and fitted to the Avrami formalism so that each one gets its own exponent, rate constant and activation energy. This page collects those parameters for two steels, a 0.32% carbon low alloy steel tempered isothermally at 350 °C, and AISI 1050 tempered at constant heating rates of 5 to 30 K per minute. For a tool steel, where the tempering response of a specific grade is what matters, the practical curves are on the tool steel tempering chart; this page is the reaction kinetics behind those curves.

The two processes a dilatometry curve shows

The steel used for the isothermal work contains 0.32% C, 1.12% Mn, 0.67% Si, 0.07% Ni, 0.02% P and 0.05% S, the balance iron. Specimens 9 mm in diameter and 100 mm long were austenitized at 900 °C for 30 minutes, quenched in water at ambient temperature, and then held at 350 °C in a dilatometer with a length resolution of 10-3 mm. Its martensite start temperature is about 365 °C, so the 350 °C hold sits just below Ms and the whole tempering sequence happens in one isothermal step.

Isothermal dilatometry curve for a low alloy steel tempered at 350 degrees C showing the two tempering processes
Fig. 5. Isothermal dilatometry record for the 0.32% carbon low alloy steel held at 350 °C. The curve starts at the beginning of the hold, t0, and shows the first process (I) ending at t end 1, after which the second process (II) continues. Reproduced from the source figure.

Two contractions appear in the length of the specimen. The first is the faster of the two and ends at t end 1, the second takes over after it and runs on to the end of the record. Each was fitted separately, using only a work interval well away from the start and the end of the reaction, and the fitted parameters are collected below.

Parameter, first processValue
Best interval in the experimental fraction transformed0.5050 to 0.7625
Work interval in the first process fraction after the final iteration0.5310 to 0.7750
N, slope of the linear fit1.58
Correlation coefficient R for N0.993
n1, Avrami exponent1.49
Correlation coefficient R for n10.992
k1, rate constant (s-1)1.2 x 10-4

Source, Table 10.b of the source chapter. Kinetic parameters for the first process on tempering of the 0.32% carbon low alloy steel, fitted by the method the source calls KTOPT.

Parameter, second processValue
Best interval in the experimental fraction transformed2.864 to 3.963
Normalized interval in the transformed fraction0.4506 to 0.6544
Work interval in the second process fraction after the final iteration0.5975 to 0.8010
N, slope of the linear fit-0.669
Correlation coefficient R for N0.998
n2, Avrami exponent0.669
Correlation coefficient R for n20.998
k2, rate constant (s-1)4.4 x 10-5

Source, Table 10.a of the source chapter. Kinetic parameters for the second process on the same steel and the same 350 °C hold.

The rate constants differ by a factor of about three, k1 at 1.2 x 10-4 per second and k2 at 4.4 x 10-5 per second, which is the quantitative statement of the first process being the faster one.

Which reaction each process is

The Avrami exponents identify the mechanisms, because the value of n is fixed by how a reaction nucleates and grows. An exponent of 1.5 corresponds to diffusion controlled growth from a fixed number of nuclei, that is, growth with a nucleation rate of zero, and that is the value returned for the first process. The first stage of tempering is the decomposition of martensite into a transition carbide and a less tetragonally distorted martensite, and the exponent is consistent with that reaction when it happens in a specimen that was already partly transformed during the quench. The study notes that the transition carbide may in part have formed before the 350 °C hold began, since Ms is about 365 °C.

An exponent of about 0.67 corresponds to precipitation on dislocations with diffusion controlled growth, for which the literature value is 0.66. The second process is the disappearance of the transition carbide and the formation of stable cementite, and cementite is reported to nucleate on dislocations and on lath and grain boundaries, which is exactly the mechanism that exponent describes. Applying the uncertainty in the measured fraction transformed to the exponents gives n1 of 1.5 with bounds of about 0.2, and n2 of 0.67 with bounds of about 0.02.

Non-isothermal tempering of AISI 1050

The second half of the work takes a different steel and a different thermal history. The steel is AISI 1050, with 0.48 to 0.55% C, 0.6 to 0.9% Mn, 0.04% P and 0.05% S, and it was heated from ambient temperature to 600 °C at constant rates of 5, 10, 15, 20 and 30 K per minute in a dilatometer with a resolution of 10-4. Plotting the rate of change of the relative length against temperature separates the two reactions into two bands, and the higher the heating rate, the further both bands move to the right.

Derivative of the relative length change against temperature at heating rates of 5 to 30 K per minute
Fig. 8. Rate of change of the relative length against temperature for AISI 1050 at the five heating rates. The two bands are the two tempering processes. Reproduced from the source figure.

The temperature of the inflection point of each band, at each heating rate, is the raw input to the rest of the calculation. An inflection point is used rather than a peak because the analysis requires the point at which the second derivative of the fit is zero.

β (K/min)Ti (K), Process ITi (K), Process II
5411.6584.3
10418.9593.5
15424.3599.1
20428.5603.1
30432.4608.8

Source, Table 11 of the source chapter. Temperatures corresponding to the inflection points at different heating rates, for the two processes on tempering of AISI 1050.

The shift is modest in absolute terms but consistent in direction. Raising the heating rate from 5 to 30 K per minute moves the first process from 411.6 to 432.4 K, a shift of about 21 K, and the second from 584.3 to 608.8 K, a shift of about 25 K. That displacement is the basis of every Kissinger-like method, because a reaction with a higher activation energy shifts further for the same change in heating rate.

Activation energy and frequency factor

Two routes to the activation energy are reported. A Kissinger-like plot uses only the temperatures at which a fixed fraction is transformed, and it returns the effective activation energy without any assumption about the reaction model. A second, non-linear regression route fits both the activation energy and the frequency factor at the same time. The results are set out below, and the two routes agree.

ProcessActivation energy E (kJ/mol)Frequency factor K0 (min-1)Confidence interval for ERegression coefficient RP value for E
I1172.99 x 1014104 to 130 (linear)0.9984.8 x 10-10
II2069.53 x 1017196 to 216 (linear)0.9999.0 x 10-14

Source, Table 12 of the source chapter. Best parameters E and K0 calculated by linear and non-linear regression for the two processes on tempering of AISI 1050. The source also prints narrow non-linear confidence intervals, 116.5 to 117.4 kJ/mol for the first process and 205.6 to 206.4 kJ/mol for the second, and gives the P value as the probability that another model fits the same data better.

The first process needs 117 kJ/mol and the second 206 kJ/mol. The second value sits in the range normally quoted for the diffusion of carbon in ferrite and for the coarsening of cementite, which is consistent with the mechanism assigned to it above. The frequency factors span three orders of magnitude between the two processes, and the two values are printed here as the source gives them.

The Avrami exponent at each heating rate

The same non-isothermal records also yield the Avrami exponent, but the authors take it from a transcendental equation rather than from the regression used above, because they consider the regression route too sensitive to the approximations behind it. The exponents and their uncertainties are given per heating rate.

First process on tempering
β (K/min)Ti (K)nδn
5411.61.20.2
10418.91.10.2
15424.31.00.2
20428.51.10.2
30432.41.00.2
Second process on tempering
β (K/min)Ti (K)nδn
5584.30.70.1
10593.50.70.1
15599.10.70.1
20603.10.70.1
30608.80.70.1

Source, Table 13 of the source chapter. The Avrami exponent for the two processes on tempering of AISI 1050 at each heating rate, with the uncertainty calculated from the uncertainty in temperature and in activation energy. The two frames are the source’s own division of one table, first process above and second process below.

The first process scatters between 1.0 and 1.2 with a stated uncertainty of 0.2 on every row, so the values are not distinguishable from one another, and the mean sits close to the 1.5 obtained isothermally once the uncertainty is taken into account. The second process is constant at 0.7 to within 0.1 at every heating rate, which is the same exponent the isothermal fit returned.

Activation energy and frequency factor at each heating rate

Fixing the exponent permits the frequency factor to be calculated for each heating rate separately, and the resulting activation energies can be compared with the single value from the Kissinger-like plot above.

β (K/min)Ti (K)n0K0 (min-1)E (kJ/mol)δE (kJ/mol)
5411.61.05.4 x 10141196
10418.91.06.87 x 10131126
15424.31.03.3 x 10141176
20428.51.02.0 x 10151247
30432.41.01.1 x 10141136

Source, Table 14 of the source chapter. Activation energies and frequency factors for the first process on tempering of AISI 1050, calculated for each heating rate from the inflection point temperature with the exponent n0 fixed at 1.0.

β (K/min)Ti (K)n0K0 x 1017 (min-1)E (kJ/mol)δE (kJ/mol)
5584.31.41.319611
10593.51.41.219610
15599.11.41.319610
20603.11.41.619710
30608.81.41.319610

Source, Table 15 of the source chapter. Activation energies and frequency factors for the second process on tempering of AISI 1050, with the exponent n0 fixed at 1.4. The frequency factor column is printed by the source in units of 1017 per minute, while the activation energy is given per mol.

The second process is the more reproducible of the two. Its activation energy is 196 to 197 kJ/mol at every heating rate, inside the error bound of about 10 kJ/mol, and the frequency factor varies by less than a factor of 1.4 across a sixfold change in heating rate. The first process scatters more, from 112 to 124 kJ/mol, which is what an exponent that is not better determined than plus or minus 0.2 would be expected to produce.

Schematic of the relative length change against temperature with the start, intermediate and end stages marked
Fig. 6. Definition of the fraction untransformed from a non-isothermal dilatometry curve. The measured relative length change at a temperature T is read against the increments at the start and the end of the reaction, so the fraction untransformed is the ratio of the two distances. Reproduced from the source figure.
Fraction untransformed against temperature for the two tempering processes at 5 K per minute
Fig. 9. Fraction untransformed against temperature for the two processes at the 5 K per minute heating rate. The points are experimental and the curves are calculated from the fitted activation energy, frequency factor and exponent. Reproduced from the source figure.

Calculating the fraction untransformed across the whole temperature range, using the fitted values, and plotting it against the experimental points is the check on the whole set. The two curves in Fig. 9 follow the points through both bands. The study notes that other nucleation and growth protocols, which produce a different value of n0, generate visible deviations between the calculated and experimental curves over the temperature range where the reactions actually develop, so the fit is sensitive to that choice.

What the numbers mean in a tempering furnace

Three practical points follow from the parameters above. The two tempering reactions are separated in temperature by roughly 170 K at these heating rates, which is why a single hold at a temperature below about 400 °C does only the first reaction, and why the second needs either a higher temperature or a long hold. The apparent activation energy of the second reaction, near 200 kJ/mol, is high enough that its rate changes sharply with temperature, so a 20 K change in the tempering set point moves it much more than it moves the first reaction. That is the quantitative reason a furnace that is 20 K off its set point still produces a usable tool when the temper is a low temperature one, and produces a large spread in hardness when the temper is being run to form cementite.

The second point concerns the measurement itself. Every number on this page was extracted from the change in length of a specimen, and the resolution of that measurement sets the noise in the fraction transformed, which in turn sets the uncertainty on the exponent. The uncertainties printed in Tables 13 to 15 are calculated from a length resolution of 10-4 and a temperature uncertainty of 0.1 K, so a laboratory reading a dilatometer with a poorer resolution would return wider bounds on the same steel rather than a different mechanism.

The third point is about what the kinetic parameters are not. They describe two reactions in a low alloy steel and in AISI 1050, not the secondary hardening that a molybdenum or vanadium bearing tool steel develops at higher tempering temperatures, which is a third reaction with its own kinetics and its own temperature range. For a tool steel, the tempering response that matters is the secondary hardening peak, and the hardness a given grade reaches at a given temperature is plotted on the tool steel tempering chart.

Source: Alloy Steel, Properties and Use, InTech, 2011.