Revisiting the Rist Diagram for Predicting Operating Conditions in Blast Furnaces With Multiple Injections
Revisiting the Rist Diagram for Predicting Operating Conditions in Blast Furnaces With Multiple Injections
Revisiting the Rist Diagram for Predicting Operating Conditions in Blast Furnaces With Multiple Injections
RESEARCH ARTICLE
1 2
Abstract
Background: The Rist diagram is useful for predicting changes in version 1
blast furnaces when the operating conditions are modified. In this 29 Nov 2021 view view
paper, we revisit this methodology to provide a general model with
additions and corrections. The reason for this is to study a new
1. Johannes Schenk , Montanuniversität
concept proposal that combines oxygen blast furnaces with Power to
Gas technology. The latter produces synthetic methane by using Leoben, Leoben, Austria
renewable electricity and CO2 to partly replace the fossil input in the
2. Eleonora Bargiacchi, IMDEA Energy
blast furnace. Carbon is thus continuously recycled in a closed loop
and geological storage is avoided. Institute, Madrid, Spain
Methods: The new model is validated with three data sets
corresponding to (1) an air-blown blast furnace without auxiliary Any reports and responses or comments on the
injections, (2) an air-blown blast furnace with pulverized coal injection article can be found at the end of the article.
and (3) an oxygen blast furnace with top gas recycling and pulverized
coal injection. The error is below 8% in all cases.
Results: Assuming a 280 tHM/h oxygen blast furnace that produces
1154 kgCO2/tHM, we can reduce the CO2 emissions between 6.1% and
7.4% by coupling a 150 MW Power to Gas plant. This produces 21.8
kg/tHM of synthetic methane that replaces 22.8 kg/tHM of coke or 30.2
kg/tHM of coal. The gross energy penalization of the CO2 avoidance is
27.1 MJ/kgCO2 when coke is replaced and 22.4 MJ/kgCO2 when coal is
replaced. Considering the energy content of the saved fossil fuel, and
the electricity no longer consumed in the air separation unit thanks to
the O2 coming from the electrolyzer, the net energy penalizations are
23.1 MJ/kgCO2 and 17.9 MJ/kgCO2, respectively.
Discussion: The proposed integration has energy penalizations
greater than conventional amine carbon capture (typically 3.7 – 4.8
MJ/kgCO2), but in return it could reduce the economic costs thanks to
diminishing the coke/coal consumption, reducing the electricity
consumption in the air separation unit, and eliminating the
Page 1 of 30
Open Research Europe Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
Keywords
Blast furnace, Rist diagram, ironmaking, operating diagram, Carbon
capture, Power to Gas, Oxyfuel combustion, CO2
Page 2 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
Page 3 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
The second major novelty of the paper is analyzing for the first
time the OBF–PtG integration under the operating diagram
methodology and, therefore, by using consistent operation data
sets. Besides, the operating lines of these blast furnaces are
obtained, which cannot be found elsewhere in literature. The
third major novelty is to provide full operation data sets for dif-
ferent blast furnace, with detailed composition of all streams and
their most relevant operating parameters (e.g., temperature of
the thermal reserve zone, heat evacuated by the staves, and the
temperature of the flame). The availability of this information
in literature is very scarce, especially for OBFs.
CO 2 + C → 2CO ∆ H = + 172.47 kJ /mol (9) can be considered that sulfur dissolves into the hot metal
by Eq.(18), and then transfers to the slag by Eq.(19)20.
H 2 O + C → H 2 + CO ∆ H = + 124.48 kJ /mol (10)
Fe + 0.5S2 → FeS (18)
Apart from the direct reduction, other relevant processes
occur in the lower zone. The burden contains various impuri- FeS + CaO + C → Fe + CaS + CO (19)
ties that will either dissolve in iron or will form part of the slag.
For example, the Al2O3, CaO and MgO oxides are not reduced Per mole of S that ends in the slag, 1 mole of CO is added to
under the blast furnace conditions and therefore they trans- the reducing gas.
fer fully into the slag. In the case of SiO2, MnO and P2O5,
they are partially reduced and dissolved in the hot metal. It 2.2 Temperature profile
can be considered that these impurities are directly reduced The blast furnace process can be divided in three different tem-
by solid carbon following Eq.(11), Eq.(12) and Eq.(13). The peratures zones (Figure 1). At the lower part, the flame tem-
final silicon, manganese and phosphorous contents in the hot perature is normally between 2000 and 2300°C (defined as
metal are much less than the equilibrium values. the temperature reached by the raceway gas when all C and
H2O have been converted to CO and H2)22. This raceway gas
SiO 2 + 2C + 3Fe → Fe 3Si + 2CO ∆ H = + 622.5 kJ /mol (11) provides heat for the direct reduction process and for the
melting of the hot metal and slag20. The hot metal exits at
MnO + C → Mn + CO ∆ H = + 286.92 kJ / mol (12) 1350–1450°C and the slag at 1500–1550°C16. The gas, which
has been cooled to about 1000°C, ascends to the middle zone.
P2 O5 ⋅ 3CaO + 5C + 6 Fe →
∆ H = + 1184 kJ / mol (13) The middle zone is a region of thermal equilibrium. In prac-
2 Fe 3 P + 3CaO + 5CO tice, a non-zero temperature difference remains between
Another component that will end up dissolved in the hot metal gas and solids, but it passes through a minimum value in a
is carbon. As in the previous case, carbon never reaches satu- region of slow heat exchange. The temperature is kept almost
ration in pig iron20. At tapping temperatures (1350–1450°C) constant around 800–1000°C10.
the carbon content may vary from 2.5 to 4.5%16. As simplifica-
Lastly, in the upper zone, the gas and the burden exchanges
tion, it is assumed that the dissolved carbon forms, austenite and
heat rapidly. The gas is cooled down from 800–900°C to
cementite, in the hot metal, according to Eq.(14) and Eq.(15),
100–200°C as it leaves the furnace top, and the burden is heated
respectively11.
from ambient temperature to 800°C while descending16,20.
C(coke) → C(austenite) ∆ H = + 34.7 kJ / mol (14) The temperature profile and the reduction zones (pre-reduction,
indirect reduction, and direct reduction) more or less coincide,
3 Fe + C(coke) → Fe 3 C ∆ H = + 6.69 kJ / mol (15) so these three zones can be used to study the blast furnace.
Also, in the lower zone, in front of the tuyeres (i.e., the race- 3. Generalized Rist diagram with multiple
way), coke burns with the oxygen of the hot blast, thus pro- injectants
viding the process with heat and CO-reducing gas. The total The Rist diagram is named so in reference to its author, who
reaction of coke in the raceway can be considered as an incom- elaborated a model for predicting changes in blast furnaces
plete combustion due to the shortage of oxygen ((Eq.(16))16. when the operating conditions are modified9–13. The model is
Actually, in the inner part of the flame, complete combustion based on the graphical representation of carbon, oxygen, and
also occurs but the CO2 ends up dissociating by Eq.(9). hydrogen balances through an operation line that depicts the
participation of these elements in the formation and utiliza-
1 (16) tion of the reducing gas9. Additionally, the diagram includes
C+ O → CO ∆ H = − 113.68 kJ / mol
2 2 an equilibrium line to delimit the maximum oxidation state of
the gas according to the Chaudron diagrams for the Fe-O-H and
In case of injecting auxiliary fuels to diminish coke consump-
Fe-O-C systems23.
tion, incomplete combustion is assumed to follow Eq.(17),
where Z denotes the ashes in the case of pulverized coal. No
water is present since it rapidly dissociates by Eq.(10) as it The construction of the Rist diagram is introduced here with
occurred for CO216. additions and corrections with respect to the original work of
Rist. The model methodology is now described for the gen-
1 eral case of multiple injectants treated separately (instead of
CH 2 aO 2 b N 2 c S 2 d Zz + O → CO + aH 2 + bO 2 + cN 2 + dS2 + zZ (17) for an overall single injection). This is especially important
2 2
because it will allow to properly calculate the heat of decom-
When the injected fuel contains sulfur (e.g., pulverized coal), position of each auxiliary fuel, as well as to specify different
this will end forming part of the slag. As simplification, it inlet temperatures for each injectant. Additionally, a detailed
Page 5 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
x v + 2 xe + ∑ ( a j + 2b j ) x j + xS i + xM n + xP + xS + xk + xd
(20)
= 1 mol of reducing gas
Page 6 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
The corresponding terms of Eq.(20) and Eq.(22) form sets of through calcination)12. Moreover, it should be noted that such
proportional numbers and, as such, they can be read on two a correspondence between the abscissa and the composition
rectangular axes as the projections of the same straight line. of the mixture does not exist in the interval 0<X<1, where the
The Figure 2 shows the straight line thus obtained, called the segments could be arranged in any order.
“operating line”9, whose slope has the units of Eq.(25).
3.3 Equation of the operating line
C + H2
µ= (25) According to the theory described above, the equation of the
Fe operating line can be written as Eq.(28). The slope, μ, is the
This corresponds to the total number of moles of reducing gas, number of moles of reducing gas required for the produc-
per unit of Fe, according to the sources of hydrogen and carbon. tion of 1 mole of Fe. The intercept, YE, represents the moles of
H2 and O coming from sources other than iron oxides that con-
3.2 Utilization of the reducing gas (Rist diagram in the tribute to the formation of the reducing gas (negative sign by
range 1<X<2) convention, Eq.(23))12.
Following the same principle, it is simple to depict the utiliza-
tion of the reducing gas mixture (Eq.(1) – Eq.(6)) by means of a Y = µ ⋅ X + YE (28)
segment representing the oxygen removed from the iron oxides.
The number of oxygen moles transferred from the iron oxides If the operating line is characterized, relevant information can
to the reducing gas is denoted by xi when referred to 1 mole be deduced from it. The slope accounts for the total reducing
of reducing gas mixture, and by yi when referred to 1 mole agent rate required (in terms of C and H2 per mole of Fe) so, if
of Fe (same convention than before). Thus, the ratio an auxiliary fuel is introduced, the decrease in the input rate
between both variables (i.e., yi /xi) is the number of moles of of coke can be computed5. The intercept stands for the hydro-
reducing gas mixture per mole of Fe produced. Since the gen and oxygen brought into the furnace (except for the O2 con-
reduction of iron oxides does not change the total number of tained in the iron ore), therefore the necessary air flow rate can
moles of the reducing gas, the ratio yi /xi is constant and equal to be calculated by subtracting the other O2 and H2 sources (mois-
μ. In other words, in all equations from Eq.(1) to Eq.(6) the gas ture, auxiliary fuels, coke and impurities)5. Also, the initial
gets oxidized without increasing or decreasing the number of oxidation state of the iron oxides introduced in the blast furnace
moles in the gas (1 mole of CO gives 1 mole of CO2, and 1 mole (YA) allows to know the final degree of oxidation of the gas leav-
of H2 gives 1 mole of H2O). Therefore they can also be read as ing the top of the furnace (XA – 1). Finally, the ratio between
projections of the same straight line of slope μ (segment in direct and indirect reduction is identified by construction. The
the range 1<X<2 of the Rist diagram, Figure 2)9. abscissa X = 1 gives the oxygen removed by direct reduc-
tion, yd (Figure 2), and then the oxygen removed by indirect
This segment in the range 1<X<2 is particularly useful because reduction is easily calculated as yi = YA – yd.
it provides information on the average oxidation state of
the reducing gas mixture (abscissa) and of the iron oxides 3.4 Equilibrium line
(ordinate)12. The abscissa can be interpreted as the reducing gas In the blast furnace, the reducing gas can never oxidize the sol-
having an average oxidation state equal to X – 1 (see Eq.(26)). ids. For this reason, the operating line in the segment 1<X<2
This means that at the abscissa X = 1, we find a reducing must necessarily remain on the left of the equilibrium line of
gas mixture composed by CO and H2, while at the abscissa the Fe-O-H-C system (Figure 2). If we were at some point at
X = 2 the gas is completely oxidized to CO2 and H2O. On the the right of the equilibrium line, the way to reach equilibrium
other hand, the Y coordinates represent the oxidation state would be displacing us upwards (i.e., providing O to the Fe)
of the iron oxides according to Eq.(27). This interpretation or leftwards (removing O from the gas), what in both cases
allow identifying the point A in the Rist diagram, whose ordi- means to oxidize the solids.
nate is the initial oxidation state of the burden (e.g., YA = 1.5
for Fe2O3)12 and whose abscissa is the final degree of The contour of the equilibrium line is delimited by five points,
oxidation of the gas leaving the top of the furnace plus 113. which we will denoted by F, W, W2, M, M2. The point W is of
special interest since it corresponds to the chemical equilib-
CO 2 + H 2 O
X =1+ (26) rium between gases and solids at the beginning of the middle
CO + CO 2 + H 2 + H 2 O zone, where pure wüstite is found if the blast furnace oper-
ates under ideal conditions. The coordinates of these points are
O
Y= (27) given in Table 1. The ordinates of the five points are easily cal-
Fe
culated as the ratio of the oxygen and iron atoms of the corre-
The necessary condition to understand the segment 1<X<2 sponding components (Eq.(27))24. The abscissae (Xw Eq.(29)
in this way is that the total number of moles of reducing gas and XM Eq.(30), which are equivalent to Eq.(26)) depend on
keeps constant and that the oxygen supplied to the gas must the molar fraction Xh that relates the hydrogen and water con-
come only from the reduction of the iron oxides (i.e., no addi- tent of the reducing gas mixture (Eq.(31)). This is used to com-
tional injections in the middle or upper zone, and no CaCO3 bine the state of oxidation at equilibrium ω for the individual
introduced, which would decompose into CaO and CO2 CO-CO2 and H2-H2O mixtures (Figure 3). It should be noted
Page 7 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
Table 1. Delimiting points of the equilibrium that the molar fraction Xh is independent of the state of oxi-
line for the Fe-O-H-C system in a Rist diagram. dation of the reducing gas (i.e., independent of the abscissa
X), and can be calculated as a function of ye, yk, ajyj and μ.
Abscissa X Ordinate Y
XW = 1 + (1 − x h ) ωWC+ x h ωWH (29)
F Iron Fe XF=XW YF=0
2
ω = a 0 + a1 T + a 2T + a 3T 3+ a4 T 4 + a 5T 5 (32)
Table 2. Coefficients for the calculation of the oxidation state at equilibrium in Fe-O-C and Fe-O-H systems
(Eq.(32)). Valid from 500 to 1100°C.
Page 8 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
known, it is easy to compute the slope and the Y-intercept of the In practice, a blast furnace does not operate under ideal condi-
operating line (Eq.(33) and Eq.(34)): tions, so a zone of pure wüstite cannot be distinguished (the
operating line no longer passes through the point W). For these
YR − YP
µ= (33) cases, it is defined the chemical efficiency of the furnace, r,
XR − XP representing the oxygen actually exchanged to the oxygen the-
oretically exchangeable (typically, r is around 0.92)5,10. The
Y − YP
YE = YR − XR R (34) coordinates of the new point R, through which the operation
X R − XP line passes, is calculated by Eq.(35) and Eq.(36) as a function
of the point W, the chemical efficiency r, and the initial
In the following subsections, it is explained how to calcu-
oxidation state of the burden YA.10.
late the points R and P. The additions and corrections to the
original work of Rist are introduced mainly during the calcu-
X R = 1 + r ( X W − 1) (35)
lation of the energy balance (explained more in detail in the
appendixes25). These are the following:
YR = YA + r ( YA − YW) (36)
▪ ll calculations are given as a function of the temperature of
A
the middle zone, TR (i.e., of the temperature of the thermal Despite there is no chemical equilibrium, the temperature is
reserve zone where thermal equilibrium exists). still nearly constant in the middle zone (thermal reserve zone
▪ ach auxiliary fuel can enter at different temperatures
E
where thermal equilibrium exists), so this temperature is used
for the calculations (TR).
because they are now treated separately, rather than as an
overall injection.
4.2 Point P: energy balance in the elaboration zone
▪ he heat associated to the direct reduction of FeO is now
T The operating line, Eq.(28), depends on different parameters
corrected taking into account the moisture of the hot blast. which are not all independent. In particular, any energy balance,
whether global or partial, imposes a relationship between them.
▪ he sensible heat of the hot metal is calculated as a function
T An option is to use the energy balance of the elaboration zone
of its composition. (i.e., middle zone plus lower zone), which follows Eq.(37).
▪ he sensible heat of the slag is calculated as a function of
T
As in the original work of Rist, we are going to work in kcal
per mole of Fe10. Moreover, the reference temperature for the
its SiO2, Al2O3, CaO and MgO content.
energy balance is chosen at TR, i.e., at the temperature of the
▪ he heat of carburization is calculated as a function of
T thermal reserve zone.
the austenite and cementite content in iron.
qc yv + qv yv + ( qiwYW − δ ) = qg ( yd − ye ) + qk yk + qe ye +
▪ he heat associated to the lack of thermal ideality is calculated
T
∑ qj yj − δ ′ + qSi ySi + qM n yM n + (37)
as a function of the burden and coke composition.
qP yP + qγ γ + f + l + p + C∆ TR
▪ he heat associated to the lack of chemical ideality is
T
calculated as a function of the burden in the thermal The input/production of energy is on the left side of the equa-
reserve zone. tion, and the output/consumption is on the right side of the
equation. Regarding the former, we have the term qcyv that
▪ he heat associated to the lack of chemical ideality now
T represents the heat released by the incomplete combustion
accounts also for the hydrogen coming from auxiliary fuels and of carbon, qvyv for the sensible heat of the air, qiwYw as the heat
for the moisture of the hot blast. released by the reduction of wüstite, and δ to account for the
▪ supplementary model is added to compute the heat of
A
lack of chemical ideality in wüstite reduction. In the right
side of the equation we have qg (yd – ye) which is the heat
decomposition of coal when injected as auxiliary fuel.
absorbed during direct reduction because of CO2 dissociation,
▪ n additional energy balance in the upper zone of the blast
A qkyk as the heat consumed due to the hydrogen entering with
furnace is added to compute the final composition of the the coke, qeye for the overall heat absorbed by the moisture of
blast furnace gas, instead of only computing the final the air, ∑ q j yj to quantify the overall heat absorbed by the
oxidation state. auxiliary injections, δ ′ accounting for the lack of chemical ide-
ality in the conversion of H2 to H2O, (qSiySi + qMnyMn + qPyP)
4.1 Point W and point R: chemical and thermal reserve for the heat absorbed by the reduction of the accompanying ele-
zones ments SiO2, MnO and P2O5, then qγγ as the heat absorbed by
In an ideal blast furnace, the reducing gas and the burden are the carburization of the iron, f denoting the heating and melt-
in chemical equilibrium after the upper zone. This is known ing of the hot metal, l for the heating and melting of the
as the chemical reserve zone, in which wüstite is the only iron slag, p for the heat removed by the staves and, lastly, (CΔTR)
oxide present. This point is denoted by W and was already accounting for the lack of thermal ideality (in case the tem-
identified during the construction of the equilibrium line perature of the gas and the solid is not the same at the thermal
(Table 1). reserve zone).
Page 9 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
A = qc + qv − eqe + eqk (1 − r ) + eqg (39) It is worth to mention that both coordinates of P depend on the
chemical efficiency of the furnace, r, through A and C. How-
B = qg (40) ever, Rist only considered the dependence on the chemi-
cal efficiency through C because he did not decompose δ ′.
Therefore, he found that Xp was independent of r10. Here we
(
C = − qiw Yw + δ + qk yk + ∑ q j yj − (1 − r ) yk + ∑ a j yj qk ) (41)
have shown that is not.
+ qS i yS i + qM n yM n + qP yP + qγ γ + f + l + p + C∆ TR
4.3 Additional results derived from the operating line
Here, the variable e denotes the moles of H2O per O moles in The relevance of characterizing the operating line comes
the air (i.e., e = ye/yv), which is a convenient notation to write from the possibility of deducing operational data such as the
the water injected with the air as a function of the air injected required reducing agent rate, the air consumption, the top gas
through the tuyeres (i.e., to include the term qeye in A). Simi- composition, the ratio between direct and indirect reduction,
larly, the term δ ′ was decomposed according to Eq.(92) to and the flame temperature.
separate the terms that depends on yv (see Appendix A.925).
It should be noted that Rist did not decompose δ ′ and there- 4.3.1 Reducing agent rate (coke consumption). The reduc-
fore he wrongly included a term that is actually dependent ing agent rate, μ, obtained from the operating line, denotes
on yv in C instead of adding it to A10. Moreover, Rist did not the number of moles of C and H2 needed inside the blast fur-
included the term +eqg that appears in A, which corrects the heat nace as reducing gas for the production of 1 mole of Fe as hot
absorbed during direct reduction (part of the direct reduction metal. When solving the Rist diagram, we assume the auxil-
takes place through Eq.(10), see Appendix A.5). iary injections to be known. Therefore, we can compute the
required amount of coke by subtracting the contributions of
Now, we can impose the relation of Eq.(38) from the energy the injections to μ. The carbon that ends up dissolved in the hot
balance to the operating line defined by Eq.(28). To do that, metal must be also taken into account (which increases the
we first rewrite Eq.(28) as Eq.(42), taking into account that required reducing agent rate), as well as the H2 that is pro-
μ = –YE + yd according to Eq.(22) and Eq.(23). duced when the H2O from the hot blast is dissociated (which
decreases the required reducing agent rate). Thus, the mass
Y = − YE (X − 1) + yd X (42) flow of coke is calculated by Eq.(49).
Page 10 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
ΩK,C ΩK,H Most authors use the CO utilization ratio (ηCO ≡ nBFG,CO2/(nBFG,CO
�K = 103 + (51) + nBFG,CO2)) as the 15th equation to complete the system of
M C M H2 equations, instead of using the energy balance in the upper
The term nHM,Fe is the number of moles of Fe in hot metal per zone. They assume to know this parameter, since in practice
ton of hot metal (Eq.(74)). they would be able to measure it at the flue gas of a real iron-
making plant. Indeed, in those cases in which the operating
4.3.2 Air flow rate. The intercept of the operation line, denoted line is characterized for a real operation where the ηCO infor-
by YE, represents the H2 and O brought into the furnace (except mation is reliable, to use this procedure instead of the energy
for the O contained in the iron ore) with negative sign by con- balance of the upper part makes the methodology easier.
vention. By subtracting the O and H2 sources other than the However, this is not a valid procedure when analyzing new oper-
hot blast, the necessary O flow rate as air can be calculated ating lines for potential blast furnace configurations. In such
(i.e., calculation of yv, Eq.(52)) case, the assumption of the value of ηCO would be arbitrary
and, almost certainly, it will lead to inconsistencies in the upper
−(YE + ∑ (a j + 2b j ) y j + ySi + yMn + yP + yS + yk ) zone (energy balance not fulfilled). Inconsistencies which are
yv = (52)
1 + 2e not detected because of not performing the corresponding
energy balance in the upper zone. For this reason, it is com-
Once we know yv, the mass of air (dry) per ton of hot metal is pletely necessary to use the energy balance of the upper zone
calculated by Eq.(53). instead of ηCO when proposing new operating lines.
yv nHM,Fe 0.79
mvd = M O2 + M (53) 4.3.4 Amount of direct and indirect reduction. By construc-
2000 0.21 N 2 tion, the abscissa X = 1 allows computing yd (Figure 2), which
is the oxygen removed from the iron oxides by direct reduction
4.3.3 Blast furnace gas composition. Once the coke and air (Eq.(55)).
flow rate are known, the input streams to the blast furnace
become completely characterized. In order to characterize also yd = µ + YE (55)
the output streams, we have to find 15 unknown variables.
These are the mass flow rate of each component in the blast Then, the oxygen removed by indirect reduction is easily
furnace gas (mBFG,CO, mBFG,CO2, mBFG,H2, mBFG,H2O, mBFG,N2), the calculated as Eq.(56).
mass flow rate of hot metal mHM (the individual mass flow
rate of each component in hot metal is calculated through its yi = YA − yd (56)
mass fraction, which is assumed known from the beginning),
and the mass flow rate of each component in the slag (mSlag,SiO2, 4.3.5 Flame temperature. The flame temperature is the tem-
mSlag,Al2O3, mSlag,CaO, mSlag,MgO, mSlag,MnO, mSlag,CaS, mSlag,P2O5, mSlag,Fe2O3, perature that the raceway gas reaches when all oxygen from
mSlag,FeO). The system of 15 equations to solve the balance com- hot blast and injections has been used for the incomplete
prises 11 mole balances for Fe, Si, Al, Ca, Mg, Mn, C, H, N, combustion of C to CO, and all water has been dissociated
S, P (the O balance was already accounted in the elabora- to CO and H2. It can be considered as a control parameter to
tion of the operating line), two mass balances for mSlag,Fe2O3, check if the studied configuration of blast furnace is reason-
mSlag,FeO (they correspond to the Fe2O3 and FeO content of the able, since the injection of auxiliary fuels drops the flame
coal ashes, which enter at the tuyeres and are not reduced), temperature as a consequence of their lower C/H ratio com-
one equation related to the final oxidation state of the blast pared with coke (flame temperature should be kept between
furnace gas (information coming from the operating line, 2000 and 2300°C for a proper operation). In such cases, blast
Eq.(54)), and one energy balance of the preparation zone (i.e., oxygen enrichment may be required to maintain a suitable flame
of the upper zone). temperature.
nBFG,CO + n BFG,H YA − YE From a theoretical point of view, this flame temperature can
2O
ηCO − H ≡ 2
= − 1 (54)
2 n BFG,CO + n BFG,CO + n BFG,H + n BFG,H µ be calculated from an energy balance in the raceways. We use
2 2 2O
Eq.(57), where we made a similar construction to that of
Eq.(37). The term qcyv stands for the combustion of C to CO
It should be noted that knowing the final oxidation state of the
by using the O from the hot blast, the term qvyv represents
gas exiting the top of the furnace (i.e., Eq.(54)) is not suffi-
the sensible heat that the dry hot blast is providing, while the
cient to compute the final composition of the top gas because
term qs,CyC is the sensible heat of the coke carbon used in com-
the gas is not in equilibrium at the upper part of the fur-
bustion (which comes at Tf from the lower zone). At the right
nace. In other words, the water–gas shift reaction ((Eq.(7))
side of the equation, we have (qer + qes)ye which is the con-
changes the BFG composition without modifying ηCO–H2,
tribution of the moisture of the hot blast (decomposition
but the extent of this reaction is unknown. For this reason,
the energy balance of the preparation zone is required. This ( )
and sensible heat), also ∑ j q jd + q js − 2b j qc y j for the aux-
energy balance is detailed in Appendix C25. iliary injections (decomposition, sensible heat and combustion
Page 11 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
of the O entering with this injections), and qs,N2yrg,N2 + qs,COyrg,CO Where Vv is the dry blast volume per ton of hot metal, Tv is the
+ qs,H2yrg,H2, as the energy required to heat the reducing gas temperature of the blast, ωv,O2 is the oxygen content in the
up to the temperature of the flame Tfl. Lastly, it is included blast, ηnat is the natural moisture in the blast, η is the total mois-
the heating from TR to Tfl of the ashes from coal injections ture in the blast, VNG is the natural gas consumption, mPC is
( )
∑ j Z j qs , Z y j . the pulverized coal consumption, ΩPC,C is the carbon content
in pulverized coal, ΩPC,H is the hydrogen content in pulver-
ized coal, ωNG,CH4 is the methane content in natural gas, and
qc yv + qv yv + qs ,C yC = ( qer + qes ) ye
Cp,v, CP,H2O and Cp,rg are the specific heats of blast, moisture and
+ ∑ ( q jd + q js − 2b j qc ) y j reducing gas. The units are provided in the nomenclature list.
j (57)
+ ∑ qs ,i yrg,i + ∑ z j qs ,Z y j 5. Model validation
i j
In order to validate the model, we are going to use three differ-
ent data sets. The first one comes from the original work of Rist10,
We can say that the flame temperature is a result indirectly
corresponding to a conventional air-blown blast furnace with-
derived from the operating line since for its computation we need
out auxiliary injections. The second one is elaborated from
to know the amount of air injected, the amount of coke, and the
the work of Babich et al.18, for an air-blown blast furnace with
amount and composition of the reducing gas (obtained from the
pulverized coal injection and O2-enriched air. The last one is
composition of the BFG). The variables yv, yC, yrg,N2, yrg,CO and
taken from the work of Sahu et al.27 for an oxygen blast furnace
yrg,H2 are calculated by Eq.(52), and Eq.(58) to Eq.(61). In
with top gas recycling. When solving the operating diagram,
the case of yC, we can compute it as the oxygen available for
we assume the composition of iron ore, coke, auxiliary injec-
combustion minus the C coming from the injections (because
tions and hot metal to be known. For iron ore and injections,
1 mole of O consumes 1 mole of C), or as the C coming from
the total mass flow is also known. For the hot blast, the mois-
coke minus the C used in reactions other than combustion. The
ture is fixed. Regarding energy balances, all the inlet and outlet
variable τj is equal to 1 when the auxiliary injection contains
temperatures are known. Finally, the operating conditions on
carbon and 0 when not.
chemical efficiency, heat removed by the staves, the proportion
of heat evacuated between the preparation and elaboration
yC = yv + ∑ (2b j − τ j ) y j
j zone, and the temperature of the thermal reserve zone are set
(58) to known values. Under this framework, the model allows to
103 Ù K,C m K
= − (γ + y Si + y Mn + y P + y S + yd ) compute the mass flows of coke, hot blast, hot metal and blast
M C nHM , Fe furnace gas. Moreover, the composition of the blast furnace
gas and slag can be calculated. Regarding operation condi-
yrg , N 2 = n BFG,N n (59) tions, the slope and intercept of the operating line, the amount
2 HM,Fe
of direct and indirect reduction, the flame temperature and the
yrg , CO = (n BFG,CO + n BFG,CO2 ) n HM,Fe (60) oxidation state of the BFG are obtained. The list of inputs and
outputs are summarized in Figure 4 together with a concep-
yrg ,H2 = (n BFG,H + n BFG,H O ) n HM,Fe tual scheme of the blast furnace. Exceptionally, for those cases
2 2
(61)
in which we are not introducing air in the blast furnace, we
assume the mass of the host blast to be known (mv = 0) and we
The calculation of the rest of the terms were already explained
calculate instead the mass flow of the injection that carries the
in Appendix A25 for the energy balance of the elaboration
main oxygen input (mj).
zone. All the q terms are tabulated in Table 6 according
to the adjustment to Eq.(113) or Eq.(114). The variable T in
these equations must be substituted by TR for qc and qer, by
Tv for qv and qes, by Tf for qs,C, by Tj for qjd and qjs, and by Tfl
for qs,i and qs,Z. Thus, the only missing variable in Eq.(57) is Tfl.
Page 12 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
5.1 Air-blown blast furnace without auxiliary injections 5.2 Air-blown blast furnace with pulverized coal
(Rist, Data set 1) injection and O2-enriched air (Babich, Data set 2)
This data set was elaborated from the original work of Rist10, This data set was elaborated from the work of Babich et al.18,
corresponding to an air-blown blast furnace without auxiliary corresponding to an air-blown blast furnace with pulver-
injections. The input mass flow of iron ore and its composition ized coal injection and O2-enriched air. We consider the O2-
was calculated through mass balances based on Rist’s results, enriched hot blast as an auxiliary injection. The proximate
assuming typical mass fractions for those oxides other than analysis of the coal is 74.1% ΩFC, 17.2% ΩVM, 7.6% ΩZ and
iron. Also, it was assumed a typical mass fraction distribution 1.2% ΩM18. Typical mass fractions for iron ore and coke were
for oxides in coke28. The outlet temperature of BFG was assumed. Moreover, since the operating parameters r, TR
assumed to be 180°C to estimate the gas composition. and θst are not provided in the work of Babich et al., typical
values are chosen.
As explained through the model description and appendixes25,
three potential errors were identified in the original model All results show a discrepancy below the 1% with respect
provided by Rist: (1) underestimation of the heat required to the data of Babich et al.18, except for the coke consump-
inside the blast furnace by not accounting for the heating of tion, which varies a 2.1%. The flame temperature calculated by
slag between TR and Tf (see Appendix A.1225), (2) overestima- Eq.(62) gives the same result than the one provided by
tion of the heat required because of wrongly computing the Babich, since he developed that formula. The error in the
term l with Cp,SiO2 in wrong units (see Appendix A.1325) and flame temperature calculated by our model is 1%. The oper-
(3) overestimation of the heat required because of considering ating line of this blast furnace is depicted in Figure 5, which
twice the heat absorbed during direct reduction for a number of no longer passes through the point W of the equilibrium line
moles equals to Ye (see Appendix A.525). For the sake of com- (chemical efficiency is 0.92). The direct reduction represents
parison, we intentionally introduced the same three errors a 30.6% and the indirect reduction a 69.4%.
in our model (third column in Table 3). This way, the model
reproduces the results of Rist with a discrepancy below 3.5%, 5.3 Oxygen blast furnace with pulverized coal injection
thus validating the consistency of the model. and top gas recycling (Sahu, Data set 3)
In order to validate the model also under oxygen regimes, a
third data set from the literature is used. The availability of
The results of the model are also presented with the three men-
complete data sets in literature is limited (probably because
tioned corrections included (fifth column in Table 3). In
of non-disclosure agreements) and the information gets even
this case, the variation with respect to the data provided by Rist
more scarce when it comes to oxygen blast furnaces. For this
is beyond 5% in some cases, without overpassing the 6%. For-
reason, we take the only full data set we found (Sahu et al.27),
tunately for Rist, the three errors counterbalanced, and reason-
despite it has some underlying inconsistencies (which are con-
able results could be achieved despite of them. Making the
sequently reflected in the results of our model). Neverthe-
comparison fair, the variation between the results of our model
less, it gives a good idea on the possibilities of the model when
with and without the corrections is in the range 1–2%. This
assessing oxygen blast furnaces. The data set corresponds to an
clearly shows how well the three errors counterbalanced.
oxygen blast furnace with pulverized coal injection and top
gas recycling (Figure 6). The proximate analysis of the coal
Looking for further validation, we compare the flame tempera- is 58% ΩFC, 27% ΩVM, 10% ΩZ and 5% ΩM.
ture computed through our model with the calculated by Eq.(62).
The latter gives 1898°C (see Appendix D25) and we obtained The results of the model are compared to the reference data in
1879°C with our model, what means we have a discrepancy Table 4. Most of the results have an error below the 8%, so it
of only 1.0%. In both cases, the flame temperature is below can be considered validated. The greatest discrepancy comes
2000°C, which is not suitable for blast furnaces. Neverthe- from the nitrogen content of the blast furnace (38% error)
less, the purpose of this data set is validation only, which can be because the nitrogen mass balance is not correct in the data
assumed achieved. The type of configuration provided by this set provided by Sahu et al.26. The total inlet mass of nitrogen is
data set does not correspond to a state-of-the-art blast furnace, 17.6 kg/tHM (from coal and top gas injection), while the total
since nowadays most blast furnaces use auxiliary injections outlet mass was reported as 26.7 kg/tHM (in the BFG). Assum-
like pulverized coal and burdens of greater quality that lead ing that the correct content of nitrogen in the BFG should
to lower slag productions. be 17.6 kg/tHM, the relative error of the result of the model
drops to a reasonable 6%. Other relevant inconsistencies are
The operating line of this blast furnace is depicted in the incorrect mass balances of Fe, Mn and P. The former
Figure 5, showing that 23% of the oxygen bonded to iron leads to an overestimation of the hot metal produced, and the
oxides is removed by direct reduction, and the remaining 77% latter to negative mass flows of MnO and P2O5 in the slag (the
through indirect reduction. The chemical efficiency was assumed inlet mass flows of MnO and P2O5 are not enough to reach the
1, so the operating line passes through the point W of the specified Mn and P content in the hot metal). Despite these
equilibrium line. issues, the model can reproduce the behavior of the oxygen
Page 13 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
Table 3. Comparison of model results with literature data sets for air-blown blast furnaces, for validation
purposes. *This column shows the results of the model when corrections are not included (see text).
Input streams (kg/tHM) Rist10 Model* Δ (%) Model Δ (%) Babich18 Model Δ (%)
Coke (25°C) 598.9 612.8 2.3 619.9 3.5 283.0 289.0 2.1
# SiO2 26.5 27.1 2.3 27.4 3.5 19.7 20.1 2.1
# Al2O3 10.8 11.0 2.3 11.2 3.5 9.1 9.3 2.1
# CaO 36.1 37.0 2.3 37.4 3.5 0.8 0.8 2.1
# MgO 5.4 5.5 2.2 5.5 3.6 0.3 0.3 2.1
Hot blast (700°C Rist / 1200°C 2125.6 2198.7 3.4 2242.5 5.5 1077.8 1086.6 0.8
Babich)
Page 14 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
Hot metal (1400°C Rist / 1500°C 1000.0 1000.0 0.0 1000.0 0.0 1000.0 1000.0 0.0
Babich)
Slag (1450°C Rist / 1550°C Babich) 1000.0 1002.8 0.3 1003.8 0.4 260.0 261.8 0.7
BFG (180°C Rist / 150°C Babich) 3000.0 3085.2 2.8 3135.1 4.5 1953.2 1966.9 0.7
Operating conditions
Oxidation state of the BFG, ηCO–H2 (-) - 0.411 - 0.410 - 0.486 0.487 0.2
blast furnace described by Sahu et al.27 within reasonable dis- one reported by Sahu et al.27 (2126°C). Although they did
crepancy limits. Moreover, the model has turn out to be a useful not explain how they calculated this temperature, they prob-
tool to identify potential inconsistencies in data sets. ably did not account for the heating of the injected top gas, while
our model does. In fact, one of the reasons of recirculating
Regarding the flame temperature, in this case we cannot use top gas is to substitute N2 as heat sink during oxygen regimes,
the Eq.(62) to compute it since this formula does not account so probably that is why Sahu et al.26 had to set the inlet tem-
for top gas injections. By using our model, we found a flame perature of the oxygen to 25°C (otherwise they would found
temperature of about 1770°C, which is a 17% lower than the excessive flame temperatures because of not accounting for
Page 15 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
Figure 5. Operating lines obtained from the Data sets 1, 2 and 3 (Table 3 and Table 4).
Figure 6. Conceptual schemes summarizing the mass flows and relative errors obtained during the validation of the model for
Data set 1, 2 and 3 (the complete data sets are presented in Table 3 and Table 4).
the recycled gas). If we calculate in our model the temperature 6. Predicting the operating line for an oxygen blast
of the flame without taking into account the sensible heat of furnace integrated with PtG
the recycled gas, we found temperatures between 2200 and Once the model has been validated, we are going to analyze
2300°C, which are closer to that provided by Sahu et al.26. a potential configuration for a blast furnace operation under
an oxygen regime, with top gas recycling, and injections of
The operating line of this blast furnace is depicted in Figure 5. pulverized coal and synthetic natural gas. The data set was
It passes through the point W of the equilibrium line because inspired in the works of Sahu et al.27 and Jin et al.29, aiming to
it is assumed chemical ideality. It can be seen that that the reproduce similar results to them, but keeping reasonable
direct reduction represents only a 0.3% of the total reduc- operating conditions (e.g., flame temperature between 2000–
tion, which is unrealistic for a real operation. Therefore, for 2300°C, direct reduction about 10–15%, slag production above
this reason and the other mentioned inconsistencies, we do not 200 kg/tHM). The proximate analysis of the coal for this data
recommend using this data set for further analyses. set is 70.7% ΩFC, 17.2% ΩVM, 10.8% ΩZ and 1.2% ΩM.
Page 16 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
Table 4. Data sets for oxygen blast furnaces with top gas recycling. The Data set 3 is taken from Sahu et al.27
and compared with the results of our model. The Data sets 4, 5 and 6 were elaborated in the present paper.
Page 18 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
The proposed oxygen blast furnace would be integrated with to recycle the top gas, which can benefit from the exothermal
a PtG plant, which renewably produces the synthetic methane. heat available from the methanation process. Furthermore,
This methane is used to reduce the consumption of coke or the electrolyzer by-produces 87 kg/tHM of O2 that can be used
coal, thus substituting a fossil fuel by a renewable fuel. The in the blast furnace, thus reducing the energy requirements
conceptual scheme of the blast furnace is the one depicted in in the air separation unit that enriches the hot blast. For
Figure 7 (Data set 5 when coke is replaced and Data set 6 these calculations, it was assumed an electrolysis efficiency
when coal is replaced). We assume a 150 MW electrolyzer and of 68%30.
a 280 tHM/h blast furnace. This means that the electrolyzer
produces 11 kg/tHM of hydrogen, which are converted to For comparison purposes, the blast furnace was also modelled
21.8 kg/tHM of synthetic methane by consuming 60 kg/tHM of without the integration of the PtG plant, and therefore with-
CO2. This carbon is continuously recycled in a closed loop, out the injection of synthetic methane (Data set 4, which is
and therefore the corresponding emissions are avoided. The considered as the base case). In all cases, the mass flow of
CO2 would come from the carbon capture stage that is used injected top gas was kept constant at 600 kg/tHM (despite its
Figure 7. Conceptual schemes summarizing the mass flows obtained for Data sets 4, 5 and 6 (see Table 4).
Page 19 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
composition changes). The three data sets are presented in These values of energy penalization for the avoidance of CO2
Table 4. In the base case (Data set 4), the total CO2 that would are clearly above the typical consumptions of conventional
be emitted after the combustion of the BFG is 1154 kgCO2/tHM, amine carbon capture, which are usually in the range
of which the 41.4% (i.e., 478 kgCO2/tHM) could be directly cap- 3.7–4.8 MJ/kgCO233,34. Nevertheless, the integration of PtG
tured since they come from the carbon capture stage. When presents the additional benefits of diminishing the coke/coal
implementing PtG, 60 kg/tHM of CO2 are recycled in closed loop, consumption in the blast furnace, reducing the electricity
which allow to reduce the coke consumption by 8.7% or alter- consumption in the air separation unit, and suppressing the
natively the coal consumption by 16.3%. In case of substituting requirement of geological storage for the avoided CO2 (it is
coke, the total emissions become 1083 kgCO2/tHM, of which the kept in a closed carbon loop thanks to PtG). Therefore, an over-
31.8% can be directly captured and stored (i.e., 345 kgCO2/tHM). all economic and energy analysis of the whole integrated steel
Thus, in overall terms, the CO2 is diminished by 71 kgCO2/tHM plant should be necessary to reach farther conclusions. Fur-
with respect to the base case by using 1.93 GJ/tHM of electric- thermore, there exist other potential PtG integrations that may
ity in the electrolyzer, so the energy penalty of the CO2 avoid- lead to lower energy penalizations, such as the utilization of
ance is 27.1 MJ/kgCO2. This value is line with the results of the BFG in the methanation stage rather than the captured CO2
Perpiñán et al.8 for similar PtG integrations in the steel indus- (which allow producing more SNG with the same amount
try (he reported 34 MJ/kgCO2). Actually, since we have cut the of H2) or the injection of the H2 in the blast furnace without
coke consumption by 22.8 kg/tHM, which, in terms of electric- considering a methanation stage.
ity, is equivalent to a reduction of 0.24 GJ/tHM (assuming a
coke heating value of 27.3 MJ/kg and a subcritical power plant Regarding the flame temperature, it is reduced about 90°C when
net efficiency of 38%31), and additionally we diminished the injecting 22 kg/tHM of synthetic methane for replacing coke.
O2 that has to be produced in the air separation unit by This is a reduction of 4.0°C per kgSNG/tHM, which is in good
77 kgO2/tHM (we need inject 10 kgO2/tHM more than in the base but agreement with the value reported by Babich et al. for natural
we have 87 kgO2/tHM available from electrolysis), which translates gas (4.5°C per kgNG/tHM16). When we replace coal by syn-
into a reduction of 0.05 MJ/tHM of electric consumption (assum- thetic natural gas, the temperature of the flame is reduced by
ing a typical ASU specific consumption of 0.61 MJ/kgO232), 3.1°C per kgSNG/tHM, which also matches the value reported by
the net energy penalization of the CO2 avoidance becomes Babich et al. for this type of replacement (3.0°C per kgNG/tHM16).
23.1 MJ/kgCO2. Alternatively, in the case of using the syn-
thetic methane to replace part of the coal, the total emissions Lastly, the percentage of direct reduction is 15% for the base
become 1068 kgCO2/tHM. Following the same calculations just case and about 12% when injecting synthetic natural gas. The
mentioned, the gross energy penalization of CO2 avoidance in decrease in this value was expected since direct reduction takes
this case would be 22.4 MJ/kgCO2 and the net energy penaliza- place through solid carbon, and the coke/coal input is dimin-
tion 17.9 MJ/kgCO2 (the heating value of coal is 29.4 MJ/kg, ished when PtG is integrated. The operating lines are presented
Eq.(130)). Table 5 summarizes these calculations. in Figure 8.
Table 5. Comparison of CO2 emissions, CO2 avoided and energy penalization between data sets 4
(oxygen blast furnace), 5 (oxygen blast furnace with PtG integration for the replacement of coke
by synthetic methane) and 6 (oxygen blast furnace with PtG integration for the replacement of
coal by synthetic methane).
Page 20 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
Figure 8. Operating lines of oxygen blast furnaces, obtained from Data sets 4, 5 and 6 (Table 4).
Page 21 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
performance of the blast furnace when coke is replaced to the CΔT sensible heat of the burden between TR – ΔTR and TR
R
performance when coal is replaced. It was assumed a 280 tHM/h (lack of thermal ideality), kcal/molFe
oxygen blast furnace producing 1154 kgCO2/tHM, coupled to a dj number of moles of S2 in injectant j per number of
PtG plant with 150 MW electrolysis power capacity that pro- moles of injectant j, molS2/molj
duces 11 kgH2/tHM of hydrogen (then converted to 21.8 kgSNG/tHM
e number of moles of H2O per number of moles of O
of synthetic methane by using pure CO2 from a carbon cap-
in the air (i.e., e = ye/yv), molH2O/molO
ture stage). The model shows fuel replacement ratios of
1.05 kg/kgSNG for coke and 1.38 kg/kgSNG for coal, which lead f sensible heat of the hot metal between TR and Tf
to CO2 emission reductions of 6.1% and 7.4%, respectively. (outlet temperature), kcal/molFe
As the electricity consumed in the PtG plant is 1.93 GJ/tHM, the gi specific Gibbs free energy of compound i, kcal/moli
gross energy penalization of the CO2 avoidance is 27.1 MJ/kgCO2 ΔG Gibbs free energy change of reaction, kcal/mol
when coke is replaced and 22.4 MJ/kgCO2 when coal is replaced.
Considering the energy content of the fossil fuel that is saved, h enthalpy, kcal/mol
and the electricity that is no longer consumed in the air Δhf enthalpy of fusion, kcal/mol
separation unit, the net energy penalizations are 23.1 MJ/kgCO2 Δch enthalpy of combustion, kcal/mol
and 17.9 MJ/kgCO2, respectively. These values are several times
Δfh enthalpy of formation, kcal/mol
greater than the specific consumption of amine carbon capture
(typically 3.7–4.8 MJ/kgCO2). However, the integration of Keq equilibrium constant of the water-gas shift reaction, -
PtG diminishes the coke/coal consumption in the blast fur- l sensible heat of the slag between TR and Tl (outlet
nace, reduces the electricity consumption in the air separation temperature), kcal/molFe
unit, and eliminates the requirement of geological storage for mj,i mass of compound i in stream j per ton of hot metal,
the avoided CO2. Therefore, a detailed economic comparison kg/tHM
between both methods should be necessary to reach firm
M molar weight, kg/kmol
conclusions. Furthermore, there exist other potential PtG inte-
grations that may lead to lower energy penalizations, such nj,i number of moles of compound i in stream j per ton
as the utilization of the BFG in the methanation stage rather of hot metal, mol/tHM
than captured CO2, what would allow producing more SNG p heat removed by the staves in the elaboration zone,
with the same amount of H2. kcal/molFe
p′ heat removed by the staves in the preparation zone,
Lastly, the present paper has provided six full data sets of dif- kcal/molFe
ferent blast furnaces operations, specifying the composition of qc heat released at TR by the reaction C (coke) + 0.5 O2 →
all streams, as well as the most relevant operating parameters CO, kcal/molO
(e.g., temperature of the thermal reserve zone, heat evacuated
by the staves, and the temperature of the flame). The availabil- qe heat required by the H2O in hot blast due to disso-
ity of this information in literature is very scarce, especially for ciation, reverse water-gas shift and sensible heat,
oxygen blast furnaces. kcal/molH2O
qer heat absorbed at TR by the reaction C + H2O → CO + H2,
Nomenclature kcal/molH2O
In this paper we used calorie as unit of energy for calculations qes enthalpy change of water between Tv and TR, kcal/molH2O
related to operating diagram in order to facilitate comparing qf,i heat required to melt and increase the temperature
results with the original work of Rist. Elsewhere, the SI unit is of compound i from TR to Tf, kcal/moli
used (Joule).
qFe sensible heat of wüstite between TR – ΔTR and TR,
0.95O
kcal/molO
Symbols
ai calculation parameter, 1/°Ci qFe C heat absorbed at TR by the reaction 3Fe + C → Fe3C,
3
kcal/molC
aj number of moles of H2 in injectant j per number of
moles of injectant j, molH2/molj qFe O sensible heat of magnetite between TR – ΔTR and TR,
3 4
kcal/molFe3O4
A calculation parameter, kcal/mol
qg heat absorbed at TR by the reaction C + CO2 → 2CO,
bj number of moles of O2 in injectant j per number of kcal/molC
moles of injectant j, molO2/molj
qim heat released at TR by the reaction 1/4 Fe3O4 + CO → 3/4
B calculation parameter, kcal/mol Fe + CO2, kcal/molO
cj number of moles of N2 in injectant j per number of qiw heat released at TR by the reaction Fe0.95O + CO → 0.95
moles of injectant j, molN2/molj Fe + CO2, kcal/molO
cp heat capacity, kcal/(mol·K) qj thermal demand by injectant j, CH2aO2bN2cS2dZz (or
C calculation parameter, kcal/mol H2aO2bN2c or O2bN2c) due to dissociation, sensible heat,
Page 22 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
reverse water-gas shift of the H2, incomplete com- T temperature, °C (unless otherwise specified)
bustion with the O2, and transfer of S to the slag, ΔTR difference of temperature between the gas and the
kcal/molj solid at the beginning of the middle zone, °C
qjd heat absorbed at Tj by the dissociation reaction of Vj volume flow of stream j consumed in the blast
injectant j, CH2aO2bN2cS2dZz → C + aH2 + bO2 + cN2 + furnace per ton of hot metal, Nm3/tHM
dS2 + zZ, kcal/molj
xd number of moles of O removed from wüstite by direct
qjs,i enthalpy change of element i from injectant j between reduction per total moles of reducing gas mixture,
Tj and TR, kcal/moli molO/(molC+molH2)
qk heat absorbed at TR by the total H2 in the furnace when xe number of moles of H2O in hot blast per total moles
considering that it is completely converted to H2O of reducing gas mixture, molH2O/(molC+molH2)
through the reverse water-gas shift reaction, kcal/molH2 xh molar fraction of hydrogen and water in the reducing
ql,i heat required to melt and increase the temperature of gas mixture, -
compound i (SiO2, Al2O3, CaO, MgO) from TR to Tl, xi number of moles of O transferred from the iron oxides
kcal/moli to the gas by indirect reduction per total moles of
qMn heat absorbed at TR by the reaction C +MnO → CO + reducing gas mixture, molO/(molC+molH2)
Mn, kcal/molO xj number of moles of injectant j (overall formula
qP heat absorbed at TR by the reaction 1/5 P2O5 ⋅ 3CaO + CH2aO2bN2cS2dZz, H2aO2bN2c or O2bN2c) per total moles
C + 6/5 Fe → 2/5 Fe3P + 3/5 CaO + CO, kcal/molO of reducing gas mixture, molj/(molC+molH2)
qrmc heat consumed by the reaction 1.5 Fe2O3 + 0.5 CO → xk number of moles of H2 in the coke per total moles
Fe3O4 + 0.5 CO2, including also the heat exchange of reducing gas mixture, molH2/(molC+molH2)
between reactants and products at different temperatures, xMn number of moles of O removed by direct reduc-
kcal/molFe3O4 tion of MnO per total moles of reducing gas mixture,
qrwc heat consumed by the reaction 0.5 Fe2O3 + 0.5 CO molO/(molC+molH2)
→ FeO + 0.5 CO2, including also the heat exchange xP number of moles of O removed by direct reduc-
between reactants and products at different temperatures, tion of P2O5 per total moles of reducing gas mixture,
kcal/molFeO molO/(molC+molH2)
qR,i sensible heat of compound i (SiO2, Al2O3, CaO, MgO, C) xS number of moles of O replaced by S in the slag per
between TR – ΔTR and TR, kcal/moli total moles of reducing gas mixture, molO/(molC+molH2)
qS heat absorbed at TR by the reaction 1/2 S2 + CaO + C → xSi number of moles of O removed by direct reduc-
CaS + CO, kcal/molS tion of SiO2 per total moles of reducing gas mixture,
qs,i sensible heat of component i between Tfl and TR, molO/(molC+molH2)
kcal/moli xv number of moles of O in hot blast per total moles
qs,nr,i sensible heat of non-reacting compound i, kcal/moli of reducing gas mixture, molO/(molC+molH2)
qSi heat absorbed at TR by the reaction C + 1/2 SiO2 + 3/2 X abscissa in the Rist diagram, (molO+molH2)/
Fe → CO + 1/2 Fe3Si, kcal/molO (molC+molH2)
qst heat removed by the staves from the blast furnace per yC number of moles of C (coke) that are consumed in
ton of hot metal, kcal/thot metal combustion in the raceways per mol of Fe produced,
molC/molFe
qv energy available from the sensible heat of the air
yd number of moles of O removed from wüstite by
between Tv and TR, kcal/molO
direct reduction per mol of Fe produced, molO/molFe
qy heat absorbed at TR by the carburization of the iron,
ye number of moles of H2O in hot blast per mol of Fe
kcal/molC
produced, molH2O/molFe
qyFe heat absorbed at TR by the reaction C(coke) →
yi number of moles of O transferred from the iron
C(austenite), kcal/molC
oxides to the gas by indirect reduction per mol of Fe
qε heat absorbed at TR by the reaction H2 + CO2 → H2O + produced, molO/molFe
CO, kcal/molH2
yj number of moles of injectant j (overall formula
q ′ε heat absorbed by the reaction H2 + CO2 → H2O + CO CH2aO2bN2cS2dZz, H2aO2bN2c or O2bN2c) per mol of Fe
taking into account the temperature of reactants and produced, molj/molFe
products, kcal/molH2(reacting)
yk number of moles of H2 in the coke per mol of Fe
r chemical efficiency of the blast furnace, - produced, molH2/molFe
R ideal gas constant, kcal/(mol K) yMn number of moles of O removed by direct reduction
Si specific entropy of compound i, kcal/(moli K) of MnO per mol of Fe produced, molO/molFe
Page 23 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
ynr,i number of moles of compound i that is traversing μ slope of the Rist diagram, i.e. , number of moles of reduc-
the preparation zone without reacting per mol of Fe ing gas per mol of Fe produced, (molC+molH2)/molFe
produced, moli/molFe τj calculation parameter that is 1 when the auxiliary
yP number of moles of O removed by direct reduction injection j contains carbon and 0 when not, -
of P2O5 per mol of Fe produced, molO/molFe ϕ percentage of H2 consumed in the preparation zone, -
yrg,i number of moles of component i in the reducing gas ω molar fraction, -
at the raceways per mol of Fe produced, moli/molFe
Ωji mass fraction of compound i in stream j, -
yrmc number of moles of Fe3O4 produced in the preparation
zone per mol of Fe produced, molFe3O4/molFe Subscripts and superscripts
yrwc number of moles of FeO produced in the preparation A initial oxidation state of the iron oxides at the inlet
zone per mol of Fe produced, molFeO/molFe of the blast furnace
yS number of moles of O replaced by S in the slag per BFG blast furnace gas
mol of Fe produced, molO/molFe C related to carbon or C-CO2 mixtures
ySi number of moles of O removed by direct reduction d dry basis (superscript), decomposition (subscript)
of SiO2 per mol of Fe produced, molO/molFe
e moisture in hot blast
yv number of moles of O in hot blast per mol of Fe
f hot metal (at the outlet)
produced, molO/molFe
FC fixed carbon
ywgs number of moles of H2 reacting in the preparation zone,
molH2/molFe Fe3C referred to the number of moles of C in hot metal
dissolved as cementite
Y ordinate in the Rist diagram, (molO+molH2)/molFe
fl flame
YE intercept of the operating line representing the moles
of H2 and O coming from sources other than iron H related to hydrogen or H2-H2O mixtures
oxides that contribute to the formation of the reduc- HM hot metal
ing gas per mol of Fe produced (negative sign IN inlet of solids at the top of the furnace
by convention), (molO+molH2)/molFe
IO iron ore
YE* terms of YE that are independent of yv, (molO+molH2)/
j injectant (overall formula CH2aj O2bj N2cj S2dj ZZJ or
molFe
stream
YM* number of moles of O per mole of Fe3O4 (i.e., YM* ),
K coke
molO/molFe3O4
slag (at the outlet)
Zj number of moles of ashes in injectant j per number
of moles of injectant j, molash/molj l magnetite or moisture
m natural moisture in blast
Greek symbols NG natural gas
γ number of moles of C dissolved in the hot metal,
nr non-reacting in the preparation zone
molC/molFe
P characteristic point of the operating line referring to
γFe C number of moles of C dissolved in the hot metal
3 the energy balance of the blast furnace
as cementite, molC/molFe
PM primary volatile matter
γγFe number of moles of C dissolved in the hot metal
as austenite, molC/molFe R characteristic point of the operating line referring to
the thermal reserve zone
δ decrease in the available heat due to the presence of
magnetite in the elaboration zone, kcal/molFe SM secondary volatile matter
δ ′ decrease in the heat absorbed by the reverse v hot blast / air
water-gas shift reaction because of the lack of chemical VM volatile matter
ideality, kcal/molFe W referring to wüstite or to the characteristic point
Δ1 calculation parameter, - of the equilibrium line in which pure wüstite is in
Δ2 calculation parameter, (molO+molH2)/molFe equilibrium with the gas
𝜖K number of moles of C and H2 in coke per kg of P ashes
coke, (molC+molH2)/molK γFe referred to the number of moles of C in hot metal that
η humidity in the air, (g/Nm )
3 are dissolved as austenite
ηCO–H oxidation state of the blast furnace gas leaving the Data availability
2
top, - Extended data
θst fraction of the total heat removed by the staves that Zenodo: Extended data for ‘Revisiting the Rist diagram for pre-
is coming from the elaboration zone (i.e., from the dicting operating conditions in blast furnaces with multiple
control volume of the energy balance), - injections’ https://doi.org/10.5281/zenodo.563758039
Page 24 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
The project contains the following extended data: Appendix D – Calculation of the flame temperature through
Appendix A – Terms of the energy balance in the elaboration Eq.(62)
zone DISIPO_OPEN DATA_1.xlsx
Appendix B – Heat of decomposition of coal (dry basis)
Data are available under the terms of the Creative Commons
Appendix C – Energy balance in the preparation zone Attribution 4.0 International license (CC-BY 4.0).
References
1. Mac Dowell N, Fennell PS, Shah N, et al.: The role of CO2 capture and
utilization in mitigating climate change. Nat Clim Chang. 2017; 7: 243–249. ironmaking by modeling with aspen plus. ISIJ Int. 2015; 55(2): 340–347.
Publisher Full Text Publisher Full Text
2. Ueckerdt F, Bauer C, Dirnaichner A, et al.: Potential and risks of hydrogen- 18. Babich A, Senk D, Solar J, et al.: Efficiency of biomass use for blast furnace
based e-fuels in climate change mitigation. Nat Clim Chang. 2021; 11: injection. ISIJ Int. 2019; 59(12): 2212–2219.
384–393. Publisher Full Text
Publisher Full Text 19. Prego AF, Sanchez AC, Babich A, et al.: Operación de hornos altos con
3. Bailera M, Lisbona P, Romeo LM, et al.: Power to Gas projects review: Lab inyección de carbón pulverizado en diferentes condiciones tecnológicas.
pilot and demo plants for storing renewable energy and CO2. Renew. Sustain Rev Metal. 2001; 37(3): 423–436.
Energy Rev. 2017; 69: 292–312. Publisher Full Text
Publisher Full Text 20. Biswas AK: Principles of blast furnace ironmaking. Theory and practice.
4. Bailera M, Lisbona P, Peña B, et al.: A review on CO2 mitigation in the Iron and Cootha Publishing house, 1981.
Steel industry through Power to X processes. J CO2 Util. 2021; 46: 101456. Reference Source
Publisher Full Text 21. Babich KTMA, Senk D, Gudenau HW: IRONMAKING. 2008; 402.
5. Hisashige S, Nakagaki T, Yamamoto T: CO2 emission reduction and exergy Reference Source
analysis of smart steelmaking system adaptive for flexible operating 22. Geerdes M, Chaigneau R, Lingiardi O, et al.: Modern Blast Furnace
conditions. ISIJ Int. 2019; 59(4): 598–606. Ironmaking. An Introduction, IOS Press, 2020.
Publisher Full Text Reference Source
6. Rosenfeld DC, Böhm H, Lindorfer J, et al.: Scenario analysis of implementing a 23. Spanlang A, Wukovits W, Weiss B: Development of a Blast Furnace Model
power-to-gas and biomass gasification system in an integrated steel plant: with Thermodynamic Process Depiction by Means of the Rist Operating
A techno-economic and environmental study. Renew Energy. 2020; 147(Part Diagram. BHM Berg- Und Hüttenmännische Monatshefte. 2020; 165: 243–247.
1): 1511–1524. Publisher Full Text
Publisher Full Text
24. Wagner D, Devisme O, Patisson F, et al.: A laboratory study of the reduction
7. Cavaliere P: Clean Ironmaking and Steelmaking Processes. Efficient of iron oxides by hydrogen. Sohn Int Symp. 2006; 2: 111–120.
Technologies for Greenhouse Emissions Abatement. 2020. Reference Source
Publisher Full Text
25. Bailera M, Nakagaki T, Kataoka R: Extended data for ‘Revisiting the Rist
8. Perpiñán J, Bailera M, Romeo LM, et al.: CO2 Recycling in the Iron and Steel diagram for predicting operating conditions in blast furnaces with
Industry via Power-to-Gas and Oxy-Fuel Combustion. Energies. 2021; 14(21): multiple injections,’. Zenodo. 2021.
7090. http://www.doi.org/10.5281/zenodo.5637580
Publisher Full Text
26. Babich A: Methodical instruction for the flame temperature calculation.
9. Rist A, Meysson N: A dual graphic representation of the blast-furnace mass Donetsk, 1988; 24.
and heat balances. JOM. 1967; 19: 50–59.
27. Sahu RK, Roy SK, Sen PK: Applicability of Top Gas Recycle Blast Furnace with
Publisher Full Text
Downstream Integration and Sequestration in an Integrated Steel Plant.
10. Rist A, Meysson N: Étude graphique de la marche du haut fourneau avec Steel Res Int. 2015; 86(5): 502–516.
vent humide et injections aux tuyères. Rev Metall. 1965; 62(11): 995–1040. Publisher Full Text
Publisher Full Text
28. Kasahara S, Inagaki Y, Ogawa M: Flow sheet model evaluation of nuclear
11. Rist A, Meysson N: Recherche graphique de la mise au mille minimale du hydrogen steelmaking processes with VHTR-IS (very high temperature
haut fourneau a faible température de vent. Rev Métallurgie. 1964; 61(2): reactor and iodine-sulfur process). ISIJ Int. 2012; 52(8): 1409–1419.
121–146. Publisher Full Text
Publisher Full Text
29. Jin P, Jiang Z, Bao C, et al.: Mathematical Modeling of the Energy
12. Meysson N, Weber J, Rist A: Représentation graphique de la décomposition Consumption and Carbon Emission for the Oxygen Blast Furnace with Top
des carbonates dans le haut fourneau. Rev Métallurgie. 1964; 61(7–8): Gas Recycling. Steel Res Int. 2016; 87(3): 320–329.
623–633. Publisher Full Text
Publisher Full Text 30. Development of Systems using Renewable Energy-derived Hydrogen
13. Kundrat DM, Miwa T, Rist A: Injections in the Iron Blast Furnace: A Graphics Development of Power-to-gas System to Synthesize Methane from
Study by Means of the Rist Operating Diagram. Metall Trans B. 1991; 22B: Renewable Hydrogen and Exhaust CO2 for Supplying via Conventiontional
363–383. Gas Grid. 2019.
Publisher Full Text Reference Source
14. Tang J, Chu M, Li F, et al.: Mathematical simulation and life cycle assessment 31. Rocha DHD, Siqueira DS, Silva RJ: Exergoenvironmental analysis for
of blast furnace operation with hydrogen injection under constant evaluating coal-fired power plants technologies. Energy. 2021; 233: 121169.
pulverized coal injection. J Clean Prod. 2021; 278: 123191. Publisher Full Text
Publisher Full Text 32. Lisbona P, Bailera M, Peña B, et al.: Integration of CO2 capture and
15. Zhan WL, Wu K, He ZJ, et al.: Estimation of Energy Consumption in COREX conversion. Adv Carbon Capture. In: M.R. Rahimpour, M. Farsi, M.A. Makarem
Process Using a Modified Rist Operating Diagram. J Iron Steel Res Int. 2015; (Eds.), Woodhead Publishing, 2020; 503–522.
22: 1078–1084. Publisher Full Text
Publisher Full Text 33. Gervasi J, Dubois L, Thomas D: Simulation of the post-combustion CO2
16. Babich A, Senk D, Gudenau H, et al.: Ironmaking Textbook. RWTH Aachen capture with Aspen HysysTM software: Study of different configurations
University, 2008. of an absorption-regeneration process for the application to cement flue
Reference Source gases. Energy Proced. 2014; 63: 1018–1028.
17. Suzuki K, Hayashi K, Kuribara K, et al.: Quantitative evaluation of CO2 Publisher Full Text
emission reduction of active carbon recycling energy system for 34. Morales-Mora MA, Pretelìn-Vergara CF, Martínez-Delgadillo SA, et al.:
Page 25 of 30
Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
Environmental assessment of a combined heat and power plant Christian Senckenberg., Sendai, 1935; 665–795.
configuration proposal with post-combustion CO2 capture for the Mexican Reference Source
oil and gas industry. Clean Technol Environ Policy. 2019; 21: 213–226. 37. Perry RH, Green DW: Perry’s Chemical Engineers’ Handbook. 7th ed.,
Publisher Full Text McGraw-Hill Education, 1997.
35. Mills KC, Karagadde S, Lee PD, et al.: Calculation of Physical Properties for Reference Source
Use in Models of Continuous Casting Process-Part 2: Steels. ISIJ Int. 2016; 38. Green DW, Southard MZ: Perry’s chemical engineers’ handbook, 9th ed.
56(2): 274–281. McGraw-Hill Education, 2019.
Publisher Full Text Reference Source
36. Saburo U: On the Specific Heat of Iron-Carbon System at High 39. Leśniak B, Słupik ł, Jakubina G: The determination of the specific heat
Temperatures, and the Heat Changes Accompanying Those of Phase. In: capacity of coal based on literature data. Chemik. 2013; 67: 560–571.
335th Rep. Res. Inst. Iron, Steel Other Met., Universitätsbibliothek Johann Reference Source
Page 26 of 30
Open Research Europe Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
Version 1
https://doi.org/10.21956/openreseurope.15400.r28319
© 2022 Bargiacchi E. This is an open access peer review report distributed under the terms of the Creative
Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
Eleonora Bargiacchi
IMDEA Energy Institute, Madrid, Spain
The paper investigates the integration of Power to Gas (PtG) technology into oxygen blast
furnaces (OBF). For this purpose, a corrected version of the Rist diagram is proposed to predict the
system operating conditions. The major novelties of the paper rely on the revisitation of the blast
operating diagram and on the integration of OBF and PtG technologies.
The paper is novel, well organized, clear, and scientifically robust. I hereby suggest some
modifications and corrections that the authors may wish to implement in future versions of the
paper:
○ In the introduction, a better contextualization of the advantages and disadvantages of OBF
with respect to other technologies would be useful to the reader. Similarly, the proposed
solution (OBF+PtG) should be better contextualized (advantages/disadvantages with respect
to other solutions). Overall, some useful references proposing advances in steelmaking
plants with respect to carbon capture and emission reduction should be added (e.g.: Arasto
et al.)
○ "Within this framework, several authors have studied the application of power to gas (PtG)
to ironmaking processes based on the reduction of iron ores with coke in a blast furnace
(BF)." The references to these authors/papers are missing.
○ After the integration of the blast furnace with the PtG plant the author states: "Thus, in
overall terms, the CO2 is diminished by 71 kgCO2/tHM with respect to the base case by using
1.93 GJ/tHM of electricity in the electrolyzer, so the energy penalty of the CO2 avoidance is
27.1 MJ/kgCO2. ". Are the impacts associated to the electricity supplied to the electrolyzer
accounted for? If not, I recommend including a discussion on the impacts of the electricity
supplied to the electrolysis process. Also, the avoided CO2 emissions and CO2 emissions
avoidance ratio should be recalculated by adding the impacts of the electrolysis (and table 5
should be modified accordingly).
Minor comment:
○ I recommend the use of the simple past instead of the present perfect (e.g.: in the
Page 27 of 30
Open Research Europe Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
"Conclusions" section).
References
1. Arasto A, Tsupari E, Kärki J, Lilja J, et al.: Oxygen blast furnace with CO 2 capture and storage at
an integrated steel mill—Part I: Technical concept analysis. International Journal of Greenhouse Gas
Control. 2014; 30: 140-147 Publisher Full Text
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and does the work have academic merit?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Yes
Are all the source data underlying the results available to ensure full reproducibility?
Yes
I confirm that I have read this submission and believe that I have an appropriate level of
expertise to confirm that it is of an acceptable scientific standard.
https://doi.org/10.21956/openreseurope.15400.r28384
© 2022 Schenk J. This is an open access peer review report distributed under the terms of the Creative Commons
Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the
original work is properly cited.
Johannes Schenk
Montanuniversität Leoben, Leoben, Austria
Generations of blast furnace operators use the Rist diagram to evaluate the operating condition of
their furnace. With this publication, the authors have shown that this diagram can also be used to
Page 28 of 30
Open Research Europe Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
evaluate new alternative process variants for the blast furnace. They have extended the original
model approach of Rist to make the application more universal. In particular, the following
extensions were made: (1) nearly any auxiliary injection composition and material is possible, (2) a
new model approach for the energy balance that takes into account all possible energy sinks and
sources, (3) for the operating line, the chemical efficiency is taken into account for the position
relative to the Wüstit point W and in the determination of the energy balance point P.
With this publication, the authors have made a significant contribution to the future application of
the Rist diagram for the blast furnace process. The published model approach is universal. A
valuable contribution is the detailed description of submodels in appendixes A to D. In the case
studies presented, new blast furnace process and integration concepts are investigated and
evaluated. These results are scientifically and technologically new. I appreciate reviewing the
paper.
I have some comments and questions that the authors might consider in the final version:
○ Equation (49) - The variable тj is not clear. Is the “+” operator after т_j correct, or should it be
“x”?
○ Final degree of oxidation of BFG, xi - Is the LOI and moisture of burden and coke
considered?
○ Calculation of flame temperature: Which temperature is assumed for the coke entering the
raceway?
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and does the work have academic merit?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Yes
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Reviewer Expertise: My area of expertise is ironmaking. I am dealing with blast furnace, direct
Page 29 of 30
Open Research Europe Open Research Europe 2021, 1:141 Last updated: 02 MAR 2023
I confirm that I have read this submission and believe that I have an appropriate level of
expertise to confirm that it is of an acceptable scientific standard.
Page 30 of 30