Numerical simulation of packed-bed reactor for oxidative coupling of methane

Abstract

A three-dimensional geometric model of the oxidative coupling of methane (OCM) packed-bed reactor loaded with Na2WO4-Mn/SiO2 particulate catalyst was set up, and an improved Stansch kinetic model was established to calculate the OCM reactions using the computational fluid dynamics method and Fluent software. The simulation conditions were completely the same with the experimental conditions that the volume velocity of the reactant was 80 mL/min under standard state, the ratio of CH4/O2 was 3, the temperature and pressure were 800 °C and 1 atm, respectively. The contour of the characteristics parameters in the catalyst bed was analyzed, such as the species mass fractions, temperature, the heat flux on side wall surface, pressure, fluid density and velocity. The results showed that the calculated values matched well with the experimental values on the conversion of CH4 and the selectivity to products (C2H6, C2H4, CO2, CO) in the reactor outlet with an error range of ±2%. The mass fractions of CH4 and O2 decreased from 0.6 and 0.4 in the catalyst bed inlet to 0.436 and 0.142 in the outlet, where the mass fractions of C2H6, C2H4, CO and CO2 were 0.035, 0.061, 0.032 and 0.106, respectively. Due to the existence of laminar boundary layer, the contours of each component bent upwards in the vicinity of the boundary layer. This OCM reaction was volume increase reaction and the total moles of products were greater than those of reactants. The flow field in the catalyst bed maintained constant temperature and pressure. The fluid density decreased gradually from 2.28 kg/m3 in the inlet of the catalyst bed to 2.22 kg/m3 in the outlet of the catalyst bed, while the velocity increased from 0.108 m/s to 0.115 m/s.