1. Introduction
Hydrogen is mostly produced today through steam reforming of natural gas [1,2]. The main reaction is assumed to be represented by steam methane reforming (SMR), an endothermic equilibrium reaction:
CH4 + H2O ⇆ CO + 3 H2 ΔH298 = +206.63 kJ/mol
Establishing together with the water gas shift (WGS) equilibrium:
CO + H2O ⇆ CO2 + H2 ΔH298 = −41.16 kJ/mol
Resulting, at least formally, in a formal “global reforming reaction” (GRR):
CH4 + 2 H2O ⇆ CO2 + 4 H2 ΔH298 = +165.47 kJ/mol
Thus producing a “syngas” containing H2, CO, CO2 and unreacted CH4.
Large size SMR reactors for hydrogen production from methane are well known and widely applied on an industrial scale. On the other hand, hydrogen is now also at the cutting edge of research in the so-called Power-to-X (P2X) framework. P2X studies, where X represents an unknown, are in progress in order to identify the best options to store or use electrical power, often produced by renewable sources, exceeding the electrical grid demand. One option currently under investigation is to generate hydrogen through water electrolysis. Subsequently, hydrogen can be employed in a methanation reaction, to produce methane to be fed into the natural gas distribution pipeline. Methanation reaction consists of the synthesis of methane from hydrogenation of COx. Formally, methanation of carbon monoxide (MCO) is the reverse of SMR (MCO = revSMR):
CO + 3 H2 ⇆ CH4 + H2O ΔH298 = −206.63 kJ/mol
While methanation of CO2 (MCO2) is formally the reverse of GRR (MCO2 = revGRR):
CO2 + 4 H2 ⇆ CH4 + 2 H2O ΔH298 = −165.47 kJ/mol
The MCO2 reaction could result from the previous conversion of CO2 into CO with a reverse water gas shift (revWGS):
CO2 + H2 ⇆ CO + H2O ΔH298 = +41.16 kJ/mol
Followed by MCO.
P2X studies, and, more in general, all the studies about new energy scenarios where the use of hydrogen is introduced as a medium for decarbonization, are drawing increasing attention towards SMR and MCO2 chemical reactors. For example, in a recent work [3], a scheme was proposed, where a reversible solid oxide fuel cell (SOFC) is coupled to SMR and MCO2 reactors, aiming at electrical grid stabilization. It must be pointed out that many of these advanced energy scenario studies address distributed generation, where the size of the SMR and MCO2 reactors is much smaller than the typical size of industrial SMR reactors. Just to give an example of the size of typical industrial SMR reactors applied in refineries, for 1 Mt/y of crude refining capacity, about 1 t/h of hydrogen is necessary, which requires about 25 tubes (typical dimensions 13 m length and 0.1 m diameter) to be produced through traditional multi-tubular SMR reactors. On the other hand, an example of a small size SMR tubular reactor to be embedded into a solid oxide fuel cell (SOFC) system is reported in [4], with the catalytic section featuring 0.4 m in length and 0.125 m in diameter.
In this context, as a first step for the optimized design of small-scale reactors for advanced energy scenario applications, two laboratory reactors are set up, one for SMR and the other for MCO2. Three different catalysts are employed and tested, two commercial and one homemade. The materials used in the laboratory plant are chosen as similar as possible to the materials to be used in small-scale reactors, and in particular typical stainless-steel piping is employed.
In parallel, a simulation model is developed, to be used as the basis to guide the future scale-up steps. The model is 1D and pseudo-homogeneous and includes mass and energy local balances. For local kinetics, the model proposed by Xu and Froment [5] is implemented, based on the assumption that the SMR, WGS and GRR reactions are reversible and can proceed simultaneously. The model is used to analyze the behavior of the laboratory reactors in detail. Indeed, from a review of the models developed in the literature for methanation reactors [6] and for SMR reactors [7,8,9], it emerges that pseudo-homogeneous models with isothermal temperature profiles are often employed for the simulation of laboratory reactors, with the purpose of identifying the correct kinetic model and/or the kinetic parameters of the selected kinetic model. On the other hand, more detailed heterogeneous models with thorough evaluation of 1D or 2D concentration and temperature profiles inside the individual catalyst particles, and in the solid and gas phase along the full-size reactor, are employed for a large scale industrial reactor design and deployment [6,10,11,12,13,14,15,16]. In this work, a model is proposed for a laboratory scale reactor, where the hypothesis of pseudo-homogeneous behavior is retained (due to the small size of the catalyst particles), but thermal effects along the reactor are analyzed in more detail by including the local energy balance. In this way, 1D temperature profiles along the laboratory reactor are evaluated, and it is found that, under relevant experimental conditions, temperature profiles can easily deviate from uniformity, due to the high enthalpy change of the reaction. Furthermore, possible development of the WGS reaction in the discharge piping is investigated, which often takes place [17,18,19]. In this contribution, instead of attempting to minimize all the temperature-related effects described above at the experimental level, it is proposed to include them in the simulation model.
2. Experimental 2.1. Materials
In this paper, an investigation of two commercial and one homemade catalyst is reported. Commercial Johnson Matthey (JM) Katalco 57-4Q and 25-4Q quadrilobe SMR catalysts were used, after gentle grinding. According to the literature [20,21,22], the 57-4Q catalyst contains 13 wt % Ni on calcium aluminate cement, probably with a Ca:Al ratio of about 1:5 [23]. The surface area was 29 m2 g−1. The 25-4Q catalyst was a slightly alkalized version, doped with small amounts (1.8 wt %) of K2O [20]. The preparation of the homemade catalysts with a 26 wt % nickel loading is reported in Section 2.2.—Homemade Catalyst Preparation.
2.2. Homemade Catalyst Preparation
A “homemade methanation catalyst” (HMMC) was prepared and used as a reference, where Siralox 5/170 (5 wt % SiO2 and 95 wt % Al2O3) was used as the support after calcination at 1073 K for 5 h. Nickel was then deposited through wet impregnation of an aqueous solution of Ni(NO3)2∙6H2O by achieving the desired Ni loading (26 wt % as wtNi × 100/wtcatalyst). The resulting powder was dried at 363 K for 5 h and calcined at 1073 K for 5 h. The composition of this catalyst was comparable both with those of pre-reforming catalysts [20,24] and of commercial COx methanation catalyst. The resulting catalyst had a surface area of 80 m2/g, and the XRD analysis shows the feature of a nearly cubic spinel of NiO. Extensive characterization data of this catalytic material were reported previously [25].
2.3. Laboratory Reactors
The MCO2 and the SMR experiments were performed in two different laboratory scale fixed-bed catalytic tubular reactors [26], respectively. The geometry was the same in both cases (Figure 1), whereas the dimensions of the two reactors were different (main data in Table 1).
In the MCO2 experimental reactor, the reactor tube was made of silica glass. The catalytic packed bed (indicated as PBR, packed bed reactor, in Figure 1) was 26 mm long, and was formed by 88.2 mg of Ni-based catalyst, crushed into particles of about 0.2 mm size and dispersed with 700 mg of silica glass particles of 0.21–0.25 mm size (corresponding to 60–70 mesh sieved). The catalytic bed was held in position by two holders made of quartz wool (H), about 5 mm long. Downstream, the gaseous reactants flow through the empty silica glass tube. This section was denominated post-reactor 1 (P-R1), and it was about 150 mm long. The silica glass tube, containing the PBR, was placed inside a tubular electric furnace with controlled temperature. At the exit of P-R1, the silica glass tube was connected to an empty stainless-steel pipe (AISI 316L, typical composition code 18/8/3, indicating 18 wt % Cr, 8 wt % Ni and 3 wt % Mo), denominated post-reactor 2 (P-R2). The P-R2, 300 mm long, was directly exposed to the ambient temperature. At the exit of P-R2 the gaseous mixture was sampled for analysis.
The SMR experimental reactor had a similar configuration, and was based on a silica glass tube containing a catalytic bed (PBR) consisting of 88.2 mg of Ni-based catalyst, crushed into particles of about 0.2 mm size, and dispersed with 440 mg of silica glass particles of 0.21–0.25 mm size. Here, the PBR was 12 mm long and, again, was held in position by two holders made of quartz wool (H), about 5 mm long. Downstream, the empty silica glass tube (P-R1) was 50 mm long. The silica glass tube, containing the PBR, was placed inside a tubular furnace with controlled temperature. P-R2 was present again (AISI 316L), 40 mm long, where the temperature was maintained at approximately 523 K using a heating ribbon. At the exit of P-R2, the gaseous mixture was sampled for analysis. For all the tubes employed and described above, the internal diameter d was 6 mm and the tube wall thickness s was 1 mm. 2.4. Catalytic Experiments All catalytic experiments were performed under steady state conditions. The MCO2 experiments were performed with feed gas composition: 6.1% CO2, 29.8% H2, N2 balance, with 81.6 mL min−1 total flow rate (at RTP, reference temperature and pressure of 273.15 K and 101.325 kPa respectively), corresponding to 6.7 × 103 h−1 gas hourly space velocity (GHSV) calculated (at RTP) on the basis of the volume of the catalytic section (PBR) of the reactor. Online-product analysis was performed through a Nicolet 6700 FT-IR spectrometer. A water condensation step was set before the IR cell. Quantitative results were obtained after a careful calibration, by looking at the frequencies 2293 cm−1 for CO2, 2170 cm−1 for CO and 1333 cm−1 for CH4.
For the SMR experiments, the feed gas composition was 5% CH4, 20% H2O and He balance. Two total flow rates were experimented: 120 mL/min (at RTP), corresponding to GHSV = 2.1 × 104 h−1 (at RTP) calculated on the basis of the volume of the catalytic section (PBR) of the reactor, and 80 mL/min (at RTP) corresponding to GHSV = 1.4 × 104 h−1 (at RTP). An Agilent 4890 gas-chromatograph, equipped with a “Molsieve 5A/Porabond Q Tandem” column and TCD and FID detectors in series, was used for gas analysis. Moreover, GC/MS (FOCUS and ISQ from Thermo-Fisher) analysis was carried out in order to have a precise identification of the compounds. The catalytic line used in the experiments is schematized in Figure 2.
3. Modeling Each laboratory reactor was simulated as a sequence of three tubular reactors: the PBR, P-R1 and P-R2. The holders (H) were not simulated. The models developed for the three tubular sections are reported and discussed below. 3.1. PBR Model
A PBR is an extremely complex random geometrical arrangement. PBR models have been developed with different degrees of detail, depending on the desired accuracy, the available information on the packed bed and the required computational efforts. Systematic reviews are given by Lemcoff et al. [27], and by Koning [28]. Among the models, the pseudo-homogeneous ones have obtained extensive application, due to their convenient mathematical formulation. In particular, when dealing with the simulation of steady-state reactors at conditions not too close to run-away, the one-dimensional pseudo-homogeneous plug-flow model is generally considered adequate [29], so it is adopted in this work. The fixed-bed catalytic reactor is simulated as one-dimensional along the reactor axis, and as a pseudo-continuum, where the solid and fluid phases are considered as interpenetrating continua, i.e., as if they coexist all over the tube. The essential feature differentiating pseudo-homogeneous from heterogeneous models [30] is that the former models do not include a detailed simulation of mass and energy transport inside the catalytic particle (internal transport limitations) or in the boundary layer around the catalytic particles (external transport limitations). In other words, temperature and concentration profiles in the catalytic particles are assumed to be uniform, with values identical to those in the bulk of the reacting fluid, which is consistent with the assumption of an effectiveness factor η equal to 1. Conversely, internal and/or external transport limitations are accounted for in heterogeneous models.
Per se, the use of small size laboratory reactors does not guarantee that the effectiveness factors are unity [31]. Instead, intraparticle transport limitations are primarily influenced by the catalyst particle size, while external transport limitations are primarily influenced by the fluid flow velocity. In this work, due to the small size of the catalyst particles employed in the experiments, the generalized Thiele modulus [29] is <1, and thus intraporous mass and energy transport resistances are neglected [29], i.e., the effectiveness factor η is assumed equal to 1 for all the reactions involved, as already reported for similar catalysts and reacting systems [5,32,33,34]. External transport limitations are also neglected [34].
The model is based on local steady-state mass and energy balances coupled to a local kinetics model.
3.1.1. Local Balance Equation
Based on the evaluation of the axial mass and heat Pe numbers (Section 3.1.3.—Evaluation of Pe), axial dispersion is neglected. The mass balance reads as follows [26]:
dFidV=ri=∑jνij rj
where V is the reactor volume, and rj is the kinetics of the jth reaction.
In the energy balance, the temperatures of gas and solid are assumed to be the same, and effective heat transport properties lump all the many heat-transfer processes. The energy balance includes the convective term of gas temperature variation along the reactor axial coordinate, the enthalpy variation of the reactions, and the heat transfer between the furnace and the PBR [26]:
∑iFi cPidTdV=∑j(−ΔHj)rj+Ua(Tfurnace−T)
where cpi is the specific heat of the ith component, T is temperature, Tfurnace is the furnace temperature, ΔHj is the enthalpy variation of reaction j, U is the global heat transfer coefficient and a is the ratio between tube external surface and volume:
a = 4/d
The mass and energy balances form a system of coupled ordinary differential equations (ODEs). Boundary conditions are composition, temperature and flow rate of the gaseous mixture at the PBR inlet. The temperature of the gaseous mixture at the reactor inlet is assumed to be equal to the furnace temperature.
3.1.2. Heat Transport Properties
The local energy balance Equation (8) is a one-dimensional model of the axial temperature profile including radial heat transfer, calculated on the basis of a single overall heat transfer coefficient U. Previous works [35,36] adopting a similar approach, pointed out that U must be calculated considering all the radial heat transfer phenomena, which, in our case, consist of three mechanisms in series: (i) convection inside the tube, (ii) conduction in the reactor wall and (iii) radiation between the furnace and the reactor tube. Radiation is considered as ideal, i.e., the external tube wall temperature is assumed identical to the furnace temperature. Thus, only the heat transfer phenomena (i) and (ii) are considered, and their resistances are lumped in one single parameter, i.e., the overall heat transfer coefficient U, which is defined as [37]:
U=(1h+skw)−1
where h is the convective heat exchange coefficient inside the tube, kw is the thermal conductivity of the tube wall, and s is its thickness. In turn, the convective heat exchange coefficient h depends on the fluid (gas) chemical–physical and transport properties, and also on the flow pattern (i.e., laminar or turbulent). The Nusselt (Nu), Reynolds (Re) and Prandtl (Pr) dimensionless numbers are defined as follows [38,39]:
Re= ρf·us·Lμf
Pr=cPf·μfkr
Nu=h·Lkr
where ρf is the gas density, us is the superficial velocity, μf is the dynamic viscosity, cPf is the gas specific heat, kr is the effective radial thermal conductivity of the packed bed and h is the convective heat exchange coefficient. L is a characteristic dimension, which, in PBRs, is the catalyst-particle diameter. The effective radial thermal conductivity kr was evaluated considering heat conduction within the fluid, between the fluid and the solid catalyst particles, and among the solid particles. The literature correlation proposed by Zehner and Schlünder [40,41] was adopted:
kr=kf·{(1−1−ε)+21−ε1−Bκ−1[B(1−κ−1)(1−Bκ−1)2ln(κB)−B−11−Bκ−1−B−12]}
where kf is the fluid (gas) thermal conductivity, ε is the bed void degree (assumed as 0.4 [42]), k is the ratio between solid and gas thermal conductivities (ks/kf) and B is a parameter defined, for a bed of spheres, as follows:
B=1.25(1−εε)1.11
The thermal conductivity of the gaseous mixture kf is calculated as the molar weighted average of the single components thermal conductivities [40]. Since the packed bed is formed by crushed catalyst highly dispersed in silica glass (Section 2.3.—Laboratory Reactors), the solid thermal conductivity ks was approximated to the value of the silica glass thermal conductivity. Additionally, the tube wall was in silica glass, and thus kw = ks = 1.6 W m−1 K−1 (from [43]).
Several correlations are proposed in the literature to evaluate Nu as a function of Re and Pr [44,45,46,47,48,49,50,51]. In this work, the correlation proposed by Singhal et al. [51] was adopted:
Nu=2.67+0.53·Re0.77·Pr0.53
Singhal et al. [51] developed Equation (16) fitting the results of numerical simulations of heat transfer in packed beds of spherical particles. Equation (16) does not include any void fraction dependency since there is a little variability in the void fraction of realistic packings of monodisperse spheres [52]. Furthermore, the limiting behavior for Re → 0 is demonstrated to be physically sound [52].
3.1.3. Evaluation of Pe
The calculation of the axial mass Péclet number Pem,ax is based on the evaluation of the expression [53]:
Pem,ax=us·lε·Dax
where l is the PBR length, and Dax is the axial mass diffusivity. Since the evaluation of Dax is not straightforwardly related to molecular diffusivity, and only a limited number of experimental investigations are reported in the literature [52], Gunn [54] developed a simple correlation between Pem,ax and Pemol, where Pemol is the Pe number evaluated from Equation (17) in which Dax is substituted with the molecular diffusivity Dmol. The axial mass Péclet number was obtained from [55]:
Pem,ax=(1τ·1Pemol+12)−1.
where the tortuosity τ was determined using the theoretical model of Lanfrey et al. [55]:
τ=ε(1−ε)4/3
The evaluation of Pem,ax for the reactors under consideration in this work gives Pem,ax ≅ 2. According to [29,56,57], this value indicates that a regime is established where convective mass transport is dominating over axial dispersion, which justifies the assumption of neglecting the axial dispersion term in the mass balances (Equation (7)).
The calculation of the axial thermal Péclet number Peth,ax is based on the evaluation of [40]:
Peth,ax = Reax ∙ Prax
where Reax and Prax are the Reynolds and Prandtl numbers evaluated in the axial direction (Equations (11) and (12), where the characteristic length L is the tube diameter d, and the thermal conductivity is the axial one kax). Several studies are proposed in the scientific literature for the evaluation of kax [58,59,60,61]. In this work, the literature correlation proposed by Votruba et al. [61] was adopted:
kax=kf·κ[0.28−0.757·log(ε)−0.057log(κ)]
For the reactors under consideration in this work, Peth,ax ranges between 2 and 40 depending on the operating conditions. According to [28], this value indicates that a regime is established where convective energy transport is dominating over axial dispersion, and justifies the assumption of neglecting the axial dispersion term in the energy balance (Equation (8)).
3.1.4. Ni-Based Catalyst Kinetics
The reaction scheme proposed by Xu and Froment [5] was adopted, based on the reaction rates of the SMR, WGS and GRR reactions. Reaction rates are reversible and of the Langmuir-Hinshelwood type:
rSMR=kSMRpH22.5·(pCH4·pH2O−pH23·pCOKp,SMR)/DEN2
rWGS=kWGSpH2·(pCO·pH2O−pH2·pCO2Kp,WGS)/DEN2
rGRR=kGRRpH23.5·(pCH4·pH2O2−pH24·pCO2Kp,GRR)/DEN2
DEN=1+KCO·pCO+KH2·pH2 +KCH4·pCH4+KH2O·pH2OpH2
with the kinetic parameters kj and adsorption parameters Ki expressed as reported below:
kj=kj, Tr·exp[−EjR·(1T−1Tr)] j=MSR, WGS, GRR
Ki=Ki, Tr·exp[−ΔHiR·(1T−1Tr)] i=CO,H2,H2O,CH4
where Tr is a reference temperature. The kinetic parameters appearing in Equations (22) to (27) were discussed in Section 3.3.—Kinetic Parameters. Their values are reported in Table 2 and Table 3.
3.2. P-R1 and P-R2 Model
The exit line from the reactor was also simulated. In this case, the empty pipes P-R1 and P-R2 were simulated by applying a steady-state, non-isothermal plug flow reactor (PFR) approach, consisting in practice of the same local mass (Equation (7)) and energy (Equation (8)) balances employed for the PBR simulation. The global heat transfer coefficient U was also evaluated in P-R1 and P-R2, using the same definition employed for the PBR reactor (Equation (10)). For Nu, Re and Pr, the definitions given in Equations (11)–(13) apply, with L being the tube hydraulic diameter. The evaluation of h was based on the evaluation of Nu. For an empty tube, with constant wall temperature and with laminar fluid flow regime as in the present case, Nu is constant, equal to 3.66 [37]. The effective thermal conductivity coincides with the thermal conductivity of the gaseous mixture kf [39].
The mass and energy balances of P-R1 and P-R2 form a system of coupled ODEs. Boundary conditions are the inlet composition, temperature and flowrate of the gaseous mixture, which coincide with those at the exit of the previous reactor section.
3.2.1. Gas Phase WGS Kinetics
Downstream from the catalytic bed, the reacting mixture flows through the post-reactor P-R1. Here, according to the literature [18,19], the quartz wool and the silica glass tube were not expected to catalyze the reactions. However, WGS was expected to occur in the gas phase, and the following kinetics was adopted [17,19]:
rWGS=k′WGS·CCO·CH2O(1−pH2·pCO2pCO·pH2O·Kp,WGS)
where:
k′WGS=k′WGS,0·exp[−E′WGSRT]
with k′WGS,0 = 7.4 × 1011 (cm3 mol−1)0.5 s−1 and E′WGS = 288.3 kJ mol−1 [17,19].
3.2.2. Metal Tube WGS Kinetics
Downstream from P-R1, the gaseous mixture flows through the metallic post-reactor P-R2. It is widely reported in the literature that the WGS reaction is catalyzed by metallic piping [17]. For the WGS kinetic expressions, it was proposed here to apply the WGS kinetic equation proposed by Xu and Froment (Equation (23)), withkWGS, Tr= 4.7 kmol bar−1 m−2 h−1, referred to the unit area of inner surface of the tube.
An alternative option follows the approach proposed by Bustamante et al. [18,19] who reported a very fast reaction rate measured with Inconel® 600 piping (typical composition 14–17 wt % Cr, 72 wt % Ni, 6–10 wt % Fe, 1 wt % Mn, 0.5 wt % Cu, 0.5 wt % Si, 0.15 wt % C and 0.015 wt % S), for which they assumed the following kinetic expression [62]:
rWGS=k″WGS·CCO(1−pH2·pCO2pCO·pH2O·Kp,WGS)
where:
k″WGS=k″WGS,0·exp[−E″WGSRT]
with k″WGS,0 = 4.26 × 106 cm2 s−1 and E″WGS = 102.4 kJ mol−1 [17,19].
3.3. Kinetic Parameters
The parameters of the local kinetic model were evaluated taking as a reference the values proposed by Xu and Froment [5] and reported in Table 2. Here, in order to adapt them to the specific catalysts under study, they were adjusted through the multiplicative coefficients ζ1 to ζ6 until satisfactory agreement was found with the experimental data. The model was formed by a system of ODEs, coupled to non-linear analytical equations, requiring numerical integration (Section 3.5.—Numerical Integration). Due to the intrinsic difficulties in implementing an automatic algorithm of minimization of the model error relative to the experimental data, the procedure adopted to perform the adjustment is a trial-and-error procedure. The full set of experimental data, collected in a wide range of operating conditions, used to perform the adjustment of the coefficients ζ1 to ζ6, was reported and discussed in a companion paper [63]. Considering the full set of results, the error between simulated and experimental data was 8.6% with the HMMC catalyst, 11.8% with the 57-4Q catalyst and 9.7% with the 25-4Q catalyst [63].
The proposed values for the coefficients ζ1 to ζ6 are reported in Table 3. ζ6, i.e., the adjusting coefficient for the H2O adsorption coefficient, was about two orders of magnitude higher for the commercial catalysts than for the laboratory-made catalyst, and this is in line with the expectations, considering that the commercial catalysts have Ca and K added in order to promote H2O adsorption. Regarding the activation energies of the SMR and WGS reactions, the adjusting coefficients ζ4 and ζ5 proposed in this work are similar to those published in [33], where values corresponding to ζ4 = 1.07 and ζ5 = 1.33 are reported. In the same work, the activation energy of the GRR reaction is was multiplied by a factor of 0.97, while in our case the multiplying coefficient was 1. The catalyst employed in [33] is a 18 wt % NiO supported on α-Al2O3 provided by Johnson Matthey Plc, tested for SMR and GRR at temperatures around 600 °C, for WGS around 300 °C. Thus, we can conclude that the values proposed here for the adjustable coefficients of the local kinetics model are consistent with those reported in the literature for similar catalysts.
3.4. Thermodynamic Equilibrium
For both the WGS and the SMR reactions, Kp,WGS and Kp,SMR were obtained through the equation below [64]:
Kp=A·exp(−EKR·T)
The values of parameters A and EK are reported in Table 4. Equation (32) provides a fit of the data of Kp provided in [65].
For GRR:
Kp,GRR = Kp,SMR Kp,WGS
3.5. Numerical Integration The model equations form a system of ODEs coupled to non-linear analytical equations, requiring numerical integration, which was performed through MATLAB® R2018a. The system of coupled ODEs (local mass and energy balances) was solved through the ode45 solver. Input data of the numerical tool were the reactor geometrical data, flow rate and composition of the reacting gas mixture and the furnace temperature. The reactor simulation, for one single operating condition, typically ran in less than 1 s on a workstation Dell Precision Tower 3420, with Intel Xeon CPU running at 3.60 GHz with 64 GB of RAM. 4. Results and Discussion
Table 5 reports the values of the main chemical-physical parameters and dimensionless numbers in the operating conditions adopted in this study. The simulated profiles of compositions and temperature along the reactor are presented and discussed in Section 4.1.—MCO2 Laboratory Reactor, and in Section 4.2.—SMR Laboratory Reactor.
4.1. MCO2 Laboratory Reactor
Figure 3 reports simulation results for the MCO2 laboratory experimental reactor, obtained with a furnace temperature of 623 K. Modeling results are reported for all the three catalysts employed in this work. Experimental data of compositions, measured at the sampling point, are also reported.
The first remark is that, with the JM catalysts, the simulated SMR and revSMR reaction rates exhibited high sensitivity towards the water content in the reactant gas mixture. In the MCO2 mode, where no H2O was present in the reactor feeding gas, the revSMR and revWGS reactions were thus very fast at the reactor inlet, especially with the 57-4Q catalyst. This is clearly visible in Figure 3, particularly in the CH4 molar fraction plot (Figure 3a), showing a steep increase at the reactor inlet, higher than with both the HMMC and the 25-4Q catalysts. Further down in the reactor, water, being a reaction product, was present, as displayed in Figure 3e. Therefore, the rate of formation of CH4 over the 57-4Q catalyst visibly slowed down, so that at the exit of the PBR the CH4 molar fraction was only slightly higher than that obtained with the HMMC catalyst. This was confirmed by the experimental results of CH4 molar fraction at the exit of the reactor: 5.1% with the 57-4Q, 5.0% with the HMMC catalyst and 2.9% with the 25-4Q catalyst.
Regarding the simulated temperature profiles (Figure 3f), the fast kinetics of the exothermal revSMR at the entrance of the PBR caused a steep temperature increase of about 10 K with the 57-4Q catalyst. A similar behavior, with a smaller temperature increase, was also found with the other two catalysts. Further down in the reactor, due to the reduced rate of the revSMR, the temperature of the PBR cooled down and at the exit of the PBR it was practically equal to the temperature of the furnace for all three simulated catalysts. Additionally in the downstream P-R1, where only gas phase WGS was considered to occur, the temperature remained equal to the furnace temperature and no appreciable variations in composition were found (results not reported in Figure 3f). In the subsequent post-reactor P-R2, due to the efficient heat exchange with the surrounding environment, the temperature quickly dropped to the ambient temperature (298 K) in a narrow length of about 7 mm. More in detail, the temperature decreased to 523 K in a short length of 0.3 mm, and then to 366 K in the subsequent 1 mm of length. Here, the simulation model, as discussed in Section 3.2.2.—Metal Tube WGS Kinetics, embedded a simulation of the revWGS reaction, catalyzed by the metallic piping. However, due to the low furnace operating temperature adopted in the MCO2 operating mode, and also to the almost instantaneous temperature dropped down to the ambient temperature occurring in the post-reactor P-R2, no variation in composition was calculated by the model in the P-R2. No results are reported in Figure 3 for the P-R1 or P-R2 simulation results, and, in Figure 3, the experimental compositions measured at the sampling point (exit of P-R2) were plotted at the exit of the PBR section.
4.2. SMR Laboratory Reactor
Figure 4 reports the simulation results obtained for the SMR experimental reactor, with furnace temperature 853 K and GHSV = 2.1 × 104 h−1. The left-hand sides of the figures show the modeling results obtained in the catalytic section of the reactor (PBR), while the right-hand sides display the P-R2 modeling results.
In this case, 5% CH4 was fed in together with 20% steam, and therefore the simulated reaction rate with the 25-4Q and 57-4Q catalysts was hindered by the presence of water, and thus was visibly slower than in HMMC operating mode. In particular, with the HMMC catalyst, thermodynamic equilibrium was established after about 2 mm from the entrance of the PBR, corresponding to an almost complete consumption of methane and to the generation of an H2 molar fraction of about 16%. Conversely, with the commercial catalysts, reaction rates were much lower and thermodynamic equilibrium was not reached within the PBR; H2 molar fractions at the exit of the PBR were 14.5% and about 12% with 57-4Q and 25-4Q respectively.
Regarding the simulated temperature profiles within the PBR section, the fast reaction rate at the entrance of the HMMC PBR was associated with a temperature drop, after which, the temperature rose to the furnace temperature. This deep temperature drop may be somehow overestimated due to the absence of the thermal axial dispersion term in the model; nevertheless, this possible inaccuracy did not influence the evaluation of the compositions at the exit of the PBR since thermodynamic equilibrium was reached quickly with the HMMC catalyst in these operating conditions. On the other hand, with the commercial catalysts, as a result of the slower kinetics, the temperature drop was smoother and milder (20 K with the 57-4Q and 15 K with the 25-4Q catalysts respectively). For all catalysts, the temperature at the exit of the PBR was equal to the furnace temperature. Like with the MCO2 reactor, also with the SMR reactor, no temperature or composition variations were evaluated by the simulation model in the post-reactor PR-1, which was not represented in Figure 4. Instead, in the post reactor PR-2, the model calculated variations not only in temperature but also in composition. This is due to the fact that, in the SMR operating mode, on the one hand, the furnace temperature was higher than in the MCO2 operating mode; on the other hand, the P-R2 was exposed to an external temperature of 523 K, higher than in the MCO2 case. As a result of the higher temperatures established in the P-R2, the simulation results indicate a non-negligible advancement of the WGS reaction. The WGS kinetics in P-R2, as already mentioned in Section 3.2.2.—Metal Tube WGS Kinetics, was based on the kinetic equation proposed by Xu and Froment (Equation (23)), withkWGS, Tr = 4.7 kmol/(bar m2 h), referred to the unit area of inner surface of the tube. Figure 4 reports the simulation results only for the first 3 mm of PR-2, where the temperature was simulated to drop to the external temperature of 523 K (Figure 4f), and where the model evaluated a slight, yet non negligible, advancement of the WGS reaction. Figure 4 shows good agreement between simulated and experimentally measured molar fractions at the sampling point (exit of P-R2).
For the sake of completeness, other literature WGS kinetic models were also tested for the simulation of P-R2. The kinetic model proposed by Bustamante et al. [18,19] (Section 3.2.1.—Gas Phase WGS Kinetics) for the gas phase WGS reaction was implemented in the simulation of the P-R2 discharge pipe, without giving any advancement in WGS conversion. A further simulation was performed employing the kinetic model proposed by Bustamante et al. [18,19] (Section 3.2.2.—Metal Tube WGS Kinetics) for the WGS reaction occurring in an Inconel® 600 pipe. This kinetic model, when implemented into the P-R2 discharge pipe, results in thermodynamic equilibrium being reached for the WGS reaction inside P-R2 (Kp,WGS = 86.7 at 523 K), exceeding by far the real WGS conversion measured experimentally at the exit of the P-R2 piping. This indicates that the WGS kinetics in the Inconel® 600 piping was faster than in the AISI 316L tubes, probably due to the higher amount of Ni (72 wt % in Inconel® 600, versus 8 wt % in the AISI 316L).
5. Conclusions The conclusions of the study reported in this paper can be summarized as follows:
- A steady-state, pseudo-homogeneous 1D non-isothermal model was proposed, which was specifically developed for the laboratory experimental reactors, and was intended to guide the development of small size SMR and MCO2 reactors in the context of P2X studies.
-
The model was based on local mass and energy balances, written in the typical form adopted for PBRs, coupled to a local kinetics. The Xu and Froment local kinetics model [5] was applied.
- Simulated results of composition and temperature profiles along the reactors show that non-uniformities were significant in the laboratory experimental reactors.
-
Progress of the WGS reaction in the discharge piping was studied, which is more easily included in the simulation model rather than prevented in the experimental apparatus [17,18,19]. The WGS reaction occurring in the metal piping downstream of the reactor was well simulated through the corresponding kinetic equation proposed by Xu and Froment [5].
Dimensions (10−3 m) | MCO2 Reactor | SMR Reactor |
---|---|---|
H length | 5 | 5 |
PBR length | 26 | 12 |
P-R1 length | 150 | 50 |
P-R2 length | 300 | 40 |
d | 6 | 6 |
s | 1 | 1 |
Parameter | Value | Notes |
---|---|---|
kSMR,Tr | 1.842 × 10−4 (kmol bar0.5 kgcat−1 h−1) | Tr = 648 K |
kWGS,Tr | 7.558 (kmol bar−1 kgcat−1 h−1) | Tr = 648 K |
kGRR,Tr | 2.193 × 10−5 (kmol bar0.5 kgcat−1 h−1) | Tr = 648 K |
ESMR | 240.1 (kJ mol−1) | |
EWGS | 67.13 (kJ mol−1) | |
EGRR | 243.9 (kJ mol−1) | |
KCH4,Tr | 0.1791 (bar−1) | Tr = 823 K |
KCO,Tr | 40.91 (bar−1) | Tr = 648 K |
KH2,Tr | 0.02960 (bar−1) | Tr = 648 K |
KH2O,Tr | 0.4152 (-) | Tr = 823 K |
ΔHCH4 | −38.28 (kJ mol−1) | |
ΔHCO | −70.65 (kJ mol−1) | |
ΔHH2 | −82.90 (kJ mol−1) | |
ΔHH2O | 88.68 (kJ mol−1) |
HMMC | 57-4Q | 25-4Q | ||
---|---|---|---|---|
ζ1⋅kSMR | ζ1= | 15 | 25 | 8 |
ζ2⋅kWGS | ζ2= | 15 | 100 | 130 |
ζ3⋅kGRR | ζ3= | 15 | 10 | 10 |
ζ4⋅ESMR | ζ4= | 1 | 1.17 | 1.22 |
ζ5⋅EWGS | ζ5= | 1 | 1.5 | 1.7 |
ζ6⋅KH2O | ζ6= | 1 | 90 | 95 |
WGS | SMR | |
---|---|---|
A | 1.412 × 10−2 (-) | 7.846 × 1012 (atm2) |
EK | 220.2 (kJ mol−1) | −37.72 (kJ mol−1) |
Variable | MCO2 | SMR | Units | ||
---|---|---|---|---|---|
T = 473 K | T = 773 K | T = 773 K | T = 1173 K | ||
GHSV @ RTP | 6.7 × 103 | 2.4 × 104 | h−1 | ||
uS | 0.082 | 0.133 | 0.200 | 0.459 | m s−1 |
μ | 2.03 × 10−5 | 2.85 × 10−5 | 3.56×10−5 | 6.49 × 10−5 | kg m−1 s−1 |
cp | 5.06 | 5.16 | 4.52 | 4.71 | kJ kg−1 K−1 |
ρ | 0.305 | 0.187 | 4.74 × 10−2 | 2.06 × 10−2 | kg m−3 |
kf | 9.89 × 10−2 | 1.42 × 10−1 | 2.42 × 10−1 | 0.464 | W m−1 K−1 |
ks | 1.57 | 1.99 | 1.99 | 3.39 | W m−1 K−1 |
kax | 0.945 | 1.25 | 1.46 | 2.57 | W m−1 K−1 |
Re = ρ uS d/μ | 7.36 | 5.243 | 1.600 | 0.878 | - |
Pr = μ cp/kf | 1.04 | 1.03 | 0.665 | 0.659 | - |
Nu | 3.43 | 2.28 | 2.91 | 2.83 | - |
U | 380 | 503 | 584 | 1019 | W m−2 K−1 |
Author Contributions
Conceptualization, P.C.; Formal analysis, P.C. and F.P.; Investigation, T.C. and G.G.; Methodology, P.C. and F.P.; Software, P.C. and F.P.; Supervision, G.B.; Validation, P.C. and G.G.; Writing-original draft, P.C. and G.B.; Writing-review and editing, F.P. and G.G. All authors have read and agreed to the published version of the manuscript.
Funding
This research received no external funding.
Acknowledgments
The authors acknowledge Agnese Buscaldi and Beatrice Antonucci for the experimental tests made during their Master and Bachelor's theses respectively.
Conflicts of Interest
The authors declare no conflict of interest.
Abbreviations
GHSV | Gas Hourly Space Velocity | |
GRR | Global Reforming Reaction | |
HMMC | Homemade Methanation Catalyst | |
JM | Johnson Matthey | |
MCO | Methanation of CO | |
MCO2 | Methanation of CO2 | |
ODE | Ordinary Differential Equation | |
PBR | Packed Bed Reactor | |
P-R | Post Reactor | |
P2X | Power-to-X | |
revGRR | reverse of Global Reforming Reaction | |
revSMR | reverse of Steam Methane Reforming | |
revWGS | reverse of Water Gas Shift | |
RTP | Reference temperature and pressure (273.15 K and 101.325 kPa) | |
SMR | Steam Methane Reforming | |
SR | Steam Reforming | |
WGS | Water Gas Shift | |
Symbols used | ||
A | pre-exponential factor in Equation (32) | (-) for WGS-(atm2) for SMR |
a | surface to volume ratio | (m-1) |
B | parameter in Equation (15) | (-) |
Ci | molar concentration of the ith component | (mol m-3) |
cpi | specific heat for the ith component | (kJ mol-1 K-1) |
D | mass diffusivity | (m2 s-1) |
d | tube diameter | (m) |
dP | particle diameter | (m) |
E,E',E" | activation energies | (kJ mol-1) |
EK | energy parameter in Equation (32) | (kJ mol-1) |
Fi | ith component molar flowrate | (mol s-1) |
h | convective heat exchange coefficient | (W m-2 K-1) |
ΔH°298K | standard reaction enthalpy | (kJ mol-1) |
ΔHi | adsorption enthalpy of component i (Equation (27)) | (kJ mol-1) |
KCH4,KCO,KH2 | adsorption constants in Equation (27) | (bar-1) |
KH2O | adsorption constant for H2O in Equation (27) | (-) |
Kp,GRR,Kp,SMR | equilibrium constants of GRR and SMR | (Pa2) |
Kp,WGS | equilibrium constant of WGS | (-) |
kGRR,kSMR | kinetic rates of GRR and SMR | (kmol bar0.5 kgcat-1 h-1) |
kWGS | kinetic rate of WGS | (kmol bar-1 kgcat-1 h-1) |
k'WGS | kinetic rate of WGS in Equation (28) | (m1.5 mol-0.5 s-1) |
k"WGS | kinetic rate of WGS in Equation (30) | (m2 s-1) |
k | thermal conductivity | (W m-1 K-1) |
L | characteristic dimension | (m) |
l | length | (m) |
Nu | Nusselt number | (-) |
pi | partial pressure of the ith component | (bar) |
Pem | mass transfer Péclet number | (-) |
Peth | thermal Péclet number | (-) |
Pr | Prandtl number | (-) |
rj | reaction rate of the jth reaction | (mol m-3 s-1) |
R | ideal gas constant | (J mol-1 K-1) |
Re | Reynolds number | (-) |
s | tube wall thickness | (m) |
Si | selectivity to product i | (-) |
T | temperature | (K) |
U | overall heat transfer coefficient | (W m-2 K-1) |
us | superficial velocity | (m s-1) |
V | volume | (m3) |
XCH4, XCO2 | methane and carbon dioxide conversion | (-) |
Yi | yield in product i | (-) |
Greek letters | ||
ε | catalyst bed void fraction | (-) |
ζi | parameters in Table 3 | (-) |
ηj | effectiveness factor for the jth reaction | (-) |
κ | ratio between solid and gas thermal conductivities | (-) |
ρ | density | (kg m-3) |
νij | stoichiometric coefficient of the ith component in the jth reaction | (-) |
τ | tortuosity | (-) |
Subscripts | ||
ax | axial | |
f | fluid | |
i | component index | |
in | inlet | |
j | reaction index | |
out | outlet | |
r | radial | |
s | solid | |
w | wall |
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2020. This work is licensed under http://creativecommons.org/licenses/by/3.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
To support the interpretation of the experimental results obtained from two laboratory-scale reactors, one working in the steam methane reforming (SMR) mode, and the other in the CO2 hydrogenation (MCO2) mode, a steady-state pseudo-homogeneous 1D non-isothermal packed-bed reactor model is developed, embedding the classical Xu and Froment local kinetics. The laboratory reactors are operated with three different catalysts, two commercial and one homemade. The simulation model makes it possible to identify and account for thermal effects occurring inside the catalytic zone of the reactor and along the exit line. The model is intended to guide the development of small size SMR and MCO2 reactors in the context of Power-to-X (P2X) studies.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer