Transient Cold Flow Simulation of Fast-Fluidized Bed Air Reactor with Hematite as an Oxygen Carrier for Chemical Looping Combustion

: Chemical looping combustion (CLC) is the most reliable carbon capture technology for curtailing CO 2 insertion into the atmosphere. This paper presents the cold ﬂow simulation results necessary to understand the hydrodynamic viability of the fast-ﬂuidized bed air reactor. Hematite is selected as an oxygen carrier due to its easy availability and active nature during the reactions. The dense discrete phase model (DDPM) approach using the commercial software Ansys Fluent is applied in the simulation. An accurate and stable solution is achieved using the second-order upwind numerical scheme. A pressure difference of 150 kPa is obtained between the outlet and inlet of the selected air reactor, which is necessary for the movement of the particle. The stable circulating rate of hematite is achieved after 28 s of particle injection inside the air reactor. The results have been validated from the experimental results taken from the literature.


Introduction
Carbon dioxide (CO 2 ) accumulation in the atmosphere, from different sources, is a serious global concern. Global mean temperatures are rising much faster and if the rate of CO 2 emission remains the same for the coming years, then the problem of global warming will increase exponentially. According to the Mauna loa observatory, Hawaii [1], which tracks the CO 2 presence in the atmosphere daily, the concentration of CO 2 in the atmosphere was determined to be 415.52 ppm in January 2021. This value is far more than 280 ppm, which is the safest level according to climate scientists.
The majority of CO 2 emissions are generated by coal-fired power plants, steel and iron industries, cement industries, and automobiles. In a developing nation like India, the backbone of the power sector is coal, which contributes to 55.65% of the total installed capacity [2]. Therefore, our primary focus should be to reduce the emissions from coal-fired thermal power plants.
Although many technologies have already been developed to reduce emissions, such as integrated gasification combined cycle (IGCC) and ultra-supercritical boilers, due to the energy penalty and high initial investment cost of these technologies, their implementation on large scale is economically not feasible. Carbon capture and storage (CCS) is one of the most reliable solutions to cut this increasing rate of CO 2 to the safest level. At least 90% of emissions from the world's largest sources can be captured by this technology [3,4]. CCS is a three-stage process in which CO 2 is captured from the sources first, and then it is transported to a suitable geological site, where it can be stored for a long time. The technologies used in CCS are well known to the world now, but their implementation on a large scale is still not achieved due to the high energy penalty associated with them. On (1) fuel reactor and (2) an air reactor, as shown in Figure 1. The fuel reactor operates as a 51 bubbling bed and the air reactor operates in a fast fluidization regime. In the CLC process, 52 air and fuel (coal) do not directly mix or react. The oxygen required for the combustion of 53 coal is provided by a metal oxide (MexOy), which is also called an oxygen carrier (OC). 54 These OCs are continuously circulated from the air reactor to the fuel reactor. In the air 55 reactor, OC is oxidized by the air, while in the fuel reactor, the OC supplies oxygen for the 56 combustion of the fuel. In this way, this process loop continues [5]. The products of this 57 combustion reaction are CO2 and steam. It is possible to isolate the steam from CO2 by 58 cooling and condensation. 59 The reactions [6] which took place inside the fuel reactor and the air reactor, 60 respectively, are 61 Fuel + nMexOy ⟶ nMexOy-1 + H2O + CO2 (1) MexOy-1+½O2⟶MexOy (2) where m and n are the numbers of moles of the reactants and the products. Since the 62 ignition of fuel in the fuel reactor takes place without air, only oxygen is utilized for the 63 combustion, which is given by the OC. The vent stream from the fuel reactor is not 64 polluted or weakened by different gases, for example, nitrogen. This gives a high-virtue 65 CO2 stream, accessible to retrieve at the fuel reactor outlet without the requirement of an 66 energy-expensive process for separation. The only energy cost associated with CLC is the 67 separation of OC from ash. 68

70
OC is the most important part of the CLC technique. The ideal OC should be highly 71 reactive. It should immediately oxidize in the air reactor, as well as reduce in the fuel 72 reactor for good combustion efficiency. It should also be stable, and show good resistance 73 towards fragmentation and attrition. In addition, it should have good fluidization 74 properties so that it does not agglomerate. The materials that can be used for OC are 75 copper, nickel, manganese and iron oxides, but these are not limited, as other elements 76 have also been developed [7]. Several numerical studies have described the reduction and 77 oxidation behavior of oxygen carriers [8][9][10][11][12][13]. 78 Various researchers have performed simulations on CLC using different OCs [14-79 21]. Deng et al. [14] developed a reaction kinetic model of a fuel reactor using a mixed OC 80 of calcium sulphate and hydrogen (CaSO4+H2). The computational result showed a low 81 fuel conversion rate with this OC. Arjmand et al. [20] evaluated the performances of 82 The reactions [6] which took place inside the fuel reactor and the air reactor, respectively, are Fuel + nMe x O y →nMe x Oy-1 + H 2 O + CO 2 (1) where m and n are the numbers of moles of the reactants and the products. Since the ignition of fuel in the fuel reactor takes place without air, only oxygen is utilized for the combustion, which is given by the OC. The vent stream from the fuel reactor is not polluted or weakened by different gases, for example, nitrogen. This gives a high-virtue CO 2 stream, accessible to retrieve at the fuel reactor outlet without the requirement of an energy-expensive process for separation. The only energy cost associated with CLC is the separation of OC from ash. OC is the most important part of the CLC technique. The ideal OC should be highly reactive. It should immediately oxidize in the air reactor, as well as reduce in the fuel reactor for good combustion efficiency. It should also be stable, and show good resistance towards fragmentation and attrition. In addition, it should have good fluidization properties so that it does not agglomerate. The materials that can be used for OC are copper, nickel, manganese and iron oxides, but these are not limited, as other elements have also been developed [7]. Several numerical studies have described the reduction and oxidation behavior of oxygen carriers [8][9][10][11][12][13].
Various researchers have performed simulations on CLC using different OCs [14][15][16][17][18][19][20][21]. Deng et al. [14] developed a reaction kinetic model of a fuel reactor using a mixed OC of calcium sulphate and hydrogen (CaSO 4 +H 2 ). The computational result showed a low fuel conversion rate with this OC. Arjmand et al. [20] evaluated the performances of different manganese ores as OCs for CLC. They found that the high char conversion rate of manganese makes it a good choice as an OC. Menon and Patnaikuni [21] simulated the fuel reactor to compare the performances of two OCs, namely, ferric oxide (Fe 2 O 3 ) and copper oxide (CuO), as regards the interaction between particle and particle, and fluid and particle, and the reaction rates. They found that CuO has a better performance when the ash agglomeration effect is not considered, and Fe 2 O 3 performed better when the ash agglomeration effect was considered. Haider et al. [22] experimentally investigated the cold flow performance of a dual circulating fluidized bed system to understand the hydrodynamic behavior of a pilot-scale advanced CO 2 capture technology (PACT) facility at Cranfield, UK. FE100 and molochite were used as the OC. They have concluded that both the risers share similar pressure and density profiles. The circulation rate of OC can be changed by altering the fluidization velocity and bed inventory. However, a stable circulation rate was achieved using pneumatic transport.
CFD modeling techniques and the fundamental equations used in the development of the CFD models have been reviewed by Singh et al. [23]. Shuai et al. [24] developed a 2D CFD model for CLC using inter-connected fluidized beds. They used an Eulerian continuum and two fluid models for both the gas phase and the solid phase. Using their model, they successfully described the hydrodynamics of gas and solid particles in the CLC process. Recently, Hamidouche et al. [25] performed a numerical simulation using the multi-phase particle-in-cell (MP-PIC) method for a biomass-fueled CLC. They found that the fluidization velocity of the air reactor greatly affects the hydrodynamic behavior. Based on their results, they also stated that the pressure balance is very important in order to minimize the gas leakage from the reactor.
From the above discussion of the literature, it can be concluded that in the development of CLC, various researchers have opted for computational approaches to study various parameters in the operation of CLC. The computational approach mainly consists of CFDbased simulation, numerical simulation, thermodynamic computation, ASPEN Plus-based simulation, the DSMC model, the direct simulation Monte Carlo method, etc., to model the process and to study the behavior resulting from variations in design and operating parameters [26]. During their research, various OCs have been investigated, but the cold flow simulation of an air reactor using hematite as the OC has not been performed. In this work, a numerical investigation using the cold flow simulation of the air reactor is performed to understand the hydrodynamic behavior of the system. The hydrodynamic study of the CLC reactors is very crucial in order to analyze the circulation of OC, the interaction among the particles, the pressure drop across the system, and the overall configuration of the system. Hematite, which is an ore of iron, is selected as OC. Iron is readily available, and its low cost makes it an attractive option for OC. In addition, its active nature of being easily oxidized and reduced during the reaction is an added benefit. The simulation results are validated with the experimental results performed by Haider et al. [22] at Cranfield, UK. For this, FE100 is considered as the OC, and the variations in static pressure at different heights of the air reactor are compared.

Materials and Methods
Ansys Fluent 17.2 was chosen for carrying out the simulation process on a six-core Intel Xenon CPU. The geometry of the air reactor is shown in Figure 2. It is a similar geometry to that used in the experimental work performed by Haider et al. [22]. The air is chosen as the primary phase and hematite as the secondary phase. Since two phases are involved in the simulation, the multiphase modeling approach is used. The Eulerian or dense discrete phase model (DDPM) approach is found suitable according to the gas-solid interaction flow regime in the granular form [27]. For each phase, the continuity and momentum equations are solved. The coupling of particles is achieved through pressure and interphase coefficients. The transient multiphase model is solved using the phasecoupled SIMPLE algorithm. It solves the momentum equation of both the phases in an isolated manner. Since the nature of these equations is highly coupled, it also considers some correction factors to solve the shared pressure and volume fraction between the two phases. The governing equations used in this simulation are described next.

Continuity Equation 140
The volume fraction is solved using the following equation [27]: 141 where Vz = volume, α = volume fraction, ρ = density and ρ z = effective density of each 142 phase. 143

Momentum Equation 144
The conservation of momentum equation for gas (fluid) is solved by the following 145 equation 146 where p = pressure of fluid, τ = shear stress of fluid, u = velocity of the fluid phase, g = 147 acceleration due to gravity and R = momentum transfer from the fluid to the solid [28]. 148 Similarly, the momentum equation used for the solid phase to define the trajectory 149 of particles is given by 150 The terms on the right-hand side in Equation (7) are forces due to gravity, interphase 151 drag, and particle-particle collisions, respectively. Here, FD is the drag coefficient given 152 by the following equation [28]: 153 where C is the coefficient of drag, dp is particle diameter, and Re represents the 154 Reynolds number, which is given by [28] 155 The interaction between the phases is defined using the Syamlal and O'Brien model 156 [29]. This model is significant in characterizing the motion of particles in a fluidized bed 157 because it corrects the terminal velocity, which is the minimum velocity required for 158 lifting the particle so that it can move outside the bed. 159

Continuity Equation
The volume fraction is solved using the following equation [27]: where V z = volume, α z = volume fraction, ρ z = density andρ z = effective density of each phase.

Momentum Equation
The conservation of momentum equation for gas (fluid) is solved by the following equation where p f = pressure of fluid, τ f = shear stress of fluid, u f = velocity of the fluid phase, g = acceleration due to gravity and R sg = momentum transfer from the fluid to the solid [28]. Similarly, the momentum equation used for the solid phase to define the trajectory of particles is given by The terms on the right-hand side in Equation (7) are forces due to gravity, interphase drag, and particle-particle collisions, respectively. Here, F D is the drag coefficient given by the following equation [28]: where C D is the coefficient of drag, d p is particle diameter, and Re p represents the Reynolds number, which is given by [28] Re The interaction between the phases is defined using the Syamlal and O'Brien model [29]. This model is significant in characterizing the motion of particles in a fluidized bed because it corrects the terminal velocity, which is the minimum velocity required for lifting the particle so that it can move outside the bed.
where V r,P is the terminal velocity correction factor, which is given by 85}. The force due to the collision of the particle is obtained by the particle pressure, which is obtained by the kinetic theory of granular flows (KTGF), given by [28]

Results
The geometry of the air reactor considered in this simulation is taken to be similar to that in the experimental work of Haider et al. [22]. The height of the riser is 7.3 m, with an interior diameter of 0.1 m, while the diameter of the inlet and outlet pipes is 0.04 m. Figure 3a,b show the mesh generated on the upper and lower parts of the computational domain. Figure 3c shows the top view of the mesh. The grid size with the mesh count of 516,332 elements is found to be optimum. At the inlet and outlet of the particles, the meshing is coarse because of the maximum number of interactions between the particles. In the current simulation, hematite is chosen as the OC with a particle diameter of 60 µm and a density of 5150 kg/m 3 . Simulation is done to understand the distribution of the pressure and volume fraction of particles at different bed heights. During simulation, the solid mass flow circulation is kept constant, considering the experimental conditions. An accurate and stable solution is achieved using the second-order upwind numerical scheme. For injection purposes, the standard parcel release method is used [30,31].
where V , is the terminal velocity correction factor, which is given by 160 where A = ∝ . and The force due to the collision of the particle is obtained by the particle pressure, which 162 is obtained by the kinetic theory of granular flows (KTGF), given by [28]

Results 164
The geometry of the air reactor considered in this simulation is taken to be similar to 165 that in the experimental work of Haider et al. [22]. Initially, the inlet velocity of the air from the bottom and the velocity of the injection of the OC is taken as zero. The injection of OC occurs at the particle inlet for 1 s only and almost 198,565 particles are released. The velocity of the air is maintained at zero for the next second so that the particle that falls due to gravitation settles down at the bottom, and its kinetic energy becomes negligible. This is the initial boundary condition of the simulation. After the bed of OC is formed, the inlet velocity of the air is increased gradually up to 2.50 m/s, which is slightly greater than the minimum fluidization velocity for the particle [22]. This velocity is supposed to lift the particle outside the reactor towards the mixture outlet. When the air velocity is increased, the particles in the settled bed experience sudden lift force, and they begin to move upward. Particles that are at the top of the bed reach the outlet of the reactor at around 15 s. The whole simulation runs for almost 45 days due to the large number of particles. The results of the simulation are monitored at different heights, i.e., 0 m (air inlet), 0.5 m, 0.75 m, 6 m, 7 m, and mixture outlet. The mass flow rate of the particle is set as 0.04 kg/s, similar to the solid circulation rate of Haider et al. [22]. Various other key modeling parameters required to run the simulation are given in Table 1.  Figure 4 shows the variation in the mass flow rate with respect to time at the outlet of the reactor. It can be concluded that at the start of the simulation, particles escape the reactor at a discontinuous rate. However, after 28 s, the particles start to escape almost at a continuous rate of 0.0007 kg/s. This is of the same order as was found in the experiment performed by Haider et al. [22]. The stable circulation of hematite is achieved, which shows it can be used as an OC. in the reactor. This is an important study because it is the pressure difference of the 207 particles at the inlet and outlet that drives them into the reactor. It is evident from the 208 graph that the pressure at the mixture outlet is zero, while the pressure inside the reactor 209 after 15 s is almost 1.5 kPa at any height. Since the height of the particle bed is much lower 210 compared to the height of the reactor, the pressure differences along the height are almost 211 the same. It is also noticeable from the figure that up to 2 s, the pressure inside the reactor 212 is also zero, because the particles are settling on the bed during this period in the reactor. 213 After the settling period is over, an instantaneous spike in pressure can be seen, the 214 magnitude of which is different at different heights. This is because when the air is injected 215 through the bed, at that moment, it experiences the resistance from the particle bed. The 216 magnitude of the pressure spike is its highest (6.2 kPa) at air inlet and its lowest (1.9 kPa) 217 at the top of the reactor for a 7 m height. The particle inside the bed begins to move 218 upward as a result of this pressure buildup. The bed stabilizes after 15 s when the fast 219 fluidization regime is fully developed. 220  Figure 5 shows the variation in static pressure (gauge) with time at different heights in the reactor. This is an important study because it is the pressure difference of the particles at the inlet and outlet that drives them into the reactor. It is evident from the graph that the pressure at the mixture outlet is zero, while the pressure inside the reactor after 15 s is almost 1.5 kPa at any height. Since the height of the particle bed is much lower compared to the height of the reactor, the pressure differences along the height are almost the same. It is also noticeable from the figure that up to 2 s, the pressure inside the reactor is also zero, because the particles are settling on the bed during this period in the reactor. After the settling period is over, an instantaneous spike in pressure can be seen, the magnitude of which is different at different heights. This is because when the air is injected through the bed, at that moment, it experiences the resistance from the particle bed. The magnitude of the pressure spike is its highest (6.2 kPa) at air inlet and its lowest (1.9 kPa) at the top of the reactor for a 7 m height. The particle inside the bed begins to move upward as a result of this pressure buildup. The bed stabilizes after 15 s when the fast fluidization regime is fully developed. in the reactor. This is an important study because it is the pressure difference of the 207 particles at the inlet and outlet that drives them into the reactor. It is evident from the 208 graph that the pressure at the mixture outlet is zero, while the pressure inside the reactor 209 after 15 s is almost 1.5 kPa at any height. Since the height of the particle bed is much lower 210 compared to the height of the reactor, the pressure differences along the height are almost 211 the same. It is also noticeable from the figure that up to 2 s, the pressure inside the reactor 212 is also zero, because the particles are settling on the bed during this period in the reactor. 213 After the settling period is over, an instantaneous spike in pressure can be seen, the 214 magnitude of which is different at different heights. This is because when the air is injected 215 through the bed, at that moment, it experiences the resistance from the particle bed. The 216 magnitude of the pressure spike is its highest (6.2 kPa) at air inlet and its lowest (1.9 kPa) 217 at the top of the reactor for a 7 m height. The particle inside the bed begins to move 218 upward as a result of this pressure buildup. The bed stabilizes after 15 s when the fast 219 fluidization regime is fully developed. 220 The particle volume fraction is an important parameter because it helps to understand the fractional distribution of the solid particle at different heights during different time intervals. The volume fraction of solid particles with respect to time at different heights is shown in Figure 6. The volume fraction of solids is much lower because of the small volume occupied by them, as compared to the volume of the reactor. However, the surge in volume fraction is noticed initially at a height of 0.5 m and 0.75 m. The maximum volume fraction is observed at a height of 0.5 m at 2 s. This happens because of the gathering of the particles at the bottom planes of the reactor, as a result of particle movement. The particle volume fraction at higher planes is not affected at that instant due to the small number of particles transferred from the bottom. A stable particle distribution is observed after 28 s when the oscillations in volume fraction are reduced.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 8 of 12 The particle volume fraction is an important parameter because it helps to 223 understand the fractional distribution of the solid particle at different heights during 224 different time intervals. The volume fraction of solid particles with respect to time at 225 different heights is shown in Figure 6. The volume fraction of solids is much lower because 226 of the small volume occupied by them, as compared to the volume of the reactor. 227 However, the surge in volume fraction is noticed initially at a height of 0.5 m and 0.75 m. 228 The maximum volume fraction is observed at a height of 0.5 m at 2 s. This happens because 229 of the gathering of the particles at the bottom planes of the reactor, as a result of particle 230 movement. The particle volume fraction at higher planes is not affected at that instant due 231 to the small number of particles transferred from the bottom. A stable particle distribution 232 is observed after 28 s when the oscillations in volume fraction are reduced. 233  The velocity of the gas phase (air) plays an important role in fluidizing the bed as well as lifting the particles from the bed to the outlet. Figure 7 shows the transient variation in air velocity along different heights of the reactor.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 8 of 12 The particle volume fraction is an important parameter because it helps to 223 understand the fractional distribution of the solid particle at different heights during 224 different time intervals. The volume fraction of solid particles with respect to time at 225 different heights is shown in Figure 6. The volume fraction of solids is much lower because 226 of the small volume occupied by them, as compared to the volume of the reactor. 227 However, the surge in volume fraction is noticed initially at a height of 0.5 m and 0.75 m. 228 The maximum volume fraction is observed at a height of 0.5 m at 2 s. This happens because 229 of the gathering of the particles at the bottom planes of the reactor, as a result of particle 230 movement. The particle volume fraction at higher planes is not affected at that instant due 231 to the small number of particles transferred from the bottom. A stable particle distribution 232 is observed after 28 s when the oscillations in volume fraction are reduced. 233 It can be observed from the graph that the initial velocity of air at the air inlet is 0 and after 2 s it is 2.5 m/s, according to the boundary condition. The velocity of air at a height of 0.25 m increases due to the internal turbulence created by the particles injected Appl. Sci. 2021, 11, 2288 9 of 12 into the reactor for 1 sec. However, this decreases as the particles get settled down at the bed. As the bed settles and the velocity of air at the inlet is increased gradually up to the terminal velocity, there is a sudden spike in the velocity, which can be seen in the graph along with the different heights of the reactor. This happens because of the turbulence created by the particle and the increasing air velocity at the inlet. However, once stable circulation is achieved, the velocity becomes constant to maintain the continuous flow of the particles to the outside. Figure 8 gives the particle velocity profile at different heights of the reactor with time. The velocity of the particles depends on the air velocity from the bottom of the reactor. Initially, the particle velocity is zero at all the heights of the reactor, but it gradually increases with the increment in the air velocity. The maximum variation in the velocity is observed at the outlet of the reactor. At the outlet of the reactor, the particles attain the maximum velocity of 13.74 m/s because of the decreased area of the outlet pipes. However, this decreases as the number of particles inside the reactor decreases. It can be observed from the graph that the initial velocity of air at the air inlet is 0 and 241 after 2 s it is 2.5 m/s, according to the boundary condition. The velocity of air at a height 242 of 0.25 m increases due to the internal turbulence created by the particles injected into the 243 reactor for 1 sec. However, this decreases as the particles get settled down at the bed. As 244 the bed settles and the velocity of air at the inlet is increased gradually up to the terminal 245 velocity, there is a sudden spike in the velocity, which can be seen in the graph along with 246 the different heights of the reactor. This happens because of the turbulence created by the 247 particle and the increasing air velocity at the inlet. However, once stable circulation is 248 achieved, the velocity becomes constant to maintain the continuous flow of the particles 249 to the outside. 250 Figure 8 gives the particle velocity profile at different heights of the reactor with time. 251 The velocity of the particles depends on the air velocity from the bottom of the reactor. 252 Initially, the particle velocity is zero at all the heights of the reactor, but it gradually 253 increases

259
For the validation, the simulation with FE100 replacing hematite as the OC has now 260 been carried out, and the values of static pressure at different heights are compared with 261 the experimental values presented in Haider et al. [22] The parameters including reactor 262 geometry, the diameter of OC, maximum fluidization velocity, etc., are taken from the 263 data available in the experimental work [22]. 264 For the validation, the simulation with FE100 replacing hematite as the OC has now been carried out, and the values of static pressure at different heights are compared with the experimental values presented in Haider et al. [22] The parameters including reactor geometry, the diameter of OC, maximum fluidization velocity, etc., are taken from the data available in the experimental work [22]. Figure 9 shows the comparison of experimental and numerical data of static pressure at various heights. It can be observed that the graphs of static pressure are almost of the same magnitude at different heights, except for the difference at the inlet. This difference occurs because of the different initial processes of the simulation and experimental work. In the experiment, the particles are injected from the narrow tube thoroughly. On the other hand, in the simulation, the particles are considered to be injected from some height of the reactor, and they are then left for 1 s to settle before the air injection. Changing the drag law can also help in minimizing the deviation of numerical results from the experimental result [32]. The decay in static pressure after the inlet is found to be consistent with the experimental results, along with the height of the reactor. Figure 10 shows the error bar graph for the difference between experimental and numerical data. It can be seen from the graph that the errors are much smaller at any height, except at the inlet. Appl. Sci. 2021, 11, x FOR PEER REVIEW 10 of 12 265 Figure 9 Comparison of numerical and experimental data of static pressure at various heights.
266 Figure 9 shows the comparison of experimental and numerical data of static pressure 267 at various heights. It can be observed that the graphs of static pressure are almost of the 268 same magnitude at different heights, except for the difference at the inlet. This difference 269 occurs because of the different initial processes of the simulation and experimental work. 270 In the experiment, the particles are injected from the narrow tube thoroughly. On the other 271 hand, in the simulation, the particles are considered to be injected from some height of the 272 reactor, and they are then left for 1 s to settle before the air injection. Changing the drag 273 law can also help in minimizing the deviation of numerical results from the experimental 274 result [32]. The decay in static pressure after the inlet is found to be consistent with the 275 experimental results, along with the height of the reactor. Figure 10 shows the error bar 276 graph for the difference between experimental and numerical data. It can be seen from the 277 graph that the errors are much smaller at any height, except at the inlet. 278

Conclusions 281
In this work, the DDPM model of Ansys Fluent 17 was used to carry out a cold flow 282 simulation of the fast-fluidized bed air reactor with hematite as the OC for CLC. A proper 283 particle distribution and mass flow rate of the circulating OC are necessary in order to 284 achieve stable circulating conditions. The mass flow rate and volume fraction of the 285   Figure 9 shows the comparison of experimental and numerical data of static pressure 267 at various heights. It can be observed that the graphs of static pressure are almost of the 268 same magnitude at different heights, except for the difference at the inlet. This difference 269 occurs because of the different initial processes of the simulation and experimental work. 270 In the experiment, the particles are injected from the narrow tube thoroughly. On the other 271 hand, in the simulation, the particles are considered to be injected from some height of the 272 reactor, and they are then left for 1 s to settle before the air injection. Changing the drag 273 law can also help in minimizing the deviation of numerical results from the experimental 274 result [32]. The decay in static pressure after the inlet is found to be consistent with the 275 experimental results, along with the height of the reactor. Figure 10 shows the error bar 276 graph for the difference between experimental and numerical data. It can be seen from the 277 graph that the errors are much smaller at any height, except at the inlet. 278

Conclusions 281
In this work, the DDPM model of Ansys Fluent 17 was used to carry out a cold flow 282 simulation of the fast-fluidized bed air reactor with hematite as the OC for CLC. A proper 283 particle distribution and mass flow rate of the circulating OC are necessary in order to 284 achieve stable circulating conditions. The mass flow rate and volume fraction of the 285 Figure 10. Error bar graph for experimental and numerical data.

Conclusions
In this work, the DDPM model of Ansys Fluent 17 was used to carry out a cold flow simulation of the fast-fluidized bed air reactor with hematite as the OC for CLC. A proper particle distribution and mass flow rate of the circulating OC are necessary in order to achieve stable circulating conditions. The mass flow rate and volume fraction of the particle stabilizes after 28 s at the particle outlet, and a constant pressure gradient of 1.5 kPa is maintained between the inlet and outlet of the reactor after 15 s of particle injection. The results are validated by comparing the variation in static pressure along the height of the reactor with the experimental work performed by Haider et al. [22] at Cranfield, UK, considering FE100 as the OC. On this basis, it can be concluded that the DDPM model is suitable for determining the hydrodynamics of a fast-fluidized bed reactor, and hematite can be opted as the OC for CLC due to its capacity for fast oxidation and reduction during the combustion process.