Approximate Closed-form Solutions for a System of Coupled Nonlinear Partial Differential Equations: Brusselator Model.

Transcription

1 ISSN print), online) International Journal of Nonlinear Science Vol.6) No.4,pp.- Approximate Closed-form Solutions for a System of Coupled Nonlinear Partial Differential Equations: Brusselator Model S.O. Ajadi, O.A. Adesina, O.A. Jejeniwa, Department of Mathematics, Obafemi Awolowo University,Ile-Ife, 5, Nigeria. Received 5 March, accepted 7 September ) Abstract: In this article, approximate analytical expressions for the solutions of a system of coupled nonlinear reaction-diffusion equations have been obtained using the the homotopy pertubation methodhpm) framework. The effectiveness of this method is elucidated by applying this procedure to the reaction diffusion Brusselator model in one and two space dimensions. For the cases considered, the study reveals that the HPM reduces computational work and converges rapidly to its closed-form solutions. Graphical demonstration of these solutions shed more lights on the behavior of the system. Keywords: HPM; Brusselator model; closed-form solutions; non-linear partial differential equation. Introduction The nonlinear partial differential equations arise in many areas of application in science and technology particularly in mechanics, chemical kinetics, and biological phenomena. The study of non-linear reaction-diffusion systems has been an active area of research and continues to be so [,,4-6, -5,]. A typical example is the well known Brusselator model, which was developed in the 97s in Brussels by Nobel laureate Ilya Prigogine among others. It is one of the simplest chemical models exhibiting a pattern forming instability called Turing instability. This class of reaction diffusion system includes some significant pattern formation equations arising from the modeling of kinetics of chemical or biochemical reactions and from the biological pattern formation theory. In other words, the Brusselator is a trimolecular chemical reaction model used to show cooperative phenomenon in chemically reacting systems. At the same time, it is one of the simplest nonlinear) models of chemical systems for which the relative concentration of the constituents can vary either in space, or in time, or can even exhibit nonlinear traveling waves. Thus, the hypothetical reaction scheme described as the Brusselator model is a system whose chemical reactions scheme are given by[,4,7,9,-] A U, R = k B + U V + D, R = k U + V U, R = k U E, R 4 = k 4 ) where U and V are spatially and temporally varying active centresradical) concentrations, whereas the concentrations of the chemicals A and B are input chemicals, D and E are output chemicals, while R i k i ) are the reaction rates for each stepi). The first reaction is known as the radical initiation stepis), the second and third are known as the radical propagation stepsps), while the fourth one is regarded as the radical termination stepts). Using the law of mass action [8,9 and ], one can derive the mathematical equations corresponding to the reaction scheme defined by Eq. ). Arising from this, we consider a general class of coupled nonlinear partial differential also known as diffusion reaction system) equations of the form U = d U + F U, V ), ) Corresponding author. address: sajadi@oauife.edu.ng V = d V + GU, V ), ) Copyright c World Academic Press, World Academic Union IJNS.../77

2 S.O. Ajadi, O.A. Adesina et.al: Approximate Closed-form Solutions for a System of Coupled Nonlinear where d, d are the diffusion coefficients, F U, V ) and GU, V ) are the reaction terms. Thus, the reaction terms in the Brusselator model in Eq. ) are given as { F U, V ) = k A k BU + k U V k 4 U, GU, V ) = k BU k U 4) V By using the variables: U u, V v, A Ak /k, B Bk /k, k 4 k, t tk and after dropping primes, give rise to { F u, v) = A B + )u + u v Gu, v) = Bu u 5) v This reduces to a system of two nonlinearly coupled reaction diffusion equations given by[],[] u = d u + A B + )u + u v, t, r), ) Ω, r = x, y, z) 6) v = d v + Bu u v, t, r), ) Ω 7) ICs : u, r) = u r), v, r) = v r) 8) BCs : ut, r ) = u, ut, r ) = u, vt, r ) = v, vt, r ) = v, t >, Ω. 9) The Brusselator model is still attracting more attentions from researchers because of its theoretical relevance[9] and relationship to many physical phenomenon Kuptsov et al.[], Karafyllis et al.[], Chowdhury et al.[]. Manaa et al.[] solved the Brusselator model numerically in one and two dimensions by using two finite differences methods. For the one-dimension, they showed that the Crank-Nicolson method is more accurate, though more rigorous than the explicit method, while for the two dimensions, the Alternating Direction ImplicitADI) is more accurate than the Alternating Direction ExplicitADE) method. Rauber and Raunger[] considered the spatial discretization of nonlinear partial differential equations PDEs) results in large systems of nonlinear ordinary differential equations ODEs) using the Brusselator equation as a characteristic example. For the parallel numerical solution of the Brusselator equation, they used an iterated Runge-Kutta method. A theoretical analysis of the resulting speed up values shows that the efficiency cannot be improved considerably. Considering that another source of stiffness is the Brusselator system through the translation of diffusion terms by divided differences into a large system of ordinary differential equations in one spatial variable, Hairer and Warner[] obtain numerical solutions by using the finite difference method. More recently, chowdhury et al[] have obtained a term HPM solutions for the reaction diffusion Brusselator model in two dimension. However, they have applied the HPM procedure to an example which does not have known Exact or closed-form solutions, which makes it impossible to measure the relative accuracy of the HPM method. The motivation for this work is to obtain more accurate and efficient method for solving system of PDEs, which has long been an active research undertaking. The effectiveness of this method is implemented to solve some coupled system of nonlinear Brusselator equations and compare with existing closed-form or other approximate methods[manaa et al.[],rauber and Runger []. It was observed that the 5 term HPM solutions obtained are very close to the closed-form solutions, where they exist. The paper has been organized as follows. To start with, the notations and basic definitions of the homotopy pertubation method has been introduced, then the effectiveness of this method has been carried out by applying to some systems of non-linear reaction diffusion equations. Finally the accuracy of this method has been elucidated by comparing with established solutions and displayed graphically and in tabular forms. Homotopy perturbation method Traditionally, solving a differential or system of differential equations simply means writing down formulae which satisfy the equations. The best form of solution is the closed-formexact) solution or formula, which is a set of algebraic functions under the operation of integration. Due to the difficulty of obtaining closed-form solutions particularly for non-linear problems, less approximate solutions in the form of power series can often be found. Unfortunately, series constructed by such method does not converge [Rapp et al. [7]]. The homotopy perturbation methodhpm) is a novel and effective method for solving a wide range of problems whose mathematical models yield differential equation or system of differential equations. The solution obtained by this method is also expressed in series form and is evidently convergent[he[8 9]]. The HPM deforms a difficult problem into an infinite set of problems which are easier to solve without compromising IJNS homepage:

3 International Journal of Nonlinear Science, Vol.6), No.4, pp. - the nonlinearity. In line with Ajadi[] and Chowdhury et al. [], we illustrate homotopy perturbation method for the nonlinear equation Au) fr) =, r Ω ) with the boundary conditions Bu, u ) =, r Ω) ) n where A is a general differential operator, B is a boundary operator, fr) is a known analytic function, and Ω is the boundary of the domain Ω. The operator A are generally divided into two parts; L and N, where L and N are linear and nonlinear parts of A respectively. Therefore, Eqn. ) may be written as We construct a homotopy vr, p) : Ω [, ] R or Lu) + Nu) fr) =. ) Hv, p) = [Lv) Lu )] + plu ) + p[nv) fr)] =, ) Hv, p) = [Lv) Lu )] + p[au) fr)] =, 4) where p [, ] is called the homotopy parameter and u is an initial approximation of ). At the two extremes p = and p =, we have Hv, ) = Lv) Lu ) = and Hv, ) = Au) fr) =. 5) In the interval < p <, then the homotopy Hv, p) deforms from Lv) Lu ) to Au) fr). Thus, the solution of Eqns. ) and ) may be expressed as v = v + pv + p v + p v + 6) Eventually, at p =, the system takes the original form of the equation and the final stage of deformation gives the desired solution. Thus taking limit u = lim v = v + v + v + 7) p Homotopy perturbation method solutions In order to solve Eqns. 6)-9) and in line with Ajadi[], we carefully choose an initial approximation and construct the homotopy[partial t-solution] as u u [ = p d u + A B + )u + u v u ], 8) Suppose that the solution of ur, t) and vr, t) are expressed in the form v v [ = p d v + Bu u v v ]. 9) { ur, t) = u r) + pu r, t) + p u r, t) + p u r, t) + p 4 u 4 r, t) + p 5 u 5 r, t), vr, t) = v r) + pv r, t) + p v r, t) + p v r, t) + p 4 v 4 r, t) + p 5 v 5 r, t), where r = x, y, z, ). Substituting Eqn.) into Eqns.8) and 9) and picking terms in p, we have p : u = d u + u v B + )u + A u p : u = d u + u u v + u v ) B + )u p : u = d u + u v + u u v + u v + u u v ) B + )u p 4 : u 4 = d u + u u v + u u v + u v + u u v + u u v ) B + )u p 5 : u 5 = d u 4 + Hu, u, u, u, u 4, u 5, v, v, v, v, v 4, v 5 ) B + )u 4 ) ) IJNS for contribution: editor@nonlinearscience.org.uk

4 S.O. Ajadi, O.A. Adesina et.al: Approximate Closed-form Solutions for a System of Coupled Nonlinear and p : v = d v u v + Bu + A v p : v = d v u u v + u v ) + Bu p : v = d v u v + u u v + u v + u u v ) + Bu p 4 : v4 = d v u u v + u u v + u v + u u v + u u v ) + Bu p 5 : v 5 = d v 4 Hu, u, u, u, u 4, u 5, v, v, v, v, v 4, v 5 ) + Bu 4 where H = u v + u u v + u u 4 v + u u v + u u v + u v + u u v + u u v.the appropriate initial and boundary conditions of the partial differential equations in Eqns.) and ) can be obtained by substituting Eqn. ) into Eqns. 8) and 9). Hence { v, t) = v, v i, t) =, u, t) = u, u i, t) = 4) v x, ) = v x), v i x, ) =, u x, ) = u x, ), u i x, ) =, i =, If we start with u r and v r and all the linear equations above can be solved, then the solutions of 8) and 9) can be obtained by setting p = in Eqn. ). Thus { ur, t) = u r) + u r, t) + u r, t) + u r, t) + u 4 r, t) + u 5 r, t) 5) vr, t) = v r) + v r, t) + v r, t) + v r, t) + v 4 r, t) + v 5 r, t) 4 Test examples In order to test the validity and accuracy of this method, it is important to apply this method to some specific examples in literature with known exact or numerical solutions. 4. One-dimensional brusselator model We start by applying the homotopy technique to an initial boundary value problem of the form subject to some initial boundary conditions: Case []: A =, B =, d = d = 5. ) u = d u x + A B + )u + u v, 6) v = d v x + Bu u v, 7) u, t) = u, t) =, v, t) = v, t) =, ux, ) = + sin πx, vx, ) =. 8) Case []: A =.6, B =., d = d = 4. In line with Eqns. 8) and 9), the homotopy is expressed as u u [ u = p d u, t) = u, t) = A, v, t) = v, t) = B A, ux, ) = A + x x), vx, ) = B A + x x). 9) v v = p x + A B + )u + u v u [ v d x + Bu u v v From Eqn. 6), the solutions of ux, t) and vx, t) are expressed in the form ], ) ]. ) ux, t) = u x) + pu x, t) + p u x, t) + p u x, t) + p 4 u 4 x, t) + p 5 u 5 x, t), vx, t) = v x) + pv x, t) + p v x, t) + p v x, t) + p 4 v 4 x, t) + p 5 v 5 x, t). ) IJNS homepage:

5 4 International Journal of Nonlinear Science, Vol.6), No.4, pp. - Substituting Eqn. ) into Eqns. 8) and 9), and the conditions u i = and v i = i ) are determined by substituting Eqn. ) into Eqns. 6) and 7) for cases I and II respectively. The analytical solutions can easily be obtained by successive integration. After each integration, we have ensured that the boundary conditions satisfy the solutions by an interactive procedure implemented on the Wolfram Mathematica 5. codeappendix A). Using these schemes, the expressions for u, u, u, u 4, u 5... and v, v, v, v 4, v 5... are easily deduced. For simplicity and space economy, we have decided to present only the expressions for u, u,v and v as follows Case u = [5 B 4d π ) sin πx + sin πx]t, u = 6 [cos 8πx ) + B + ) cos 6πx + 8B d π ) cos 4πx) + 66B 88 8B + d π 64Bd π 8d π 4 ) sin πx]t, u = u 4 = u 5 = v = [B 6) sin πx sin πx]t, v = 6 [cos 8πx ) + B + ) cos 6πx + 8B + 96d π + 9d π ) cos 4πx) +74B 8 8B + 9d π + Bd π + 9d π Bd π ) sin πx]t, v = v 4 = v 5 = ) Case 4 u = t A [A xx ) A x x ) Bx x ) + Axx )x 4 x + x + B )], u = t A [A 6 xx ) A 5 x x ) B x x ) + A 4 xx ) +ABd xx ) + x x ) x 4 x + x + B )) +A 4d xx ) + x x ) x 4 4x + x + B 4)) A d 5x 4 x + 8x x) + xx )x 8 4x 7 + 6x 6 4x 5 + 6B )x 4 +6 B)x + B )x + B ) ))], u = u 4 = u 5 = v = t A [A xx ) A + B)x x ) + Axx )x 4 x + x + B)], v = t A [A 6 xx ) A 5 x x ) B x x ) +ABd xx ) + x x ) x 4 6x + x + B )) A B xx ) + B6x 6 8x 5 + 8x 4 6x x x) + xx )d xx ) +x x ) x 4 x + x ))) +A 4d xx ) + x x ) x 4 4x + x + B ))], v = v 4 = v 5 = 4) The final solutions are obtained by taking limit as p, ux, t) = u x) + u x, t) + u x, t) + u x, t) + u 4 x, t) + u 5 x, t), vx, t) = v x) + v x, t) + v x, t) + v x, t) + v 4 x, t) + v 5 x, t). 5) 4. Two-dimensional model The procedure has been applied to some known examples in two-dimensions in the form u ) = d u x + u y + A B + )u + u v, IJNS for contribution: editor@nonlinearscience.org.uk

6 S.O. Ajadi, O.A. Adesina et.al: Approximate Closed-form Solutions for a System of Coupled Nonlinear 5 subject to some initial boundary conditions: v ) = d v x + v y + Bu u v, 6) Case 5 []: A =, B =.4, d = d = ] u n = v n =, ux, y, ) = + y, vx, y, ) = + 5x. 7) Case 6 []: A =, B =, d = d = 4.] u n, v n) = x + y) exp t/ x y), expt/ + x + y)), ux, y, ) = exp x y), vx, y, ) = expx + y). 8) In two dimensions, we define the homotopy in line with Eqns. 8) and 9) as u u [ ) u = p d x + u y + A B + )u + u v u ], 9) v v [ ) = p v d x + v y + Bu u v v ]. 4) Suppose that the solutions of ux, y, t) and vx, y, t) are expressed in the form ux, y, t) = u x, y) + i= pi u i x, y, t) = u x, y) + pu x, y, t) + p u x, y, t) + vx, y, t) = v x, y) + i= pi v i x, y, t) = v x, y) + pv x, y, t) + p v x, y, t) +. Substituting Eqn. 4) into Eqns. 8) and 9), we have ) p : u = d u x + u y + u v B + )u + A u ) p : u = d u x + u y + u u v + u v ) B + )u ) p : u = d u x + u y + u v + u u v + u v + u u v ) B + )u p 4 : u 4 = d u x + u y ) + u u v + u u v + u v + u u v + u u v ) B + )u p 5 : u 5 = d u 4 x + u 4 y ) + Hu, u, u, u, u 4, u 5, v, v, v, v, v 4, v 5 ) B + )u 4 4) 4) and ) p : v = d v x + v y + u v Bu v ) p : v = d v x + v y + u u v + u v ) Bu ) p : v = d v x + v y + u v + u u v + u v + u u v ) Bu p 4 : v 4 = d v x + v y ) + u u v + u u v + u v + u u v + u u v ) Bu p 5 : v 5 = d v 4 x + v 4 y ) + Hu, u, u, u, u 4, u 5, v, v, v, v, v 4, v 5 ) Bu 4 4) Based on the initial conditionsu and v ) and parameters in cases I and II, the solutions of Eqns. 4) and 4) can easily be obtained by successive integration using the initial-boundary conditions in Eqns. 6) and 7) and implemented on the Wolfram Mathematica 5. Code[Appendix A]: IJNS homepage:

7 6 International Journal of Nonlinear Science, Vol.6), No.4, pp. - Case 7 u = [A + B) + y) + + 5x) + y) ]t, u =! { d + 5x) + + y) [B + y) + 5x) + y) ] + B)[A + B) + y) + + 5x) + y) ] + + 5x) + y)[a + B) + y) + + 5x) + y) ]}t, u = u 4 = u 5 = v = [B + y) + 5x) + y) ]t, v =! { d + 5x) + y) [B + y) + 5x) + y) ] +B[A + B) + y) + + 5x) + y) ] + 5x) + y)[a + B) + y) + + 5x) + y) ]}t, v = v 4 = v 5 = 44) Case 8 u = A B d ) exp x y)t, u =! {A B) + [d B + d + B d B + d ] exp x y) +B ) exp x y) + d d expx + y)}t, u = u 4 = u 5 = v = [B ) exp x y) + d expx + y)]t, v =! {AB ) + [B 4d B + d B 4d + d B] exp x y) + B) exp x y) + 4d expx + y)}t, v = v 4 = v 5 = The final solutions are obtained by taking limit as p, ux, y, t) = i= u ix, y) ti i! = u x, y) + u x, y)t + u x, y)t /! + u x, y)t /! + vx, y, t) = i= v ix, y) ti i! = v x, y) + v x, y)t + v x, y)t /! + v x, y)t /! + 45) 44) 5 Conclusion The choice of the non-negative initial datau and v ) and the positive constants A, B, d and d in literature[, 6] are for easy comparison. However, in applying the homotopy perturbation method, it is important to state that the rate of convergence of the homotopy perturbation method is greatly influenced by the choice of initial data[liao [7]]; which is the power of the method. In all the cases considered, we observed that the HPM method provides a better and more accurate results than other approximate methods. For case, Figs. 4 shed more lights on the behaviour of the system. Figs. and 4 6 show the variation and emergence of oscillatory pattern for ux, t) and vx, t) respectively with t for some fixed x, while Figs. 7 9 and show variation and emergence of oscillatoty pattern for ux, t) and vx, t) respectively with t for some fixed x. In Figs. - and Figs. 4 6, it can be observed that oscillatory pattern starts in the vicinity of t = and becomes more pronounced as t increases for ux, t) and vx, t) respectively. The D plots in Figs. 4 has not only confirm known results [], they also confirm that the oscillatory behaviour of the system is noticed at t >. In case, Figs. 5 6 are the ux, t) and vx, t) profiles for t =. The D plots for case are shown in Figs. 7. It can be observed that as t increases from t = to t =., the frequency of oscillations of vx, t) is more feasible than that of ux, t). Case 4 is an example that provides us with the basis of comparing HPM solutions with Exact solutions []. Figs. and Figs. 4 are representations of ux, t) and vx, t) profiles respectively. All these examples reveal that HPM is highly comparable with exact methods and a better improvement above other approximate methods. IJNS for contribution: editor@nonlinearscience.org.uk

8 S.O. Ajadi, O.A. Adesina et.al: Approximate Closed-form Solutions for a System of Coupled Nonlinear 7 U.5 t.5 t Figure : Plot of ux, t) vs t for x =. Figure : Plot of ux, t) vs t for x = U V x,t.5 t.5 t Figure : Plot of ux, t) vs t for x =.. Figure 4: Plot of vx, t) vs t for x =.. V x,t V x,t.5 t.5 t Figure 5: Plot of vx, t) vs t for x =. Figure 6: Plot of vx, t) vs t for x = x x Figure 7: Plot of ux, t) vs x for t =.. Figure 8: Plot of ux, t) vs x for t =. IJNS homepage:

9 8 International Journal of Nonlinear Science, Vol.6), No.4, pp x x Figure 9: Plot of ux, t) vs x for t =.. Figure : Plot of vx, t) vs x for t = x x Figure : Plot of ux, t) vs x for t =. Figure : Plot of vx, t) vs x for t =.. U x,t - - t x 5 V x,t -5 5 t x Figure : D plot of ux, t) vs x and t. Figure 4: D plot of vx, t) vs x and t. IJNS for contribution: editor@nonlinearscience.org.uk

10 S.O. Ajadi, O.A. Adesina et.al: Approximate Closed-form Solutions for a System of Coupled Nonlinear 9.8 u x,t.7.6 v x,t.4 Figure 5: D plot of ux, t) vs x and t. Figure 6: D plot of vx, t) vs x and t. 6 u x,t 4 x y 6 v x,t 4 x y Figure 7: D plot of ux, t) for t = Figure 8: D plot of vx, t) for t =. 6 u x,t 4 x y 6 v x,t 4 x y Figure 9: D plot of ux, t) for t =.. Figure : D plot of vx, t) for t =.. IJNS homepage:

11 International Journal of Nonlinear Science, Vol.6), No.4, pp. -.8 u x,y u x,y Figure : D plot of ux, t) for t = [HPM]. Figure : D plot of ux, t) for t = [EXACT]. 8 v x,y v x,y Figure : D plot of vx, t) for t = [HPM]. Figure 4: D plot of vx, t) for t = [EXACT]. References [] S. O. Ajadi and M. Zuilino, Approximate analytical solutions of reaction-diffusion equations with exponential source term: Homotopy perturbation methodhpm), Applied Mathematics Letters 4), [] S. O. Ajadi and P. Vegulla, Instabilities In A Reaction Diffusion Model: Spartially Homogeneous And Distributed System, Applied Mathematics E-Notes, ),6-46. [] E. Hairer and G. Wanner996): Solving Ordinary Differential Equation II. H. Stiff and Differential-Algebraic Problems. Springer Series in Comput. Math., vol. 4, nd edition. [4] S. L. Hollis, R. H. Martin, M. Pierre, Global existence and boundedness in reaction-diffusion systems, SIAM J, Math. Anal, 8),987), [5] D. Lesnic, The decomposition method for Cauchy reaction-diffusion problems, Applied Mathematics Letters 7),4-48. [6] M. Ganjiani and H. Ganjiani, Solution of coupled system of non-linear differential equations using homotopy analysis method, Nonlinear Dynamics9)56), [7] S. Liao, Comparison between homotopy analysis method and homotopy perturbation method, Applied Mathematics and Computation 695) [8] J. H. He, Homotopy perturbation method: A new nonlinear analytical technique, Applied Mathematics and Computation 5),7-79. [9] J. H. He Homotopy perturbation method for solving boundary value problems, Physics Letters A 56), [] A. Yildirim, Application of He s homotopy perturbation method for solving the Cauchy reaction-diffusion problem, Comput Math. Appl.579),6-68. [] A. Yildirim and H. Kocak, Homotopy perturbation method for solving space-time fractional advectiondispersion equation, Advances in Water Resources, 9),7-76. [] T. Rauber and G. Runger, Aspects of a distributed solution of the Brusselator equation, Proc. of the first Aizu IJNS for contribution: editor@nonlinearscience.org.uk

12 S.O. Ajadi, O.A. Adesina et.al: Approximate Closed-form Solutions for a System of Coupled Nonlinear International Symposium on Parallel Algorithms and Architecture Synthesis, 995),4-. [] Yuncheng You, Global Dynamics of the Brusselator Equations, Dynamics of PDE, Vol.4, No., 67-96, 7. [4] A. Yildirim and Y. Gulkanat, Analytical approach to fractional Zakharov-Kuzuetsov equations by He s homotopy perturbation method, Communication in Theoretical Physics, 5),5-. [5] S. Momani and A. Yildirim, Analytical approximate solutions of tne fractional convection-diffusion equation with nonlinear source term by He s homotopy perturbation method, International Journal of Computer Mathematics, 57), [6] Y. Moriakawa, T. Yamaguchi and T. Amemiya, Analytical Approach for Oscillation Properties of Soft Materials, Forma, 5, 4956,. [7] P. E. Rapp, T. I. Schmah and A. I. Mees, Models of knowing and the investigation of dynamical system, Physica D 999), -49. [8] G. C. Wake, J. G. Graham-Eagle and B. F. Gray, Oscillating chemical reactions: the well stirred and spatially distributed cases, Lect. Appl. Math., 4986), -55. [9] M. P. McDowell and J. M. Powers, Mathematical Modeling of the Brusselator, Department of Aerospace and Mechanical Engineering, University of Notre Dame, Notre Dame, Indiana8) AME [] M. S. H. Chowdhury, T. H. Hassan and S. Mawa, A New Application of Homotopy Perturbation Method to the reaction-diffusion Brusselator Model, Procedia Social and Behavioral Sciences 8), [] S. A. Manna, R. K. Saeed and F. H. Easif, Numerical Solution of Brusselator Model by finite Difference Method, Journal of Applied Sciences Research, 6) ), [] I. Karafyllis, P. D. Christofides and P. Daoutidis, Dynamics of a reaction-diffusion system with Brusselator kinetics under feedback control, PHYSICAL REVIEW E 59)999). [] P. V. Kuptsov, S. P. Kuznetsov and E. Mosekilde, Particle in the Brusselator model with flow, Physica D 6), Appendix A The Wolfram Mathematica program below provides the platform for obtaining the analytical solutions of the differential equations 4) and 5) subject to the initial-boundary conditions 6) and 7). By using the evaluated expression, then using the plot and plotd graphic options, it is easy to obtain the and dimensional plots. A := A; B := B; d := d ; d := d ; x := x; t := t; y := y; u := u x); v x); u = Evaluate[d D[u, x, x] + u v B + ) u + A] v = Evaluate[d D[v, x, x] u v + B u ] u = Evaluate[d D[u, x, x] + u u v + u v B + ) u ) B + ) u ] v = Evaluate[d D[v, x, x] u u v + u v ) + B u ] u = Evaluate[d D[u, x, x] + u v + u u v + u v + 4 u u v ) B + ) u ] v = Evaluate[d D[v, x, x] u v + u u v + u v + 4 u u v ) B u ] u 4 = Evaluate[d D[u, x, x] + 6 u u v + v u v + 6 u + v +6 u u v + 6 u u v ) B + ) u ] v 4 = Evaluate[d D[v, x, x] 6 u u v + u u v + 6 u + v +6 u u v + 6 u u v ) B u ] u 5 = Evaluate[d D[u 4, x, x] + 6 u v + 8 u u v + u u 4 v +4 u u v + 8 u u v + u v + u u v + 8 u u v ) B + ) u 4 ] v 5 = Evaluate[d D[v 4, x, x] + 6 u v + 8 u u v + u u 4 v +4 u u v + 8 u u v + u v + u u v + 8 u u v ) B u 4 ] Appendix B The Wolfram Mathematica program below provides the platform for obtaining the analytical solutions of the differential equations 4) and 5) subject to the initial-boundary conditions 6) and 7). Taking advantage of the evaluated IJNS homepage:

13 International Journal of Nonlinear Science, Vol.6), No.4, pp. - expression and the plotd graphic options on Mathematica, we obtain the and dimensonal plots. A := A; B := B; d := d ; d := d ; x := x; t := t; y := y; u := u x); v x); u = Evaluate[d D[u, x, x] + d D[u, y, y] + u v B + ) u + A] v = Evaluate[d D[v, x, x] + d D[v, y, y] u v + B u ] u = Evaluate[d D[u, x, x] + d D[v, y, y] + u u v + u v B + ) u ) B + ) u ] v = Evaluate[d D[v, x, x] + d D[v, y, y] u u v + u v ) + B u ] u = Evaluate[d D[u, x, x] + d D[u, y, y] + u v + u u v + u v + 4 u u v ) B + ) u ] v = Evaluate[d D[v, x, x] + d D[v, y, y] u v + u u v + u v + 4 u u v ) B u ] u 4 = Evaluate[d D[u, x, x] + d D[u, y, y] + 6 u u v + v u v + 6 u + v +6 u u v + 6 u u v ) B + ) u ] v 4 = Evaluate[d D[v, x, x] + d D[v, y, y] 6 u u v + u u v + 6 u + v +6 u u v + 6 u u v ) B u ] u 5 = Evaluate[d D[u 4, x, x] + d D[u 4, y, y] + 6 u v + 8 u u v + u u 4 v +4 u u v + 8 u u v + u v + u u v + 8 u u v ) B + ) u 4 ] v 5 = Evaluate[d D[v 4, x, x] + d D[v 4, y, y] + 6 u v + 8 u u v + u u 4 v +4 u u v + 8 u u v + u v + u u v + 8 u u v ) B u 4 ] IJNS for contribution: editor@nonlinearscience.org.uk