Isothermal and Non-Isothermal Reduction Behaviors of Iron Ore Compacts in Pure Hydrogen Atmosphere and Kinetic Analysis

This study examines the isothermal and non-isothermal reduction behaviors of iron ore compacts in a pure hydrogen atmosphere and compares the results obtained during the reduction process by CO. The different phases accompanying the reduction reactions were identified using X-ray diffraction (XRD) and its morphology was microscopically examined. In isothermal experiments, temperature plays a significant role in the reduction process. At any given temperature, the reduction rate during the initial stages is higher than that during the final stages. The reduction rate in H2 atmosphere was faster than in CO gas. The comparison of activation energy values suggested that reduction with H2 is more efficient than with CO. At the same temperature, the time required to achieve a certain degree of reduction was lower when using H2 gas than CO atmosphere. In non-isothermal tests, the heating rate has a significant effect on the reduction rate and reduction extent. At the same heating rate, the degree of reduction was higher in H2 atmosphere than in CO gas. Based on experimental data, the parameters of reaction kinetics were deduced by application of model-free and model-fitting methods. The reduction in H2 atmosphere was controlled by nucleation model (Avrami-Erofeev model), while the CO reduction reaction was controlled by gas diffusion.


Introduction
Iron and steelmaking sector is one of the most important sectors due its great impact on the global growth, economy, and development. In recent years, the steel production rate has increased sharply [1]. By 2050, steel demand is expected to increase to 1.5 times higher than the current levels in order to meet the needs of a growing population [2,3]. About 70% of the total steel production relies directly on inputs from coke and coal [3,4]. The CO 2 emission from iron and steelmaking was 2.3 billion tons in 2007, while by 2050 it is expected to reach 3.0 billion tons [5]. Nowadays, the main challenges for the steel industry are energy consumption and environmental pollution. Therefore, the growing of the steel sector requires serious attention in order to establish processes that are economically viable and environmentally friendly [6].
Direct reduction (DR) process is one of the promising alternative ironmaking processes to effectively reduce CO 2 emissions. Therefore, recently great efforts have been paid to develop DR processes. The natural gas and non-coking coal are main reductants and sources of energy used in these DR processes [7]. In Sweden, three companies, steel manufacturer SSAB, mining company LKAB, and energy company Vattenfall, are exploring the use of hydrogen in steel production processes. This technology is known as Hydrogen Breakthrough Ironmaking Technology (HYBRIT) and aims to totally replace the coke with hydrogen gas [8,9]. In view of increasing demand for direct reduced iron in the manufacturing of steel, a good deal of emphasis is being given to promote the study of DR process [10]. Recently, attention has been paid on studying the kinetics and mechanisms of iron oxide reduction in the different gas mixture such as H 2 , CO, and CO/ H 2 gas mixtures.
The effect of H 2 -H 2 O on the iron ore pellet reduction was isothermally investigated at 700-1100°C [11]. The influence of a mixture of CO-CO 2 -H 2 -H 2 O-N 2 reducing gas was compared with that of a mixture of reducing gas containing only CO-CO 2 -N 2 . The results indicate that the reduction degree increases with increasing the temperature in both reducing gas atmospheres. Spreitzer et al. [12] indicated that, by increasing H 2 content in the reducing atmosphere, the reduction degree of pellets improved, owing to improving the diffusion rate due to hydrogen's small molecular size.
El-Geassy et al. studied the reduction mechanism of iron oxide in H 2 /CO mixtures in a temperature range of 800 to 1100°C [13]. They concluded that the reduction with either CO or CO/H 2 gas mixtures is controlled by mixed mechanism (combination of chemical reaction and gaseous diffusion) with a greater contribution of the gaseous diffusion as the CO content increases in the gas mixture. El-Geassy investigated the influence of doping with CaO and/or MgO on the stepwise reduction of hematite compacts [14]. They found that, during the first and second step (Fe 2 O 3 -Fe 3 O 4 and Fe 3 O 4 -FeO), the reduction rate increased in its early stages when any doping oxide was used. At later stage (FeO-Fe), the presence of MgO impeded the reduction process. El-Geassy et al. [15,16] studied the volume change of iron oxide compacts and the influence of gas composition on reduction behavior at 800 to 1100°C. They found that increased CO content in the gas mixture had a significant influence on swelling and maximum swelling (224%) was observed for samples reduced at 900°C. Decreasing CO content in the gas mixture led to a decreased in swelling by 24% in pure H 2 atmospheres.
The reduction of iron ore with gas takes place in either two or three stages, depending on the applied temperature. At temperature > 570°C, hematite is first reduced to magnetite then to wustite in a non-stoichiometric phase (Fe 1-y O) and finally to metallic iron. At temperatures ≤ 570°C, hematite is reduced to magnetite and then directly transformed into metallic iron [17]. Daniel et al. [18] study the reduction kinetics of hematite iron ore fines by hydrogen gas using a laboratory fluidized bed reactor in a temperature range 873-1073 K. They concluded that the reduction takes place in three steps: Fe 2 O 3 -Fe 3 O 4 , Fe 3 O 4 -FeO, and finally FeO-Fe. Lin et al. [19] were using temperature-programmed reduction (TPR) methods to study the reduction mechanism of iron oxide with hydrogen. They found that reduction takes place in two steps: the first step (Fe 2 O 3 -Fe 3 O 4 ) was confirmed by the unimolecular model, while the second step (Fe 3 O 4 -Fe) matched a nucleation twodimensional model according to the Avrami-Erofeev model.
Piotrowski et al. [20] investigated the effect of gas composition on the kinetics of iron oxide reduction by hydrogen. It was found that the reaction rate increases with both temperature and the hydrogen content. The reduction was initially proved to be a surface-controlled process. Once a thin layer of lower oxides (magnetite and wustite) formed on the surface, then the reduction mechanism shifts to control by gas diffusion. Bonalde et al. [21] investigated the reduction kinetics of hematite pellets at 850°C. The reduction experiments were carried out using H 2 , CO, and gas mixture. Reduction rates with H 2 were better than with CO. Skorianz [22] studied the reduction behavior of different iron ores with different porosities in a fluidized bed reactor. The reducibility of ores was demonstrated to be in this order: limonitic with the highest porosity hematite then magnetite with low porosity.
During the industrial DR process, the reduction of iron ore happens at variable temperatures (non-isothermal) rather than at a constant temperature (isothermal). Therefore, it is vital to understand the kinetic analysis of non-isothermal conditions, which would be a better choice for non-blast furnace ironmaking, i.e., the DR process and hydrogen metallurgy. The latter in recent years becomes quite popular in the global steel industry and in academic field.
Non-isothermal reduction reaction of Fe 2 O 3 nanopowder by hydrogen was studied at various heating rates [23]. It was found that the reduction process takes place in the sequence of Fe 2 O 3 to Fe 3 O 4 and Fe 3 O 4 to Fe. The apparent activation energy Ea for kinetic analysis was determined using the Kissinger-Akahira-Sunose (KAS) method. During the reduction step (Fe 3 O 4 to Fe), the reduction rate was retarded by sintering of oxide particles leading to an abrupt increase of activation energy. As the volume fraction of the metallic iron increased, the reaction rate accelerated, which can be attributed to the catalytic role of the reduced Fe particles facilitating reduction of the adjacent oxide particles. Jie Dang et al. [24] investigate the reduction of titanomagnetite (TTM) powders by hydrogen gas under a non-isothermal condition by using thermogravimetry (TG). It was found that the reduction of TTM proceeded via a dual-reaction mechanism. The first one was the reduction of TTM to wustite and ilmenite, while the second reaction was the reduction of wustite and ilmenite to metallic iron and titanium dioxide. The calculated apparent activation energies were 90-98 and 115-132 kJ/mol for the first and second reactions, respectively.
Tiernan et al. [25] investigate the reduction kinetics of iron oxide catalysts under non-isothermal conditions. Both temperature-programmed reductions (TPR) and constant rate temperature-programmed reduction (CR-TPR) were utilized. Fe 2 O 3 was found to reduce to Fe in two steps: Fe 2 O 3 -Fe 3 O 4 and Fe 3 O 4 -Fe. The rate-controlling mechanism was described by the random formation and growth of nuclei. CR-TPR profiles allowed the measurement of apparent activation energies as a function of the reduction extent. For the first step, Ea was found to be 96 kJ/mol. For the second step, Ea decreased from 69 to 59 kJ/mol as the reduction proceeded.
The present study aims to investigate the isothermal and non-isothermal reduction kinetics of iron ore fines in a hydrogen atmosphere, comparing the experimental results of H 2 with those of CO. The study will predict the kinetics and correlate these predictions with the morphological structure to deduce the reduction mechanisms.

Materials
The iron ore fines used in this study were collected from a DR plant (Al-Ezz El-Dekheila Steel Company, Alexandria) in the form of a waste material produced by the transportation and handling of pellets. The sample was characterized using X-ray diffraction (XRD; PW 1730 with Cu radiation, λ = 1.5Ǻ, at 40 kV and 30 mA). Figure 1 shows that hematite was the main mineral phase in the sample.
The chemical composition of the sample was identified using X-ray fluorescence (XRF advanced Axios, Netherlands); the results are listed in Table 1. Iron (III) oxide (Fe 2 O 3 ) was the main oxide in the sample accounting for about 96 wt.%.
Size analysis was conducted for representative sample of iron ore fines to determine the size distribution. A laboratory "Wedag" Ro-tap sieve shaker was employed for dry screening of sample. Shaking was performed at 60 Hz for 20 min, after which the weight of oversize on sieves was recorded to calculate weight passed and/or retained. The particle size distribution in the iron ore fine sample is shown in Table 2, indicating that more than 60% of the sample consisted of particles in the < 100 μm size range.
The structure of the sample was examined using both a scanning electron microscope (SEM JEOL, JSM-5410, 10Kv) and reflected light microscope (RLM, Zeiss Axio Scope A1) as shown in Fig. 2a and b, respectively. The SEM image shows that the structure contains fragments of different shapes and sizes which came out of pellet disintegration and abrasion during pellet transportation and handling. The RLM photomicrograph indicates that the internal structure includes a relatively dense matrix with more or less homogeneously distributed pores.

Compact Preparation
In order to ensure that the compacts were of uniform size and shape, iron ore fines were compressed in a stainless steel cylindrical mold. The compacts were created by mixing iron ore fines with 6% distilled water as a binder. Then, equal weights (2 g) of the mixture were compressed at 10 atm (10 kg/cm 2 ) using a hydraulic press. Compacts with a diameter of 10 mm and a height of 12 mm were dried at 110°C and kept in a dissector for subsequent analysis.
The total reducible oxygen content (w 0 ) of iron ore fines was estimated by preliminary reduction experiment at 1100°C using H 2 gas until the weight of the sample reaches constant (Eq. 1). Then, the extent of the reduction at time t (R t %) was calculated using Eq. 2.
where w i is the weight of the sample at zero time, w t is the  weight of the sample at time t, and w 0 is the total oxygen content.

Reduction System and Procedure
Thermogravimetry (TG) was used to measure both isothermal and non-isothermal reactions. The reduction apparatus is illustrated in Fig. 3. The main parts of the reduction system include a vertical tube furnace (F), fitting with an alumina reaction tube (A) and an automatic sensitive balance (B) at the top. The data from the balance (O) is connected to a recording unit (R) for continuous monitoring of the weight loss of samples during the experiments. The control parts (K) are two thermocouples T1 and T2. The temperatures of the furnace and sample are measured with two different thermocouples, one in the hot zone (T1) and the other behind the sample (T2). A gas supply system and gas purification unit (E) is used to remove moisture and other contaminants from the gases.

Isothermal Experimental Procedure
In experiments, the furnace is heated up to a target temperature and maintained at this temperature for around 5 min. Pure argon (Ar) gas is then introduced into the furnace at a flow rate of 1.0 l/min. After that, the sample is placed in a basket which hangs from the balance arm by a wire; then, the basket is positioned in the hot zone of the furnace. Under these conditions, Ar gas is replaced by 1.0 l/min of H 2 . During the experiment, the weight loss of the sample due to O 2 removal is continuously recorded until the weight remains constant. At the end of the experiment, H 2 gas is replaced by Ar gas and the furnace is cooled back down to room temperature. The reduced sample is taken out and kept in a desiccator for later characterization.

Non-Isothermal Experimental Procedure
In these tests, the furnace is first flushed with Ar at a flow rate of 0.2 l/min and held in this condition for 5 min. Then, the Ar gas switched to H 2 , which introduced at a flow rate of 1.0 l/min. The sample is heated in the reducing gas from room temperature to 1000°C at different heating rates (5, 10, 15, and 20°C/min). At the end of the experiment, the H 2 is replaced by Ar gas and the furnace is cooled back down to room temperature. The reduced sample is taken out and kept in a desiccator for later characterization.

Reduction Behavior
The reduction curves of Fe 2 O 3 compacts by H 2 gas at temperatures 800 to 1100°C are illustrated in Fig. 4. For any given reduction temperature, the rate of reduction is highest in the early stages and gradually decreases until the end of the reduction reaction. At low temperatures, the reduction takes longer and the reduction time decreases as reduction  In the case of reduction using CO at temperatures 80--900°C , samples did not complete reduction and stopped at different extents of reduction, depending on the applied temperature [26]. Moreover, during experiments, an increase in the weight of the compacts was observed. This phenomenon previously observed elsewhere [27] indicates that the rate of carbon participation as a result of the thermal decomposition of carbon monoxide (2CO = C + CO 2 ) is greater than the rate of oxygen removal.
The microstructure of the completely reduced compacts with H 2 at 950 and 1100°C is given in Fig. 5a and b, respectively. Figure 5a illustrates that the structure includes elongated grains of metallic iron and coalesces to each other with separated small macropores. As the reduction temperature increases, the matrix structure became relatively dense due to the coalescence and connection of metallic grains as shown in Fig. 5b.

Kinetics and Mechanism of Iron Ore Reduction
According to previous literature [12], iron ore reduction is supposed to take place from the surface and progress to the Fig. 3 Schematic diagram for reduction apparatus of TGA system Fig. 4 Reduction behaviors of iron ore compacts at 800-1100°C by pure H 2 core and the diameter of the unreacted core becomes smaller and smaller until zero. The reduction mechanisms of iron ore during both the initial and final stages can be deduced from the correlations between the apparent activation energy values, the gas-solid chemical reaction, and the microstructure of the reduced samples in the reduction procedure [26,[28][29][30][31]. The kinetics of iron ore reduction can be described using a shrinking unreacted core model to elucidate the reduction reaction. In the beginning, iron ore particles are assumed to have 100% hematite and have an initial diameter. With the increase in the reduction reaction, a product layer of iron is created and the unreacted core of pellet particle is decreased to a smaller diameter. A layer of reactant gas exists at the initial diameter of the iron ore particle. A chemical reduction reaction can be divided into three stages: external diffusion, internal diffusion, and chemical reaction. The chemical reaction is expressed by the Arrhenius equation in Eq. 3: where K r is the reaction rate constant, K o is the frequency factor, Ea is the apparent activation energy, R is the gas constant, and T is the absolute temperature. The rate of reaction (K) is calculated using the isothermal curves given in Fig. 4. The relationships between lnK and (1/T) are plotted during both the initial stages (5-15%) and the final stages of reduction (85-95%) as shown in Fig. 6. The computed apparent activation energy values for the reduction by H 2 and CO are listed in Table 3. The comparison of activation energy values indicated that reduction with H 2 is more efficient than with CO.
At the initial stages, the Ea value indicates that the reduction process seems to be controlled by a mixed mechanism of gaseous diffusion and interfacial chemical reaction. During the final stages, the Ea value reveals that the reaction is controlled by interfacial chemical reaction and this matched with other studies [17,26,32].
The mathematical models that are typically adopted for predicting kinetics mechanisms in heterogeneous gas-solid reactions and the grain model developed by Szekely et al. [28] are the most accurate one for the similar experimental conditions, in which long cylinder samples are produced from compacting fines of solid materials. Gaseous diffusion, interfacial chemical reaction, and mixed control mathematical formulae have been tested in this study to evaluate the control stage in the overall reduction process. The three mathematical kinetic equations are as follows [28,32]: (a) For gaseous diffusion mechanism: (b) For chemical reaction mechanism: (c) For mixed control mechanism: where K and X express the rate of reaction and the fractional reduction degree, respectively. In order to confirm the reduction mechanism prediction from the activation energy values, the above chemical reaction   [26] formulations were tested against experimental results as presented in Figs. 7 and 8. During the initial stages, linear correlations were provided from the application of the mixed control mechanism (Eq. 6) as shown in Fig. 7. During the final stages, linear correlations were obtained when applying an interfacial chemical reaction (Eq. 5) as illustrated in Fig. 8. These are consistent with the mechanisms suggested by the activation energy values. The morphological structures of the reduced samples at 1000°C during the initial and final stages are shown in Fig.  9a, b, c, d, respectively. Figure 9 a shows the formation of a relatively porous structure which facilitates the access of gas diffusion. Figure 9b shows the formation of the outer/middle interface (red line) of the compact, which is comprised of ironwustite phase diffusion in the wustite matrix. This indicates that the reducing gas is easily diffused through the outer shell and consequently, iron is nucleated and grows up on the surface of the wustite particles allowing the development of the iron-wustite intergrowth layer. Resistance to gas diffusion decreases near the surface of wustite grains and the chemical reaction resistance has a major effect as a controlling mechanism. Therefore, mixed mechanism combines gas diffusion and interfacial chemical reaction controlling the reaction. Figure 9c illustrates that the grains of metallic iron are coalescing to each other with small separate macropores. As the reduction proceeds, large macropores are observed due to the coalescing of metallic iron grains as shown Fig. 9d. Accordingly, gas diffusion becomes easier through relatively porous structures, so the reaction is controlled by chemical reaction mechanism.

Reduction Behavior
Non-isothermal experiments were carried out at temperatures ranging from room temperature up to 1000°C in H 2 atmosphere using different heating rates (5, 10, 15, and 20°C/min). The results are shown in Fig. 10. The heating rate has a considerable effect on the reduction process. The degree of reduction decreases as the heating rate increases. Increasing heating rates leads to shorter time for the sample at any particular temperature. At the same heating rate, the degree of reduction is higher in H 2 than in CO atmosphere [26].
For any given reduction curve presented in Fig. 10, the reduction process can be categorized into three steps: i. Incubation step, from the beginning up to ≈ 400°C, during which no reduction reaction occurs. ii. Induction step, from temperature > 400 up to 800°C, where the reduction rate greatly increases with the rise in temperature and depending on the applied heating rate. iii. Deceleration step, over 800°C, at which the reduction rate decreases until the end of experiment.

Kinetics and Mechanism
There are many methods available for evaluating nonisothermal kinetic parameters [33,34]. This study employs two mathematical analysis methods to determine the kinetic parameters: a model-free isoconversional analysis and model fitting [35,36]. The overall reaction rate can be generally described by [34] Under non-isothermal condition, when the sample is heated at a constant heating rate, the explicit temporal in Eq. 7 is eliminated through the trivial transformation where β is the heating rate. Upon integration, Eq. (8) gives Model-free isoconversional methods estimate the kinetic parameters without modelistic assumptions [35][36][37]. The basic assumption of this method is that the reaction rate at a given conversion degree is only a function of temperature [38]. A typical example of this kind of approach is the Kissinger-Akahira-Sunose (KAS Eq. 10) method, which has been widely used to analyze the reaction of kinetic parameters as a function of the conversion degree. ln β Figure 11 shows the KAS plots of ln β/T 2 against 1/T for different conversion degree values. The apparent activation energy values listed in Table 4 were computed using the slopes of the lines drawn. With the increase in reduction during the process, the Ea value becomes smaller and smaller from 92.57 to 20.19. In order to predict the most probable reduction mechanism, model-fitting methods were employed. The model offering the best statistical fit was chosen as the reaction model that was used to determine the kinetic parameters. Coats and Redfern (CR) have provided one of the most popular modelfitting methods [39].
ln g a ð Þ where T is the mean experimental temperature and g(a) is a function that mathematically describes the reaction kinetics. Several models tabulated elsewhere [40] for the integral form g(α) were tested according to CR method. The plot of ln [g (α)/T 2 ] versus 1/T will be linear with a high correlation coefficient. Figure 12 reveals that the Avrami-Erofeev model provides a linear trend with a correlation coefficient (R) higher than 0.99, indicating that the reduction is controlled by the Avrami-Erofeev model.
The calculated activation energy values (Ea) at different conversion degrees (α) are presented in Fig. 13. Higher activation energy values were obtained at earlier stages and decrease as the reduction extents increased. The reduction starts at a lower temperature and requires high energy (Ea) to proceed. The decreasing in Ea values can be attributed to the effect of temperature, which has a considerable effect on the reduction process.
The morphological changes in the completely reduced samples at 20°C/min and 5°C/min analyzed using an SEM microscope are shown in Fig. 14a and b, respectively. Figure 14a illustrates that the microstructure of iron ore includes significant porosity and the formation of metallic iron whiskers. Furthermore, the formation of dense structures of   Fig. 13 Calculation Ea against conversion degree of iron ore compact reduced by H 2 Fig. 12 Coats-Redfern plot of reduction of iron ore compacts with H 2 gas metallic iron grain results from the connections between the metallic iron grains as shown in Fig. 14b. The microstructure of fully reduced samples at different heating rates (20,15,10, and 5°C/min) examined by a reflected light microscope is presented in Fig. 15a-d. It can be observed that the heating rate has an obvious effect on the structure of the iron grains. Figure 15a (20°C/min) shows a structure that includes bright white grains of metallic iron scattering all over the matrix. Figure 15b (15°C/min) illustrates that the grains of metallic iron connected together forming a relatively dense matrix with small separate macropores. By decreasing the heating rate as shown in Fig.  15c (10°C/min), the matrix includes large grains of metallic iron which contain micropores. As shown in Fig. 15d (5°C /min), denser structure can be observed attributable to sintering and the connection of metallic iron grains. These above results indicate that the microstructure of samples becomes denser as the heating rate decreases. This is most likely attributable to the relatively longer reaction time, which allows more time for growth of the metallic iron grains. Figure 16 shows the XRD patterns of non-isothermal reduction products at different heating rates and at 1000°C. The intensity of the metallic iron peaks increases as heating rates decrease, because of the longer retention time of the sample, and consequently enhances the recrystallization of the iron phase, which is evident from Fig. 15. Understanding the kinetic analysis of non-isothermal conditions is helpful in the DR process and in hydrogen metallurgy. The latter has become important topic in the global steel industry and in academic fields in recent years. It provides promising solution that could be very effective in decreasing the CO 2 emissions and also metallurgical coke utilization in iron and steel industry by optimization of hydrogen-rich materials such as natural gas in the ironmaking processes. Such solutions require intensive investigations in order to understand and precisely evaluate their influence on iron and steel industry and CO 2 mitigation. This study investigated isothermal and non-isothermal reductions of iron ore fines in a pure H 2 atmosphere. In the experiments, total weight loss was continuously recorded using TG and the kinetic mechanisms were studied. The experimental results and kinetic analyses suggest the following conclusions: 1. Temperature plays a significant role in the isothermal reduction process. The reduction rate during the initial stages is higher than the reduction rate during the final stages. 2. The rate-controlling mechanism of the reduction process was predicted using an activation energy calculation, through applying gas-solid mathematical models, and microstructure examination. During the initial stages, the reduction reaction was controlled by a mixed mechanism of interfacial chemical reaction and gas diffusion. During the final stages, the interfacial chemical reaction was the rate-determining step. 3. The heating rate has a significant effect on the reduction rate and reduction conversion for non-isothermal reductions. With identical heating rates, the degree of reduction was higher in H 2 than in CO atmosphere. 4. The kinetic mechanisms of the reductions could be predicted by applying model-free and model-fitting methods. The wide range of activation energy values (92.57 to 20.19) in the model-free approach was mainly a result of temperature differences. When applying the Coats-Redfern method, the reduction reaction was controlled by the Avrami-Erofeev model.

Compliance with Ethical Standards
Conflict of Interest The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.