Mathematical modeling of self-oscillations in ethane oxidation over nickel

Abstract

The methodology of constructing a phenomenological model for complex heterogeneous catalytic reactions is described in detail. The proposed approach is applicable to development of mathematical models describing the onset of self-oscillations in hydrocarbon oxidation on the transition metal surface. The approach is based on construction of a microkinetic scheme taking into account the formation of main reaction products and intermediates, on estimation of the heat of reaction, activation energy, and preexponential factor for elementary steps and includes development and a subsequent analysis of the corresponding mathematical model. Catalytic reactions are considered in the ideal adsorption layer approximation without taking into account the relationship between coverages and spatial coordinates. Accordingly, the mathematical model is an independent system of ordinary differential equations. This methodology is used to develop a point (lumped) model for ethane oxidation over nickel, which is based on a 36-step microkinetic scheme taking into account the oxidation and reduction of nickel and the formation of total (CO₂ and H₂O) or partial (CO and H₂) ethane oxidation products, as well as the dehydrogenation of ethane into ethylene. The proposed model predicts the onset of self-oscillations in this reaction at atmospheric pressure in the temperature range from 850 to 1400 K. The kinetic oscillations are caused by the cyclic oxidation and reduction of nickel. The self-oscillations of the reaction rate are accompanied by oscillations of the catalyst temperature. The results of modeling are compared with experimental data.