A mathematical model of the oscillatory regimes of CO oxidation over plantinum-group metal catalysts are discussed. The model is based on nonstationary diffusion equation containing a nonlinear term related to Michaelis-Menten kinetics of the enzymatic reaction. This paper presents the analytical and numerical solution of the system of non-linear differential equations. Here the Homotopy perturbation method (HPM) is used to find out the analytical expressions of the concentration of CO molecules, O atom and oxide oxygen respectively. A comparison of the analytical approximation and numerical simulation is also presented. A good agreement between theoretical and numerical results is observed.
Rate auto-oscillations in a heterogeneous catalytic reaction were first discussed more than two decades ago [1-3]. Some heterogeneous systems which show oscillatory behaviour on catalysts with various structures in a wide pressure range are presently known. One of the most extensively studied heterogeneous catalytic reactions exhibiting oscillatory dynamics is CO oxidation on platinum catalysts.
Various mathematical models are used in detailed analysis of the mechanism of the rate oscillations in catalytic CO oxidation. These models are based on a set of nonlinear ordinary differential equations [4-12]. Microscopic stochastic models by using the Monte Carlo method [13-18] one of the most interesting in theoretical investigation of the dynamics of fluctuating reaction systems. These stochastic models are based on detailed information concerning the elementary steps of the reaction, the structure of the catalyst surface, and the mobility of species in the adsorption layer.
Therefore, the theoretical investigation and explanation of complex dynamic phenomena on a catalyst surface cannot be comprehensive. It is necessary to develop systems of consistent mathematical models describing the evolution of reaction systems on different spatial scales. This would allow the advantages of different classes of mathematical models.
Recently Kurkina and Semendyaeva [
The oscillatory dynamics of the thickness shear model (TSM) kinetic network will be described using Lattice gas model. In this model any lattice site may be free (*) or occupied by an adsorbed carbon monoxide molecule COads, an adsorbed oxygen atom Oads or an adsorbed oxide oxygen atom. The thickness shear mode (TSM) kinetic network includes the Langmuir-Hinshelwood mechanism [
and the formation and removal of surface oxide,
Furthermore, we consider the migration of adsorbed species via a vacant-site mechanism:
Here, and are molecules in the gas phase. Two-site process occurs on adjacent lattice sites. The rate of a two-site process is defined for a pair of sites, and the rate of a one-site process is defined for one site. The surface oxide formation step (4) is viewed as a one-site process, and step (5) is a two-site process. It is assumed that oxide oxygen retains its reactivity (see step (5)) and markedly decreases the activity of the catalyst in reactant adsorption from the gas phase.
The basic kinetic equations in a multidimensional system will be solved by approximate methods. In Monte Carlo stochastic modelling, state trajectories are constructed for the reaction system in the state space. Stochastic models describe the evolution of a selected lattice fragment at the atomic level. According to the above kinetic network, the variation of, and concentrations are described by a set of following nonlinear differential equations [
The initial conditions are as follows:
Here, , , and are the concentrations of adsorbed CO molecules and O atoms and oxide oxygen, respectively, and are the rate constants of elementary steps. The normalization conditions are as follows:
To simplify, we can assume that, and The above Eqs.9-11 become
Now the boundary condition (12) (or initial conditions) becomes
where
Recently, many authors have applied the HPM to various problems by demonstrating the efficiency of the HPM in handling non-linear structures and solving various physics and engineering problems [16-19]. This method is a combination in topology and classic perturbation techniques. Ji Huan He used the HPM to solve the light hill equation [
where
The Eqs.18-20 represent the new closed form of analytical expression of concentrations of CO molecules, O atom and oxide oxygen for all values of rate constant and time.
For the case of steady state, the Eqs.14-16 become
Solving the above non-linear algebraic equation by using Maple software, we can obtain the concentration of CO molecule, O atom and oxide oxygen as follows:
In order to test accuracy of this method the non-linear differential Eqs.14-16 for the boundary conditions (Eq.17) are solved by numerical methods. The function pdex4 in Matlab software was used to solve these equations. It is a function of solving the initial boundary value problems for non-linear differential equations. The Matlab program is also given in Appendix B. The numerical results are also compared with our analytical results in Figures 1(a)-(d). A satisfactory agreement is noted here.
Figures 1(a)-(d) represent the non steady state concentration of CO molecules u(t) or, O atom v(t) or and oxide oxygen w(t) or for all values of rate constant and time. From the figure it is evident that the value of the concentration of CO molecule decrease from its initial value of concentration and reaches the steady state value zero when t = 2 sec for all values of rate constant. From the figures it is also inferred that the value of concentration of O atom rises initially and reaches the maximum value when t » 0.5 sec and then decreases gradually. The concentration of O atom attains the steady state value when t = 10 sec. From this figure it is also observed that the concentration of oxide oxygen always increases when time increases and reaches the steady state value 1 when t = 5 sec for all values of rate constant. From this figure it is to conclude that for all values of time and rate constant.
The theory of the dynamics of catalysts CO oxidation in the frame work of the Thickness shear mode kinetic is described. Approximate analytical solutions to the system of non-linear reaction equations are presented using Homotopy perturbation method. A simple, straight forward and a new method of estimating the concentrations of CO molecules, O atom and Oxide oxygen are derived. This analytical result will be useful to know the behaviour of the reaction system. A good agreement with the numerical simulation data is also noted. The extension of
this method with more modelling and simulation procedure to the some of the non-linear reaction diffusion mechanism in biosensor [
This work is supported by the University Grants Commission (UGC) (Ref. No. F. No. 39 - 58/2010(SR)) and Council of Scientific and Industrial Research (CSIR) (No.01 (2442)/10/EMR-II), Government of India. The authors are thankful to the Secretary, the Principal, The Madura College, Madurai, India for their constant encouragement.
In this appendix, we indicate how Eqs.14-16 may be solved using HPM. To find the solution of Eqs.14-16, we first construct a homotopy as follows [16-19]:
where is an embedding parameter. According to HPM, we can first use the parameter p as a “small parameter” and assume that the solution of Eqs.14-16 can be written as a power series in p.
where are zero-order solution (initial approximation) and are the first order approximate solution. When p = 0, we get zero-order solution (initial approximate solution or solution of linear terms). Setting p = 1, results in the approximate solution of Eqs.14-16. Substituting Eqs.A4-A6 into Eqs.A1-A3 and comparing the coefficients of like powers of p we obtain the following differential equations.
and
The initial conditions in Equation (17) becomes
Solving the equations using the initial conditions Equation (A14), we obtain the following results.
According to the HPM, we can conclude that
where
After putting the Eqs.A15-A20 in the Eqs.A21-A24 we get the Eqs.18-20 in the text.
function main1 options = odeset ('RelTol',1e-6,'Stats','on');
%initial conditions T=100;
Xo = [0.1; 0.5; 0.01];
tspan = [0,10];
tic
[t,X] = ode45(@TestFunction,tspan,Xo,options);
toc figure hold on plot(t, X(:,1))
plot(t, X(:,2))
plot(t, X(:,3),'.')
legend('x1','x2','x3')
ylabel('x')
xlabel('t')
return function [dx_dt]= TestFunction(t,x)
a=0.9,b=0.002,k2=0.5,k4=0.03,k5=0.2,k3=1,c=0.1,d=0.01,O=0.5dx_dt(1) =a*(1-x(1)-x(2)-x(3))-b*x(1)-4*k3*x(1)*x(2)-4*k5*x(1)*x(3);
dx_dt(2) =4*k2*(1-x(1)-x(2)-x(3))^2-4*k3*x(1)*x(2)-k4*x(2);
dx_dt(3)=k4*x(2)-4*k5*x(1)*x(3);
dx_dt = dx_dt';
return