Analytical expression of concentrations of adsorbed CO molecules , O atoms and oxide oxygen

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.


INTRODUCTION
Rate auto-oscillations in a heterogeneous catalytic reaction were first discussed more than two decades ago [1][2][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][5][6][7][8][9][10][11][12].Microscopic stochastic models by using the Monte Carlo method [13][14][15][16][17][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 mobil-ity 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 [3] described the oscillatory dynamical system which includes a Monte Carlo simulation and adsorption layer approximation.To our knowledge no analytical expressions of the concentration of CO molecules, O atom and oxide oxygen have been derived from the system.The purpose of the communication is to derive the analytical expression of , , and by solving the system of non-linear differential equation using Homotopy perturbation method.

MATHEMATICAL FORMULATION OF BOUNDARY VALUE PROBLEM
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 CO ads , an adsorbed oxygen atom O ads or an adsorbed oxide oxygen atom ads .The thickness shear mode (TSM) kinetic network includes the Langmuir-Hinshelwood mechanism [3].
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.

STOCHASTIC AND DETERMINISTIC MODELS
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 ads , ads and ads concentrations are described by a set of following nonlinear differential equations [4]: To simplify, we can assume that , and The above Eqs.9-11become Now the boundary condition ( 12) (or initial conditions) becomes

ANALYTICAL EXPRESSIONS OF CONCENTRATIONS USING HOMOTOPY PERTURBATION METHOD [HPM]
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][17][18][19].This method is a combination in topology and classic perturbation techniques.Ji Huan He used the HPM to solve the light hill equation [20], the Duffing equation [21] and the Blasius equation [22].The idea has been used to solve non-linear boundary value problems, integral equations and many other problems [23][24][25][26][27].The HPM is unique in its applicability, accuracy and efficiency.The HPM uses the imbedding parameter p as a small parameter and only a few iterations are needed to search for an asymptotic solution.Using this method (see Appendix A), we can obtain the following solution to Eqs.14-16 for the given boundary conditions (Eq.17).
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.

ANALYTICAL EXPRESSION OF CONCENTRATION FOR STEADY STATE CONDITION
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:

NUMERICAL SOLUTION
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 O 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 CO for all values of time and rate constant.

CONCLUSION
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 Copyright © 2013 SciRes.
Openly accessible at http://www.scirp.org/journal/ns/this method with more modelling and simulation procedure to the some of the non-linear reaction diffusion mechanism in biosensor [28] seems possible.
( ) ( ) ( ) ( ) 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.