Original Research ARTICLE
Influence of the Catalyst Particle Size on the Aqueous Phase Reforming of n-Butanol Over Rh/ZrO2
- 1Department of Chemical and Metallurgical Engineering, Aalto University, Espoo, Finland
- 2Catalytic Processes and Materials, Department of Science and Technology, MESA+ Institute for Nanotechnology, University of Twente, Enschede, Netherlands
- 3Department of Applied Science and Technology, Politecnico di Torino, Turin, Italy
Butanol is a by-product obtained from biomass that can be valorized through aqueous phase reforming. Rh/ZrO2 catalysts were prepared and characterized, varying the size of the support particles. The results showed a relatively mild effect of internal mass transport on butanol conversion. However, the influence of internal transport limitations on the product distribution was much stronger, promoting consecutive reactions, i.e., dehydrogenation, hydrogenolysis, and reforming of propane and ethane. Hydrogen consuming reactions, i.e., hydrogenolysis, were more strongly enhanced than hydrogen producing reactions due to internal concentration gradients. Large support particles deactivated faster, attributed to high concentrations of butyraldehyde inside the catalyst particles, enhancing deposit formation via aldol condensation reactions. Consequently, also the local butyric acid concentration was high, decreasing the local pH, enhancing Rh leaching. The influence of internal transfer limitation on product distribution and stability is discussed based on a reaction scheme with three main stages, i.e., (1) formation of liquid intermediates via dehydrogenation, (2) formation of gas via decarbonylation/decarboxylation reactions, and (3) hydrocarbon hydrogenolysis/reforming/dehydrogenation.
The environmental issues and the depletion of conventional sources of energy demand development of alternative and sustainable technologies. Among the several possibilities, biomass is seen as a strategic feedstock for the production of renewable energy and materials. One of the possible products of biomass exploitation is hydrogen, through thermochemical or biological routes (Balat and Kirtay, 2010).
In the last years, a considerable effort has been put on the production of hydrogen from oxygenated hydrocarbons, e.g., via aqueous phase reforming (APR) (Cortright et al., 2002). The Dumesic research group demonstrated that hydrogen can be produced from alcohols in water in the condensed phase, with a noteworthy energetic advantage compared to the conventional steam reforming as evaporation of water is circumvented (Davda et al., 2005).
APR is a promising strategy for valorization of aqueous side-streams. Among the possible reactants for this process, oxygenates with 1 to 1 O to C ratios are preferred for H2 production, via reforming and subsequent water gas shift reaction (Cortright et al., 2002; Shabaker et al., 2003a; Davda et al., 2005). Indeed, methanol is the most investigated mono-alcohol, thanks to its optimal carbon/oxygen ratio.
However, little attention has been paid so far to APR of butanol, despite its strategic importance as it can be produced from biomass via fermentation of sugar cane (Kumar et al., 2017). Also, aqueous waste streams of e.g., flash pyrolysis contain butanol, making it an interesting model compound for study. As a matter of fact, butanol has been studied for hydrogen production vis supercritical water reforming for its representativeness of oxygenates present in the bio-oil aqueous phase (Gutiérrez Ortiz et al., 2016; Gutiérrez Ortiz and Campanario, 2018).
Roy et al. investigated for the first time APR of butanol over Ni-based catalysts supported on ceria or alumina (Roy et al., 2011). Successively, the same research group used these catalysts also in harsher conditions, i.e., for steam reforming of butanol, enlarging the range of operating conditions as reaction temperatures, pressure, concentration, and flow rate of the feed (Roy et al., 2014).
In previous work we investigated the steam reforming of butanol over a Rh/ZrO2 catalyst (Harju et al., 2015, 2016). This catalytic system showed promising results in terms of hydrogen productivity and stability, since coke formation is slower compared to other supports. APR of higher alcohols results in significant formation of hydrocarbons on more conventional metal catalysts (Roy et al., 2011, 2014; Lobo et al., 2012) because of limited C-C cleavage. Rh however is known to be active for C-C cleavage (Sinfelt, 1973; Bond et al., 1996) and low temperature steam reforming (Kolb et al., 2004; Halabi et al., 2010). It has been recently showed that, among several metals (Pd, Ru, Re, Ir, and Cr) Rh was the best promoter in a bimetallic Pt-based catalyst for the APR of glycerol (Larimi and Khorasheh, 2019). Furthermore, ZrO2 is known as one of the few oxides capable of resisting the harsh hydrothermal conditions of APR (Elliott et al., 1993), unlike more commonly used supports (De Vlieger et al., 2012). Therefore, this study explores the use of Rh/ZrO2 for the valorization of butanol in APR conditions.
APR is a three phase (G-L-S) system and therefore issues related to mass transfer limitations may arise. Hydrogen mass transfer is of paramount importance in APR as reported by Neira D'Angelo et al. (2013, 2014a,b) reporting that microchannel reactors enhance mass transfer, increasing the hydrogen yield by suppressing sequential reactions consuming H2.
The goal of this work is to determine the effect of internal mass transfer on aqueous phase reforming of butanol over Rh/ZrO2 catalyst by varying the dimension of the catalyst particles using a reactor design that increase the sensitivity for internal mass transfer. For example, in contrast to previous works (Neira D'Angelo et al., 2013, 2014a), the reactor was not flushed with inert gas, so that the concentration of gas products outside the catalyst particles was relatively high, retarding diffusion out of gaseous products. The obtained results allow not only to discuss the effect of internal mass transport on product distribution and deactivation, but also to propose a reaction scheme.
Materials and Methods
The catalysts were prepared by vacuum assisted dry impregnation method followed by calcination at 850°C, described in more detail in earlier work (Harju et al., 2015). To obtain the different catalyst particle sizes, the ZrO2 support (monoclinic ZrO2, MEL Chemicals) was crushed and sieved to the desired particle sizes prior to impregnation. The support particles sizes remained unchanged during catalyst preparation.
The Rh content of both fresh and spent catalyst was determined by X-Ray Fluorescence performed with Malvern Panalytical Axios mAX 3 kW. The 7-point BET surface area of both fresh and spent catalyst was performed with Coulter Omnisorp 100 CX, using methods described elsewhere (Kaila et al., 2007, 2008). Pore volume of both fresh and spent catalyst were determined during the same N2 physisorption measurement as the BET surface area, using the NLDFT method, assuming spherical and cylindrical pores. Unfortunately, the equipment used for BET measurement did not allow experiments with small particles and only the samples with the largest particle size could be measured. The Rh dispersion on both fresh and spent catalyst was determined using pulsed H2 chemisorption in Ar (Chemisorb 2750, Micromeritics). The Rh particle size distribution on spent catalyst was determined by scanning tunneling electron microscopy (STEM) using Jeol 2200FS equipped with spherical aberration corrector. Elemental C, H, and O analysis of the deposits on spent catalyst was done using methods described in previous work (Harju et al., 2016). The elemental analysis was also done for fresh catalyst to determine the contribution of hydroxyl and carbonate groups of atmospheric origin.
A schematic representation of the set-up used for APR experiments is reported in Figure 1. The reactant solution (5 wt% n-butanol in de-ionized water) was pressurized and fed by a setup of dual ISCO D-pumps and pre-heated up to the desired reaction temperature (220°C) before the catalytic bed. The catalytic bed was 7 mm in diameter. The catalyst loading in the reactor was varied (1.5, 1.0 and 0.1 g for the 250–420, 60–100, and 40–60 μm catalysts, respectively, resulting is acceptable bed height), keeping the LHSV constant at 150 h−1 for all the experiments by adjusting the volumetric flow rate. The set temperature was reached about 30 min after starting heating, and it was considered in the results as the zero time. Afterwards the effluent was cooled down to room temperature by a water-cooled tubular heat exchanger. The pressure was kept constant at ~35 bar by a manual back-pressure regulator (TESCOM, model 26-1764-24-090) following the heat exchanger. Prior to phase separation at ambient pressure in a separation vessel, nitrogen sweep gas was introduced to the stream to aid in gas product purging and to keep a steady minimum flow present to improve accuracy of both gas phase analysis and flow rate measurement. The gaseous products were periodically sampled from the top of the de-mister section of the phase separation vessel for analysis in a Varian micro-GC, while the gas flowrate was measured by a gas mass flow meter. The liquid flowrate was integrally measured by weighing the eluent collection vessel mounted on a scale.
Figure 1. Schematic of the experimental set-up. 1. Piston pumps. 2. Pre-heater. 3. Heater. 4. Catalyst section. 5. Heat exchanger. 6. Back-pressure valve. 7. Gas-liquid separation vessel. 8. Varian Micro-GC. 9. Eluent collection vessel. 10. Gas mass-flow meter.
The time on stream (TOS) was 150 min for the 40–60 and 60–100 μm particles and 210 min for the 250–420 μm particles. After cooling the system overnight, the catalyst was removed from the reactor and dried at 100°C overnight.
Product Analysis and Calculation
The gas products were analyzed with a Varian micro-GC equipped with two columns with a TCD detector to determine the yields. H2, O2, N2, CH4, and CO were analyzed with a Molsieve 5A column (argon carrier, temperature column 100°C), while CO2, C2H4, C2H6, C3H6, and C3H8 were analyzed over a PoraPLOT U column (helium carrier, temperature column 85°C).
The liquid phase was sampled every 30 min and analyzed with a Shimadzu HPLC (Prominence) equipped with an Aminex HPX-87H column and a refractive index detector (RID). The flow rate of the mobile phase (5 mM H2SO4 aqueous solution) was fixed at 0.6 ml/min and the working temperature at 30°C. External calibration curves were used for the quantification. Liquid flow rates were measured according to the variation of the mass in the collecting vessel, assuming density of the liquid product equal to water at room temperature. Gas flow rates were measured with a Brooks gas flow meter; as it was calibrated for nitrogen, the actual value was obtained using conversion factors, knowing the gas composition. The gas flowrate out of the reactor was not constant, coming out in larger bursts, occasionally causing scatter in the micro-GC results and gas mass flow measurement. This problem occurs when gas bubbles form upstream of the backpressure regulator in case of high conversion and high gas-yields. These errors were mitigated by averaging over 30 min periods.
The liquid product yields (C1 based) were calculated according to Equation (1), where YC1, i is the C1 based yield of any carbon containing component i, Fi is the molar flow rate (mol min−1) of component i, NCi is the carbon number of component i. The gas product yields (C1 based) were measured according to Equation (2) as the difference between butanol conversion and liquid yields, based on the fact that the carbon mass balance closed within 10%. The gas selectivity was defined as the ratio between the gas product yield and the butanol conversion. Finally, the selectivity toward a specific component i was defined as the ratio between its molar flow and moles of reacted butanol (Equation 4).
Please note that the maximum selectivity depends on the component, i.e., 4 for C1 compounds, 2 for C2 compounds, 1 for C4 compounds as well as C3 compounds, considering that a butanol molecule can deliver max 1 C3 molecule, and 12 for H2.
XRF, Physisorption and Chemisorption
Table 1 shows the results of N2 physisorption of both fresh and spent 250–420 μm catalyst. Surface area calculated with the seven-point BET method and the NLDFT model are in agreement, showing a small increase during the experiment. The total pore volume is also increased during the experiment and the pore size distribution shifts to slightly larger pores. The catalyst is mesoporous. Surprisingly, application in APR causes a minor but significant increase in surface area and pore volume, despite the fact that ZrO2 is reported to be stable under APR conditions (Elliott et al., 1993). Also, no significant changes in the morphology of ZrO2 were observed during preliminary aging experiments in pure water at 200°C for 12 h (not shown).
Table 2 shows the Rh loading in the fresh and spent catalysts, as measured by XRF. The Rh leaching increased significantly with increasing dimension of the catalyst particles. The reduction in Rh content is caused by leaching. Please note that any apparent decrease of Rh loading due to catalyst mass being increased by deposit formation is in all cases much smaller than the observed effects, therefore leaching is the main cause of the decreasing Rh loading.
Table 3 shows the Rh dispersions and metal surface areas in fresh and spent catalysts as measured by H2 chemisorption as well as the average Rh crystallite size calculated from both Rh surface area and STEM images, shown in Figure 2. The fresh catalysts show very good dispersion and small Rh crystallite size. During the reaction, the available Rh surface area is reduced by factor of typically 2, without any significant effect of the zirconia particle. The increase in metal particle sizes according chemisorption are in reasonable agreement with the observed particle sizes observed with STEM, presented below.
Table 3. H2 chemisorption and particle size determined with STEM; the typical experimental relative error in dispersion and surface area is 10%.
STEM Analysis of Spent Catalyst
Figure 2 shows dark field STEM images of the spent 40–60 and 250–420 μm catalysts. The Rh crystallites in the fresh catalysts were oxidized, making it impossible to see them with STEM without a reductive pretreatment. Rh in the spent catalysts is found in mostly round and sometimes also oval particles. The particle size distribution is quite narrow, as seen in Figure 3 for the 40–60 μm catalyst. Rh particles were considerably harder to locate in the larger catalyst particles, leading to a somewhat limited data set (30 particles as compared to 173 on the 40–60 μm catalyst). However, the distribution in 250–420 μm catalyst (Supplementary Figure A-1) reveals no significant difference with the distribution shown in Figure 3.
Figure 4 shows the C, H, and O elemental analysis results of spent catalysts; the amount of deposits was calculated by subtracting the amount of C, H, and O detected on fresh catalyst from the amounts detected on spent catalyst. The average rate of deposition (in mg h−1) is calculated by dividing the amount of deposits by the TOS in order to take into account small differences in the TOS, as also presented in Figure 4. Formation of deposits clearly increases with increasing particle size, particularly the amounts of carbon and oxygen as well as the total mass of deposit (mg ) and the average rate of deposition.
Figure 4. Elemental composition and amounts of deposits as well as averaged rate of formation of deposits in spent catalysts; contamination on fresh catalysts was subtracted.
Conversion and Products Yields
Figure 5 shows conversion and carbon-based selectivity to gaseous product over time. The conversion over the 40–60 and 60–100 μm catalysts is similar, whereas the conversion over the 250–420 μm catalyst is consistently 30–40% lower. Conversion on all catalysts declined at a similar rate over time. Selectivity to gaseous products remained high on 40–60 and 60–100 μm catalysts, declining slowly over time from ca. 99 to 97% in 150 min. On 250–420 μm catalyst, gas selectivity starts similarly high, declining severely over time. The conversion on bare support was too low to enable reliable quantification (Supplementary Figure A-3A).
Figure 6 shows selectivities to butyraldehyde and butyric acid, the main products in liquid phase, over time. In addition, traces of acetic acid and ethanol were detected. On bare support (Supplementary Figure A-3B), only butyraldehyde is observed in very low amounts, decreasing with TOS. Clearly, butyraldehyde selectivity is much higher compared to the butyric acid selectivity. The liquid product selectivities over 40–60 and 60–100 μm catalysts are very similar, increasing slowly with time on stream. In contrast, the 250–420 μm catalyst produces much more butyraldehyde and butyric acid, both increasing significantly with time on stream. The molar ratio of aldehyde to acid also increases over time, particularly on the 250–420 μm catalyst, from an initial value 4 to about 14 by the end of the experiment.
Figure 6. Selectivity of products in liquid phase; (A) Butyraldehyde and (B) Butyric acid. Note the difference scale between (A) and (B).
Figures 7, 8 show the selectivities to products in the gas phase, i.e., H2, CO2, and CH4 (Figure 7) and C2−3 hydrocarbons (Figure 8). The selectivity to CO was always very low, below 0.05 mol mol−1 (Supplementary Figure A-2). On the bare support formation of H2, CO2, and C3 was detected at such low concentration that quantification was not possible. H2 selectivities (Figure 7A) were similar on all catalysts, showing no effect of the catalyst particle size. CO2 selectivities (Figure 7B) decreased with decreasing size in the order 250–420 > 60–100 > 40–60. The CO2 selectivity decreased over time on the 250–420 μm catalyst, whereas such a trend is less clear on the smaller catalysts. Formation of methane (Figure 7C) increased with decreasing particle size, opposite to CO2. In fact, the 250–420 μm catalyst produced hardly any methane. Also selectivity to C2 compounds (Figures 8A,B) increased with decreasing catalyst particle size, similar to methane. Furthermore, ethane (Figure 8B) was the main C2 product over the 40–60 and 60–100 μm catalysts, while the 250–420 μm catalyst produced more ethylene than ethane (Figure 8A). Propane (Figure 8D) was the main C3 product on all catalysts, despite the relatively large scatter in the data. Furthermore, large catalyst particles (250–420 μm) produce less C3 compounds than the 40–60 and 60–100 μm catalysts, showing similar C3 selectivities.
Mass Transfer Criteria
The significance of external mass transfer was evaluated using the Mears' criterion (Equation 5; Mears, 1971):
The detailed calculation is presented in the Supplementary Material. The highest value obtained for the most active and largest catalyst particles was 2.3*10−3, leading to the conclusion that external mass transfer limitations can be disregarded.
The significance of internal mass transfer was estimated using the Weisz-Prater criterion (Equation 6; Weisz and Prater, 1954). Assuming first order reaction, the criterion becomes:
In which r'A is the observed reaction rate (kmolbutanol s−1), R is the catalyst particle radius (m), ρc is the density of the catalyst (kg m−3), CAs is the concentration of butanol on the external catalyst surface (kmol m−3) and De is the effective diffusion coefficient (m2 s−1), calculated according to Equation (7):
in which DAB is the binary diffusion coefficient of butanol in water (m2 s−1), ϕp is the catalyst pellet porosity and τ is the tortuosity factor. Catalyst density ρc and porosity ϕp were estimated using the bulk density of non-porous ZrO2 and the measured pore volume of the catalyst. As the external mass transfer was considered not limiting, the concentration at the external catalyst surface was assumed to be equal to the bulk concentration, CAs ≈ CAb. Since the exact value of the tortuosity factor was unknown, the value of CWP was calculated as a function of catalyst particle size with values of τ ranging from 1 to 9, shown in Figure 9, using the highest observed reaction rate at time zero. As clearly inferred from the graph, internal diffusion is not limiting for the two smaller particle sizes, but may be limiting for the largest particle size, depending on the value of the tortuosity factor. Therefore, the role of internal diffusion required further experimental evaluation, as will be discussed later. In order to discuss the effect of internal transport on the product distribution, it is necessary to discuss first the main reactions contributing to the overall conversion.
Catalyst Performance; Comparison to Literature
APR of alcohols and polyalcohols has been widely investigated, in particular with C-to-O ratio of 1:1 (Shabaker et al., 2003a; Kirilin et al., 2010, 2014). Ethanol, propanol or butanol were typically less studied; the presence of the alkyl groups makes complete reforming less facile and formation of alkanes is reported.
Lobo et al. (2012) studied APR of n-propanol over Pt catalysts (5 wt.% propanol solution, 250°C temperature, 69 bar pressure), reporting reaction rates per gram catalyst about one order of magnitude higher. Unfortunately, the Pt loading is not reported. Ethane and carbon dioxide are the main products, with a ratio close to 1, whereas propanal found in the liquid phase.
Godina et al. studied the reforming of alcohols with 3 carbon atoms using Pt-based catalyst on polymer based spherical activated carbons (Godina et al., 2018). Even in this case, exclusively propanal and propionic acid are reported in the liquid phase, whereas ethane was the most formed alkane. On the other hand, methane was not detected, in contrast to our results in the case of the smaller catalyst particles. The presence of propanal and propionic acid from 1-propanol APR was also confirmed by Wawretz and co-workers (Wawrzetz et al., 2010) with an alumina supported Pt catalyst.
Coronado et al. studied the APR of ethanol and propanol with nickel based catalysts on ceria-zirconia supports (Coronado et al., 2018, and literature cited therein), reporting propanal as main liquid product, together with a small amount of propionic acid. On the other hand, methane was found in the gas phase, which was attributed to methanation of CO and CO2.
Pipitone et al. (2019) studied the aqueous phase reforming of butanol with a Pt-based catalyst in a batch system. The produced gas phase contained hydrogen carbon-dioxide and propane in a 2-1-1 molar ratio, in accordance with the previous works on C3 alcohols (Godina et al., 2018). On the other hand, the nature of the active sites seems to have an effect on the product distribution in the gas phase. Indeed, Roy et al. (2011, 2014) studied APR of n-butanol over alumina and ceria supported Ni catalysts (20 wt.% Ni loading, 5 wt.% butanol solution, 185–215°C and 10–31 bar), resulting in reaction rates an order of magnitude lower. The reported product distribution on Ni is similar to the present work, with the exception of somewhat higher selectivity to C2 hydrocarbons, which the authors attributed to Fischer-Tropsch reactions. The main reaction pathway was proposed to proceed through consecutive dehydrogenation and decarbonylation steps, very similar to the findings of Lobo et al. for APR of propanol with Pt catalysts (Lobo et al., 2012). Tishchenko coupling of propionaldehyde followed by hydrolysis and propionic acid decarboxylation was observed as a side reaction pathway on Pt. Moreover, it was suggested that propane was reformed to hydrogen and carbon monoxide, as likely occurred also in the reaction conditions reported during the present investigation. This result may be explained by the higher activity of nickel toward C-C bond breaking compared to platinum (Sinfelt, 1973).
Reaction Scheme for Rh Catalyst
As no reaction scheme has been reported so far for APR of butanol over Rh catalysts, a proposal for such a scheme is described below (Figure 10), based on the observations described above and inspired by literature on APR of butanol on Ni catalysts (Roy et al., 2011, 2014) and of 1-propanol on Pt catalysts (Lobo et al., 2012). Please note this scheme is required for discussing the effect of mass transfer limitations on the product distribution. The scheme is divided in three stages: (I) initial conversion, (II) initial gas formation, and (III) hydrocarbon reactions and reforming.
Stage I is the initial conversion of butanol to liquid intermediates via reactions 1 and 2. Reaction 1, dehydrogenation of butanol into butyraldehyde, mostly takes place on Rh. Figure 6A shows a significant yield of butyraldehyde over the Rh catalyst, whereas the yield over plain zirconia (Supplementary Figure A-3B) is <0.5 mol% of butanol. This agrees with literature, where dehydrogenation of alcohols is mostly reported on metals (Roy et al., 2011, 2014; Lobo et al., 2012), although it is also reported on zirconia (Sabatier and Mailhe, 1910; Shinohara et al., 1997). Interestingly, the butyraldehyde selectivity (Figure 6A) increases over time and the selectivity to gas products (Figure 5B) declines simultaneously. This is especially clear over the 250–420 μm catalyst, indicating that the reactions forming gasses deactivate more than the dehydrogenation reaction. Furthermore, no C3 oxygenates were observed in the liquid phase. Therefore, reactions involving C-C cleavage in butanol are not significant and reaction via butyraldehyde dominates, in agreement with suggestions in literature for Ni and Pt (Roy et al., 2011, 2014; Lobo et al., 2012). Reaction 2 is the formation of butyric acid, most likely via Tishchenko esterification on Rh, supported by the fact that butyric acid selectivity closely follows the selectivity for butyraldehyde (Figures 6A,B). The intermediate butyl-butyrate was not observed, as would be expected considering that in aqueous medium the equilibrium strongly favors hydrolysis of butyl-butyrate.
Stage II is the formation of COx according reactions 3 and 4, respectively, decarbonylation and decarboxylation. The claim that the C-C bond is cleaved via decarbonylation and/or decarboxylation is supported by the observations that propane and C4 oxygenate products are the dominant products, as well as by the absence of any C3 oxygenates and the virtual absence of C2 oxygenates. Decarbonylation of butyraldehyde (reaction 3) is the main reaction forming gaseous products, as much more butyraldehyde than butyric acid is formed (Figure 6). The reaction mainly takes place on the metal, but can also be observed on the bare support (Supplementary Figure A-3A), in agreement with literature (Bie et al., 2013; Miao et al., 2016). Decarboxylation of butyric acid (reaction 4) requires the presence of an active metal according literature (Bie et al., 2013; Miao et al., 2016).
Stage III consists of the reactions of hydrocarbons (reforming and (de-)hydrogenation) and (R-)WGS (reactions 5–10). The water-gas-shift reaction (WGS, reaction 5) converts virtually all CO to CO2 and H2. The aqueous environment is expected to strongly drive WGS toward the products (Shabaker et al., 2003b), causing extremely low CO yields (Supplementary Figure A-2).
Propane produced in reactions 3 and 4 can be dehydrogenated over Rh (reaction 7) as indicated by the formation of propylene in Figure 8C. The shorter C1−2 hydrocarbons are likely formed via propane hydrogenolysis (reaction 8). Although propane hydrogenolysis has not been reported in aqueous phase, it is known to be readily catalyzed by Rh in this temperature range in gas phase (Bond et al., 1996).
Methanation (reaction 6) has been reported on Rh in APR (Davda et al., 2003); however, methanation is slow at the low temperature of operation (Mutz et al., 2015). We speculate that the dominant pathway to methane is hydrogenolysis of propane based on the observation that the methane selectivity (Figure 7C) is very similar to the C2 selectivity (sum of ethane and ethylene in Figures 8A,B), consistent with formation of methane, ethane and ethylene according reactions 8 and 9. Interestingly, selectivity of both methane and C2 increase when operating with smaller catalyst particles.
Reforming of hydrocarbons (an aggregate reaction 10) probably contributes, possibly via the activation mechanism discussed above, based on three observations. First, on all catalysts the selectivity to CO2 (Figure 7B) is higher than 1 mol/mol butanol, indicating that next to reactions 3 and 4 an additional pathway to CO2 must exist. Second, H2 formation via exclusively the sequence dehydrogenation-decarbonylation-WGS (reactions 1, 3, 5) would result in selectivity to hydrogen of 2 mol per mole butanol converted. Figure 7A shows that the H2 selectivity varies between 3 and 6, providing clear evidence that an additional pathway contributes significantly. Third, the combined selectivity of C2−3 products (Figure 8) is <1 mol per mol butanol converted, clearly indicating that C2−3 products are being consumed in a consecutive reaction, most likely reforming. In addition to experimental observations, Rh is more active for steam reforming at low temperature compared to e.g., Pt or Ni (Kolb et al., 2004; Halabi et al., 2010), supporting the suggestion that reforming of alkanes contributes at the temperatures used in the present work. Further research would be needed to confirm this hypothesis.
Effect of Mass Transfer
Internal mass transfer can affect the local concentration of reactants as well as intermediate reaction products, possibly affecting the consecutive reactions and therefore the final product distribution. The data shows several clear trends as discussed below.
Butanol conversion is limited by diffusion in large catalyst particles. Figure 5A shows that conversion of butanol over the 250–420 μm catalyst is lower than over the smaller catalyst particles, whereas the conversion is similar over both of the smaller particle sizes. This indicates that diffusion limits the reaction on large catalyst particles, in agreement with the estimated values of Weisz-Prater criterion in Figure 9. Possibly, subtle differences in the metal dispersion, metal loading and metal surface area (Table 3) cause the relative high activity of the 60–100 μm fraction. On the other hand, difference in performance between the 40–60 μm fraction and the 250–420 μm fraction are clearly due to mass transfer effects.
Butyraldehyde coupling is enhanced on large particles as the selectivity to butyric-acid is higher on the larger particles (Figure 6B). The Tishchenko coupling reaction (reaction 2) is 2nd order (Anderson and Peters, 1960) and thus rates are strongly influenced by the butyraldehyde concentration. Slow diffusion increases the local concentration of butyraldehyde in the center of the catalyst particles, thus resulting in high butyric-acid yield.
Conversion of hydrocarbons is enhanced on large catalyst particles, including dehydrogenation (reaction 7), hydrogenolysis (reaction 8), and reforming (reactions 6, 10). This is based on three observations. First, the selectivity to methane (Figure 7C) and C2−3 hydrocarbons (Figure 8) is the lowest on large catalyst particles, suggesting more reforming, attributed to sluggish diffusion of dissolved C1, 2, 3 hydrocarbons. Second, the ratio of olefins to alkanes is at the same time higher on the large catalyst particles, indicative for more dehydrogenation. Third, CO2 selectivity (Figure 7B) is the highest over the large catalyst particles. Furthermore, the CO2 selectivity is well above 1 mol per mol butanol converted, especially on large catalyst particles (Figure 7B), also indicating reforming activity as discussed above. The products of reactions 3 and 4 (i.e., CO, CO2, and propane) are all gasses and diffusion out of the catalyst particles is controlled by diffusion of these molecules dissolved in water. However, the H2 selectivity decreases with catalyst particle size (Figure 7A) despite reforming activity of large particles, indicating that longer diffusion length increases consecutive reactions consuming hydrogen more than consecutive reactions producing hydrogen.
This argument is analogous to the mechanism proposed by Neira D'Angelo et al. (2013), accounting for decreasing H2 yield in case of mass transfer limitation because of the lower selectivity, while the conversion is not strongly affected. The difference in the operation of the reactors should be noted. In our case, we do not co-feed any gas, and formation of gas phase can occur only via nucleation of oversaturated solutions. In the case of Neira D'Angelo et al. (2013) as well as most APR studies in continuous operation, inert gas is added to the reactant stream and the reactor operated in trickle-phase mode, resulting in stripping of gaseous components in the reactor. This will enhance internal diffusion by decreasing the concentration at the external surface of the catalyst particles, i.e., in the bulk of the liquid. Therefore, our experiments are more sensitive for effects of internal diffusion of intermediate products.
The particle size affects not only the liquid and gaseous products distribution, but also strongly influences the stability of the catalyst. As reported in Figure 5B, the gas selectivity declined slowly for the 40–60 and 60–100 μm catalysts, while decreasing significantly for the 250–420 μm catalyst. This can be understood considering the high butyraldehyde concentration present inside the large catalyst particles. The aldehyde is a precursor of deposits via aldol-condensation reactions enhanced by the acid-base properties of the zirconia support, causing deactivation of the catalyst (Takanabe et al., 2006; Koichumanova et al., 2018). This hypothesis is also supported by the fact that larger particles form more deposits (Figure 4) and the high oxygen and hydrogen content of the deposits agrees with formation via butyraldehyde condensation, similar to results with the same catalyst system in gas phase reforming (Harju et al., 2016). Furthermore, it is also clear that larger support particles suffer more from Rh leaching than small particles (Table 2) which is in line with the fact that the formation of butyric acid is enhanced, resulting in more acidic conditions inside the larger catalyst particles, enhancing leaching.
Finally, it should be noted that description of internal mass transfer according Thiele modulus and Weisz-Prater criterion can be criticized for the case of APR. These models assume molecular diffusion of dissolved species in pores filled with water. However, APR produces molecules (H2, CO, CO2, CH4, and other light hydrocarbons) that are forming a new phase during the reaction, i.e., gas phase. The critical question is whether bubbles form inside or outside the catalyst particles and at this time this question remains unanswered. Theoretical work of Datsevich (2003a; 2003b; 2004; 2005, Oehmichen et al., 2010) has shown that formation of gas bubbles can influence transfer inside particles dramatically, whereas we have recently demonstrated in a microfluidic device, mimicking pores in a catalyst support, that both retardation as well as enhancement of transport on the pore is possible (Espinosa et al., 2018). Very recently, a simulation using APR of glycerol as model reaction showed the effect of bubbles on the kinetics and transport phenomena in a 2D system (Ripken et al., 2019). This clearly needs further research to decide if such effects experimentally contribute in practical 3D porous catalysts.
The influence of internal mass transfer in Rh/ZrO2 catalysts for APR of 1-butanol was studied by varying the support particle size. Larger support particles cause a minor but significant decrease in activity, which can be attributed to internal mass transfer limitations. The effect on selectivity and stability is much stronger though.
A reaction scheme is proposed in order to discuss the effects of internal mass transfer on product distribution and stability. The reaction scheme includes three main stages i.e., (1) formation of liquid intermediates via dehydrogenation, (2) formation of gas via decarbonylation/decarboxylation reactions and (3) hydrocarbon hydrogenolysis/reforming/dehydrogenation. The main pathway to hydrogen involves stages 1 and 2, while methane, ethane, and ethylene are formed via hydrogenolysis of propane.
Internal mass transport limitation of butyraldehyde enhances formation of butyric acid. Slow diffusion of methane, ethane, and propane dissolved in water promotes consecutive reactions, i.e., dehydrogenation, hydrogenolysis, and reforming. Furthermore, large support particles deactivate faster, attributed to high concentrations of butyraldehyde inside the catalyst particles, enhancing deposit formation via aldol condensation reactions, as well as to local high acidity caused by butyric acid.
Data Availability Statement
The datasets generated for this study are available on request to the corresponding author.
HH and LL: conceptualization, methodology, validation, formal analysis, investigation, project administ, and ration. LL: resources, supervision, and funding acquisition. HH, GP, and LL: data curation and writing—review and editing. HH: writing—original draft preparation.
The authors declare that this study received funding from TEKES. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors thank Ing. B. Geerdink and K. Altena-Schildkamp for technical assistance, Dr. Kaisa Vikla for performing part of the experiments and Ellen de Wit for the preliminary aging experiments. The XRF analyses were made with the assistance of Dr. Esko Ahvenniemi. This work made use of the Aalto University Nanomicroscopy Center (Aalto-NMC) premises. Dr. Hua Jiang is acknowledged for performing the microscopy.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2020.00017/full#supplementary-material
Bie, Y., Gutierrez, A., Viljava, T. R., Kanervo, J. M., and Lehtonen, J. (2013). Hydrodeoxygenation of methyl heptanoate over noble metal catalysts: catalyst screening and reaction network. Ind. Eng. Chem. Res. 52, 11544–11551. doi: 10.1021/ie4012485
Coronado, I., Pitínová, M., Karinen, R., Reinikainen, M., Puurunen, R. L., and Lehtonen, J. (2018). Aqueous-phase reforming of Fischer-Tropsch alcohols over nickel-based catalysts to produce hydrogen: Product distribution and reaction pathways. Appl. Catal. A Gen. 567, 112–121. doi: 10.1016/j.apcata.2018.09.013
Datsevich, L. B. (2003b). Some theoretical aspects of catalyst behaviour in a catalyst particle at liquid (liquid–gas) reactions with gas production: oscillation motion in the catalyst pores. Appl. Catal. A Gen. 247, 101–111. doi: 10.1016/S0926-860X(03)00091-7
Datsevich, L. B. (2004). Oscillation theory: Part 1. Temperature difference between the center of a catalyst particle and its surface: contradiction to Thiele/Zeldovich model. Appl. Catal. A Gen. 262, 149–153. doi: 10.1016/j.apcata.2003.11.023
Davda, R. R., Shabaker, J. W., Huber, G. W., Cortright, R. D., and Dumesic, J. A. (2003). Aqueous-phase reforming of ethylene glycol on silica-supported metal catalysts. Appl. Catal. B Environ. 43, 13–26. doi: 10.1016/S0926-3373(02)00277-1
Davda, R. R., Shabaker, J. W., Huber, G. W., Cortright, R. D., and Dumesic, J. A. (2005). A review of catalytic issues and process conditions for renewable hydrogen and alkanes by aqueous-phase reforming of oxygenated hydrocarbons over supported metal catalysts. Appl. Catal. B Environ. 56, 171–186. doi: 10.1016/j.apcatb.2004.04.027
De Vlieger, D. J. M., Mojet, B. L., Lefferts, L., and Seshan, K. (2012). Aqueous phase reforming of ethylene glycol - role of intermediates in catalyst performance. J. Catal. 292, 239–245. doi: 10.1016/j.jcat.2012.05.019
Elliott, D. C., Sealock, L. J., and Baker, E. G. (1993). Chemical processing in high-pressure aqueous environments. 2. Development of catalysts for gasification. Ind. Eng. Chem. Res. 32, 1542–1548. doi: 10.1021/ie00020a002
Godina, L. I., Tokarev, A. V., Simakova, I. L., Mäki-Arvela, P., Kortesmäki, E., Gläsel, J., et al. (2018). Aqueous-phase reforming of alcohols with three carbon atoms on carbon-supported Pt. Catal. Today 301, 78–89. doi: 10.1016/j.cattod.2017.03.042
Gutiérrez Ortiz, F. J., and Campanario, F. J. (2018). Hydrogen production from supercritical water reforming of acetic acid, acetol, 1-butanol and glucose over Ni-based catalyst. J. Supercrit. Fluids 138, 259–270. doi: 10.1016/j.supflu.2018.04.023
Gutiérrez Ortiz, F. J., Campanario, F. J., and Ollero, P. (2016). Supercritical water reforming of model compounds of bio-oil aqueous phase: acetic acid, acetol, butanol and glucose. Chem. Eng. J. 298, 243–258. doi: 10.1016/j.cej.2016.04.002
Halabi, M. H., De Croon, M. H. J. M., Van Der Schaaf, J., Cobden, P. D., and Schouten, J. C. (2010). Low temperature catalytic methane steam reforming over ceria-zirconia supported rhodium. Appl. Catal. A Gen. 389, 68–79. doi: 10.1016/j.apcata.2010.09.004
Harju, H., Lehtonen, J., and Lefferts, L. (2016). Steam reforming of n-butanol over Rh/ZrO2 catalyst: role of 1-butene and butyraldehyde. Appl. Catal. B Environ. 182, 33–46. doi: 10.1016/j.apcatb.2015.09.009
Kaila, R. K., Gutiérrez, A., Korhonen, S. T., and Krause, A. O. I. (2007). Autothermal reforming of n-dodecane, toluene, and their mixture on mono- and bimetallic noble metal zirconia catalysts. Catal. Lett. 115, 70–78. doi: 10.1007/s10562-007-9075-z
Kaila, R. K., Gutiérrez, A., Slioor, R., Kemell, M., Leskelä, M., and Krause, A. O. I. (2008). Zirconia-supported bimetallic RhPt catalysts: characterization and testing in autothermal reforming of simulated gasoline. Appl. Catal. B Environ. 84, 223–232. doi: 10.1016/j.apcatb.2008.04.001
Kirilin, A. V., Tokarev, A. V., Manyar, H., Hardacre, C., Salmi, T., Mikkola, J. P., et al. (2014). Aqueous phase reforming of xylitol over Pt-Re bimetallic catalyst: effect of the re addition. Catal. Today 223, 97–107. doi: 10.1016/j.cattod.2013.09.020
Kirilin, A. V., Tokarev, A. V., Murzina, E. V., Kustov, L. M., Mikkola, J. P., and Murzin, D. Y. (2010). Reaction products and transformations of intermediates in the aqueous-phase reforming of sorbitol. ChemSusChem 3, 708–718. doi: 10.1002/cssc.200900254
Koichumanova, K., Vikla, A. K. K., Cortese, R., Ferrante, F., Seshan, K., Duca, D., et al. (2018). In situ ATR-IR studies in aqueous phase reforming of hydroxyacetone on Pt/ZrO2 and Pt/AlO(OH) catalysts: The role of aldol condensation. Appl. Catal. B Environ. 232, 454–463. doi: 10.1016/j.apcatb.2018.03.090
Kolb, G., Zapf, R., Hessel, V., and Löwe, H. (2004). Propane steam reforming in micro-channels-results from catalyst screening and optimisation. Appl. Catal. A Gen. 277, 155–166. doi: 10.1016/j.apcata.2004.09.007
Larimi, A., and Khorasheh, F. (2019). Renewable hydrogen production over Pt/Al2O3 nano-catalysts: effect of M-promoting (M=Pd, Rh, Re, Ru, Ir, Cr). Int. J. Hydrogen Energy 44, 8243–8251. doi: 10.1016/j.ijhydene.2019.01.251
Lobo, R., Marshall, C. L., Dietrich, P. J., Ribeiro, F. H., Akatay, C., Stach, E. A., et al. (2012). Understanding the chemistry of H2 production for 1-propanol reforming: pathway and support modification effects. ACS Catal. 2, 2316–2326. doi: 10.1021/cs300405s
Miao, C., Marin-Flores, O., Davidson, S. D., Li, T., Dong, T., Gao, D., et al. (2016). Hydrothermal catalytic deoxygenation of palmitic acid over nickel catalyst. Fuel 166, 302–308. doi: 10.1016/j.fuel.2015.10.120
Mutz, B., Carvalho, H. W. P., Mangold, S., Kleist, W., and Grunwaldt, J. D. (2015). Methanation of CO2: structural response of a Ni-based catalyst under fluctuating reaction conditions unraveled by operando spectroscopy. J. Catal. 327, 48–53. doi: 10.1016/j.jcat.2015.04.006
Neira D'Angelo, M. F., Ordomsky, V., Van Der Schaaf, J., Schouten, J. C., and Nijhuis, T. A. (2013). Aqueous phase reforming in a microchannel reactor: the effect of mass transfer on hydrogen selectivity. Catal. Sci. Technol. 3, 2834–2842. doi: 10.1039/c3cy00577a
Neira D'Angelo, M. F., Ordomsky, V., Van Der Schaaf, J., Schouten, J. C., and Nijhuis, T. A. (2014a). Continuous hydrogen stripping during aqueous phase reforming of sorbitol in a washcoated microchannel reactor with a Pt-Ru bimetallic catalyst. Int. J. Hydrogen Energy 39, 18069–18076. doi: 10.1016/j.ijhydene.2014.02.167
Neira D'Angelo, M. F., Schouten, J. C., Van Der Schaaf, J., and Nijhuis, T. A. (2014b). Three-phase reactor model for the aqueous phase reforming of ethylene glycol. Ind. Eng. Chem. Res. 53, 13892–13902. doi: 10.1021/ie5007382
Oehmichen, T., Datsevich, L. B., and Jess, A. (2010). Influence of bubble evolution on the effective kinetics of heterogeneously catalyzed gas/liquid reactions. Part I: reactions with gaseous products. Chem. Eng. Technol. 33, 911–920. doi: 10.1002/ceat.200900624
Pipitone, G., Zoppi, G., Ansaloni, S., Bocchini, S., Deorsola, F. A., Pirone, R., et al. (2019). Towards the sustainable hydrogen production by catalytic conversion of C-laden biorefinery aqueous streams. Chem. Eng. J. 377:120677. doi: 10.1016/j.cej.2018.12.137
Ripken, R. M., Wood, J. A., Gardeniers, J. G. E., and Le Gac, S. (2019). Aqueous-phase reforming in a microreactor: the role of surface bubbles. Chem. Eng. Technol. 42, 1–9. doi: 10.1002/ceat.201900142
Roy, B., Sullivan, H., and Leclerc, C. A. (2014). Effect of variable conditions on steam reforming and aqueous phase reforming of n-butanol over Ni/CeO2 and Ni/Al2O3 catalysts. J. Power Sources 267, 280–287. doi: 10.1016/j.jpowsour.2014.05.090
Shabaker, J. W., Davda, R. R., Huber, G. W., Cortright, R. D., and Dumesic, J. A. (2003a). Aqueous-phase reforming of methanol and ethylene glycol over alumina-supported platinum catalysts. J. Catal. 215, 344–352. doi: 10.1016/S0021-9517(03)00032-0
Shabaker, J. W., Huber, G. W., Davda, R. R., Cortright, R. D., and Dumesic, J. A. (2003b). Aqueous-phase reforming of ethylene glycol over supported platinum catalysts. Catal. Letters. 88, 1–8. doi: 10.1023/A:1023538917186
Shinohara, Y., Nakajima, T., Suzuki, S., Mishima, S., and Ishikawa, H. (1997). A computational chemical investigation of the dehydration and dehydrogenation of ethanol on oxide catalysts. J. Chem. Softw. 4:89. doi: 10.2477/jchemsoft.4.89
Takanabe, K., Aika, K., Ken-ichi, A., Seshan, K., and Lefferts, L. (2006). Catalyst deactivation during steam reforming of acetic acid over Pt/ZrO2. Chem. Eng. J. 120, 133–137. doi: 10.1016/j.cej.2006.04.001
Wawrzetz, A., Peng, B., Hrabar, A., Jentys, A., Lemonidou, A. A., and Lercher, J. A. (2010). Towards understanding the bifunctional hydrodeoxygenation and aqueous phase reforming of glycerol. J. Catal. 269, 411–420. doi: 10.1016/j.jcat.2009.11.027
Keywords: aqueous phase reforming, hydrogen, mass transfer, reaction pathway, rhodium
Citation: Harju H, Pipitone G and Lefferts L (2020) Influence of the Catalyst Particle Size on the Aqueous Phase Reforming of n-Butanol Over Rh/ZrO2. Front. Chem. 8:17. doi: 10.3389/fchem.2020.00017
Received: 29 October 2019; Accepted: 08 January 2020;
Published: 28 January 2020.
Edited by:Dmitry Yu. Murzin, Åbo Akademi University, Finland
Reviewed by:Benjaram M. Reddy, Indian Institute of Chemical Technology (CSIR), India
Atte Aho, Åbo Akademi University, Finland
Copyright © 2020 Harju, Pipitone and Lefferts. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Leon Lefferts, email@example.com