2. Thermodynamic Characteristics of Chromium and Carbon Oxidation
The greatest losses of chromium occur during the melting of the scrap material in the EAF in the stainless steels production process. When we talk about losses, we mean the oxidation of chromium, whereby the formed chromium oxide binds and dissolve into slag and is difficult to reduce it. The source of oxygen during the melting process can be an oxidized, moist scrap material, air intrusion through the door of the EAF, excessive oxidation of the steel melt, or last but not least, an incorrect ratio between natural gas and oxygen [
3,
4,
5,
18,
19]. The reaction of chromium oxidation can be written as Equation (1):
A low carbon content should be reached when producing stainless steels, which can be most easily achieved by oxidation of the dissolved carbon. Oxidation of the dissolved carbon in liquid steel melt can be written as Equation (2):
From the given Equations, we can conclude that both carbon and chromium are oxidized during melting, and which reaction will preferentially take place depends on the activity and conditions in the furnace. Equations (1) and (2) can be combined and written as follows:
Now the equilibrium between chromium and carbon oxidation can be calculated.
Table 1 shows the values for the standard enthalpy (Δ
H°), entropy (Δ
S°), and Gibbs free energy (Δ
G°) at 1873 K for Equation [
3] according to different literature sources [
1,
4,
9,
14,
15,
16,
17,
20,
21]. When calculating Gibbs free energies, the initial and final states of elements and compounds are important, the reactions must be accurately described. Equation [
3] shows that carbon and chromium were/are dissolved in the liquid iron solution, therefore the Gibbs free energy of elements dissolution (Δ
Gdis) should be considered. According to some authors [
4,
9,
16,
20], only standard enthalpies and entropies of formation were considered in calculating the Gibbs free energy. Wang et al. [
15], Ma et al. [
17], and Heikkinen et al. [
14] on the other hand include in their calculations the energy used to dissolve C and Cr into the iron solution (Δ
Hdis, Δ
Sdis). When considering the required energy to dissolve C and Cr in an iron solution, the values from Sigworth et al. [
22] were taken, it was necessary to consider the standard state used (Raoultain, Henrian, and 1 wt. % standard state).
In
Table 1, calculated values based on given data from literature sources are marked with
1,2,3. When calculating the Gibbs free energies, Wang et al. [
15] and Ma et al. [
17] took into account the energy required to dissolve solids C and Cr and considered the 1 wt. % standard state in the calculations. Heikkinen et al. [
14] also took into account the dissolution energy of carbon and chromium when calculating the Gibbs free energy, except that they considered the dissolution energy of chromium in the liquid state, and performed the calculations according to the Henrian standard state.
Table 2 gives the values for the enthalpy (Δ
Hdis), entropy (Δ
Sdis), and calculated Gibbs free energy (Δ
Gdis) for the infinite dilution of C and Cr in liquid iron from different states at a temperature of 1873 K.
When calculating the Gibbs free energy (Δ
G), the reaction constant (
K) must also be taken into account, so we can write the expanded expression for calculating the Gibbs free energy:
Where
T and
R represent the temperature and the gas constant,
K3 is the constant of Equation [
3],
pCO represents partial pressure of carbon monoxide,
aCr,
aC,
are activities of chromium, carbon, and chromium oxide, respectively.
In multicomponent melts, as in alloy steels, the activity coefficient of solute
i is affected by the alloying elements for which Wagner’s formalism can be written. The activity of elements can be calculated using different standard states if we take the 1 wt. % standard state, Wagner’s formalism is written with the following expressions [
23]:
Where,
,
are the activity and the activity coefficient of element
i in the iron solution, using 1 wt. % standard state,
and
are the first-order self-interaction parameter of the element
i and first-order interaction parameter of element
j on element
i in dilute liquid iron, respectively. [wt. %
i] and [wt. %
j] is the weight percent of the dissolved element
i and
j, respectively. The values of interaction parameters
,
,
and
also vary in different references and are given in
Table 3 [
5,
19,
22,
24,
25,
26].
Figure 1 gives the equilibrium solubility of chromium and carbon in a liquid iron considering Gibbs free energies taken from the literature (
Table 1). In addition to the Gibbs free energies listed in
Table 1, the dissolution energies (Δ
Gdis) of the elements in the liquid iron melt were used to calculate the equilibrium solubility (from the solid state to dissolved in liquid iron, taking into account 1 wt. % standard state). As the summation of the Gibbs free energy increases, the equilibrium solubility of chromium in the liquid iron melt decreases, thus, for example, the equilibrium concentration of chromium at 0.8 wt. % C and 1600 °C according to Carboni et al. [
16] is 22.35 wt. % Cr and according to Heikkinen et al. [
14] is 16.85 wt. % Cr. The values for enthalpy and entropy were taken from the HSC 8.0 software database for further calculations. The curve obtained from these values shows an approximate mean value. With data from HSC 8.0, the equilibrium value of chromium is 19.60 wt. % Cr at 0.8 wt. % C and 1600 °C.
Figure 2 shows the equilibrium concentration curves of C and Cr using the Gibbs free energy from the HSC 8.0 database at 1600 °C and different values of the interaction parameters
,
,
and
. When choosing the interaction parameters, there are also discrepancies between the authors, as can be seen from
Table 3. The result of using different interaction parameters for Equation [
3] are curves with different slopes, as can be seen from
Figure 2. The equilibrium solubility of chromium is the highest when interaction parameters according to Ghosh et al. [
5], Seetharaman et al. [
19] and JSPS (The Japan Society for the Promotion of Science) [
25] were used. The lowest equilibrium solubilities are given by calculations using interaction parameters according to Total Materia [
26]. For further calculations, the interaction parameters according to Sigworth et al. [
22], whose curve in
Figure 2 shows an approximate median value.
In metallurgy, it is practical to calculate activity according to the 1 wt. % standard state, as the concentrations of elements in solutions are usually given in weight percentages. Using the 1 wt. % standard state is much more practical from the point of view of engineering practice since no conversion and adjustment of parameters is required. However, for more accurate calculations it is necessary to calculate with mole fractions (
X). Wagner [
23] and later Lupis and Elliott [
27] proposed the following notation for calculating activity coefficients:
Where
, are the activity and the activity coefficient of element
i, using Henrian standard state,
is the zero-order activity coefficient,
is the first-order interaction parameter of the influence of element
j on element
i for calculation taking into account mole fractions (
X),
Xi,
Xj are the mole fractions of element
i and
j in the solution, respectively. The first-order interaction parameter
can be calculated from the interaction parameter of the 1 wt. % standard state [
22,
27]:
Where Mj is the molar mass of element j, M1 is the molar mass of the solvent element of the solution, and is the first-order interaction parameter of the 1 wt. % standard state.
Expression (9) describes well the state of diluted solutions which is described by Henry’s law, but the activity and activity coefficient deviate as the concentration of alloying elements in the solution increases. Therefore, Pelton and Bale [
28,
29], Ma et al. [
17] proposed modification of the interaction parameter formalism and added the solvent activity coefficient (
) to Equation [
9] and thus enabled a more accurate calculation of the activity at higher concentrations. Pelton and Bale [
28,
29] wrote the following expressions:
Figure 3 shows the equilibrium solubility curves of chromium and carbon in iron solution, taking into account 1 wt. % standard state, Henrian standard state, and modified Henrian standard state for calculating the activity coefficient of the elements. The results were calculated according to the given expressions and then converted to percent by weight. When calculating with the 1 wt. % standard state, the solubility of chromium is higher than in the case of calculation with mole fractions. When calculating with mole fractions, there is no significant difference in the equilibrium solubility of Cr and C at low concentrations (up to 0.2 wt. % C), but with an increase in the carbon concentration, the equilibrium solubility of chromium increases faster, taking into account the solvent activity coefficient. Direct links between the solubility of carbon and chromium using the ThermoCalc program were unable to be calculated. The data is protected, and such a display is not possible.
Figure 4a,b show the dependence of carbon and chromium activity on the molar fraction of the element in the iron solution with 15 wt. % Cr or 1 wt. % C at 1600 °C. Activities were calculated according to Raoultain, Henrian, and 1 wt. % standard state. In the calculation according to Pelton and Bale [
28,
29] and Ma et al. [
17], it involves modified calculations based on Henry’s standard state, taking into account the activity of the solvent. The diagrams also show the activity of the C and Cr in the solution determined by the ThermoCalc program using the TCFE12 database. The activity of carbon in the melt is challenging to determine and poses a problem, as indicated by the results obtained using the ThermoCalc program. Determining activity according to the standard state of 1 wt. % is exceptionally imprecise; activities sharply increase even at low concentrations of both C and Cr. Calculations of activity using other models yield values that are comparable to each other.
Based on
Figure 3 and
Figure 4, it can be observed that as the solubility of C decreases, the solubility of Cr also decreases. In other words, this implies that during oxidation, in addition to the combustion of carbon, there is also the combustion of chromium. Stainless steels typically contain a low concentration of carbon, necessitating intensive oxidation, leading to significant chromium losses. These chromium losses can be reduced by selecting an appropriate oxidation temperature and configuration, as outlined in the following. Calculations of the equilibrium solubility of chromium and carbon as a function of temperature, the activity of chromium oxide, the partial pressure of carbon monoxide, and the activity of dissolved elements in the stainless steel melt were made to show the influence of individual parameters on the solubility of chromium and carbon in the steel melt. Based on the results of the calculations the process window for managing the production of stainless steel melt with as little chromium loss as possible and greater process and economic efficiency can be determined. In the following calculations, the enthalpy and entropy values from the HSC 8.0 database were used, and the interaction parameters
and the dissolution energies of the elements were calculated from the data taken according to Sigworth et al. [
22].
Figure 5 shows the equilibrium solubility curves of chromium and carbon between temperatures 1550 and 1800 °C. As the temperature increases, the equilibrium concentration of carbon decreases. From a thermodynamic point of view, the influence of temperature is logical, since at higher temperatures gaseous phases (CO) are more stable, while metal oxides (Cr
2O
3) are less and less stable and break down into lower oxidation states. A large difference appears in the equilibrium values for the lowest and highest calculated temperature at 15 wt. % Cr, which amounts to 1.21 wt. % C at 1550 °C and 0.28 wt. % C at 1800 °C. Decarburization of steels with an elevated chromium concentration in the EAF is advisable to conduct at an elevated temperature to minimize chromium losses. When melting scrap material, which already contains a high percentage of chromium, precautions should be taken to prevent oxidation.
Figure 6 shows the curves of equilibrium concentrations of C and Cr for different activities of chromium oxide. As the proportion of chromium oxide increases in the slag, its activity increases until it begins to precipitate in solid form and its activity is equal to 1. The equilibrium solubility of chromium in liquid iron is lower at lower activities of chromium oxide, which can be concluded that in the initial stages of melt oxidation, due to the lower the activity of chromium oxide, the more intensive oxidation of chromium takes place. At the smallest calculated activity of chromium oxide,
= 0.01, the equilibrium value of carbon at 15 wt. % Cr and 1600 °C is equal to 2.24 wt. % C, on the other hand the concentration of carbon at the activity of chromium oxide
= 1 is equal to 0.91 wt. % C.
Figure 7 shows the influence of the partial pressure of carbon monoxide (CO) on the equilibrium solubility of C and Cr. By constantly removing the gaseous product from the system, which is achieved by lowering the pressure, the direction of the reaction is favoured towards the formation of the gaseous phase. In our case, by lowering the partial pressure of CO, the carbon oxidation is accelerated and the solubility of chromium in the liquid iron is increased. In the EAF, it is not possible to influence the gas pressure above the level of the melt, but the procedure is used in subsequent steps when the steel melt is processed in VOD and AOD converters.
Figure 7 shows that the equilibrium solubility of chromium increases drastically with decreasing
pCO, which is why the oxidation of carbon in the further steps of processing the steel melt is extremely fast and efficient.
From the notation of Wagner’s formalism, it is clear that by mixing alloying elements into the solution, the activity of the selected elements is directly influenced. Silicon has a positive influence on the activity of carbon
= 0.08 and a negative influence on the activity of chromium
= –0.0043 [
22]. Increasing the concentration of silicon in the solution on the one hand increases the activity of carbon, while on the other hand, the activity of chromium is decreased, therefore, increasing the concentration of silicon in the iron solution increases the solubility of chromium in the melt and directs the balance of Equation (3) towards the oxidation of carbon.
Figure 8 shows the effect of silicon concentration on the equilibrium solubility of C and Cr in the iron solution. At the maximum calculated concentration of silicon, i.e., 2 wt. % Si, and 15 wt. % Cr at 1600 °C, the equilibrium solubility of carbon is 0.72 wt. % C, at 0 wt. % Si and the equilibrium value of carbon is 0.91 wt. % C at 15 wt. % Cr.
The value of the interaction parameter of oxygen to carbon and oxygen to chromium is comparable and negative in both cases (
= –0.34 and
= –0.14 [
22]). Therefore, the addition of oxygen does not have much effect on changing the equilibrium solubility of carbon and chromium in the iron solution.
Figure 9 shows the calculated equilibrium solubility curves of chromium and oxygen at different oxygen concentrations in the solution.
In addition to the interaction parameters and activity of the elements in the solution, the oxidation reactions of other elements and their Gibbs free energies should also be discussed.
Figure 10 shows the curves of Gibbs free energies as a function of temperature for the oxidation of various alloying elements. The diagram in
Figure 10 indicates that the oxidation of carbon increases with temperature, as confirmed by the diagram in
Figure 5. Decarburization of the stainless steel melt is advisable to carry out at high temperatures. The Gibbs free energies for the oxidation of silicon and aluminium are lower compared to the energy for the oxidation of chromium, meaning that the oxidation of chromium from the melt is reduced in the presence of dissolved Si and Al.
Figure 10 shows only the dependence of temperature on the Gibbs free energy of reactions, so the results only indicate which elements have a greater affinity for oxygen. For more accurate results, it is mandatory to take into account the activities of the dissolved elements in the solution.
Based on the presented calculations, we can drag the guidelines for more efficient production of melting steels with an elevated chromium concentration. In the initial phase, at low temperatures and low activities of chromium oxide in the slag, there is a high probability of chromium oxidation. Oxidation can be reduced by protecting the scrap with the addition of carbon (anthracite), silicon, and aluminium, and it is necessary to add an appropriate amount of CaO and MgO to obtain slag with suitable properties, namely slag basicity, viscosity, and reactivity with refractory materials [
2,
3,
10,
30]. Secondly, it is necessary to limit the access of oxygen, so it is recommended that melting does not take place with the help of natural gas burners. Large losses of chromium in the initial phase of melting can also be avoided by carefully stocking the baskets, if possible, by stocking the initial baskets of the scrap with the lowest possible concentration of chromium. As can be seen from
Figure 5, as the temperature of the melt increases, the solubility of chromium increases, so it makes sense that the decarburization of the steel melt takes place at temperatures as high as possible (normally higher than 1650 °C), and the carbon concentration at the end of decarburization should not be too low, because the oxidation of chromium also takes place in parallel. Decarburization is more effective at a lower partial pressure of CO, which can only be ensured in VOD and AOD aggregates. The kinetics of the decarburization of the steel melt can be dynamically controlled by off-gas analysis. Echterhof et al. [
30] presented a method in which they estimated the kinetics of melt decarbonization based on the proportion of CO
2 in off-gases. According to the authors, in the first stage, the proportion of CO
2 in the off-gases increases with the amount of injected oxygen, which can be connected to the oxidation of silicon from the melt and the increasing trend of carbon oxidation. In the second stage, intensive decarbonization takes place, and the amount of CO
2 in off-gases reaches a plateau. The authors associate the drop of CO
2 fraction in the off-gases with the reduced kinetics of decarburization and increased activity of chromium in stainless steel melt, and therefore, the beginning of more intense chromium oxidation, the oxygen injection should be stopped at this point. Straightforward and practical methods were implemented to reduce the chromium losses in steel plants. Additionally, using decarburization control during EAF melting, the oxygen consumption and production time were reduced.
Based on the thermodynamic data from the HSC 8.0 database and the parameters listed in Sigworth et al. [
22], a model was created to calculate the decarburization of the stainless steel melt, also partially presented in our previous work [
31]. The model defines the oxidation of elements based on the amount of oxygen injected into the solution and the Gibbs free energies of dissolved elements.
Figure 11a–c show the results of the amount of oxidized C, Cr, and Si during the oxidation of the solution with 15 wt. % Cr and 2 wt. % C.
Figure 11a shows the amount of oxidized carbon and chromium for the solution at the oxidation start temperature of 1500 °C. In the first stage, oxidation of carbon takes place, and above 3.5 kg of injected oxygen per ton of steel, oxidation of chromium also occurs. At 20 kg of injected oxygen per ton of steel, 7 kg of chromium, and 13 kg of carbon per ton of steel oxidized.
Figure 11b shows the amount of chromium and carbon oxidized at the oxidation start temperature of 1650 °C. Due to the higher solubility of chromium, as a consequence of the higher temperature of the system, carbon oxidizes for a longer time, and oxidation of chromium occurs at over 16 kg of injected oxygen per ton of steel. At 20 kg of injected oxygen per ton of steel, it oxidizes 2.5 kg of chromium and 13.5 kg of carbon per ton of steel. In
Figure 11c, the steel melt at 1500 °C contained 0.5 wt. % Si. The results show that the addition of silicon protects the melt well from chromium oxidation, despite the low temperature, chromium oxidation occurs above 17 kg of injected oxygen per ton of steel. In this case, silicon oxidation takes place in the first stage, followed by the second, more intense phase of carbon oxidation. Above 17 kg of injected oxygen per ton of steel, chromium also oxidizes. At 20 kg of injected oxygen per ton of steel, it oxidizes 0.3 kg of chromium, 12 kg of carbon, and 3 kg of silicon. The model results approximate the thermodynamic oxidation calculations of Echterhof et al. [
30]. The calculated model confirms that a high melt temperature is necessary for effective decarburization with low chromium losses, and secondly, the model shows that it is possible to protect chromium losses at low temperatures by adding silicon.
Slag also plays an important role in reducing chromium losses in the production of stainless steels. In general, slag protects the melt from oxidation, loss of heat and material, and additionally, slag enables better efficiency of electricity and increases the productivity of the process. With the presence of foaming slag, the productivity and energy efficiency of the process increase even more. Due to the formation of chromium oxides, the formation of foaming slag in the production of stainless steels can be difficult. At low concentrations of chromium oxides (up to 8 wt. %), there are no significant differences, above 13 wt. % of chromium oxides, the viscosity of the slag decreases at the expense of an increase in the concentration of solid particles. With an increased concentration, the foaming index decreases and makes the formation of foaming slag more difficult. Another obstacle to the formation of foaming slag is the low concentration of FeO in the slag, as a result of which the kinetics of the formation of gaseous CO is low, and the amount or height of the foaming slag is correspondingly lower [
2,
18,
32,
33,
34].
The decarburization phase of the melt is usually followed by the reduction of the slag with silicon and aluminium, with the aim of reducing the chromium oxides in the slag. Slag reduction is effective in the case of intensive mixing of melt and slag, which usually occurs during tapping both, steel and slag into the ladle [
9,
10,
11,
12,
13,
16].