Format And Type Fonts CHEMICAL ENGINEERING TRANSACTIONS VOL. 39, 2014 A publication of The Italian Association of Chemical Engineering www.aidic.it/cet Guest Editors: Petar Sabev Varbanov, Jiří Jaromír Klemeš, Peng Yen Liew, Jun Yow Yong Copyright © 2014, AIDIC Servizi S.r.l., ISBN 978-88-95608-30-3; ISSN 2283-9216 DOI: 10.3303/CET1439161 Please cite this article as: Silva J.D., 2014, Dynamic simulation of the steam reforming of methane for the production of hydrogen in a catalytic fixed bed membrane reactor, Chemical Engineering Transactions, 39, 961-966 DOI:10.3303/CET1439161 961 Dynamic Simulation of the Steam Reforming of Methane for the Production of Hydrogen in a Catalytic Fixed Bed Membrane Reactor Jornandes D. Silva Polytechnic School - UPE, Laboratory of Environmental and Energetic Technology; Rua Benfica - 455, Madalena, Recife – PE, Brazil, Cep: 50750-470 jornandesdias@poli.br This paper presents an one-dimensional isothermal modelling for a CFBMR (H2 permeable). The work has been developed to simulate the SRM to produce H2. The dynamic simulation for molar flow rates of H2 was analyzed in the permeation side and shell side zones at different temperatures at the exit of the proposed system. The model has allowed the validation for the conversion of CH4 by comparing of optimized values. Additionally, the molar flow rates for chemical species (CH4, H2O, H2, CO and CO2) were analyzed along of the CFBMR (shell side) at t = 14 s, a comparison between permeation side and shell side for H2 and the molar flow rates of H2 in the permeation side and shell side zones at different temperatures. 1. Introduction Hydrogen is widely used as an important chemical feedstock and it is also regarded as a promising future clean fuel. Currently, many chemical, petrochemical and petroleum refining industries and well perceived as a pollution-free primary energy carrier for future transportation fuel as well as electricity generation. The steam reforming of methane (SRM) is the primary hydrogen production route on industrial scale. Methane (CH4) is the main component of conventional natural gas (NG) which can also be obtained from associated gas, stranded NG and shale gas. It was previously determined that three reactions are predominant (Xu and Froment, 1989), the (SRM), water gas shift (WGS) and overall steam reforming of methane (OSRM). - Reaction of the SRM; 10 K298224 molkJ8205HH3COOHCH   .;  (1) - Reaction of the WGS; 10 K298222 molkJ141HHCOOHCO   .;  (2) - Reaction of the OSRM. 10 K2982224 molkJ9164HH4COOH2CH   .;  (3) The steam reforming process of methane (SRPM) for the hydrogen production can be studied using various reactor configurations such as fixed be reactors, fluidized bed reactors and membrane reactors (fixed bed membrane reactors or fluidized bed membrane reactors). The use of a permeable membrane reactor offers the potential for a compact unit in combining both the reaction and separation in a simple unit. Several studies have been carried out over SRPM for separating H2 from a reaction mixture using Pd-based membrane reactors (Marín et al., 2012; Sadooghi and Rauch, 2013). The most of these researches have been experimentally and numerically investigated to produce H2 at high temperature. The main focus of the 962 experimental studies was to develop the robust catalysts which are resistant for carbon deposition. While the numerical works have been reported to simulate the hydrogen production using Pd-based membrane and robust catalysts, respectively. The catalytic fixed bed membrane reactors (CFBMRs) can significantly improve the production of thermodynamically limited reactions. Hydrogen-selective membranes permit hydrogen produced by reactions to escape under the influence of the difference in the hydrogen partial pressures to the permeate side where an inert gas such as nitrogen can be used to purge hydrogen. This type reactor is promising for thermodynamically limited reactions such as steam reforming. In this CFBMRs, the effect of the intraphase transport by considering diffusion inside the catalyst pellets on both sides of the coupled CFBMR. The Fickian diffusion model has been used for the reactants for assuming the diffusion process through a stagnant bulk phase with the diffusivities of components in the reacting mixture along the radial direction inside catalyst pellets. The commercial hydrogen separation by metallic membranes is mainly focused on palladium alloys. Palladium (Pd) membranes have been proposed for application in hydrogen separation for many years due to their exclusive selectivity for hydrogen. Pd is known to be highly soluble for hydrogen and it can tolerate hydrocarbon containing streams. Pd-based membranes also put forward the possibility of hydrogen separation at high temperature aiming at combining a separation process with a chemical reaction in one unit. Membrane application in the chemical, petrochemical and petroleum industries for separation and reaction processes were recently reviewed. The objective of this study is to theoretically investigate the SRM to produce H2 (permeation side and shell side) in the CFBMR at different temperatures at the exit of the system in study. 2. Kinetic mechanism The kinetic models treat the mechanism, rates of chemical reactions and the resulting concentrations of species at some point in time and space of a CFBMR. These models are generally related to specific processes, providing important considerations for relating to the mechanisms of reactions and ways to increase the rate of a reaction or process. The process discussed here has involved the main reactions of SRM. In normal operation, the reaction (1) is the SRM, the reaction (2) shows water-gas shift and the reaction (3) presents the overall SRM. Reactions (1)-(3) are reversible (equilibrium). Thermodynamically, reactions (1) and (3) are strongly endothermic. On the other hand, the reaction (2) is exothermic. Components of reactions (1), (2) and (3) are defined like methane (CH4), water (H2O), carbon monoxide (CO), hydrogen (H2) and carbon dioxide (CO2). The stoichiometric coefficients of components were presented in Table 1. Table 1: Stoichiometric coefficients ij for components of reactions Reactions CH4 H2O CO H2 CO2 1 -1 -1 +1 +3 0 2 0 -1 -1 +1 +1 3 -1 -2 0 +4 +1 2.1 Mathematical modelling on intrinsic kinetic process The kinetic mathematical model applied here considers only three homogeneous reactions (1), (2) and (3) according to Langmuir-Hinshelwood concept (Xu and Froment, 1989).  2 SRM CO 3 H OHCH 52 H SRM SRM DEN K PP PP P k R 2 24 2             . ;  2 WGS COH OHCO H WGS WGS DEN K PP PP P k R 22 2 2           (4)  2 OSRM CO 4 H2 OHCH 53 H OSRM OSRM DEN K PP PP P k R 22 24 2            . ; 2 22 4422 H OHOH CHCHHHCOCO P PK PKPKPK1DEN  (5) Where Rj (j = SRM, WGS and OSRM; DEN is denominator in Eq(4) and (5)) are the rates of reactions (1)-(3). The reaction rate coefficients (kj, j = SRM, WGS and OSRM) and the adsorption constants of gases ( 4CH K , 963 OH2 K , 2H K and COK ) have an Arrhenius type dependence with temperature. On the other hand, Pi (i = CH4, H2O, CO, H2 and CO2) are the partial pressures of the respective components within the catalyst and DEN is a dimensionless parameter. The Arrhenius expressions for the reaction rate coefficients, adsorption constants for the components CH4, H2O, H2 and CO and the equilibrium constants are listed in Table 2. Table 2: Kinetic, adsorption and equilibrium parameters for the reactions involved in steam reforming process of methane (Xu and Froment, 1989) Parameters of Eqs(4) to (5) Units Rate coefficients          s 16 SRM T K87928 10x4909k . exp. ;          s 16 OSRM T K33629 10x2902k . exp. 1150 a hkgkPkmol .          s 4 WGS T K38074 10x3904k . exp. 111 a hkgkPkmol  Adsorption coefficient constants           s 4 CH T K204604 10x6506K 4 . exp. ;          s 3 OH T K3510666 10x7701K 2 . exp. 1 akP            s 11 H T K139971 10x1206K 2 . exp. ;           s 7 CO T K718497 10x2308K . exp. 1 akP  Equilibrium constants          1130 T K83026 7610266K s SRM . . exp. ;          0634 T K4400 K s WGS .exp 2 akP                           0634 T K4400 1130 T K83026 7610266K ss OSRM .exp. . exp. 2 akP The net rates for each component (consumption and formation) type of reactions (1), (2) and (3) have been obtained by using the following equation (Oliveira and Silva, 2012).    n 1j jiji Rr  (6) Where ri are the net rates of that species (i) in reactions (Rj), ij are the stoichiometric coefficients of species and Rj are the reaction (1), (2) and (3) rates. Effectiveness factors i are used to account for the intraparticle transport limitation. Thus, the Eq(6) was rewritten as given below.    n 1j jijji Rr  (7) The Eq. (7) can be used to obtain rCH4, rH2O, rH2, rCO and rCO2. The details of these net rates have given by Oliveira and Silva, (2013). 3. Mathematical modelling A conceptual schematic of the integrated CFBMR was developed to investigate the SRM according to the Figure 1. The inner is the tube is the supported membrane and the outer one is the non-permeable shell. The catalyst (Ni/-Al2O3) is packed in the shell side where the reactants (methane and steam) are introduced. The sweep gas (nitrogen or steam) is fed into the tube side in the co-current flow mode. In this system, three reforming reactions occur in the non-permeable shell by producing light gases (H2, CO and CO2). The mathematical modelling and computer simulation for the CFBMR are in continuous development aiming to improve the knowledge of the phenomenological processes for this simple unit (combining both the reaction and separation). According to its applications, the CFBMRs are applied in the chemical, petrochemical and oil refining industries. For this work, an one-dimensional mathematical model without axial dispersion taking into account diffusional limitations in the solid porous network was developed. This model has been projected like 964 isothermal dynamic pseudo-homogeneous model. The mass balance equations are given by the following partial differential equation system. Figure 1: Schematic representation of the CFBMR for the SRM 3.1 Balance equations for the shell side zone The differential balances of the molar flow rates (Fi) for the model components i ( i = CH4, H2O, H2, CO and CO2) in the reaction side (shell side) are given as follows.   0Jthus,Hiwhen;JR2r1LR z F d4 Q4 t F g u s Hi2 s iShell 3 1j 5 1i iijjbs 2 Shell i 2 shell gis 2             (8) - The initial and boundary conditions for the Equation (8) are given as: productsforFandtstanreacforFFzallforFF tiitiiti 0;; 00,00,0   (9) 0tforFF Enti0zi   ;., ; 0tfor0 z F Lz i     ; (10) 3.2 Balance equation for the permeation side zone The mass balance equation for the molar flow rate (FH2) of H2 in the permeation side zone (inner tube of CFBMR) is written as follows:                                       50 t H 50 s H 12H 2 tubeinner gHs 22 22 PP A T 38767 10x00611 z F d4 Q4 t F g u ... exp.   (11) - The initial and boundary conditions for the Eq(11) are given as: zallfor0F 0t H 2 ;  ; 0tforFF EntH 0z H 22   ;., ; 0tfor0 z F Lz H 2     ; (12) 4. Solution of the model equations using the technique of the Laplace transform The Laplace transform methodology was applied for changing the set of the partial differential equations (PDFs) into a set of ordinary differential equations (EDOs) (Silva and Oliveira, 2013). Then, the new set of EDOs has been integrated by a subroutine based on the Runge-Kutta Gill with automatic step size and double precision to ensure accuracy (Oliveira and Silva, 2013). 5. Results and discussion The mathematical modelling was developed to analyze the molar flow rates of chemical species (CH4, H2O, H2, CO and CO2) on the shell side zone as well as the molar flow rate of H2 into the permeation side zone. The proposed model for this work has been used to relate the evolution of the molar flow rates over space (z) 965 along the CFBMR (shell side zone and permeation side zone). For this simulation, the computational code was fed with parameters shown in Table 3. Table 3: Parameters for the modelling involved in steam reforming of methane Categories Properties Numerical values References Operating conditions Shell side feed inlet temp., Tin,ss(K) 800 Ji et al., (2003) Per. side feed inlet temp., Tin,ps(K) 800 Ji et al., (2003) Shell side pressure Pss, (bar) 14 Ji et al., (2003) Permeation side pressure, (bar) 14 Ji et al., (2003) Feed composition Inlet flow rate of CH4, FCH4 (kmol s -1 ) 2.5x10 -6 Ji et al., (2003) Inlet flow rate of H2O, FH2O (kmol s -1 ) 3.5x10 -6 Ji et al., (2003) Inlet flow rate of H2, FH2 (kmol s -1 ) 0 Ji et al., (2003) Inlet flow rate of CO, FCO (kmol s -1 ) 0 Ji et al., (2003) Inlet flow rate of CO2, FCO2 (kmol s -1 ) 0 Ji et al., (2003) CH4/H2O ratio 0.667 Ji et al., (2003) Inlet partial pres. of CH4, pCH4 (bar) 5.46 Pantoleontos et al.,(2012) Inlet partial pres. of H2O, pH2O (bar) 18.34 Pantoleontos et al.,(2012) Inlet partial pres. of H2, pH2 (bar) 0.68 Pantoleontos et al.,(2012) Inlet partial pres. of CO, pco (bar) 0.00 Pantoleontos et al.,(2012) Inlet partial pres. of CO2, pCO2 (bar) 0.31 Pantoleontos et al.,(2012) Gas solid Porosity of catalyst bed, b (dimen.) 0.43 Ji et al., (2003) Catalyst density, B (kg m -3 ) 2,100 Ji et al., (2003) Effect. Factors (i) for Eqs. 1, 2, 3 0.07, 0.07 and 0.57 Xu and Froment, (1989) Shell side radius, Rshell (m) 3.00x10 -3 Ji et al., (2003) Thickness of dense membrane, (m) 5.0x10 -5 Ji et al., (2003) (a) (b) 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 F i ( k m o l s -1 ) Length of the CFBMR (m) at t=14sec H 2 , H 2 O, CO CO 2 , CH 4 (c) 0 6 12 18 24 30 0.0 0.2 0.4 0.6 F H 2 (k m o l s -1 ) Permeation side zone Shell side zone z = 1m Time (sec) 1070K 1070K (d) 0 5 10 15 20 0.0 0.2 0.4 0.6 0.8 1.0 950K 950K 1150K1250K 1150K 1250K F H 2 (k m o l s -1 ) H 2 psz , H 2 psz and H 2 psz H 2 ssz , H 2 ssz and H 2 ssz Time (sec) 966 Figure 2: (a) Validation of the model, (b) Profiles of chemical species (CH4, H2O, H2, CO and CO2) along the CFBMR, (c) Evolutions of molar flow rates in the permeation side and shell side zones and (d) Evolutions of molar flow rates in the permeation side and shell side zones at different temperatures Figure 2(a) shows a validation before simulating the CFBMR, the main program which is used to simulate the conventional CFBMR has been checked with data of the literature. The simulation results for conversion of CH4 and data for different authors were validated by comparison. Figure 2(b) reports the profiles of CH4, H2O, H2, CO and CO2 versus the space variable along the CFBMR (shell side zone). Figure 2(c) shows the evolution of the molar flow rates in the permeation side and shell side zones versus the time variable at the exit. Figure 2(d) presents the evolutions of the molar flow rates of H2 at different temperatures (1,250K, 1,150K and 950K) in the permeation side and shell side zones. 6. Conclusions Conducted in the context of reform of the line of chemical components, this research will resort numerical methodology in order to carry out the development process in isothermal catalytic membrane reactor dynamics of fixed bed. In conditions allowed for this research, the method of Runge-Kutta Gill was used to predict the model components (CH4, H2O, H2, CO and CO2). The development of computer code to process and analyze the behavior of the variables in this research allowed the lead the following conclusions: - The validation confirmed by comparing which the results of this research and the results obtained by different authors are good agreement; - Consumption and production have shown the profiles of chemical species (CH4, H2O, H2, CO and CO2) along the CFBMR; - The molar flow rates of H2 at 1,070K shown an evolution with the time. Once, The amount of H2 is most in the permeation side zone than in the shell side zone; - The molar flow rates of H2 at different temperatures (1,250K, 1,150K and 950K) have shown an evolution with the time. Once, the amount of H2 is most in the permeation side zone than in the shell side zone. References Ji P., Van Der Kooi H.J., De Swaan Arons J., 2003, Simulation and thermodynamic of an integrated process with H2 membrane CPO reactor for pure H2 production, Chemical Engineering Science, 58(17), 3901- 3911. Marín P., Patiño Y., Díez F.V., Ordóñez S., 2012, Modelling of hydrogen perm-selective membrane reactors for catalytic methane steam reforming, International Journal of Hydrogen Energy, 37, 18433-18445. Oliveira C.C.B., Silva J.D., 2012, Dynamic modelling of the gasification region of a bubbling fluidized bed gasifier, Chemical Engineering Transactions, 29, 841-846. Oliveira C.C.B., Silva J. D., 2013, Mathematical modelling of the steam reforming of toluene for fuel gas production in a fixed bed catalytic reactor, Chemical Engineering Transactions, 35, 307-312. Pantoleontos G., Kikkinides E.S., Georgiadis M.C., 2012, A heterogeneous dynamic model for the simulation and optimisation of the steam methane reforming reactor, 37, 16346 - 16358. Sadooghi P., Rauch R., 2013, Pseudo heterogeneous modeling of catalytic methane reforming process in a fixed reactor, Journal of Natural Gas Science and Engineering, 11, 46-51. Silva J.D., Oliveira C.C.B., 2013, Mathematical modelling for the adsorption process CO2 in nanopores of catalytic particles in a fixed bed reactor using numeral Laplace transform, Chemical Engineering Transactions, 35, 829-835. Xu J, Froment G.F., 1989, Methane steam reforming, methanation and water-gas shift: I. Intrinsic kinetics, AIChE Journal, 35, 88-96.