Study of Methane Oxidation in a Monolith Reactor: Comparison between a 1D and a 3D Modelling Approach Andrea Celani 1,2 , Sam Wilkinson 1 , Li Liu 1 , Hugh Stitt 1* 1 School of Chemical Engineering, University of Birmingham, B15 2TT 2Johnson Matthey Technology Centre, Billingham, TS23 1LB *Corresponding author: Hugh.Stitt@matthey.com Highlights Study of methane oxidation reaction in a monolithic reactor using a computational approach Comparison of results for a CFD-based model with a conventional 1D+1D model Good agreement between the two approaches 1. Introduction CFD based modelling including diffusion and reaction within a meshed porous solid to represent a heterogeneous catalyst has seen a rapid growth in interest over the last decade. Early attempts to do this essentially represented the reaction simply as a heat source terms [1]. Developments in commercial and open source CFD sub-models and methods allowed meshing of the porous solid and the development of functions to model species diffusion and reaction within the porous phase. This approach, sometimes referred to as “hi-fidelity” reactor modelling, is now being widely published [2][3][4][5]. There remains though a shortage of rigorous model validation and verification. A significant number of published validations apply only a single “global” result to the validation (e.g. temperature, reaction conversion etc.) [6] which does not address the validity of the complex substructure of the overall model [7]. 2. Methods This presentation will report on a detailed comparison of results for a CFD-based model of a catalyst coated channel (a single square 1mm channel of a washcoated monolith with length of 1 cm) with those from a conventional chemical reaction engineering (CRE) 1D+1D model. Methane oxidation, at low concentrations, was used as an example reaction and experimental data for methane oxidation over a 1 wt.% Pd/Al2O3 catalyst were used to derive a Langmuir Hinshelwood kinetic model. The same kinetic expression was used in both the CFD and CRE reactor models. Steady state conditions were assumed and the reaction was modelled at different temperature (400 ºC, 425 ºC and 450 ºC) and different feed compositions (dry feed and wet feed). The CRE 1D + 1D model was developed in Athena Visual Studio while the CFD model, with the catalyst meshing, diffusion and reaction aspects based largely on that reported by Partopour & Dixon [8] was compiled in OpenFOAM. For the most part this made use of existing verified folders available within OpenFOAM (e.g. reactingFoam solver) although additional coding for diffusion models and complicated kinetic models was required. 3. Results and discussion Figure.1 shows the comparison between the two approaches for the methane oxidation reaction conducted at 400 ºC and dry feed conditions.