1. Introduction
Currently, greenhouse gases (chiefly CO
2 and CH
4), which can lead to global warming and energy shortages, are major challenges for the world. It is shocking that the CO
2 composition in the atmosphere has reached the level of 400 ppm, which has elevated the temperature of the earth by 1.5 °C [
1]. The potentiality of methane to cause global warming is 25 times higher than that of CO
2 [
2]. Thus, it is necessary to come up with a strategy that alleviates the carbon footprint of contemporary civilizations and tackles these challenges. Dry reforming of methane (DRM) is regarded among the existing strategies owing to its ability to consume CO
2 and CH
4 so as to generate syngas (CO+H
2), an environmentally safe and cheap product [
3,
4,
5]. At present, large-scale technologies for manufacturing syngas have been performed. The emerging industrial processes comprise lab-scale studies of DRM, partial oxidation, and methane steam reforming [
6]. Compounds such as those produced by hydroformylation are directly manufactured from syngas. The H
2/CO ratio can be attuned to producing dimethyl ether, methanol, oil, or fuel via the Fischer–Tropsch process. DRM is usually accompanied by a reverse water–gas shift reaction.
Noble metals have good performance in the reforming process; however, they are ruled out because of their scarcity and high cost [
7]. Transition metals, such as Ni and Co, are universally considered to be suitable catalysts for DRM because of their low cost, but deactivation due to carbon deposition hampers the industrial application of the process. The stability and activity can be enhanced by the use of proper support. Mourhly et al. used mesoporous silica as support for Ni-based catalysts in DRM [
8]. The results showed that the catalysts exhibited long-term stability with low coke formation predominantly as carbon nanotubes. Arif et al. [
9] employed Ni/ZrO
2 and Ni/CaO catalysts to generate H
2 through CO
2 reforming of glycerol. The performance and activity of the catalysts usually depend on several factors, which include the diameter, support type, promoters, acid–base features, and oxidation numbers of the metals [
10]. For Ni particles to be used as catalysts for DRM at 500 to 600 °C and 5.0 bar, the most effective size is around 2–3 nm [
11]. Computational performance studies using a single-site Ni/MgO catalyst in DRM revealed that the binding between the Ni and MgO had too low of a synergistic effect because of the limited number of active sites. However, increasing the Ni over the MgO enhanced the binding and promoted the removal of carbon [
12]. Al-Fatesh et al. investigated the DRM using an Ni catalyst supported by a combination of La
2O
3 and ZrO
2 [
13]. It was found that the Ni dispersion over the lanthana–zirconia support stabilized the tetragonal ZrO
2 phase. The support had additional amorphous La
2O
2CO
3 species and format species, which indicated a broad interaction of CO
2 over the catalyst surface. Numerous studies have evaluated the function of support in the performance of Ni catalysts. For instance, Ni/ZnAl
2O
4 exhibited the highest resistance to carbon deposition and metal sintering. Supports of alumina and silica oxidize Ni during the process, while ZrO
2 stabilizes Ni particles, thus increasing the efficiency of methane conversion into fuel. The stability is due to the water accumulation behavior of ZrO
2, which favors the formation of hydroxyl for the greater efficacy of SRM and DRM. CeO
2 and ZrO
2 are ideal supports for Ni catalysts in DRM because of their high storage capacity of oxygen, great resistance to coking, and redox properties. These features give superior catalytic properties compared with the conventional supports used so far, such as aluminum oxide and MgAl
2O
3. Doping Ni with other metals as promoters improves the activity, efficiency, stability, and durability of the catalysts. As a model, the combination of Ni with Fe to form the alloy Ni-Fe showed a methane conversion of 97.5% and exhibited selectivity of 92.9% for carbon monoxide at a temperature of 900 °C [
14]. Wysocka and colleagues examined the impact of a K promoter on Ni-based catalysts for DRM [
15]. Their results revealed that K promotion of nickel altered the nickel distribution and metal–support interactions and decelerated carbon deposition while enhancing the sorption of carbon dioxide. Abd Ghani et al. found that Nb promoters improved the catalytic activity of the Ni-based catalyst because of the improvement in the morphology and the electronic adjustment of the catalyst [
16]. The Nb promoter increased CO
2 conversion. This study aims to assess the performance activities in DRM with Ni/ZrO
2 catalysts, including aspects of the support modifiers MgO and Y
2O
3; different loading of the best modifiers; the effects of Sc, Ga, and Sr promoters; and, finally, the sensitivity test of the optimum non-promoted catalyst in terms the temperature and flow rate. Relevant characterization methods were performed to better comprehend the different features of the catalytic system.
3. Results
For verification of the performed preparations, a sample (5Ni10YZr) was tested using EDX characterization; in
Figure S1, the images display the peaks for O
2, Ni, Y, and Zr, and therefore, the elemental analysis was found to be consistent with the compositions of the tested catalyst. The quantity aspects of the elements closely resembled those of the given sample. Zirconium dioxide has three crystalline phases: a monoclinic phase below 1170 °C, a tetragonal phase between 1170 °C and 2370 °C, and a cubic phase above 2370 °C [
17,
18]. ZrO
2 exhibits thermal stability, oxygen vacancies, and electrical conductivity, and its properties can be improved by adding cations of lanthanide and alkaline earth, such as Mg
2+, Ca
2+, La
3+, or Y
3+, into the ZrO
2 lattice because it has a lower valence than Zr
4+, leading to an overall negative charge [
19]. As shown in
Figure 1, the 5Ni5YZr catalyst improved the conversion of methane and carbon dioxide. While the 5Ni5MgZr reduced the activity performance with respect to pristine catalyst (5NiZr), the CH
4 and CO
2 conversions decreased in the order 5Ni5YZr > 5NiZr > 5Ni5MgZr. The observed variations of activities with the addition of Mg and Y could be related to the metal support interactions, which, in the case of Y, led to a decrease in CH
4 and CO
2 activation energies and hence promoted the conversions. By contrast, in Mg, the opposite phenomenon was observed, in which the metal and the modified support interactions increased the energy activations of CH
4 and CO
2, and the reactivity thus decreased.
Recent studies have shown the impact of yttrium oxide (Y
2O
3) on the stabilization of ZrO
2 [
20].To determine the optimum Y
2O
3 loading, five different wt.% loading conditions (x = 0.0, 5.0, 10.0, 15.0, and 20.0 wt.%) were applied to synthesize the 5Ni xYZr. The essential features of the fresh catalysts were examined using nitrogen adsorption–desorption isotherms.
Figure S2 displays the outcomes of the isotherms. According to the IUPAC sorting of isotherms, the isotherms fell under the classification of type V, with an H1-type hysteresis loop that resulted from the evaporation at high relative pressures and the capillary condensation [
21]. The relative pressure (P/Po) augmented in the span of 0.80–1.00. In addition, the low specific surface areas of all catalysts resulted from the low specific volume N
2 adsorbed, occurring in the extent of 4.5–8.0 cm
3/g.
Table 1 outlines the textural properties. The pristine ZrO
2 displayed a BET value of 21.0 m
2/g, which diminished to 16.0 m
2/g when NiO was impregnated, an indication of proper loading. Nevertheless, the pore size expanded from 37.0 nm for meso-ZrO
2 to 42.0 nm after NiO impregnation because of the evolution of H
2O, NO
2, HNO
3, and O
2 gases and the decomposition of nickel (II) nitrate hexahydrate [
22] throughout the course of the calcination stages. The BET value of the catalysts increased when Y
2O
3 was put in, and this was related to a reduction in the NiO diameter and the creation of O
2 vacancy [
23]. The 5Ni15YZr was characterized by the lowest NiO crystal size and a high surface area.
The hydrogen reduction temperature is an important standard for the catalytic oxidation properties of metal oxide support. H
2-TPR measurements were adopted to identify the effects of the Y
2O
3 content on the reduction temperatures, and the profiles of all prepared catalysts are shown in
Figure 2. The reduction temperature range may be classified into three zones—low, medium, and high—and each zone reflects the interaction and dispersion of the active metal with the support [
24]. Mostly, all catalyst samples display a similar classification of bimodal at a medium temperature range, but with different patterns.
Figure 2 shows the absence of peaks in the low zone (free NiO). The absence of free NiO species suggests that some interaction occurred between NiO and Y
2O
3 or ZrO
2 species. In addition, the weakly bound nickel species, such as nickel oxide (NiO), were easier to reduce but also more easily sintered, resulting in larger particles that promoted carbon formation. At 600 °C and above, there were no reduction peaks, denoting the absence of NiO species. In particular, 5NiZr had a distinctive reduction peak at 434°. On the other hand, adding the Y
2O
3 to the ZrO
2 significantly altered the range of the reduction temperature. Whenever the loading of Y
2O
3 was less than 10 wt.%, the synergy between the support and the active metal grew, consequently shifting the peaks to elevated temperatures. Conversely, if the Y
2O
3 loading exceeded 10 wt.%, synergy effects of the support and the active metal appeared at moderate temperatures. The temperature of Ni reduction varied with the different loading conditions of Y
2O
3, which proves the differences in the activity and catalytic performance of each catalyst.
Because of the acidity of carbon dioxide, it is easily adsorbed on surfaces of basic sites. The basicity measurements of catalysts are essential for the regulation of carbon formation. The excessive basicity of the catalyst is detrimental to catalytic activity, as it stimulates a higher extent of CO
2 dissociation (CO
2 → C + O
2), thus deactivating the catalyst. The phenomenon intensifies when the Boudouard reaction (2CO → C + CO
2) starts to occur at high temperatures because of the enriched composition of CO upon methane reforming activity, thus resulting in a higher quantity of coke deposited on the catalyst’s surface. The basicity of the prepared catalysts was characterized by CO
2-TPD. Moreover, the CO
2 desorption peak area and temperature span regulate the basicity power. The feeble acid sites (surface hydroxyl) are physisorbed at 50–200 °C, while the medium-strength basic sites (surface oxygen anion) are chemisorbed at 200–400 °C, and the strong and very strong basic sites (bulk oxygen anion/oxygen vacancy) interact with oxygen vacancies at 400–650 °C and >650 °C, respectively [
25,
26]. The CO
2-TPD profiles for all prepared catalysts are illustrated in
Figure 3.
The CO
2-TPD profiles show that all prepared catalysts presented the CO
2 desorption peaks in two major ranges: the first desorption peak, which occurred at lower temperatures of 50–200 °C (weak), and the second desorption peak, which occurred at higher temperatures of 200–400 °C (medium) because of the strength of the basic sites. An approximate highest temperature of 665 °C was recorded for the CO
2 desorption peaks in the very strong range with the exclusion of 5Ni15YZr and 5Ni10YZr. These desorption peaks are ascribed to strong basic sites [
27]. Built on these outputs, the basicity sequence was as follows: 5Ni20YZr < 5NiZr < 5Ni5YZr < 5Ni10YZr < 5Ni15YZr. Catalytic basicity has a strong effect on the DRM reaction. However, there are other factors that influence catalytic performance, including surface area, pore size, active metal particle size, and dispersion [
27].
Figure 4 shows the X-ray powder diffraction “XRPD” images of synthesized catalysts in 2θ° = 10–80°. The XRPD image depicts the ZrO
2 phase. The images are similar, but the intensity decreased as the Y
2O
3 amount increased because of the reduction in the NiO diameter, and the Y
2O
3 peaks appeared clearly when their loading exceeded 10 wt.%. The peaks of low loading Y
2O
3 were not visible, as the particles of Y
2O
3 were incorporated into the pores of the zirconia. The diameter of the NiO is shown in
Table 1. MDI Jade
® software (version 6.5, Materials data Inc., Newtown Square, PA, USA) was used to inspect the diffraction data. The ZrO
2 monoclinic phase was identified at 24.170, 28.140, 31.470, 34.250, 38.580, 40.890, 44.870, 49.240, 54.000, 55.400, 57.250, 58.000, 60.000, 62.700, 65.000, 69.000, 71.300, 75.000, and 79.000° (JCPDS 81-1314) [
28], whereas NiO diffraction peaks were identified at 35.340 and 45.500° (JCPDS 65-5745) [
29]. The peaks of high loading Y
2O
3 were identified at 29.200, 31.000, and 50.200° (JCPDS 89–5592) [
30]. The crystalline size calculations of 5NiZr, 5Ni5YZr, 5Ni10YZr, 5Ni15YZr, and 5Ni20YZr catalysts using the Debye−Scherrer method are tabulated in
Table 1. The 5NiZr catalyst presented the highest NiO particle size and the lowest surface area, which negatively affected the activity and stability performance, whereas the 5Ni15YZr catalysts with high surface area yielded the smallest NiO size and the best performance.
Furthermore, the 5Ni15YZr catalyst had the highest conversion of methane and carbon dioxide when compared with the other prepared catalysts. It can be seen in
Figure 5 that CH
4 and CO
2 conversion decreased in the order 5Ni15YZr > 5Ni10YZr > 5Ni5YZr > 5Ni20YZr > 5NiZr. Thus, on the basis of the catalytic performance of methane and carbon dioxide conversions and the different characterization methods of BET, TPR, CO
2-TPD, XRPD, TAG, and TEM, it was shown that the 15 wt.% Y
2O
3 loading was optimal. Several studies [
31,
32,
33,
34] have investigated related systems in which the effect of the incorporation of alkaline earth and lanthanide metal oxides into the active metals, a process that develops oxygen vacancies, is observed. The blockage of the catalyst pores, feed flow, and inactivity results mainly from coke formation.
Figure 6 shows the TGA analysis of the catalysts used. TGA was performed on the spent catalysts to determine the amount of carbon deposited on the catalysts. The weight loss below was caused by the removal of the deposited carbon. The weight drop began at 470 °C because of the burning of the coke, which inhibited the catalytic activity. The percent weight drop corresponded to the extent of the carbon formed. However, the addition of Y
2O
3 to 5NiZr for loading of less than 15.0 wt.% influenced the carbon deposition insignificantly. Alternatively, the 5Ni20YZr sample generated fair coke deposition. This could be related to the fact that the high loading of Y
2O
3 covers the active Ni catalyst excessively. Thus, the reactivity of the active Ni was lowered, and consequently, the amount of carbon deposited decreased. The weight % drop sequence was as follows: 5Ni20YZr < 5NiZr < 5Ni15YZr < 5Ni10YZr < 5Ni5YZr.
The TEM for 5Ni15YZr was displayed to give a better visualization of the surface morphology and carbon deposition on the surface of the catalyst. The fresh sample had better distribution and a reduced active metal diameter. By contrast, the sample used showed slim sheets of multiwall carbon nanotubes.
Figure 7 shows the TEM micrographs of the fresh and used samples of 5Ni15YZr.
The performance of the 5Ni15YZr catalysts was further examined by combining the Ni with promoting elements, such as Cs, Ga, and Sr; 3 wt.% of the promoter was used.
Figure 8 depicts the catalytic conversion of CH
4 and CO
2 against the time on stream. The outcome showed that the first conversions of CH
4 and CO
2 for 5Ni15YZr, 5Ni153CsYZr, 5Ni3Ga15YZr, and 5Ni3Sr15YZr were 56.3%, 57.6%, 61.5%, and 66.5%, respectively, for CH
4 and 64.3%, 72.5%, 77.9%, and 79.7%, respectively, for CO
2. The Sr, Ga, and Cs promoters improved the percent conversion of CH
4 by 18%, 9%, and 2%, respectively, over the non-promoted 5Ni15YZr. This may be ascribed to the formation of alloyed structures with strongly modified electronic properties and the provision of extra active sites for reaction. Thus, the 5Ni3Sr15YZr catalyst outperformed the remaining catalysts and gave rise to the largest (66.5%) mean conversion of CH
4.
The thermogravimetric analysis of the used promoted catalysts obtained after 440 min of reaction was carried out.
Figure 9 exhibits the percent weight loss profiles. The percent weight loss is in accordance with the type of promoter. Thus, the 5Ni15SrYZr catalyst produces 9.7% weight loss, which is the highest value, while the 5Ni15CsYZr catalyst generated a minimum weight loss of 2.6% as a result of carbon deposition. For the promoted catalysts, the higher the activity of the sample, the greater the amount of carbon formed, which was in line with the scale of the reaction. The promoted samples displayed slight carbon formation in contrast to the non-promoted sample (5Ni15YZr), providing a 25.60% weight drop. The promoted catalysts formed stable ZrO
2 phases and additional mixed phases of cubic zirconium promoter oxide. Therefore, promoters enhanced the stability of the catalyst. This could be associated with the fact that the elevated loading of Y
2O
3 covered the active Ni catalyst excessively. Thus, the reactivity of the active Ni was lowered, and consequently, the amount of carbon deposited decreased. The percent weight drop sequence of the catalysts was as follows: 5Ni20YZr < 5NiZr < 5Ni15YZr < 5Ni10YZr < 5Ni5YZr.
Figure 10 shows the TEM micrographs of the fresh and used samples of the promoted 5Ni15YZr at 100 nm magnification. For a wider comprehension of carbon formation and structure, the TEM of the promoted catalysts was typically carried out to elucidate the surface structure and carbon deposition on the surface of the catalyst. The fresh samples of the promoted catalysts (A, B, and C) had better scatter and insignificant active metal diameters. By contrast, the promoted samples used (a’, b’, and c’) exhibited accumulated particles over the surface and formation of the carbon nanotubes.
Figure 11 shows the Raman spectra of the promoted catalysts used. The addition of the promoter altered the structure of the pristine catalyst 5Ni15YZr because the promoter incorporation resulted in the formation of stable ZrO
2 phases and other mixed phases of cubic zirconium promoter oxide. The promoted catalysts displayed three bands with Raman shifts of 816–857 cm
−1, 966–979 cm
−1 and 1112–1115 cm
−1 for the 5NiGa15YZr and 5NiCs15YZr catalysts, which shifted to lower values of 654.0 cm
−1, 885.0 cm
−1, and 1006.0 cm
−1 for NiSr15YZr catalyst. 5Ni15YZr showed a similar peak at 612.0 cm
−1. The peaks for the promoted samples could be associated with the instigated disorders in the graphite (d), the ideal vibration of the graphite layers (g), and the overtone of the d band (d1) [
35]. Commonly, the ratio intensities of the d and g bands denote the amount of graphite relative to carbonaceous materials [
36]. The recorded fractions were 0.840, 0.880, and 0.740 for 5NiGa15YZr, 5NiCs15YZr, and 5NiSr15YZr, respectively, denoting a dominance of graphitization as the ratios are less than unity. The best catalyst, 5NiSr15YZr, presented the lowest ratio.
The catalysts 5Ni15YZr and 5NiSr15YZr gave the best results among all catalysts for the non-promoted and promoted catalysts, respectively. As a result, these catalysts became targets for further investigation in a sensitivity analysis. Typically, in sensitivity analysis experiments, the effects of the flow rate and reaction temperature are studied.
Figure 12 presents the effect of the reaction temperature on the CH
4 and CO
2 conversions of both 5Ni15YZr and 5NiSr15YZr catalysts. As expected, the reaction temperature was directly proportional to the conversion of reactants owing to the highly endothermic feature of the DRM reaction. The effect of reaction temperatures of 500–650 °C was very weak, but the effect was strong and clear when the temperature reached 700 °C; there was a small difference between the reaction temperature results of 700 °C and 750 °C, and the highest effect of the reaction temperature was obtained at 800 °C.
Figure 13 shows the effects of two different flow rates on the CH
4 and CO
2 conversions of both the 5Ni15YZr and 5NiSr15YZr catalysts. The lower flow rate provided better activity performances of CH
4 and CO
2. This could be related to the fact that the low flow rate causes a longer residence time of the reactants on the catalyst.
Figure 13 also shows that the 5NiSr15YZr catalyst maintained better performance than the 5Ni15YZr catalyst for the different flow rates.
Figure 14 displays the TGA analysis of the spent 5Ni15YZr catalyst at different (CH
4/CO) feed flow rates. The percent weight loss for the lower feed flow rate was higher, as expected, because of the high conversion of the reactant and the long interaction time between the feed and the catalyst.
The addition of Y
2O
3 up to 15 wt.% was found to enhance the basicity and optimize the Ni size (13.8 nm). High basicity induced improved CO
2 interaction, whereas size-optimized Ni conveyed controlled CH
4 dissociation. Controlled CH
4 dissociation and sequential carbon deposit oxidation by CO
2 over the 5Ni15YZr catalyst system led to the highest catalytic activity (>50% CH
4 conversion and >60% CO
2 conversion) toward dry reforming of methane. Furthermore, the addition of the alkaline earth metal Sr may have induced additional basicity over the surface. In our previous research, Sr addition was found to upgrade the basicity of the catalyst [
37,
38], which may induce additional CO
2 adsorption over the catalyst’s surface. The reactions of Sr(OH)
2 with CO
2 and SrH
2 with CO
2 have already been established [
39,
40], and the increased interaction of CO
2 over an Sr-modified surface can be accepted. Overall, 15 wt.% Y and 15 wt.% Sr addition over the zirconia-supported Ni catalyst induced proper basicity, which is needed for optimal catalytic performance and minimal coke deposition at a 70 mL/min CH
4: CO
2 flow rate