Modeling of Porous Catalyst Pellets: Comparison of Diffusion Flux Models for Steam Methane Reforming

P. L. Mills [1], L. Olabanji [1], A. Nagaraj [1], A. Nanduri [1],
[1] Department of Chemical and Natural Gas Engineering, Texas A&M University - Kingsville, Kingsville, TX, USA
发布日期 2015

Introduction Steam reforming of natural gas has been the most common method for producing synthesis gas (CO + H2) for the production of H2, MeOH and NH3 for over half a century [1]. Production of H2 in the USA using steam methane reforming (SMR) is ca. 9 MM tons/yr. Global production of NH3 using H2 derived from SMR was ca. 109 MM/yr in 2004. In the chemical industry, the SMR process is carried out in tubular furnaces where a Ni-supported catalyst is placed in tubes that are typically made from Ni-Cr alloys. The primary reactions that occur include methane reforming, water-gas-shift, and the methanation reactions.

CH4 + H2O ⇌ CO + 3H2 ΔH298 = 206.2 kJ/mol CO + H2O ⇌ CO2 + H2 ΔH298 = -41.5 kJ/mol CH4 + 2H2O ⇌ CO2 + 4H2 ΔH298 = 164.7 kJ/mol

The catalyst pellets can vary in both size and geometry, and the compositions is usually Ni-on-Al2O3 [1, 4]. Some typical catalyst shapes are illustrated in Figure 1. Comparison of diffusion flux models to describe species transport-kinetic interactions for SMR has been reported only for a spherical catalyst shape[2]. Dixon and coworkers [5] employed CFD to compare the effect of various catalyst geometries on heat transfer, pressure drop and methane conversion. However, a detailed performance comparison for various catalyst shapes using different flux models has not been examined.
The primary objective of this study is to describe a general modeling framework for various catalyst shapes using different diffusion flux models. Use of Comsol Multiphysics A 1D steady state model for a spherical catalyst shape is developed in COMSOL Multiphysics® software using Transport of Dilute Species and Heat conduction interfaces. The transport-kinetics interactions in the porous catalyst are analysed using various diffusion flux models such as Wilke, Wilke-Bosanquet, Stefan-Maxwell, and Dusty gas models. The kinetic equations from Xu and Froment are used to describe the reaction rates [3].
The steady state results using the various diffusion flux models are shown (Figure 2). The species mole fractions obtained using both the mass and mole formulations of the pellet model will be compared for each catalyst shape and for various diffusion flux models.