Coffee Pulp Gasification for Syngas Obtention and Methane Production Simulation Using Ni Catalysts Supported on Al2O3 and ZrO2 in a Packed Bed Reactor

The methanation of CO2 is of great interest in power-to-gas systems and contributes to the mitigation of climate change through carbon dioxide capture and the subsequent production of high-added-value products. This study investigated CO2 methanation with three Ni catalysts supported on Al2O3 and ZrO2, which were simulated using a mathematical model of a packed bed reactor designed based on their chemical kinetics reported in the literature. The simulated reactive system was fed with syngas obtained from residual coffee pulp obtained after a solvent phytochemical extraction process under several gasification conditions. The results reflect a high degree of influence of the catalyst support, preparation method, and syngas composition on CO2 and H2 conversions and CH4 selectivity. For all the syngas compositions, the Ni/ZrO2 catalysts showed the best values for CO2 conversion and H2 conversion for the Ni/Al2O3 catalyst except in gasification at 700 °C and using the Ni/ZrO2p catalyst.


Introduction
Nowadays, the methanation of CO 2 , also known as the Sabatier reaction, is of great interest in power-to-gas systems since it produces a substitute natural gas (SNG) and enables the valorization of CO 2 [1].Furthermore, the resulting methane is a combustible gas and does not present many difficulties in its storage [1,2].As CO 2 is a cause and promoter of climate change, it is highly desirable to develop new processes that capture and/or reuse it [1][2][3] and thus avoid the use of fossil fuels, since their combustion generates more greenhouse gas (GHG) emissions [4,5].Additionally, being a fuel that can result from the reuse of CO 2 , in principle a zero-or low-carbon footprint can be established once it has been burned, as well as biofuels [6].
Moreover, synthesis gas (syngas) feeds, the raw material used for CO 2 methanation, can be obtained from the gasification of renewable sources such as agro-industrial residues and municipal organic waste, among others [7,8].This study investigated syngas obtained from coffee pulp, which is conventionally used as fertilizer, animal feed, and/or low-quality fuel and, if not properly treated, may cause environmental damage [9,10].Therefore, valorizing coffee pulp promotes green chemical synthesis and follows the principles of a circular economy.It also contributes to obtaining high-added-value products in order to base the economy and consumption patterns on more sustainable activities [11].However, to incorporate the biorefinery concept, the coffee pulp considered is a byproduct after a phytochemical extraction process with solvents to take advantage of all its potential to generate high-added-value products [8].
Moreover, for the formation of methane from CO 2 , there must be enough H 2 [3].Therefore, this study aimed to determine the composition of syngas from coffee pulp gasification using steam as the gasifying agent and varying the steam-biomass ratio (S/B) and operation temperature.Then, it considered simulation of the catalytic conversion of syngas into methane in a packed bed reactor.MATLAB © R2023a software was used to solve the mathematical model.Nickel-based catalysts were extensively studied for this reaction due to their high performance/cost ratio [5,[12][13][14], with Ni/Al 2 O 3 as the most widely used at the industrial scale [1,12,13].For this reason, in the present work the selected catalyst was based on Ni supported in Al 2 O 3 and was compared with ZrO 2 as a support, since this latter oxide shows high activity and selectivity, a high thermal stability, and a low carbon deposition rate in CO 2 methanation [15,16].In addition, the considered ZrO 2 catalysts were prepared using two different methodologies to evaluate the influence of these preparation methods [16].Based on chemical kinetic models, a commercial Ni/Al 2 O 3 and two Ni/ZrO 2 catalysts were compared in terms of CO 2 and H 2 conversion and CH 4 selectivity.The evaluation of the reaction system with all compounds in a packed bed reactor simulation for CO 2 methanation allowed for a more comprehensive understanding of the process, identification of potential interactions and by-products, and reactor design optimization.

Syngas from Biomass as a Raw Material for the Methanation Catalytic Process
Gasification is the thermo-chemical conversion of a carbonaceous fuel, and it is characterized as an endothermic process, meaning it requires a heat source [17].Therefore, it is performed at high temperatures typically ranging from 500 to 1400 • C and can take place under atmospheric or elevated pressures, reaching up to 33 bar [18].Gasifiers traditionally operate within distinct configurations, including fixed-bed, fluidized-bed, and moving-bed systems [19].This process involves a partial fuel oxidation by using an oxidizing agent, which could be oxygen, air, steam, or mixtures of these [20,21].Although using air as a gasifying agent is cheap, the syngas dilution by the N 2 presence can lead to a reduction in its high heating value as well as a decrease in the overall efficiency of the gasification process [22].Utilizing steam as the gasifying agent leads to the formation of a syngas characterized by its elevated calorific value, typically ranging from 10 to 15 MJ N/m 3 , and a hydrogen-rich composition [21].Additionally, Ref. [23] highlights the outcomes of syngas generation through biomass gasification with steam.This work demonstrates a yield enhancing of both hydrogen and carbon dioxide while also presenting a notably higher calorific value in comparison to gasification using oxygen or air as gasifying agents.Therefore, steam was selected as gasifying agent in the present study.
The result of gasification is a fuel gas known as syngas.The main components of syngas are carbon monoxide (CO), hydrogen (H 2 ), carbon dioxide (CO 2 ), steam (H 2 O (g) ), methane (CH 4 ), nitrogen (N 2 ) if the oxidation agent is air or a mixture of it with other agents, some hydrocarbons in a very low quantity, and contaminants such as carbon particles, tar, and ash [17,24].In the existing literature, the typical syngas composition derived from agro-industrial waste exhibits a hydrogen content falling within the range of 40% to 50%, carbon monoxide at 9.97% to 12.38%, carbon dioxide levels ranging between 25.04% and 26.50%, and consistent methane levels below 0.5% [25].However, the final syngas composition depends on operational parameters, the type of biomass, and the gasifier configuration [24].
Syngas can be used as a raw material to heat, generate electricity, or synthesize high-added-value chemical and fuel products through several conversion routes, such as methanol or synthetic fuel production [8,11].Therefore, biomass gasification has been considered a viable option for the conversion/utilization of a variety of feedstocks such as vegetable waste, agro-waste, industrial waste, kitchen waste, food waste, and agricultural waste, and even as the key to a successful substitution for petroleum derivatives [26].The conditions for this obtention are described in the Materials and Methods section.

Catalytic Methane Production
In CO 2 methanation, catalysts are needed to achieve high reaction rates, high conversions during CO 2 hydrogenation, and high selectivity toward methane formation [1,14].The more reported active metal phases used in this reaction are Ni, Ru, Rh, and Co [1,5,12].Nickel-based catalysts have been subjected to comprehensive investigation under several reaction conditions due to their relatively low cost and comparatively high catalytic activity [5].In addition, CO 2 methanation reaction is affected by the nature of the catalyst support since it plays an important role in the dispersion of metallic sites, CO 2 adsorption and activation, and metal-support interaction [5,16,27,28].The most common supports for this reaction are based on metal oxides such as Al 2 O 3 , TiO 2 , SiO 2 , ZrO 2 , and CeO 2 [5].Al 2 O 3 -and ZrO 2 -supported nickel catalysts were both reported to be active for CO 2 methanation, showing a better performance than those that used ZrO 2 as a support [15,28].The chemical reactions that can be involved in the methanation process are shown in Equations ( 1)-(3): where Equation ( 1) indicates the chemical reaction of CO 2 methanation; Equation ( 2) is the reverse water gas shift (RWGS) reaction; and Equation ( 3) is CO methanation.The carbon dioxide methanation reaction is exothermic and favored at high pressures and low temperatures [1].The RWS reaction interferes with the selectivity toward methane production.In theory, CO 2 methanation is more favorable with an H 2 /CO 2 ratio equal to or greater than the stoichiometric ratio (4:1) [14].The catalysts employed in the methanation process are highly susceptible to the presence of impurities in the stream.These impurities can lead to catalyst poisoning and deactivation, thereby diminishing the catalytic performance and overall efficiency [4,29].The impurities include chlorine and sulfur compounds, ammonia, tars, and particulate matter [29].The formation of carbon deposits on nickel-based catalysts during the CO methanation reaction has been extensively investigated and can be a significant problem, but it does not pose a concern during the CO 2 methanation reaction [4].Therefore, it is better to use a syngas with higher CO 2 content than CO as a carbon source.Furthermore, the H 2 S presence in the syngas should be analyzed if a real operation is considered.

Ni/ZrO 2 Catalyst
In the present work, the kinetic information was taken from [16].In the latter, the two Ni/ZrO 2 catalysts were tested for methane production by using carbon dioxide hydrogenation.Catalysts were prepared using the incipient wetness impregnation method.According to [16], the impregnated sample was calcinated to obtain NiO/ZrO 2 -C at 500 • C for 3 h, while NiO/ZrO 2 -P was synthesized using dielectric barrier discharge (DBD) plasma for 1 h.Subsequently, the oxidized metal catalysts were subjected to hydrogen reduction at elevated temperatures to form Ni/ZrO 2 -C and Ni/ZrO 2 -P.These catalysts were reduced in situ with pure H 2 (20 mL/min) at 500 • C for 1 h.
The chemical kinetics proposed in that study did not include a model such as Langmuir-Hinshelwood-Hougen-Watson kinetics, which explicitly considers mass transfer in chemical kinetics.However, simple power-law kinetics for direct reactions can be used for technical applications [30].According to experimental data reported in [16], the CH 4 selectivity for both prepared ZrO 2 catalysts is close to 100%.Therefore, it only considers Equation (1) to develop the kinetic model.Then, the kinetic model of Ni/ZrO 2 catalysts follows a power-law model as shown in Equation (4).

−r
where −r CO2 is the kinetic rate for CO 2 ; A is a pre-exponential factor related to the chemical kinetics in (L/g•h); R is the ideal constant of gases (8.314 × 10 −3 kJ/mol); R is the ideal constant of gases (0.082 atm•L/K•mol); Ea is the activation energy in kJ/mol; T is the reaction temperature in K, and P CO2, and P H2 are the partial pressure of the CO 2 and H 2 in atm, respectively.Table 1 shows the values of the chemical kinetic model of the two Ni/ZrO 2 catalysts.
where −r i is the kinetic rate for each reaction; Ea is the activation energy in kJ/mol; Pi is the partial pressure of component i in bar; Ki is the adsorption constant of component i; and Keq, i is the chemical equilibrium constant for reaction i.Tables 2 and 3 show the adsorption and kinetic parameters, respectively, of this catalyst model.

Biomass Chemical Composition
The chemical composition of the residual coffee pulp was expressed by its elemental analysis.Table 4 shows the C, N, H, O, and S content in the sample on a dry basis.The coffee pulp showed a low hydrogen content, which recommended using a gasifying agent, such as steam, capable of supplying an increment to the syngas H 2 concentration [31].Furthermore, the sulfur content was very low or negligible.However, the obtained syngas showed low quantities of H 2 S, which indicated that in the elemental analysis, the sulfur was below the detection limit.Therefore, it is possible the sulfur content was below the limit detection.Elemental analysis of the coffee pulp showed a similar CHNO ratio than other residual biomasses reported in [32].This study suggests that the elemental analysis composition is attributed to presence of lignin and cellulose in the biomass, which is a typical feature of agro-industrial waste.

Syngas Description
We performed the gasification of the coffee pulp at two temperatures and two steam/biomass (S/B) ratios.The compositions of the obtained syngas in %v/v are shown in Table 5. S/B ratios were used because steam was the gasifying agent selected in this work.Assuming an ideal gas behavior for syngas, the composition in %v/v is equivalent to the molar fraction.According to the results described in Table 5, the formation of CO 2 was favored under conditions where the S/B ratio was 0.5 and 1.0 at 700 • C.However, there was no significant variation in temperature that allowed establishing a better operating condition for this gasification parameter.However, at a gasification temperature of 700 • C, a slightly higher amount of CO 2 was achieved for an S/B ratio of 0.5.On the other hand, higher amounts of H 2 were achieved, as for CO 2 , for both S/B ratios at 700 • C. Additionally, upon analyzing the considered reactions, CO could contribute to methane formation according to reaction B. Therefore, according to the data in Table 5, better production was achieved for both S/B ratios and a temperature of 800 • C.However, at this temperature, the production of CO 2 was not high as at 700 • C. The H 2 /CO 2 ratios were greater than the stochiometric ratio for reaction 1 (4:1) when a gasification temperature of 800 • C for both S/B ratios was implemented.Furthermore, with an S/B of 0.5, a better value for this relation was achieved.
According to [21], raising the gasification temperature leads to a reduction of chemical species that can poison and enhance catalyst deactivation.However, this was true for our own gasification results when an S/B of 1.0 was used, because at an S/B ratio of 0.5, the behavior was contrary to what was established by these references.Therefore, the S/B ratio could affect the formation of H 2 S since there was a clear variation in the formation of this poison agent when comparing both gasification temperatures.To obtain synthesis gases from coffee pulp with the most adequate composition to avoid the possible acceleration of catalyst poisoning and deactivation, a temperature of 800 • C and an S/B of 1.0 should be selected because they presented the lowest value for H 2 S formation in the syngas.

Catalytic Simulation of Methanation
Once we obtained the syngas compositions after gasification, we analyzed the CO 2 conversion rate and H 2 conversion throughout the weight catalytic bed (w) of each one of the syngas compositions and catalysts while considering the kinetics described above.The results are shown in Figures 1 and 2. The nomenclature used in all figures with a legend contain the catalyst type followed by the S/B ratio and gasification temperature.
According to [21], raising the gasification temperature leads to a reduction of chemical species that can poison and enhance catalyst deactivation.However, this was true for our own gasification results when an S/B of 1.0 was used, because at an S/B ratio of 0.5, the behavior was contrary to what was established by these references.Therefore, the S/B ratio could affect the formation of H2S since there was a clear variation in the formation of this poison agent when comparing both gasification temperatures.To obtain synthesis gases from coffee pulp with the most adequate composition to avoid the possible acceleration of catalyst poisoning and deactivation, a temperature of 800 °C and an S/B of 1.0 should be selected because they presented the lowest value for H2S formation in the syngas.

Catalytic Simulation of Methanation
Once we obtained the syngas compositions after gasification, we analyzed the CO2 conversion rate and H2 conversion throughout the weight catalytic bed (w) of each one of the syngas compositions and catalysts while considering the kinetics described above.The results are shown in Figures 1 and 2. The nomenclature used in all figures with a legend contain the catalyst type followed by the S/B ratio and gasification temperature.At 800 °C gasification for both S/B ratios, the feed syngas streams contained a lower concentration of CO2, resulting in a decreased reactant-to-catalyst ratio; in other words, it had a lower space velocity.This condition led to higher conversion rates for these gasification process parameters.To obtain an accurate comparison, it was necessary to evaluate the results obtained at 700 °C and those obtained at 800 °C separately for all simulated catalysts.As shown in Figures 1 and 2, Ni/ZrO2-P obtained better results than Ni/ZrO2-C in terms of CO2 conversion and H2 conversion.In accordance with [16], this behavior was due to a greater number of active sites, which could be provided by highly dispersed Ni particles on Ni/ZrO2-P, leading to a higher reaction rate for CO2 methanation.The surface At 800 • C gasification for both S/B ratios, the feed syngas streams contained a lower concentration of CO 2 , resulting in a decreased reactant-to-catalyst ratio; in other words, it had a lower space velocity.This condition led to higher conversion rates for these gasification process parameters.To obtain an accurate comparison, it was necessary to evaluate the results obtained at 700 • C and those obtained at 800 • C separately for all simulated catalysts.As shown in Figures 1 and 2, Ni/ZrO 2 -P obtained better results than Ni/ZrO 2 -C in terms of CO 2 conversion and H 2 conversion.In accordance with [16], this behavior was due to a greater number of active sites, which could be provided by highly dispersed Ni particles on Ni/ZrO 2 -P, leading to a higher reaction rate for CO 2 methanation.The surface Ni concentration of Ni/ZrO 2 -P, estimated at Ni/Zr = 0.18, was higher than that for the Ni/ZrO 2 -C catalyst (Ni/Zr = 0.13).Therefore, this indicated a better Ni dispersion on the support induced by plasma decomposition.In addition, it has been proposed that this method of preparation facilitates the partial reduction of ZrO 2 to create more oxygen vacancies, thus improving CO 2 activation [16].For these ZrO 2 supported catalysts, at 700 • C it is better to use an S/B ratio of 1.0 to obtain high conversions for CO 2 .However, for the same temperature, the H 2 conversion is best when the Ni/ZrO 2 -C catalyst uses a syngas for an S/B of 1.0 and for Ni/ZrO 2 -P, an S/B of 0.5.The highest H 2 conversion was produced at 800 • C with an S/B ratio of 1.0 and CO 2 conversion at an S/B of 0.5.
Ni concentration of Ni/ZrO2-P, estimated at Ni/Zr = 0.18, was higher than that for the Ni/ZrO2-C catalyst (Ni/Zr = 0.13).Therefore, this indicated a better Ni dispersion on the support induced by plasma decomposition.In addition, it has been proposed that this method of preparation facilitates the partial reduction of ZrO2 to create more oxygen vacancies, thus improving CO2 activation [16].For these ZrO2 supported catalysts, at 700 °C it is better to use an S/B ratio of 1.0 to obtain high conversions for CO2.However, for the same temperature, the H2 conversion is best when the Ni/ZrO2-C catalyst uses a syngas for an S/B of 1.0 and for Ni/ZrO2-P, an S/B of 0.5.The highest H2 conversion was produced at 800 °C with an S/B ratio of 1.0 and CO2 conversion at an S/B of 0.5.Regarding the results for the commercial Ni/Al2O3 catalyst in Figure 1c, it achieved the highest value of CO2 conversion with a syngas obtained using an S/B ratio of 1.0 for both temperatures in the gasification process.According to Figure 2c, the H2 conversion behavior of this support was opposite for the mentioned S/B gasification condition.
Figures 3 and 4 show the comparisons of CO2 and H2 conversions, respectively, obtained with the different catalysts at the reactor outlet using all the syngas compositions.A positive and negative value indicate a stronger influence of the variable in the first and second position in the comparison, respectively.For instance, in the comparison between ZrO2p and Al2O3, positive values indicate a higher conversion of ZrO2p, while negative values indicate a greater conversion of Al2O3.
Figure 3 shows that the highest values of CO2 conversion were achieved with ZrO2.This result is in accordance with [15], which reported that the catalytic system exhibited a higher activity with ZrO2 supports than with Al2O3.ZrO2p shows better conversions than ZrO2c for all syngas compositions.Additionally, it could be attributed to a pre-exponential factor, and the reaction order was greater for Ni/ZrO2p than for Ni/ZrO2c.Therefore, the CO2 methanation must be faster for the Ni/ZrO2p catalyst.On the other hand, Figure 4 shows that H2 conversion had the most prevalent negative values for Ni/ZrO2C compared to Ni/Al2O3.Regarding the results for the commercial Ni/Al 2 O 3 catalyst in Figure 1c, it achieved the highest value of CO 2 conversion with a syngas obtained using an S/B ratio of 1.0 for both temperatures in the gasification process.According to Figure 2c, the H 2 conversion behavior of this support was opposite for the mentioned S/B gasification condition.
Figures 3 and 4 show the comparisons of CO 2 and H 2 conversions, respectively, obtained with the different catalysts at the reactor outlet using all the syngas compositions.Figure 3 shows that the highest values of CO 2 conversion were achieved with ZrO 2 .This result is in accordance with [15], which reported that the catalytic system exhibited a higher activity with ZrO 2 supports than with Al 2 O 3 .ZrO 2 p shows better conversions than ZrO 2 c for all syngas compositions.Additionally, it could be attributed to a pre-exponential factor, and the reaction order was greater for Ni/ZrO 2 p than for Ni/ZrO 2 c.Therefore, the CO 2 methanation must be faster for the Ni/ZrO 2 p catalyst.On the other hand, Figure 4 shows that H 2 conversion had the most prevalent negative values for Ni/ZrO 2C compared to Ni/Al 2 O 3 .Since the power-law kinetic model for the Ni/ZrO2 catalysts only included Equation (1) for CO2 methanation, it was not possible to establish the CO selectivity.However, chemical kinetic model for the commercial Ni/Al2O3 catalyst considered all reactions.Therefore, this catalyst allowed us to calculate the selectivity for CO, to determine the performance of gasification conditions in Sabatier process to increase the CH4 yield, and to minimize the by-product generation.Figure 5 illustrates the CH4 and CO selectivity for the commercial Ni/Al2O3 catalyst when considering all the gasification conditions.
Figure 5 shows a better behavior of the Ni/Al2O3 catalyst for CO selectivity at 800 °C when an S/B of 0.5 was applied.At 700 °C, the same effect is shown for an S/B of 1.0.Furthermore, with these gasifying conditions, the CH4 selectivity was lower.Despite achieving an H2/CO2 ratio of less than 4 under these gasification conditions, the   Since the power-law kinetic model for the Ni/ZrO2 catalysts only included Equation (1) for CO2 methanation, it was not possible to establish the CO selectivity.However, chemical kinetic model for the commercial Ni/Al2O3 catalyst considered all reactions.Therefore, this catalyst allowed us to calculate the selectivity for CO, to determine the performance of gasification conditions in Sabatier process to increase the CH4 yield, and to minimize the by-product generation.Figure 5 illustrates the CH4 and CO selectivity for the commercial Ni/Al2O3 catalyst when considering all the gasification conditions.
Figure 5 shows a better behavior of the Ni/Al2O3 catalyst for CO selectivity at 800 °C when an S/B of 0.5 was applied.At 700 °C, the same effect is shown for an S/B of 1.0.Furthermore, with these gasifying conditions, the CH4 selectivity was lower.Despite achieving an H2/CO2 ratio of less than 4 under these gasification conditions, the Since the power-law kinetic model for the Ni/ZrO 2 catalysts only included Equation (1) for CO 2 methanation, it was not possible to establish the CO selectivity.However, chemical kinetic model for the commercial Ni/Al 2 O 3 catalyst considered all reactions.Therefore, this catalyst allowed us to calculate the selectivity for CO, to determine the performance of gasification conditions in Sabatier process to increase the CH 4 yield, and to minimize the by-product generation.Figure 5 illustrates the CH 4 and CO selectivity for the commercial Ni/Al 2 O 3 catalyst when considering all the gasification conditions.
Figure 5 shows a better behavior of the Ni/Al 2 O 3 catalyst for CO selectivity at 800 • C when an S/B of 0.5 was applied.At 700 • C, the same effect is shown for an S/B of 1.0.Furthermore, with these gasifying conditions, the CH 4 selectivity was lower.Despite achieving an H 2 /CO 2 ratio of less than 4 under these gasification conditions, the concentration of CO 2 was higher than that of CO.This predominance of CO 2 concentration plays a crucial role in methane (CH 4 ) production.According to the adsorption parameters of the Al 2 O 3 catalyst given in Table 2, it was observed that CO 2 adsorption was stronger compared to CO.However, even though the kinetic constant for CO methanation was greater and the activation energy was lower than that for CO 2 methanation, this result implies that the rate-limiting step in methane production using the Ni/Al 2 O 3 catalyst was the adsorption of reactants.According to the results, a stronger influence of the S/B ratio was observed for CH 4 selectivity for both employed gasification temperatures.concentration of CO2 was higher than that of CO.This predominance of CO2 concentration plays a crucial role in methane (CH4) production.According to the adsorption parameters of the Al2O3 catalyst given in Table 2, it was observed that CO2 adsorption was stronger compared to CO.However, even though the kinetic constant for CO methanation was greater and the activation energy was lower than that for CO2 methanation, this result implies that the rate-limiting step in methane production using the Ni/Al2O3 catalyst was the adsorption of reactants.According to the results, a stronger influence of the S/B ratio was observed for CH4 selectivity for both employed gasification temperatures.

Materials and Methods
This section starts with a detailed description of the biomass elemental analysis and gasification procedure for syngas production from coffee pulp after a solvent phytochemical extraction to obtain chlorogenic acid.These residual coffee pulps were taken as a solid filtered after the extraction process using ethanol/H2O (70/30 %V/V) at ambient temperature, atmospheric pressure, and a coffee pulp-to-solvent ratio of 1 to 4.Then, a theorical kinetic study is considered for methane production using Ni-based catalysts on ZrO2 and Al2O3 supports on a packed bed reactor.It presents a mathematical model that includes the assumptions and all the equations that are necessary to simulate a packed bed reactor considering a heterogenous catalyst for the reactive system.Finally, it details the CO2 conversion, H2 conversion, and CH4 selectivity of the model constructed here when it was simulated with each catalyst using its respective chemical kinetic parameters found in the cited papers.

Elemental Analysis
The elemental analysis of the coffee pulp sample was determined using a LECO CHNS analyzer (Truspec micro model).The procedure was developed following the ASTM D-5373-08 method.The CHN analysis was conducted at 1050 °C and the sulfur analysis at 1350 °C, both in a helium atmosphere.The results are reported on a dry basis.Since the sulfur content was negligible, the oxygen content was calculated according to the difference [32].

Materials and Methods
This section starts with a detailed description of the biomass elemental analysis and gasification procedure for syngas production from coffee pulp after a solvent phytochemical extraction to obtain chlorogenic acid.These residual coffee pulps were taken as a solid filtered after the extraction process using ethanol/H 2 O (70/30 %v/v) at ambient temperature, atmospheric pressure, and a coffee pulp-to-solvent ratio of 1 to 4.Then, a theorical kinetic study is considered for methane production using Ni-based catalysts on ZrO 2 and Al 2 O 3 supports on a packed bed reactor.It presents a mathematical model that includes the assumptions and all the equations that are necessary to simulate a packed bed reactor considering a heterogenous catalyst for the reactive system.Finally, it details the CO 2 conversion, H 2 conversion, and CH 4 selectivity of the model constructed here when it was simulated with each catalyst using its respective chemical kinetic parameters found in the cited papers.

Elemental Analysis
The elemental analysis of the coffee pulp sample was determined using a LECO CHNS analyzer (Truspec micro model).The procedure was developed following the ASTM D-5373-08 method.The CHN analysis was conducted at 1050 • C and the sulfur analysis at 1350 • C, both in a helium atmosphere.The results are reported on a dry basis.Since the sulfur content was negligible, the oxygen content was calculated according to the difference [32].

Gasification Experimental System
The residual coffee pulp was subjected to gasification using a horizontal furnace reactor.The reactor comprised a quartz tube with an external diameter of 3.5 cm (internal diameter of 3.0 cm) and a length of 50 cm.The quartz tube was inserted into the annular space of the 40 cm long horizontal furnace.The furnace was heated by electric resistances, and its temperature was controlled through a proportional-integral-derivative (PID) loop, enabling continuous monitoring of the reactor temperature.To facilitate the gasification process, the sample was loaded into a quartz sample holder that was 2.0 cm in diameter and 13.5 cm in length within the reactor.Inside the sample holder, a K-type thermocouple was inserted to accurately measure the temperature of the sample bed.Scheme 1 shows a gasifier diagram employed in this study.

Gasification Experimental System
The residual coffee pulp was subjected to gasification using a horizontal furnace reactor.The reactor comprised a quartz tube with an external diameter of 3.5 cm (internal diameter of 3.0 cm) and a length of 50 cm.The quartz tube was inserted into the annular space of the 40 cm long horizontal furnace.The furnace was heated by electric resistances, and its temperature was controlled through a proportional-integral-derivative (PID) loop, enabling continuous monitoring of the reactor temperature.To facilitate the gasification process, the sample was loaded into a quartz sample holder that was 2.0 cm in diameter and 13.5 cm in length within the reactor.Inside the sample holder, a K-type thermocouple was inserted to accurately measure the temperature of the sample bed.Scheme 1 shows a gasifier diagram employed in this study.The gasification experiments were conducted at two different temperatures: 700 °C and 800 °C.Additionally, two steam-to-biomass (S/B) ratios were tested, namely 0.5 and 1.0, to evaluate their effects on the gasification process.For each S/B ratio and temperature, seven experiments were conducted with variations in the gasification time.The first sample was analyzed at 6 min of the gasification process, the second was at 10 min, and subsequent analyses were conducted every 10 min until reaching 60 min of gasification.After that, the arithmetic average of these seven measurements was taken as the result for the syngas composition.
The produced syngas was analyzed using an Agilent Micro GC model 3000.This analytical instrument was equipped with two thermal conductivity detectors (TCD).One TCD utilized a 10 m × 0.32 mm 5A molecular sieve column with argon as the carrier gas, while the other TCD employed an 8 m × 0.32 mm column with helium as the carrier gas.

Mathematical Model for the Methanation Catalytic Packed Bed Reactor Simulation
This subsection discusses the simulation of a packed bed reactor operating in steadystate conditions.In the simulation, the reactor considered a catalyst weight of 100 mg in a packed bed configuration, 1 atm of pressure, a total inlet molar flux of 0.1 mol/min, ambient pressure, and 400 °C.This temperature was selected since in previous studies, at the optimal temperature range of 350 to 450 °C, Ni/Al2O3 catalyst systems synthesized through conventional methods typically exhibited CO2 conversion rates ranging from 50% to 80% and CH4 selectivity exceeding 90% [14].Another reason to choose this temperature The gasification experiments were conducted at two different temperatures: 700 • C and 800 • C. Additionally, two steam-to-biomass (S/B) ratios were tested, namely 0.5 and 1.0, to evaluate their effects on the gasification process.For each S/B ratio and temperature, seven experiments were conducted with variations in the gasification time.The first sample was analyzed at 6 min of the gasification process, the second was at 10 min, and subsequent analyses were conducted every 10 min until reaching 60 min of gasification.After that, the arithmetic average of these seven measurements was taken as the result for the syngas composition.
The produced syngas was analyzed using an Agilent Micro GC model 3000.This analytical instrument was equipped with two thermal conductivity detectors (TCD).One TCD utilized a 10 m × 0.32 mm 5A molecular sieve column with argon as the carrier gas, while the other TCD employed an 8 m × 0.32 mm column with helium as the carrier gas.

Mathematical Model for the Methanation Catalytic Packed Bed Reactor Simulation
This subsection discusses the simulation of a packed bed reactor operating in steadystate conditions.In the simulation, the reactor considered a catalyst weight of 100 mg in a packed bed configuration, 1 atm of pressure, a total inlet molar flux of 0.1 mol/min, ambient pressure, and 400 • C.This temperature was selected since in previous studies, at the optimal temperature range of 350 to 450 • C, Ni/Al 2 O 3 catalyst systems synthesized through conventional methods typically exhibited CO 2 conversion rates ranging from 50% to 80% and CH 4 selectivity exceeding 90% [14].Another reason to choose this temperature was because the reaction rate is relatively slow at low temperatures.However, CH 4 selectivity is higher at these operating temperatures [15].Therefore, these conditions were maintained constant in all the simulations to be able to compare the behavior of the catalysts using the syngas compositions described above.The following considerations and simplifications proposed in [33,34] were implemented for this type of operation and reactor:

•
Negligible radial diffusion: concentration and temperature profiles were assumed to be constants, which led to a one-dimensional model.

•
Temperature and pressure profiles in the catalyst were assumed to be constants (homogeneous catalytic particles).

•
As in [13], the mechanisms related to catalyst deactivation, such as sulfur poisoning or carbon formation via the Boudouard reaction, were not taken into consideration or disregarded in the present study.
As the process was carried out in a packed bed reactor, the resulting model had to be adjusted to the design equation of this type of reactor, which can be consulted in [13].Equations ( 8)-( 12) are the system of ordinary differential equations (ODEs) used here for the mass balance and reactor design.They express each reactive species involved in the Sabatier process.
dF H2 dW = r H2 (10) dF CH4 dW = r CH4 (11) dF H2O dW = r H2O (12) where Fi denotes the molar flux of species i; r i is the chemical kinetics of species i; and W is the catalyst weight inside the heterogeneous reactor in grams.Since the system had three chemical reactions in parallel; Equation (1) CO 2 methanation, Equation (2) RWGS, and Equation (3) CO methanation), it needed to consider the species that appeared in more than one of these reactions because the chemical kinetics of each of the said reactions were coupled and given in parallel as well.Equations ( 13)- (17) show the global or total reactions for each species involved in the process described above, where each r i is expressed in mol/min•g cat .
r CO = r RWGS − r CO−met (17) Catalyst performance was analyzed here by measuring the CO 2 and H 2 conversion and the selectivity of CH 4 and CO.CO 2 and H 2 conversion are described in Equations ( 18) and ( 19), and the selectivity of CO to CH 4 are described in Equations ( 20) and (21).
where X CO2 and X H2 are the CO 2 and H 2 conversion inside the reactor and S CO2/CH4 and S CO2/CO are the selectivity of CO 2 to transform it into CH 4 and CO, respectively; and F i_in and F i_out are the molar flux of components according to the subindex in mol/min.Most of the kinetic expressions studied and analyzed here depend on the partial pressures of the chemical species involved, and the system of differential equations was based on molar balance.Therefore, Equation ( 22) must be used to relate the partial pressures to the molar fluxes of each species.
where P i is the partial pressure of species i inside the reactor; F i is the molar flux of component i (mol/min); ∑ n i=1 F i is the summation of all molar fluxes of the species involved (mol/min); and P T is the total system pressure.

Conclusions
This paper presented a coffee pulp gasification after a phytochemical extraction process with varying S/B ratios (0.5 and 1.0) and temperatures (700 • C and 800 • C).The resulting syngas was considered as a raw material from methane production through the Sabatier process.The adequate syngas according to the H 2 /CO 2 ratio was achieved at 800 • C and both S/B ratios because these conditions exceeded the stoichiometric ratio (4:1) for the CO 2 methanation reaction.Additionally, the temperature of 800 • C reduced the formation of H 2 S, which favored avoiding the acceleration of catalyst poisoning.However, in the present study, this effect was not considered in the catalyst simulation.
Then, a comparative analysis of three Ni catalysts in terms of CO 2 methanation performance and behavior was considered.Such analysis was based on a mathematical model that simulated their chemical kinetics reported in their respective references using MATLAB © R2023a.The chemical kinetics of the two simulated Ni/ZrO 2 catalysts were taken from [16], in which they were prepared by adopting two different methods: one was obtained using dielectric barrier discharge (DBD) plasma, and the other by calcination.In turn, the kinetics of the Ni/Al 2 O 3 commercial catalyst were taken from [1].The present study established CO 2 conversion, CH 4 selectivity, and H 2 conversion in CO 2 methanation depending on the catalyst preparation method, catalyst support, raw material composition, and gasification process conditions.
According to the results, Ni/ZrO 2 -P had a better catalytic behavior than Ni/ZrO 2 -C, which corroborates that the catalyst preparation method influences catalyst activity.Also, the catalyst support affected CO 2 conversion and H 2 conversion.ZrO 2 -supported catalysts showed better CO 2 conversion, but the Ni/Al 2 O 3 commercial catalyst had a higher H 2 conversion.However, the Ni/ZrO 2 p catalyst presented a better H 2 conversion than the commercial Ni/Al 2 O 3 at 700 • C and both S/B ratios.
Regarding the syngas composition, Ni/Al 2 O 3 was more active for CO 2 conversion and selective for CH 4 when an S/B ratio of 1.0 for both temperatures was employed.In contrast, the H 2 conversion was better for an S/B of 0.5.
Finally, the model presented in this study has the potential to be applied to evaluate the CO 2 methanation process using different types of syngas derived from various biomass sources and operating conditions during the gasification process.By utilizing the same catalysts studied in this paper, researchers can conduct theoretical exploratory research to quickly estimate the impact of different variables on CO 2 conversion, H 2 conversion, and CH 4 selectivity.This approach offers the advantage of saving both time and resources by providing a valuable method to identify the variables in the gasification process that promote favorable CO 2 conversion, CH 4 selectivity, and H 2 conversion.
A positive and negative value indicate a stronger influence of the variable in the first and second position in the comparison, respectively.For instance, in the comparison between ZrO 2 p and Al 2 O 3 , positive values indicate a higher conversion of ZrO 2 p, while negative values indicate a greater conversion of Al 2 O 3 .

Figure 3 .
Figure 3.Comparison in percentage (%) of CO2 conversion for the catalysts under study.

Figure 4 .
Figure 4. Comparison in percentage (%) of H2 conversion for the catalysts under study.

Figure 3 .
Figure 3.Comparison in percentage (%) of CO 2 conversion for the catalysts under study.

Figure 3 .
Figure 3.Comparison in percentage (%) of CO2 conversion for the catalysts under study.

Figure 4 .
Figure 4. Comparison in percentage (%) of H2 conversion for the catalysts under study.

Figure 4 .
Figure 4. Comparison in percentage (%) of H 2 conversion for the catalysts under study.

Figure 5 .
Figure 5. Selectivity for commercial Ni/Al2O3 catalyst considering all the gasification conditions (S.B0.5T700 means a syngas produced using an S/B ratio of 0.5 and a gasification temperature of 700 °C).

Figure 5 .
Figure 5. Selectivity for commercial Ni/Al 2 O 3 catalyst considering all the gasification conditions (S.B0.5T700 means a syngas produced using an S/B ratio of 0.5 and a gasification temperature of 700 • C).

Table 5 .
Composition of syngas (%v/v) obtained from coffee pulp waste at different temperatures and steam/biomass (S/B) ratios.