Subscriber access provided by University of Groningen
Article
Reduction kinetics of ilmenite ore for pressurized chemical looping combustion of simulated natural gas Yewen Tan, Firas Nor Ridha, Dennis Y. Lu, and Robin W Hughes Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.7b02648 • Publication Date (Web): 28 Nov 2017 Downloaded from http://pubs.acs.org on November 29, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
1
Reduction kinetics of ilmenite ore for pressurized
2
chemical looping combustion of simulated natural
3
gas
4
Yewen Tan*, Firas N. Ridha, Dennis Y. Lu, Robin W. Hughes
5
Natural Resources Canada, CanmetENERGY-Ottawa
6
1 Haanel Drive, Ottawa, Ontario, Canada K1A 1M1
7
KEYWORDS
8
chemical looping combustion; ilmenite; elevated pressure; thermogravimetric analysis; natural
9
gas
10
ABSTRACT
11
Reduction kinetics of ilmenite ore as an oxygen carrier for the chemical looping combustion of a
12
simulated natural gas mixture under elevated pressure was studied using a pressurized
13
thermogravimetric analyzer (PTGA). The fuel gas is a mixture of hydrocarbon, carbon dioxide
14
and nitrogen to simulate an actual combustion environment. The oxidation phase of the
15
experiments was carried out in air. Effects of temperature (1023 – 1223 K), total pressure (0.6 –
16
1.6 MPa), fuel partial pressure (0.126 – 0.34 MPa) as well as CO2 partial pressure *
Corresponding author,
[email protected] 1 ACS Paragon Plus Environment
Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 2 of 42
17
(pFuel:pCO2=0.5–1) were studied. The results showed that the presence of small amounts of
18
ethane and propane clearly led to a higher ilmenite reactivity at temperatures below 1123 K, but
19
this effect became less significant as temperature increased and completely disappeared above
20
1123 K. The results also showed that increasing CO2 partial pressure had little effect on ilmenite
21
conversion rate, though it did have some slightly negative influence on ilmenite oxygen carrying
22
capacity that was especially noticeable at lower total pressure. A higher fuel partial pressure
23
appeared to have a slightly negative impact on ilmenite oxygen carrying capacity, especially at
24
higher temperature. A kinetic model, based on a phase-boundary controlled mechanism with
25
contracting sphere, was developed by incorporating the total pressure, fuel and CO2 partial
26
pressures and the temperature and it was able to satisfactorily reproduce most of the test results
27
with a conversion ratio of up to 70%. This model predicted that the ilmenite conversion rate had
28
a strong positive correlation with the temperature and fuel partial pressure, and a relatively
29
weaker negative correlation with the total pressure and CO2 partial pressures. Overall conversion
30
rate will increase when total pressure increases, which justify the pressurized chemical looping
31
combustion technology.
32
1. Introduction
33
Chemical-looping combustion (CLC) is a technology for the conversion of fossil fuels with
34
inherent separation of CO2 and efficient use of energy. In a typical CLC setup, two reactors are
35
used. One of them is a fuel reactor and the other is an air reactor. A solid, called an oxygen
36
carrier (OC), is used as an oxidant for the fuel in the fuel reactor. A fuel, such as natural gas,
37
reacts with the OC in the fuel reactor, producing a flue gas composed mostly of CO2 and H2O,
38
thus making CO2 capture very effective. The reduced OC is then sent back to the air reactor
39
where it reacts with air for regeneration before being sent back to the fuel reactor. As a result, 2 ACS Paragon Plus Environment
Page 3 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
40
CLC does not require an air separation unit as the oxy-fuel combustion does nor does it need
41
expensive post-combustion CO2 capture units. According to the literature, improved efficiency
42
can be obtained by employing a CLC-integrated power generation system.1-3 CLC processes can
43
be also used for hydrogen and steam production. A comprehensive review of the CLC
44
technology and systems by Adánez et al. can be found in the literature.4
45
Recently, pressurized CLC (PCLC) technology has attracted the attention of some research
46
groups.5-10 The major advantages of pressurized CLC are: high pressure CO2 will be produced
47
from the fuel reactor thereby reducing the power and associated capital required to compress
48
CO2 to pipeline pressure (~13.8 MPa); meanwhile the latent heat of condensation of flue gas
49
moisture can be more effectively integrated into the steam cycle due to high dew point
50
temperature at high pressure; and the pressurized system provides improved hydrodynamics in
51
the fluidized bed resulting in better mixing of reactants. Other noted advantages of PCLC include
52
a more compact system for transport to remote areas, and the improved reaction rate with higher
53
reactant partial pressures.
54
CanmetENERGY-Ottawa has been studying CLC technology, including pressurized CLC, for
55
the past several years.11-14 The works have focused on using a commercial ilmenite ore as oxygen
56
carrier. Ilmenite ore is mainly composed of FeTiO3 (FeO·TiO2). It is an attractive oxygen carrier
57
due to its wide availability and low cost. There have been numerous studies on its characteristics
58
as an oxygen carrier.13, 15-19 Most of these studies were conducted using methane, CO or coal. For
59
natural gas or syngas as fuel, very few studies were done with ilmenite.20, 21 A few were done
60
with iron ore such as in the study of Wang et al.22
61
Natural gas is mainly composed of methane with some higher hydrocarbons. According to
62
Union Gas, a typical natural gas can be composed of 87-97% CH4, 1.5-7% of C2H6, 0.1-1.5%
3 ACS Paragon Plus Environment
Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 4 of 42
63
C3H8, 0.01-0.3% iso-butane and other minor components. Literature shows that these small
64
amounts of higher hydrocarbons can have significant impact on the rate of reaction of the
65
mixture compared to pure CH4 alone.23-25 It is hypothesized then that the presence of these higher
66
hydrocarbons may have an impact on the reactivity of ilmenite as oxygen carrier in chemical
67
looping combustion.
68
To the best of the authors’ knowledge, no prior research has been done to study the reduction
69
kinetics of ilmenite using natural gas or simulated natural gas, whether at atmospheric or
70
elevated pressures. The intent of this work is to fill this knowledge gap at elevated pressure. To
71
this end, the effects of total pressure, fuel and CO2 partial pressures and temperature were
72
studied and a preliminary model was derived based on the experimental data. The conversion
73
ratio, conversion rate and oxygen transport capacity of ilmenite ore were evaluated.
74
2. Materials and methods
75
Pressurized thermogravimetric analyzer (PTGA). The PTGA used in this work was
76
manufactured by Linseis GmbH and is described elsewhere.13 This facility operates from
77
atmospheric pressure up to 10.1 MPa at temperatures of up to 1873 K. The reaction gases,
78
provided by the gas cylinders, are mixed with steam before entering the PTGA from the top. The
79
balance itself is purged with nitrogen during the operation. The PTGA is connected to a data
80
acquisition system controlled by a computer. The test data are collected at 1 s intervals. A
81
schematic of the facility is given in Figure 1.
82
Figure 1. A schematic diagram of the PTGA
4 ACS Paragon Plus Environment
Page 5 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
83 84 85
Energy & Fuels
Sample preparation and fuel mixture. The fuel mixture used for this work was provided by Linde with a certified gas composition of 95.02% CH4, 3.96% C2H6 and 1.02% C3H8. The oxygen carrier (OC) used in this work is UKTO ilmenite ore provided by Rio Tinto
86
Canada. It is an upgraded roasted ilmenite ore (concentrate produced from raw massive ilmenite
87
ore after removing most of its gangue material by gravity separation and sulphur by roasting at
88
1073 K. 11 The material was crushed and sieved to 106-212 µm and was then calcined in a muffle
89
furnace for 2 hours at 1173 K before it was used in the PTGA tests. Its composition, analyzed by
90
X-ray fluorescence, is presented in Table 1.
91
Table 1. UKTO ilmenite ore composition (wt.%)
92
Test procedure. The test procedure has been described in detail in a previous publication13
93
and is only briefly described here. For each test, 95-100 mg of OC sample was evaluated and
94
each test consisted of 8 redox cycles. Our previous study using methane has shown that, after an
95
activation period of 3 redox cycles under elevated pressures, the reactivity of ilmenite stayed
96
constant for at least 40 cycles. Nitrogen was used to pressurize the PTGA to the test pressure and
97
as purge gas. The total flow rate was 0.5 normal L/min. Previous work13 has shown that this flow
98
rate was enough to eliminate mass transfer effects.
99
Data repeatability of the PTGA was very good as can be seen from the previous study.
100
3. Reactions of ilmenite with hydrocarbons
101
As stated above, most of the previous studies focused on using pure methane as fuel. The
102
main reactions between methane and ilmenite ore are those involving pseudobrookite and
103
hematite and can be described by the following overall reactions:18
104
4Fe2TiO5+4TiO2+CH4 → 8FeTiO3+CO2+2H2O
(R1)
105
12Fe2O3+CH4 → 8Fe3O4+CO2+2H2O
(R2)
5 ACS Paragon Plus Environment
Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
106
Page 6 of 42
The oxidation stage can be described by the following reactions:18
107
4FeTiO3+O2 → 2Fe2TiO5+2TiO2
(R3)
108
4Fe3O4+O2 → 6Fe2O3
(R4)
109
With the higher hydrocarbons under consideration here, we can write the reactions as below.
110
7Fe2TiO5+7TiO2+C2H6 → 14FeTiO3+2CO2+3H2O
(R5)
111
21Fe2O3+C2H6 → 14Fe3O4+2CO2+3H2O
(R6)
112
10Fe2TiO5+10TiO2+C3H8 → 20FeTiO3+3CO2+4H2O
(R7)
113
30Fe2O3+C3H8 → 20Fe3O4+3CO2+4H2O
(R8)
114
However, it should be noted that ilmenite ore reactions in CLC are in fact more complex than
115
the above reactions. This is because iron in the ilmenite ore can exist in several different oxide
116
forms and combine with titanium in various ways.26
117
In the literature, the reactivity of oxygen carriers is typically assessed using the conversion
118
ratio Xr, which is expressed as equation (1):
119
ܺ =
120 121 122
ି
(1)
ିೝ
Here, mo is the mass of the ilmenite ore sample after oxidation; mt is the mass at time t and mr is the sample mass at the end of the reduction stage. Conversion rate Ẋr is also used to assess the reactivity of oxygen carriers and is calculated
123
using the equation (2):
124
ௗ ܺሶ = ௗ௧ೝ
125
where t is time in seconds.
(2)
126
Since the conversion ratio, Xr is a normalized term from which the conversion rate Ẋr can be
127
determined, it is important to note that it is possible to misinterpret OC performance test results
128
because an oxygen carrier with a very low oxygen carrying capacity may appear to have a faster 6 ACS Paragon Plus Environment
Page 7 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
129
reaction rate than one with a high oxygen carrying capacity. To compensate for this, oxygen
130
carrying capacity Ro is used in combination with Ẋr in evaluating oxygen carrier performance. A
131
better oxygen carrier would be the one with a higher Ẋr as well as a higher Ro. The oxygen
132
transport capacity Ro, in percentage, can be calculated using equation (3):
133
ܴ =
134
ିೝ బ
× 100%
(3)
In this paper, ilmenite reactivity is often presented with a weight change (%) vs. time curve as
135
we feel that this can lead to easier data interpretation and comparison. The weight change C in
136
percentage, is calculated using equation (4):
137
=ܥ
138
ି బ
× 100%
(4)
With this formula, the weight change is negative during the reduction stage and positive
139
during the oxidation stage. Equation (4) shows that the absolute value of C is the same as that of
140
oxygen carrying capacity at the end of the reduction stage.
141
Carbon deposition is a common concern in chemical looping combustion with hydrocarbons.
142
Our previous study with methane has shown that, under the test conditions in Table 2, carbon
143
deposition was not observed.13 The fuel used in this work was a mixture of methane, ethane and
144
propane. In theory, ethane and propane are more likely to lead to carbon deposition because of
145
their higher carbon content. However, due to their low concentrations in the mixture, carbon
146
deposition was not observed in any of the tests presented here.
147
4. Test results and discussion
148
Table 2 summarized the test conditions where total pressure and partial pressures are given as
149
absolute pressure. Nitrogen was used as balance gas. Carbon dioxide was used to simulate a
150
realistic operating environment. In commercial operations, it may not be practical to use pure
151
natural gas to fluidize the fuel reactor, so a recycled flue gas (mainly CO2) or a mixture of CO2 7 ACS Paragon Plus Environment
Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 8 of 42
152
and steam may be used. For this reason, the partial pressure of CO2 varies for some tests in order
153
to evaluate its effects on oxygen carrier reactivity. The test data from the 8th redox cycle was
154
used in results evaluation.
155 156
Table 2. Test conditions for PTGA experiments Under elevated pressure conditions, gas dispersion in the gas supply system of the
157
experimental device becomes pronounced and results in the sample being exposed to a gas
158
mixture of varying composition when passing from one stage to the next (i.e. oxidation to
159
reduction and vice versa). This effect of gas dispersion can be seen in Figure 2 showing oxygen
160
carrier reduction for tests T915a and T505a. Gas dispersion effect becomes more pronounced
161
when total pressure increases. This phenomenon has been observed also by other researchers.7
162
Increasing gas flow rate can reduce the gas dispersion effect but cannot eliminate it. Our tests
163
showed that increasing gas flow rate from 0.5 to 1 normal L/min had no influence on ilmenite
164
reactivity and only marginally reduced the gas dispersion effect.
165
Figure 2. Ilmenite weight change (T=1223 K, pHC=0.34 MPa, pFuel:pCO2=1) for tests T915a
166
(dashed line) and T505a (solid line)
167
In order to assess the conversion ratio Xr and the conversion rate Ẋr across different total
168
pressures, the data in period of pronounced gas dispersion effect was removed from the data sets.
169
This was achieved by assessing the minimum value of the derivative
170
after that point in determining the reaction kinetics. On the other hand, the criteria
171
used in assessing Ro since the oxygen carrier was reacting with fuel while the fuel partial
172
pressure was still rising during the gas dispersion period.
173 174
ௗమ ௗ௧ మ
and then selecting data ௗ ௗ௧
< 0 was
Effect of temperature and total pressure. Figure 3 provides plots of ilmenite conversion ratio, Xr, for the oxidation of the CH4-C2H6-C3H8 mixture as a function of time at temperatures 8 ACS Paragon Plus Environment
Page 9 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
175
from 1023 to 1223 K, at total pressure of 0.9 MPa and fuel-CO2 partial pressure ratio
176
pFuel:pCO2=1, where pFuel= pCH4+ pC2H6 + pC3H8.. As such, the curves represent the
177
conversion ratio at each temperature as a function of time.
178
Figure 3. Conversion ratio at 0.9 MPa for pFuel:pCO2=1. Tests T511a, T531a, T513a, T504a,
179
T503a and T505a
180
Figure 3 shows that, up to 70% conversion, the rate of reduction of ilmenite increases with
181
increasing temperature from 1023 to 1173 K. However, increasing temperature from 1173 K to
182
1223 K did not increase the conversion rate. In fact, increasing temperature beyond 1173 K led
183
to a decrease in oxygen carrying capacity, Ro, by ~6%, as can be seen from Figure 4.
184
Figure 4. Oxygen carrying capacity at 0.9 MPa for pFuel:pCO2=1. Tests T511a, T531a, T513a,
185
T504a, T503a and T505a
186
Even though this paper focuses on the reduction of ilmenite, it is interesting to see what
187
happens at 1223 K during the oxidation stage as it can help us to interpret our observations
188
during the reduction period. The negative effect of high temperature on ilmenite reduction
189
conversion rate and capacity, i.e. the effect noted in the previous paragraph at 1223 K, was
190
observed in a more pronounced way during the re-oxidation of the sample where the oxidation of
191
the reduced ilmenite was considerably slower at 1223 K as shown in Figure 5. At the 120 s mark,
192
the total weight change (%) of the reduced ilmenite sample at 1223 K was comparable to that at
193
1023 K, even though its oxidation rate was comparable to that at 1173 K in the first 10 s. It
194
should be noted that, over an extended period of oxidation, the ilmenite ore was still fully
195
oxidized at 1223 K with a weight gain of ~4.3%.
196
Figure 5. Weight change at 0.9 MPa for oxidation stage. Tests T511a, T531a, T513a, T504a,
197
T503a and T505a
9 ACS Paragon Plus Environment
Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
198
Page 10 of 42
SEM images for the reduced samples of tests T503b (T=1173 K) and T505b (T=1223 K) are
199
presented in Figure 6. These images show that the surface morphology of the oxygen carrier
200
underwent significant changes when temperature rose from 1173 to 1223 K. At 1173 K, ilmenite
201
particles preserved their structure with no major cracks or defects apparent on the surface, as
202
shown in image 6a. At higher magnification, as shown in image 6b, the surface seems to exhibit
203
a granular morphology with grain size ranges from 1 to 5 µm and pores apparent on the surface.
204
When the temperature was increased to 1223 K, the surface morphology of the particle
205
experienced significant changes. Large agglomerates, with size of ~500 µm, were formed as
206
shown in image 6c. The reason for agglomerates formation is the development of fused regions
207
on the surface of particles that tend to form at points of contact. Fused regions are shown in
208
image 6d that presents a magnified image of the surface of particles in image 6c. Furthermore, a
209
closer examination of this image reveals the disappearance of grain boundaries as they merged
210
into each other. Such agglomeration was noted to take place during and after oxidation of deeply
211
reduced particles.
212
The above discussion is consistent with our previous study from this group that showed that
213
larger grains were formed with increasing reaction temperature.11 The same study also showed
214
that total pressure did not have any noticeable effect on surface morphology.
215 216
Figure 6. SEM images for test T503b (T=1173 K, top) and T505b (T=1223 K, bottom) b
The effect of total pressure at 0.9 MPa (T505a) and 1.6 MPa (T915a) can be seen from Figure
217
2. It shows that, as has been previously demonstrated in the literature,7, 12 increasing total
218
pressure while maintaining constant fuel partial pressure led to a clear decrease in the rate of
219
weight loss of ilmenite ore. This figure also shows that no significant effect on the oxygen
220
carrying capacity of ilmenite ore was observed with the increase of total pressure as the final c
d 10
ACS Paragon Plus Environment
Page 11 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
221
weight loss was the same at both pressures. It is important to point out here that, even though
222
total pressure by itself is a negative factor in terms of reaction kinetics, with increasing total
223
pressure in practical operation, fuel partial pressure can be increased, which will lead to a higher
224
reaction rate as shown later in the modeling section of this paper. This, combined with other
225
benefits associated with high-pressure systems, makes PCLC an overall attractive option.
226
Comparison with pure methane. Tan et al.25 showed that, in a jet-stirred reactor,
227
combustion kinetics of methane could be significantly accelerated when a small amount of
228
ethane and propane was added. It was thus reasonable to hypothesize that, when small amounts
229
of ethane and propane were present, the kinetics of the hydrocarbon mixture could also be
230
accelerated. In the TGA studies for chemical looping combustion scenarios, this would manifest
231
itself as improved reactivity of the oxygen carrier. This is shown to be true in Figure 7. At 1023
232
K, the reactivity of ilmenite for the oxidation of the hydrocarbon mixture was significantly
233
higher than for pure methane, and its oxygen carrying capacity was also much higher for the
234
mixture. For example, after 50 sec, ilmenite weight loss with the simulated natural gas mixture
235
was ~2% compared to ~1.1% for pure methane. Not surprisingly, the reactivity of ilmenite was
236
the highest for pure ethane. At 50 sec, the ilmenite weight loss with ethane was at ~3.7%. The
237
final ilmenite weight loss was ~3.2% for the simulated natural gas mixture and ~2.3% for pure
238
methane. A much deeper reduction, over 5% weight loss, was obtained when pure ethane was
239
used.
240
Figure 7. Ilmenite weight loss for pure methane, the simulated natural gas mixture and pure
241
ethane. Tests T511a, T512 and T628. T=1023 K, P=0.9 MPa, pHC 0.34 MPa, pFuel:pCO2=1
11 ACS Paragon Plus Environment
Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
242
Page 12 of 42
As temperature rose, the difference between the simulated natural gas mixture and pure
243
methane became gradually smaller and at 1173 K, no difference between the two was observed
244
as can be seen in Figure 8.
245
Figure 8. Ilmenite weight loss for pure methane and simulated natural gas mixture. Tests T503a
246
(solid line) and T922 (dotted line). T=1173 K, P=0.9 MPa, pHC=0.34 MPa, pFuel:pCO2=1
247
The first step in high temperature alkane oxidation is unimolecular breakdown of the parent
248
molecule. The rate constant for methane decomposition is considerably smaller than that of
249
ethane and propane. Because of this, in a methane-ethane-propane mixture, the initial radical
250
production will be exclusively from the breakdown of ethane and propane at a rate orders of
251
magnitudes faster than that of methane.
252
According to Tan et al.,25 when ethane and propane are added to methane, the higher
253
hydrocarbons start to decompose first,
254
C2H6 ↔ CH3 + CH3
(R9)
255
C3H8 ↔ CH3 + C2H5
(R10)
256
Methyl radicals then mainly react with molecular oxygen to produce OH radicals while ethyl
257
radicals mainly decompose and become major H-atom source,
258
CH3 + O2 ↔ CH2O + OH
(R11)
259
C2H5 ↔ C2H4 + H
(R12)
260
H + O2 ↔ OH + O
(R13)
261
Finally, H and OH radicals react with methane, increasing the radical pool,
262
OH + CH4 ↔ H2O + CH3
(R14)
263
H + CH4 ↔ H2 + CH3
(R15)
12 ACS Paragon Plus Environment
Page 13 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
264
Energy & Fuels
As temperature rises, methane becomes increasingly less stable and starts to decompose more
265
easily, so that the difference between the simulated natural gas mixture and pure methane
266
becomes increasingly smaller and eventually disappears.
267
While the above mechanism was well established for combustion of hydrocarbons in gaseous
268
environment, it is likely also valid for chemical looping combustion and can be invoked to
269
explain the increased ilmenite reactivity for the combustion of simulated natural gas mixtures
270
compared to pure methane. In this case, due to the nature of gas-solid reactions, reactions
271
involving oxygen must be rewritten. For example, R11 can be rewritten as
272
CH3 + 2FeTiO5 + 2TiO2 ↔ CH2O + OH + 4FeTiO3
(R11a)
273
CH3 + 6Fe2O3 ↔ CH2O + OH + 4Fe3O4
(R11b)
274
For R13, it is more complicated as an extra O radical is produced when H reacts with the
275
ilmenite:
276
H + FeTiO5 + TiO2 ↔ OH + 2O + FeTiO3
277 278
(R13a)
On the other hand, one less O radical is produced by reaction R13b: H + 3Fe2O3 ↔ OH + 2Fe3O4
(R13b)
279
The exact elementary reaction mechanism of hydrocarbons reacting with ilmenite oxygen
280
carrier is unknown at this time. This will require more fundamental research. However, the lack
281
of understanding of this mechanism does not affect this work, as we are only interested in the
282
overall kinetics of reactions between ilmenite and hydrocarbons.
283
The negative effect of high pressure on hydrocarbon oxidation is consistent with our previous
284
findings12, 13 as well as findings of other researchers7 and it was explained before using other gas-
285
solid systems. From a kinetics perspective, it is well known that radical recombination reactions
286
become increasingly important as system pressure increases. Reactions such as H + O2 + M ↔
13 ACS Paragon Plus Environment
Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 14 of 42
287
HO2 + M are increasingly favored in the recombination direction, thus reducing concentrations
288
of the hyperactive H atoms, leading to slower and steadier reactions. Dagaut and Cathonnet
289
argued27 that, under high pressure conditions, the oxidation route of ethyl radical oxidation
290
(producing C2H3 + O2) started to compete with the thermal decomposition route (forming H2).
291
Thus, the net production of H atoms decreases reducing the rate of the main hydrogen production
292
channel. Considering the presence of ethane in the fuel mixture, this is a plausible explanation.
293
The above discussion shows that the radical formation theories based on gas phase chemical
294
kinetics could be used to explain some of the observations seen in this study and that it could be
295
interesting to pursue research in this area to understand interactions between the gas phase fuel
296
and the solid phase ilmenite.
297
Effect of CO2 partial pressure. Most of the practical chemical looping systems have been
298
based on dual fluidized bed concepts. In such configurations, a fluidizing gas is needed to keep
299
the solids looping between two reactors to maximize combustion efficiency. For the fuel reactor,
300
in order to produce a flue gas composed of only CO2 and H2O to avoid any need for flue gas
301
separation, the options for fluidizing gas will be limited to either fuel gas or recycled flue gas
302
(mostly CO2) or a mixture of the two. It is therefore necessary to understand the effect of CO2
303
partial pressure on ilmenite reactivity.
304
To this end, several tests were done with different fuel-CO2 partial pressure ratios:
305
pFuel:pCO2 =1 and pFuel:pCO2 = 0.5, where pFuel = pCH4+ pC2H6 + pC3H8. The results are
306
shown in Figures 9-11.
307
Figure 9. Weight loss curves at T=1173 K, P=0.6 MPa, pHC=0.126 MPa. T920a, pFuel:pCO2=1
308
(solid line) and T920b, pFuel:pCO2=0.5 (dotted line)
14 ACS Paragon Plus Environment
Page 15 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
309
Figure 10. Weight loss curves at T=1173 K, P=0.9 MPa, pHC=0.126 MPa. T923a,
310
pFuel:pCO2=1 (solid line) and T923b, pFuel:pCO2=0.5 (dotted line)
311
Figure 11. Weight change at T=1173 K, P=1.6 MPa, pHC=0.126 MPa. T825a, pFuel:pCO2=1
312
(solid line) and T825b, pFuel:pCO2=0.5 (dotted line)
313
These figures, obtained for tests done at 1173 K, show that increasing CO2 partial pressure
314
had a slightly negative effect on the conversion rate of ilmenite from 0.6 to 1.6 MPa. Increasing
315
CO2 partial pressure also led to lower oxygen carrying capacity. At 0.6 MPa at 140 s, doubling
316
CO2 partial pressure led to a decrease of ~13.7% in ilmenite oxygen carrying capacity; at 0.9
317
MPa at 200 s, doubling CO2 partial pressure led to a decrease in oxygen carrying capacity by
318
~7.8% and at 1.6 MPa at 150 s, doubling CO2 partial pressure led to a decrease of about ~5%. It
319
can thus be seen that the effect of CO2 on ilmenite oxygen carrying capacity diminished as total
320
pressure increased. These results show that CO2 partial pressure is a factor that should be taken
321
into account for fuel reactor design, especially at lower total pressure conditions. The effect of
322
the observed CO2 partial pressure can be explained by the fact that, as CO2 is a product of
323
reactions R1, R2, R5-R8, when CO2 partial pressure increased, the equilibrium of these reactions
324
was favored toward their left. At higher total pressures, the increased fuel partial pressure was
325
able to compensate for the negative effects of pCO2 as the model shows in a later section.
326
Effect of fuel partial pressure. The effect of hydrocarbon partial pressure is shown in
327
Figures 12-14. As expected, higher fuel partial pressure led to higher rate of ilmenite reduction.
328
However, it is interesting to note that oxygen carrying capacity was slightly higher at lower fuel
329
partial pressure, especially as temperature increased. For example, at 1123 K, ilmenite weight
330
loss was ~9% higher at fuel partial pressure of 0.126 MPa than at 0.34 MPa; at 1173 K, the
331
weight loss was ~10% higher at fuel partial pressure of 0.126 MPa than at 0.34 MPa. At 1223 K,
15 ACS Paragon Plus Environment
Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 16 of 42
332
the difference was ~13%. Probably, at higher fuel partial pressure, the reaction proceeded so fast
333
that the active sites of the oxygen carrier surface were quickly covered so that further reaction
334
was hindered. With higher temperature, the reaction proceeds even faster so this effect became
335
more pronounced.
336
Figure 12. Weight loss curves for tests T927a and T504a. T=1123 K, P=0.9 MPa, pFuel:pCO2=1.
337
pHC=0.34 MPa (solid line) and pHC=0.126 MPa (dotted line)
338
Figure 13. Weight loss curves for tests T923a and T503a. T=1173 K, P=0.9 MPa, pFuel:pCO2=1.
339
pHC=0.34 MPa (solid line) and pHC=0.126 MPa (dotted line)
340
Figure 14. Weight loss curves for tests T926a and T505a. T=1223 K, P=0.9 MPa, pFuel:pCO2=1.
341
pHC=0.34 MPa (solid line) and pHC=0.126 MPa (dotted line)
342
5. Kinetic modeling
343
The approach for kinetic modeling of the experimental data was described in detail in a
344
previous publication.13 Essentially, the Hancock and Sharp method28 was used to establish
345
kinetic models for ilmenite ore reduction with methane. This approach has also been successfully
346
used by several research groups to model ilmenite ore redox at atmospheric pressure
347
conditions.29, 30
348 349
Following Hancock and Sharp’s method, ln(-ln(1-Xr)) against ln(t) is plotted, where t is the time of reaction in seconds, and a linear equation is fitted to the data.
350
In Hancock and Sharp’s model, the phase-boundary controlled mechanism with contracting
351
sphere, which is the adopted model for this work, could be described with the following equation:
352
݇ × = ݐ1 − (1 − ܺ )య
భ
(5)
16 ACS Paragon Plus Environment
Page 17 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
353
where t is reaction time (s), Xr is ilmenite ore’s conversion ratio and k is a variable determined by
354
total pressure, fuel partial pressure, pCO2 and temperature. In the Arrhenius form, k can be
355
expressed as
356
݇ = (ݔ݁ܣ− ோ்)
357
where A is the pre-exponential factor (s-1), E the activation energy (J mol-1 K-1), T the reaction
358
temperature (K) and R the universal gas constant (J mol-1 K-1).
359 360 361 362
ா
(6)
The conversion ratio, Xr, can then be calculated as ܺ = 1 − (1 − ݇ × )ݐଷ
(7)
By plotting 1-(1-Xr)1/3 against time t using experiments data, it is possible to derive k from the slope of the resulting line as shown in Figure 15.
363
Figure 15. Data (test T504a) fitting to obtain k. T=1123 K, P=0.9 MPa, pHC=0.34 MPa,
364
pFuel:pCO2=1
365
From equation (6), the activation energy can be calculated using the experimental data. In the
366
temperature range from 1023 – 1173 K, at total pressure of 0.9 MPa and partial pressure pCH4 of
367
0.34 MPa, the activation energy was 69 kJ mol-1 K-1 for pure methane reduction and 56 kJ mol-1
368
K-1 for the hydrocarbon mixture reduction.
369
Compared to our previous work,13 the activation energy for pure methane is quite a bit higher
370
in this work, 69 kJ mol-1 K-1 instead of 28 kJ mol-1 K-1. This is because we extended the
371
temperature range from the previous 1123 K-1173 K to 1023 K-1173 K with a lower temperature
372
boundary. This indicates the significant effect of temperature on ilmenite activation energy,
373
which in fact can be expected as literature has already shown that the final products of ilmenite
374
oxidation depend greatly on the temperature. It is generally acknowledged that complete
17 ACS Paragon Plus Environment
Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 18 of 42
375
oxidation to pseudobrookite (Fe2TiO5) occurs above 1123 K.18 Below this temperature, a mixture
376
of different oxidized compounds can exist in the product.
377
With the data from this work, in the temperature range of 1123 – 1173 K where complete
378
oxidation of ilmenite can be expected, the activation energy obtained for the hydrocarbon
379
mixture was 31 kJ mol-1 K-1, close to the 28 kJ mol-1 K-1 obtained for pure methane. On the other
380
hand, in the lower temperature range of 1023 – 1123 K, the activation energy for the
381
hydrocarbon mixture reduction was 67 kJ mol-1 K-1 while that for pure methane reduction was
382
much higher at 90 kJ mol-1 K-1.
383
In order to easily evaluate the effect of various parameters studied in this work on the kinetics
384
of ilmenite for hydrocarbon reduction, k was fitted using the experimental results obtained under
385
various conditions and the following expression for k was obtained:
386
݇ = 5.263 × 10ିଶହ × ܲି.ସଶହ × ݈݁ݑܨ.ଽଷହ × ܱܥଶ ି.ସ × ܶ .ଷ
387
(8)
Here, the pressure is in MPa(a) and the temperature T in K; P is total pressure, pCO2 is CO2
388
partial pressure and pFuel is fuel partial pressure. This expression applies to a total pressure
389
range between 0.9 and 1.6 MPa(a). Equation (8) shows that temperature T has by far the most
390
significant and positive influence on conversion rate and that k has an almost linear relationship
391
with fuel partial pressure. It also shows that the negative impact of pCO2 on k was even stronger
392
than the total pressure P. Further, equation (8) reveals that the negative effects of increasing P
393
and pCO2 could be compensated for by the positive effect of increasing pFuel. For example, for
394
a reaction gas mixture of 20 vol. % fuel and 20 vol. % CO2, when total pressure P increases from
395
0.9 to 1.6 MPa(a), k will decrease by ~21% due to increased P and ~23% due to increased pCO2,
396
but will increase by 70% due to increased pFuel. Overall conversion rate k will increase when P
18 ACS Paragon Plus Environment
Page 19 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
397
increases, as can be seen in Figure 16. This justifies our research in improving reaction rates by
398
pressurizing chemical looping combustion.
399 400
Figure 16. Individual effects on k by total pressure P, fuel partial pressure pFuel and CO2 partial pressure pCO2, assuming 20 vol.% fuel and 20 vol.% CO2 in the reaction gas mixture
401
Using equations (7) and (8), most of the experimental results can be satisfactorily reproduced
402
with up to 70% conversion, for some as high as 90% conversion, as shown in Figure 17. The
403
model does have a tendency to overestimate the negative effect of pCO2 by underestimating the
404
conversion ratio at a higher pCO2 (T919a) when conversion ratio was above 50%. Overall, the
405
model proposed in equation (8) has a strong confidence level for temperature, as well as for total
406
pressure and fuel partial pressure, but a relatively lower confidence level for CO2 partial pressure.
407
Figure 17. Comparison between model predictions (solid lines) and experimental data (symbols).
408
o: T825a; ∆: T919a; +: T503a; *: T927a
409
6. Conclusions
410
Reduction kinetics of a commercial natural ilmenite ore were studied using a pressurized
411
thermogravimetric analyzer with a simulated natural gas as fuel. The results showed that,
412
compared to reaction with pure methane, the reduction kinetics of ilmenite ore were considerably
413
improved with the hydrocarbon mixture for temperatures below 1123 K, but no significant
414
difference was observed for temperatures above 1123 K. The higher ilmenite reactivity at lower
415
temperature was attributed to the presence of ethane and propane in the hydrocarbon mixture
416
because these higher hydrocarbons are easier to decompose than methane. The results also
417
showed that increasing CO2 partial pressure had a minor negative effect on ilmenite reduction
418
kinetics but a more pronounced negative effect on its oxygen carrying capacity and this effect
419
diminishes at higher total pressure. While increasing hydrocarbon partial pressure was beneficial
420
to ilmenite reduction kinetics, it had slight negative impact on its oxygen carrying capacity, 19 ACS Paragon Plus Environment
Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 20 of 42
421
especially at higher temperature. The conclusion of this paper would draw is that the presence of
422
small amounts of higher hydrocarbons in the natural gas can have significant effect on the
423
reduction kinetics of an oxygen carrier at lower temperature so that they should be taken into
424
account when designing fuel reactors for such conditions. Finally, a kinetic model was developed
425
as a function of the total pressure (0.9-1.6 MPa(a)), fuel and CO2 partial pressures and the
426
temperature. The model showed that the ilmenite conversion rate had a strong positive
427
correlation with the temperature and fuel partial pressure, and a relatively weaker negative
428
correlation with the total pressure and CO2 partial pressures.
429
430
The authors wish to acknowledge funding from Natural Resources Canada’s Program of
431 432
Acknowledgments
Energy and Research Development.
References
433
(1) Ishida M, Zheng D, Akehata T. Evaluation of a chemical-looping combustion power-generation system by
434
graphic energy analysis. Energy 1987;12(2):147-154.
435
(2) Anheden M, Svedberg G. Exergy analysis of chemical-looping combustion systems. Energy Conversion &
436
Management 1998;39(16-18):1967–1980.
437
(3) Brandvoll O, Bolland O. Inherent CO2 capture using chemical looping combustion in a natural gas fired power
438
cycle. J. Eng. Gas Turbine Power 2004; 126(2):316–321.
439
(4) Adánez J, Abad A, García-Labiano F, Gayán P, de Diego LF. Progress in chemical-looping combustion and
440
reforming technologies. Prog. Energy & Comb. Sci. 2012; 38(2):215-282.
441
(5) Xiao R, Song Q, Song M, Lu Z, Zhang S, Shen L. Pressurized chemical-looping combustion of coal with an iron
442
ore-based oxygen carrier, Comb. & Flame 2010: 157(6):1140-1153.
443
(6) Xiao R, Chen L, Saha C, Zhang S, Bhattacharya S. Pressurized chemical-looping combustion of coal using an
444
iron ore as oxygen carrier in a pilot-scale unit. Int. J. of Greenhouse Gas Control 2012; 10:363-373.
445
(7) García-Labiano F, Adánez J, de Diego LF, Gayán P, Abad A. Effect of pressure on the behavior of copper-, iron-,
446
and nickel-based oxygen carriers for chemical looping combustion. Energy Fuels 2006; 20(1):26-33.
20 ACS Paragon Plus Environment
Page 21 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
447
(8) Xiao R, Song Q, Zhang S, Zheng W, Yang Y. Pressurized chemical-looping combustion of Chinese bituminous
448
coal: cyclic performance and characterization of iron ore-based oxygen carrier. Energy Fuels 2010; 24(2):1449-1463.
449
(9) Zhang S, Saha C, Yang Y, Bhattacharya S, Xiao R. Use of Fe2O3-containing industrial wastes as the oxygen
450
carrier for chemical-looping combustion of coal: effects of pressure and cycles. Energy & Fuels 2011; 25(10):4357-
451
4366.
452
(10) Zhang S, Xiao R, Zheng, W. Investigation of chemical-looping combustion of coal with iron ore as oxygen
453
carrier under pressurized condition. In: CFB-11: Proceedings of the 11th International Conference on Fluidized Bed
454
Technology 2014:855-860.
455
(11) Ridha FN, Duchesne MA, Lu, X, Lu DY, Filippou D, Hughes RW. Characterization of an ilmenite ore for
456
pressurized chemical looping combustion. Applied Energy 2016; 163:323-333.
457
(12) Lu X, Rahman RA, Lu DY, Ridha FN, Duchesne MA, Tan Y, Hughes RW. Pressurized chemical looping
458
combustion with CO: reduction reactivity and oxygen-transport capacity of ilmenite ore, Applied Energy 2016; 184:
459
132-139.
460
(13) Tan Y, Ridha FN, Duchesne MA, Lu DY, Hughes RW. Reduction kinetics of ilmenite ore as an oxugen carrier
461
for pressurized chemical looping combustion of methane. Energy&fuels, DOI:10.102 1/acs.energyfuels.7b01038.
462
(14) Ashrafi O, Navarri P, Hughes R, Lu D. Heat recovery optimization in a steam-assisted gavity drainage (SAGD)
463
plant. Energy 2016; 111:981-990.
464
(15) Miller DD, Siriwardane R, Poston J. Fluidized-bed and fixed-bed reactor testing of methane chemical looping
465
combustion with MgO-promoted hematite, Appl. Energy 2015; 146:111-121.
466
(16) Cuadrat A, Abad A, Adánez J, de Diego LF, García-Labiano F, Gayán P. Behavior of ilmenite as oxygen
467
carrier in chemical looping combustion. Fuel Processing Technology 2012; 94(1):101-112.
468
(17) Campos DC, Belkouch J, Hazi M, Ould-Dris A. Reactivity investigation on iron-titanium oxides for a moving
469
bed chemical looping combustion implementation. Adv. Chem. Eng. & Sci. 2013; 3(1):47-56.
470
(18) Abad A, Adánez J, Cuadrat A, García-Labiano F, Gayán P, de Diego LF. Kinetics of redox reactions of
471
ilmenite for chemical-looping combustion. Chem. Eng. Sci. 2011; 66(4):689-702.
472
(19) Bao J, Li Z, Cai N. Reduction kinetics of foreign-ion-promoted ilmenite using carbon monoxide (CO) for
473
chemical looping combustion. I&EC Research 2013; 52(31):10646-10655.
21 ACS Paragon Plus Environment
Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 22 of 42
474
(20) Pröll T, Mayer K, Bolhàr-Nordenkampf J, Kolbitsch P, Mattisson T, Lyngfelt A, Hofbauer H. Natural minerals
475
as oxygen carriers for chemical looping combustion in a dual circulating fluidized bed system. Energy Procedia
476
2009; 1(1):27-34.
477
(21) Bidwe AR, Hawthorne FM, Charitos A, Schuster A, Scheffknecht. Use of ilmenite as an oxygen carrier in
478
chemical looping combustion-batch and continuous dual fluidized bed investigation. Energy Procedia 2011; 4:433-
479
440.
480
(22) Wang B, Yan R, Lee DH, Liang DT, Zheng Y, Zhao H, Zheng C. Thermodynamic investigation of carbon
481
deposition and sulfur evolution in chemical looping combustion with syngas. Energy & Fuels. 2008;22:1012-1020.
482
(23) Lifshitz A, Scheller K, Burcat A, Skinner B. Shock-tube investigation of ignition in methane-oxygen-argon
483
mixtures. Combust. Flame 1971; 16(3):311-321.
484
(24) Crossley RW, Dorko EA, Scheller K, Burcat A. The effect of higher alkanes on the ignition of methane-
485
oxygen-argon mixtures in shock waves. Combust. Flame 1972; 19(3):373-378.
486
(25) Tan Y., Dagaut P., Cathonnet M, Boettner J.C. Natural gas and blends oxidation and ignition: experiments and
487
modeling. In: Symp. (Int.) on Combustion 1994; 25(1):1563-1569.
488
(26) Den Hoed P, Luckos A. Oxidation and reduction of iron-titanium oxides in chemical looping combustion: A
489
phase-chemical description. Oil Gas Sci. Technol. 2011; 66(2):249-263.
490
(27) Dagaut P, Cathonnet M. Kinetics of ethane oxidation in a high pressure jet-stirred reactor: experimental results.
491
J. Chim. Phys. 1990:87:1173-1185.
492
(28) Hancock JD, Sharp JH. Method of comparing solid-state kinetic data and its application to the decomposition of
493
kaolinite, brucite, and BaCO3. J. Am. Ceram. Soc. 1972; 55(2):74-77.
494
(29) Nasr S, Plucknett KP. Kinetics of iron ore reduction by methane for chemical looping combustion. Energy
495
Fuels 2014; 28(2):1387-1395.
496
(30) Jeong MH, Lee DH, Bae JW. Reduction and oxidation kinetics of different phases of iron oxides. Int. J.
497
Hydrogen Energy 2015; 40(6):2613-2620.
498
22 ACS Paragon Plus Environment
Page 23 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
499
Table 1. UKTO ilmenite ore composition (wt.%)
500
Al2O3
Fe2O3
TiO2
MgO
Others
0.36
60.90
36.50
1.72
0.52
501 502
23 ACS Paragon Plus Environment
Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 24 of 42
Table 2. Test conditions for PTGA experiments
503
Total
Partial
Partial
Partial
Partial
Partial
Pres.
Pres.
Pres.
Pres.
Pres.
Pres.
Test
P,
pCH4,
pC2H6,
pC3H8,
pFuel,
pCO2,
ID
MPa
MPa
MPa
MPa
MPa
MPa
TemperatureK
T920a
0.6
0.120
0.005
0.001
0.126
0.126
1173
T920b
0.6
0.120
0.005
0.001
0.126
0.252
1173
T923a
0.9
0.120
0.005
0.001
0.126
0.126
1173
T923b
0.9
0.120
0.005
0.001
0.126
0.252
1173
T926a
0.9
0.120
0.005
0.001
0.126
0.126
1223
T927a
0.9
0.120
0.005
0.001
0.126
0.126
1123
T922
0.9
0.340
0.000
0.000
0.340
0.340
1173
T503a
0.9
0.320
0.013
0.003
0.336
0.336
1173
T503b
0.9
0.320
0.013
0.003
0.336
0.560
1173
T504a
0.9
0.320
0.013
0.003
0.336
0.336
1123
T505a
0.9
0.320
0.013
0.003
0.336
0.336
1223
T505b
0.9
0.320
0.013
0.003
0.336
0.560
1223
T511a
0.9
0.320
0.013
0.003
0.336
0.336
1023
T512
0.9
0.340
0.000
0.000
0.340
0.340
1023
T513a
0.9
0.340
0.000
0.000
0.340
0.340
1073
T531a
0.9
0.320
0.013
0.003
0.336
0.336
1048
T628
0.9
0.000
0.34
0.00
0.340
0.340
1023
T825a
1.6
0.320
0.013
0.003
0.336
0.336
1173
24 ACS Paragon Plus Environment
Page 25 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
T825b
1.6
0.320
0.013
0.003
0.336
0.672
1173
T915a
1.6
0.320
0.013
0.003
0.336
0.336
1223
T919a
1.6
0.320
0.013
0.003
0.336
0.336
1123
504 505
25 ACS Paragon Plus Environment
Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47
Page 26 of 42
506
CO2
Air
CH4
N2
507 508
Figure 1. A schematic diagram of the PTGA
509
26 ACS Paragon Plus Environment
Page 27 of 42
0 -0.5 -1 weight change, %
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
-1.5 -2 -2.5 -3 -3.5 -4 -4.5 0
20
40
60
80
100
Time, s
510 511
Figure 2. Ilmenite weight change (T=1223 K, pHC=0.34 MPa, pFuel:pCO2=1) for tests T915a
512
(dashed line) and T505a (solid line)
513
27 ACS Paragon Plus Environment
Energy & Fuels
1 0.9 0.8 Conversion ratio
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 28 of 42
0.7
1023 K
0.6
1048 K
0.5
1073 K
0.4
1123 K
0.3 0.2
1173 K
0.1
1223 K
0 0
50
100
Time, s
514 515
Figure 3. Conversion ratio at 0.9 MPa for pFuel:pCO2=1. Tests T511a, T531a, T513a, T504a,
516
T503a and T505a
517
28 ACS Paragon Plus Environment
Page 29 of 42
5 4.5 Oxygen carrying capacity, %
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
4 3.5 3 2.5 2 1.5 1 0.5 0 1023 1048 1073 1123 1173 1223 Temperature, K
518 519
Figure 4. Oxygen carrying capacity at 0.9 MPa for pFuel:pCO2=1. Tests T511a, T531a, T513a,
520
T504a, T503a and T505a
521
29 ACS Paragon Plus Environment
Energy & Fuels
5 4.5 4 weight change, %
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 30 of 42
3.5 3 2.5
1023K 1048K 1073K 1123K 1173K 1223K
2 1.5 1 0.5 0 0
50
100
Time, s
522 523
Figure 5. Weight change at 0.9 MPa for oxidation stage. Tests T511a, T531a, T513a, T504a,
524
T503a and T505a
525
30 ACS Paragon Plus Environment
Page 31 of 42 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
a
b
c
d
526
527 528
Figure 6. SEM images for test T503b (T=1173 K, top) and T505b (T=1223 K, bottom)
529
31 ACS Paragon Plus Environment
Energy & Fuels
0 -1
Weight change, %
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 32 of 42
-2 -3 CH4
-4
Mixture -5 C2H6 -6 0
50
100
150
200
Time, s
530 531
Figure 7. Ilmenite weight loss for pure methane, the simulated natural gas mixture and pure
532
ethane. Tests T511a, T512 and T628. T=1023 K, P=0.9 MPa, pHC 0.34 MPa, pFuel:pCO2=1
533 534
32 ACS Paragon Plus Environment
Page 33 of 42
0 -0.5 -1 weight change, %
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
-1.5 -2 -2.5
mixture
-3
CH4
-3.5 -4 -4.5 -5 0
20
40
60
80
100
Time, s
535 536
Figure 8. Ilmenite weight loss for pure methane and simulated natural gas mixture. Tests T503a
537
(solid line) and T922 (dotted line). T=1173 K, P=0.9 MPa, pHC=0.34 MPa, pFuel:pCO2=1
538 539
33 ACS Paragon Plus Environment
Energy & Fuels
0 -0.5 -1 weight change, %
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 34 of 42
-1.5 -2 -2.5 -3 -3.5 -4 -4.5 -5 0
50
100
150
Time, s
540 541
Figure 9. Weight loss curves at T=1173 K, P=0.6 MPa, pHC=0.126 MPa. T920a, pFuel:pCO2=1
542
(solid line) and T920b, pFuel:pCO2=0.5 (dotted line)
543
34 ACS Paragon Plus Environment
Page 35 of 42
0 -1 Weight change, %
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
-2 -3 -4 -5 -6 0
50
100
150
200
Time, s
544 545
Figure 10. Weight loss curves at T=1173 K, P=0.9 MPa, pHC=0.126 MPa. T923a,
546
pFuel:pCO2=1 (solid line) and T923b, pFuel:pCO2=0.5 (dotted line)
547
35 ACS Paragon Plus Environment
Energy & Fuels
0 -0.5 -1 weight change, %
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 36 of 42
-1.5 -2 -2.5 -3 -3.5 -4 -4.5 -5 0
50
100
150
Time, s
548 549
Figure 11. Weight change at T=1173 K, P=1.6 MPa, pHC=0.126 MPa. T825a, pFuel:pCO2=1
550
(solid line) and T825b, pFuel:pCO2=0.5 (dotted line)
551
36 ACS Paragon Plus Environment
Page 37 of 42
0 -0.5 -1 weight loss, %
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
-1.5 -2 -2.5 -3 -3.5 -4 -4.5 -5 0
50
100
150
200
Time, s
552 553
Figure 12. Weight loss curves for tests T927a and T504a. T=1123 K, P=0.9 MPa, pFuel:pCO2=1.
554
pHC=0.34 MPa (solid line) and pHC=0.126 MPa (dotted line)
555
37 ACS Paragon Plus Environment
Energy & Fuels
0 -0.5 -1 weight loss, %
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 38 of 42
-1.5 -2 -2.5 -3 -3.5 -4 -4.5 -5 0
50
100
150
Time, s
556 557
Figure 13. Weight loss curves for tests T923a and T503a. T=1173 K, P=0.9 MPa, pFuel:pCO2=1.
558
pHC=0.34 MPa (solid line) and pHC=0.126 MPa (dotted line)
559
38 ACS Paragon Plus Environment
Page 39 of 42
0 -0.5 -1 weight loss, %
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
-1.5 -2 -2.5 -3 -3.5 -4 -4.5 -5 0
50
100 Time, s
560 561
Figure 14. Weight loss curves for tests T926a and T505a. T=1223 K, P=0.9 MPa, pFuel:pCO2=1.
562
pHC=0.34 MPa (solid line) and pHC=0.126 MPa (dotted line)
563
39 ACS Paragon Plus Environment
Energy & Fuels
0.4 0.35 0.3 1-(1-Xr)1/3
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 40 of 42
0.25 0.2 y = 0.0112x R² = 0.9948
0.15 0.1 0.05 0 0
10
20
30
40
Time, s
564 565
Figure 15. Data (test T504a) fitting to obtain k. T=1123 K, P=0.9 MPa, pHC=0.34 MPa,
566
pFuel:pCO2=1
567
40 ACS Paragon Plus Environment
Page 41 of 42
80 P
k value change, %
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Energy & Fuels
60
pFuel
40
pCO2
20 0 0.9
1.1
1.3
1.5
-20 -40
Total pressure, MPa(a)
568 569 570
Figure 16. Individual effects on k by total pressure P, fuel partial pressure pFuel and CO2 partial pressure pCO2, assuming 20 vol.% fuel and 20 vol.% CO2 in the reaction gas mixture
571
41 ACS Paragon Plus Environment
Energy & Fuels
1 0.9 0.8 0.7 Conversion
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 42 of 42
0.6 0.5 0.4 0.3 0.2 0.1 0 0
50
100
Time, sec
572 573
Figure 17. Comparison between model predictions (solid lines) and experimental data (symbols).
574
o: T825a; ∆: T919a; +: T503a; *: T927a
575
42 ACS Paragon Plus Environment