Abstract

A better understanding of the reaction mechanism and kinetics of dry reforming of methane (DRM) remains challenging, necessitating additional research to develop robust catalytic systems with high catalytic performance, low cost, and high stability. Herein, we prepared a zirconia-alumina-supported Ni-Fe catalyst and used it for DRM. Different partial pressures and temperatures are used to test the dry reforming of methane reaction as a detailed kinetic study. The optimal reaction conditions for DRM catalysis are 800°C reaction temperature, 43.42 kPa CO2 partial pressure, and 57.9 kPa CH4 partial pressure. At these optimal reaction conditions, the catalyst shows a 0.436 kPa2 equilibrium constant, a 0.7725 /gCat/h rate of CH4 consumption, a 0.00651 /m2/h arial rate of CH4 consumption, a 1.6515 /gCat/h rate of H2 formation, a 1.4386 molCO/gCat/h rate of CO formation. This study’s findings will inspire the cost-effective production of robust catalytic systems and a better understanding of the DRM reaction’s kinetics.

1. Introduction

Dry reforming of methane (DRM) has received worldwide interest in terms of its ability to reduce the concentration of greenhouse gases (CH4 and CO2) and its efficiency in producing an important synthetic feedstock known as “syngas” (CO + H2) (reaction 1). Dry reforming of methane is a highly endothermic reaction, and it is feasible between 600 and 800°C reaction temperature [1]. The reverse water gas shift reaction (RWGS) is the most significant competitive reaction (reaction 2) with DRM. As RWGS consumes H2, the presence of this reaction may affect the H2 yield of the reaction. It is the responsibility of the DRM research community to optimize the reaction parameters, such as the partial pressure of each feed gas component and reaction temperature, to maximize the catalytic activity of the DRM reaction while minimizing the other parallel side reactions (such as RWGS).

The DRM reaction proceeds in two steps: the dissociation of CH4 over catalytically active sites followed by the oxidation of dissociated CHx species by CO2. Slow dissociation of CH4 results in low activity, whereas delayed oxidation leads to coke deposition over catalytic active sites, ultimately affecting the catalytic activity of the catalyst used. This necessitates that DRM scientists address both aspects, namely how to increase the catalytic active site and provide instant oxidation of carbon deposits. Metal centres such as Ni, Co, Rh, Pd, and Pt are catalytically active sites for CH4 dissociation, with Ni having numerous advantages over the others. Ni offers low-cost preparation and has 25 times more CH4 interaction energy than Co and less CH4 dissociation energy than Pd and Pt [2, 3]. However, Ni sintering at high temperatures remains problematic, leading to the size growth of metallic Ni to the point of inactivation. For CH4 decomposition over silica-supported Fe, the catalyst was inactive, whereas redox metal oxides such as ZrO2-supported Fe catalysts were active [4, 5]. The redox property of ZrO2 enables the support to release instant lattice oxygen during the surface reaction, thus leaving a vacancy behind. Furthermore, ZrO2 is supposed to enhance CO2 dissociation, forming oxygen [6, 7]. ZrO2 may be channel oxygen flow for a surface oxidation reaction under DRM conditions of high reaction temperatures. Overall, it can be stated that Fe dispersed over ZrO2 is responsible for CH4 dissociation, and “CO2 along with ZrO2” channelizes the oxygen flow for instant oxidation of carbon deposits during the DRM reaction. ZrO2-supported Fe catalysts achieved ∼20% CH4 conversion up to 4 h TOS [4, 5].

Interestingly, alumina-silicate (Al/Si = 80 : 20) supported Fe catalyst and alumina-supported Ni catalyst showed <5% CH4 conversion after 4 h, >40% CH4 conversion after 5 h, and >80% CH4 conversion up to 5 h, respectively [4]. With increasing alumina content, acidity has increased, and high and stable CH4 conversion has been observed. Apart from this, Fe dispersion over basic support MgO was also found to be moderately active for CH4 decomposition [4, 8] (∼60% CH4 conversion, ∼50% H2 yield).

Fe dispersion over Al2O3 was highly active for CH4 decomposition, but this catalyst system is inactive in DRM. In the presence of CO2, the dissociation of CH4 may be inhibited by an iron-based catalyst, as Fe is oxidized into FeO [9], which is incapable of dissociating CH4. At the same time, in the presence of CO2, the Ni-supported catalyst was found to be quite active for CH4 decomposition. However, it is rapidly rendered inactive due to coke deposition. Optimizing Ni and Fe content was beneficial regarding catalytic activity and stability. Ni-Fe alloy is stable at extremely high temperatures; therefore, catalytically active metallic Ni is maintained under DRM conditions [10]. Al2O3-supported NiFe (3 : 1) catalyst had a specific alloy composition (Ni3Fe), resulting in improved catalytic activity (13% CH4 conversion after 3 h TOS) than Al2O3-supported Ni catalyst (8% CH4 conversion after 3 h) at 600°C [11]. In low-temperature DRM, evaporation-induced self-assembly (EISA)-prepared NiFe catalysts supported on alumina have garnered interest. It showed 26.6% CH4 conversion, 37.8% CO2 conversion, and 0.67 H2/CO ratio at 550°C [12]. Li et al. prepared a mesoporous alumina-supported NiFe (Ni/Fe = 10/7) catalyst using the EISA method and found that the catalyst was deactivated after 24 hours due to the dealloying of FeNi3 [13]. Gunduz-Meric et al. prepared a coke and sinter-resistant “Ni-iron core (ratio 4 : 1) and silica sphere” catalyst [14]. Silica shell formed SiC and protects carbon decomposition in the presence of catalytically active Ni, as well as limiting Ni dispersion. It showed more than 70% CH4 conversion and a 0.70 H2/CO ratio for up to three cycles. Fe-modified MgAl2O4-supported Ni catalysts and MgAl2O4 supported NiFe catalysts were also tested for DRM [10, 15]. At Ni/Fe = 1.4, loaded over MgAl2O4 support, the role of iron was found to be crucial in the decoking at the metal centre [10]. Here, some of the metallic Fe is oxidized to FeOx by CO2. Besides, the lattice oxygen of FeOx is superior for decoking, producing CO and metallic Fe, where the latter restores the original Fe-Ni alloy [16]. The second pathway for coke oxidation involves the dissociation of CO2 over metallic Ni into CO and O, followed by subsequent coke oxidation by surface oxygen [9]. Ni-Al (3 : 1) supported on Mg(Al)O demonstrated >1.5 molCH4 molmetal−1 s−1 CH4 conversion for up to 30 h time on stream. The effective decoking at the metal centres was due to the reaction of FeO with surface carbon [17]. Mayenite (Ca12Al14O33) support Ni catalyst is prone to deactivation in DRM, but adding 2 wt.% Fe [18] effectively suppressed coke. The support may facilitate the transfer of oxygen species from FeOx to Ni sites, thereby promoting carbon deposit oxidation.

We expect that by incorporating ZrO2 into Al2O3, lattice oxygen endowing capacity will be enhanced and that by introducing Fe with Ni, effective decoking will occur. Herein, we prepare a zirconia-alumina-supported Ni-Fe catalyst. Different partial pressures and temperatures are used to test the dry reforming of methane in a detailed kinetic study. To the authors’ best knowledge, this is the first detailed study to optimize the partial pressure of each feed gas component over a 700–800°C reaction temperature to maximize the DRM catalytic activity and minimize side reactions like RWGS.

2. Experimental

2.1. Catalyst Preparation

The catalyst 5Ni2Fe/ZrAl is synthesized by impregnating the required amounts of Ni(NO3)2.6H2O (99%; Alfa Aesar) aqueous solution (equivalent to 5 wt% NiO) and Fe(NO3)3.9H2O (403.99 g/mol; 99.99%; Alfa Aesar) (equivalent to 2 wt% of Fe2O3) with commercially available 10 wt% zirconium oxide–90 wt% alumina support. The solution was heated under stirring until a slurry was formed. It was further dried at 120°C and calcined at 700°C with a heating rate of 3°C/min for 5 hrs. The catalyst is abbreviated as 5Ni2Fe/ZrAl.

2.2. Catalyst Reaction

DRM reaction is carried out in a stainless tubular reactor (diameter 0.91 cm and length 30 cm) made by PID Eng. & Tech. Micro. Activity company. 100 mg catalyst is packed in the reactor and reduced under reductive pretreatment with H2 (flow rate 20 ml/min) at 600°C for 1 h. The gas feed is composed of CH4, CO2 and N2, which are allowed to pass through the catalyst bed at different flow rates (total flow rate = 70 ml/min) at three different reaction temperatures (700, 750, and 800°C). The products are analyzed by a gas chromatograph equipped with a TCD detector.

The expressions for mole fraction, partial pressure, and specific feed rate of each gas are shown in supporting information S1. In this manuscript, we have studied different activity terms (shown below) at a different partial pressure of gas feed at 700°C, 750°C, and 800°C reaction temperatures. conversion, , and H2/CO ratio are expressed by equations (4)–(6). The specific feed rate of a gas is defined as the flow rate of a gas per gram weight of catalyst. The rate of CH4 consumption and rate of CO2 consumption are shown by equations (6) and (7) , respectively. The details of expressions for rate of H2 formation, rate of CO formation, and rate of H2O formation are derived in supporting information S2 and the final expressions are shown in equations (9)–(11). The arial rate of gas consumption is defined as the rate of gas consumption per unit surface area per gram weight of catalyst. The expression for the arial rate of CH4 consumption and the arial rate of CO2 consumption are shown in equations (10) and (11), respectively.

3. Results and Discussion

3.1. Characterization Results

The catalyst had a 118.6 m2/g BET surface area and type IV adsorption/desorption isotherms with an H1 hysteresis loop (Figure 1(a)). The sharp inflection between 0.6 and 0.75 relative pressure regions in the isotherm indicates capillary condensation, an indication of uniformity of pore distribution in the mesoporous material [19, 20]. Pore size distribution over the catalyst surface is shown by the dV/d log W vs. W plot (where V is volume and W is the pore width) in Figure 1(a) (inset). The majority of the pores on the surface of the catalyst are 27 nm in size. Based on the BJH pore size measurement, the average pore size is determined to be 15.5 nm. The fresh 5Ni2Fe/ZrAl catalyst has a thermally stable tetragonal ZrO2 phase (at 2θ = 30, 50, 60°; JCPDS reference number: 00-024-1164), rhombohedral Al2O3 phase (at 2θ = 36.8, 45.5, 60, 66°; JCPDS reference number: 01-077-2135), and cubic NiAl2O4 phases (at 2θ = 30, 36.8, 45.5, 50, 60, 66°; JCPDS reference number: 01-071-0965) (Figure 1(b)). Previously, metallic Ni derived from thermally stable NiAl2O4 (during reduction) was claimed prominent active site for high catalytic activity in DRM [2123]. The IR spectra reveals the presence of physically adsorbed CO2 at approximately 2349 cm−1 [24], carbonate species at 1384 cm−1, and formate at 2850 and 2925 cm−1 [7] (Figures 1(c) and 1(d)). The IR peak at 1631 cm−1 and 3444 cm−1 indicates the bending and stretching vibrations of O-H, respectively [25]. Figures 1(e) and 1(f) show an HR-TEM image of a fresh and spent catalyst. The spent catalyst has carbon nanotubes with a variable diameter.

3.2. Catalytic Activity Result and Discussion

The catalyst 5Ni2Fe/ZrAl has uniform mesopores with an average size of 15.5 nm. It possessed thermally stable support (as tetragonal ZrO2 phase and rhombohedral Al2O3 phase), metallic Ni (derived from thermally stable NiAl2O4 upon reduction), catalytic active sites, and various types of “CO2 interacting species (physically adsorbed CO2, carbonate or formate)” over the catalyst surface [23, 24, 26]. The role of Fe-Ni was previously claimed in coke suppression [18]. A uniform mesoporous, thermally stable, metallic Ni dispersed (derived from NiAl2O4 after reduction), and CO2-interacting catalyst surface seems efficient for the DRM and RWGS reaction. Table 1 shows the catalytic activity results of the 5Ni2Fe/ZrAl catalyst in terms of an equilibrium constant, CH4 conversion, CO2 conversion, and H2/CO ratio. At low CH4 partial pressure and high CO2 partial pressure, CH4 conversion is very high (Table 1, Entry 5-6, 14-15, and 23-24) due to the instant oxidation of the substrate (CH4) by a large number of oxidants (CO2). Similarly, at low CO2 partial pressure and high CH4 partial pressure, CO2 conversion is very high (Table 1, Entry 1-2, 10-11, 19-20) due to the instant utilization of the oxidating agent (CO2) by a large number of substrates (CH4). From a catalytic standpoint, the “mole of CH4 or CO2 consumption per gram mass of catalyst per hour,” “mole of CH4 or CO2 consumption per gram mass of catalyst per hour per surface area,” and “mole of product formation per gram mass of catalyst per hour” are exact presentations of catalytic activity. Thus, furthermore, Table 1 contains catalytic activity data regarding the rate of CH4 consumption, areal rate of CH4 consumption, the rate of CO2 consumption, the areal rate of CO2 consumption, the rate of H2 formation, the rate of CO formation, and the rate of H2O formation at a different partial pressures of feed gas during 700, 750, and 800°C reaction temperatures.

In the dry reforming of methane, CH4 is the substrate and CO2 is the oxidant. Upon increasing the partial pressure of CH4 (substrate) or CO2 (oxidant), the rate of CH4 and CO2 consumption increases at reaction temperatures of 700, 750, and 800°C (Figures 2(a)2(d)). This suggests that the presence of an increasing amount of substrate (CH4) over a fixed oxidant (CO2) or the presence of an increasing amount of oxidant (CO2) over a fixed substrate (CH4) gives rise to more collision, and a higher rate of conversion at a given temperature. At 43.43 kPa partial pressure of CH4 and 14.48 kPa partial pressure of CO2, the H2/CO ratio is found to be 1.4, 1.86, and 1.2 at 700, 750, and 800°C, respectively (Table 1 Entry 1, 10, and 19). The high H2/CO ratio at these partial pressures is due to the availability of a high concentration of CH4 (which is primarily responsible for H2 generation) as well as a low concentration of CO2 (which is mainly responsible for CO generation) over the catalyst surface.

At constant CO2 partial pressure (43.425 kPa) and increasing CH4 partial pressure (from 14.475 to 57.9 kPa) at 800°C, the rate of CH4 consumption is significantly increasing (Table 1 Entry 23–26, Figure 2 A). At constant 43.425 kPa CH4 partial pressure at 800°C, the rate of CH4 consumption increases sharply (0.2623 /gcat/h to 0.4104 /gcat/h) up to 29 kPa CO2 partial pressure, after which it remains constant (Table 1 Entry 19–22, Figure 2(c)). This finding needs to be explained in more detail. Upon doubling the partial pressure of CO2 (from 28.95 kPa to 57.9 kPa), the rate of CH4 consumption remains constant at 800°C reaction temperature, whereas the rate of CO2 consumption doubles (0.4490 /gcat/h at 28.97 kPa CO2 partial pressure to 0.8981 /gcat/h at 58 kPa CO2 partial pressure) (Figure 2(d)). It indicates that as CO2 partial pressure increases from 29 kPa to 58 kPa, CO2 remains converted but does not oxidize CH4 (as in conventional DRM reaction). In this partial pressure range, it may oxidize the H2 or carbon deposit on the catalyst surface. Under the same conditions, the rate of hydrogen formation decreases by more than 50% (0.7821 /gcat/h at 28.97 kPa CO2 partial pressure to 0.3330 /gcat/h at 58 kPa CO2 partial pressure) (Figure 3(c)), the rate of CO formation increases by up to 1.5 times (0.8594 molCO/gcat/h at 28.97 kPa CO2 partial pressure to 1.3085 molCO/gcat/h at 58 kPa CO2 partial pressure) (Figure 3(d)) and rate of H2O formation increases by up to 12 times (0.0387 /gcat/h at 28.97 kPa CO2 partial pressure to 0.4877 /gcat/h at 58 kPa CO2 partial pressure) (Figure 3(f)). This observation suggests that the reverse water gas shift reaction is accelerating in the 29 kPa to 58 kPa CO2 partial pressure range at 800°C reaction temperature, whereas the dry reforming of methane is just continuing at the same rate.

At 750°C reaction temperature and constant CH4 partial pressure (43.42 kPa), with a rise of CO2 partial pressure from 14.475 kPa to 28.95 kPa, the rate of CH4 formation constantly increases (from 0.3408 /gcat/h to 0.3772 /gcat/h), but the rate of H2 formation decreases from 0.7893 /gcat/h to 0.7358 /gcat/h, and the rate of water formation becomes significant. This observation indicates the presence of a reverse water gas shift reaction (Table 1 Entry 10–12, Figures 2(c), 3(c), and 3(f)). However, on the further increase of CO2 partial pressure up to 57.9 kPa; a rise of rate CH4 consumption (0.3882 /gcat/h to 0.4840 /gcat/h), a rise in the rate of CO2 consumption (0.4289 /gcat/h to 0.5239 /gcat/h), a rise in the rate of H2 formation (0.7358 /gcat/h to 0.9281 /gcat/h), and a rise in the rate of CO formation (0.8171 molCO/gcat/h to 1.0081 molCO/gcat/h) were achieved without affecting the rate of H2O formation much (Table 1 Entry 12-13, Figures 2(c), 3(c), and 3(f)). Thus, RWGS is competent up to 28.95 kPa CO2 partial pressure, but at 57.9 kPa CO2 partial pressure, RWGS product formation rates are not significantly affected, whereas DRM product formation rates are significantly affected. At constant CO2 partial pressure (43.425 kPa) and increasing CH4 partial pressure (from 14.475 kPa to 57.9 kPa) at 750°C, the rate of CH4 consumption is increased to about ∼2.5 times (0.2411 /gcat/h to 0.6013 /gcat/h), rate of CO2 consumption is increased to about ∼2.5 times (0.2645 /gcat/h to 0.6444 /gcat/h), the rate of H2 formation is again increased to 2.5 times (0.4588 /gcat/h to 1.1594 /gcat/h), and the rate of CO formation is again increased to about ∼2.5 times (0.5056 molCO/gcat/h to 1.2456 molCO/gcat/h) (Table 1 Entry 14–17, Figures 2(a), 2(b), 3(a), and 3(b)). No such progressive correlation with the rate of H2O formation is found (Figure 3(e)), but it remains significant. It demonstrates that at 750°C reaction temperature, with constant CO2 partial pressure and increasing CH4 partial pressure, the dry reforming of methane reaction progressed continuously.

At a low reaction temperature of 700°C, the least activity is noticed. On constant CH4 partial pressure (43.42 kPa) and increased partial pressure of CO2 from 14.47 kPa to 57.9 kPa, the rate of H2 production decreases continuously (0.6638 /gcat/h to 0.5020 /gcat/h), the rate of CO formation increases continuously (0.4902 molCO/gcat/h to 0.7765 molCO/gcat/h), and the rate of H2O formation is found to a maximum of 0.1372 /gcat/h at 57.9 kPa CO2 partial pressure (Table 1 Entry 1–4; Figures 3(c)3(f)). It indicates that at a constant CH4 partial pressure and a rising partial pressure of CO2 at 700°C, the RWGS product formation rate influences the DRM product formation rate. However, at constant CO2 partial pressure (43.42 kPa) and rise of CH4 partial pressure (14.47 kPa to 57.9 kPa), the rate of H2 formation and rate of CO formation are increased, but no such correlation is found with the rate of water formation (Table 1 Entry 5–8; Figures 3(a), 3(b), and 3(e)). It indicates that a higher rate of H2 and CO formation may be caused by a pronounced DRM reaction in which a portion of H2 participates in the RWGS reaction, but this does not affect the high rate of H2 production.

Comparing the catalytic activity results at all temperatures reveals that at 14.47 kPa CO2 partial pressure, 43.42 kPa CH4 partial pressure, and 43.42 kPa N2 partial pressure, the RWGS reaction does not occur at 700°C, 750°C, and 800°C reaction temperatures. In these instances, DRM activity is, however, low. At 800°C reaction temperature, the 43.42 kPa CO2 partial pressure and 57.9 kPa CH4 partial pressure (with diluent N2), the equilibrium constant of DRM is maximum (0.436 kPa2). In this reaction condition, the RWGS reaction is not exiting, and the rate of CH4 consumption (0.7725 /gCat/h), the arial rate of CH4 consumption (0.00651 /m2/h), the rate of H2 formation (1.6515 /gCat/h), and the CO formation (1.4386 molCO/gCat/h) are the maximum among all tested conditions. In the mean of the maximum rate of CO2 consumption, 57.9 kPa CO2 partial pressure, 43.42 kPa CH4 partial pressure and 800°C reaction temperature are found appropriate. The equal partial pressure of CH4 and CO2 (50.66 kPa) at 750°C, 0.6145 /gCat/h rate of CH4 conversion, 1.0795 /gCat/h rate of H2 formation, and 1.3785 molCO/gCat/h CO formation were noticed.

The apparent activation energy for CH4 dissociation, CO2 dissociation, H2 formation, and CO formation is shown in Figure 4. It is interesting to note down that with an increasing flow rate of CH4 from 10 ml/min to 40 ml/min and a fixed CO2 flow rate of 30 ml/min, the apparent activation energies for CH4 dissociation, CO2 dissociation, and H2 formation are increasing (Figure 4). Upon increasing the flow rate of CO2 and fixing CH4 flow rate at 30 ml/min, no such correlation is found (Figure S3). It indicates that the C-H dissociation of CH4 is a rate-determining step (Figure 4 and Figure S3). The Mears criterion (for external diffusion concerning CH4 and CO2) and Weisz-Prater criterion (for internal diffusion for CH4 and CO2) values for the 5Ni2Fe/ZrAl catalyst were found < 0.1 and < 1 (in the Supporting Information S4), respectively, in every case [27]. The absence of both external and internal mass transfer limitations is found in the 5Ni2Fe/ZrAl catalyst at various gas feed rates.

4. Conclusion

The thermally stable catalytic active sites NiAl2O4 as well as the CO2-interacting mesoporous surface of the 5Ni2Fe/ZrAl catalyst were found to be effective in the DRM reaction and a competing RWGS reaction. C-H dissociation of CH4 is the rate-determining step. Upon different flow rates of gas feed over the 5Ni2Fe/ZrAl catalyst, external and internal mass transfer limitations are absent. At 800°C reaction temperature, the constant 43.43 kPa CH4 partial pressure, and 28.95 to 57.9 kPa CO2 partial pressure range, the reverse water gas shift reaction is accelerated over a zirconia-alumina-supported Ni-Fe catalyst. Again, at 750°C reaction temperature, at a constant 43.43 kPa CH4 partial pressure and 28.95 kPa CO2 partial pressure, RWGS is noticed. At 750°C, 43.425 kPa CO2 partial pressure and 14.475 kPa to 57.9 kPa CH4 partial pressure, the rate of DRM reaction is increased to about 2.5 times. At a reaction temperature of 700°C, the catalyst’s performance is diminished. The optimal reaction conditions for DRM catalysis are 800°C reaction temperature, 43.42 kPa CO2 partial pressure, and 57.9 kPa CH4 partial pressure. At these optimal reaction conditions, the catalyst shows a 0.436 kPa2 equilibrium constant, a 0.7725 /gCat/h rate of CH4 consumption, a 0.00651 /m2/h arial rate of CH4 consumption, a 1.6515 /gCat/h rate of H2 formation, a 1.4386 molCO/gCat/h rate of CO formation. On increasing the flow rate of CH4, the apparent activation energy for CH4 dissociation, CO2 dissociation and H2 formation are increasing which indicates that C-H dissociation of CH4 is a rate-determining step in DRM over a 5Ni2Fe/ZrAl catalyst. The mass transfer limitation is absent over this catalyst.

Data Availability

No underlying data were collected or produced in this study.

Disclosure

The views and opinions expressed in this paper do not necessarily reflect those of the European Commission or the Special EU Programmes Body (SEUPB).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to extend their sincere appreciation to Researchers Supporting Project number (RSP2023R368), King Saud University, Riyadh, Saudi Arabia. KA and RK acknowledge Indus University, Ahmedabad, India, for supporting research. Ahmed I. Osman wishes to acknowledge the support of The Bryden Centre Project (Project ID VA5048), which was awarded by The European Union’s INTERREG VA Programme, managed by the Special EU Programmes Body (SEUPB), with match funding provided by the Department for the Economy in Northern Ireland and the Department of Business, Enterprise and Innovation in the Republic of Ireland. The authors would also like to thank Charlie Farrell for proofreading the manuscript.

Supplementary Materials

The expression for mole fraction, partial pressure, and specific feed rate of each gas is shown in Supporting Information S1. The detail expressions for rate of H2 formation and rate of CO formation are shown in Supporting Information S2. Fig. S3 shows influence of the reaction temperature on the CH4 consumption rate, CO2 consumption rate, H2 formation rate and CO formation rate upon fixing CH4 flow rate. Supporting Information S4 shows calculation of mass transfer limitation. (Supplementary Materials)