Abstract
The objective of this paper is to investigate the effectiveness and performance of optimal homotopy asymptotic method in solving a system of nonlinear partial differential equations. Since mathematical modeling of certain chemical reaction-diffusion experiments leads to Brusselator equations, it is worth demanding a new technique to solve such a system. We construct a new efficient recurrent relation to solve nonlinear Brusselator system of equations. It is observed that the method is easy to implement and quite valuable for handling nonlinear system of partial differential equations and yielding excellent results at minimum computational cost. Analytical solutions of Brusselator system are presented to demonstrate the viability and practical usefulness of the method. The results reveal that the method is explicit, effective, and easy to use.
1. Introduction
The Brusselator model, the nonlinear system of partial differential equations, arises in the modeling of certain chemical reaction-diffusion processes. This Brusselator reaction-diffusion model plays a substantial role in the study of cooperative processes of chemical kinetics. This system occurs in a large number of physical problems. It arises in the creation of ozone by atomic oxygen through a triple collision, in enzymatic reactions, and in plasma and laser physics in numerous coupling among models [1]. A pair of variables are involved in dealing with these chemical reactions; intermediates with input and output chemicals, whose concentrations are likely to be controlled during the reaction process, are substantial under quite genuine conditions and are discussed by Nicolis and Prigogine in [2, 3]. This model has been revealed as the trimolecular model [4].
The two-dimensional nonlinear reaction-diffusion Brusselator system issubject to the initial condition:where are unknown functions representing the dimensionless concentrations of two reactants, , and denote spatial and temporal independent variables, respectively, and are constant concentrations of the two reactants, represent the diffusion coefficient, and and are known functions. It is evident that, for small values of diffusion coefficient , steady state solution of Brusselator system converges to the equilibrium point if . During the last few years, the researchers have keen interest in the existence of solution of the Brusselator reaction model when [5–7]. In this paper, we have made a successful attempt to find the solution of such types of Brusselator system.
Numerical methods need large size of computational works and generally the consequence of round-off error causes loss of precision in the results for system of nonlinear partial differential equations. Analytical methods mostly used for solving these equations are very restricted and can be used in very special cases. Therefore, an optimal technique is required to resolve such circumstances.
The OHAM was recently introduced by Marinca and Herişanu [8–11]. In series of papers, authors [12–16] have applied this method effectively to validate the solutions of currently important problems. Iqbal et al. [12–15] used OHAM for linear and nonlinear differential equations and time dependent partial differential equations. Similarly, Idrees et al. [16] validated OHAM for nonlinear phenomena of fluid mechanics.
In this presentation, we have extended OHAM formulation for system of partial differential equations. Particularly, the extended formulation is demonstrated by illustrative example of nonlinear Brusselator system partial differential equations.
2. OHAM Formulation for Brusselator System
The optimal homotopy asymptotic method (OHAM) is given in [8–16]; we formulate this method for fractional order partial differential equations in the following steps.(a)Write the governing partial differential equation system as subject to the initial condition: where , are unknown functions; , and denote spatial and temporal independent variables, respectively.(b)Construct an optimal homotopy for system of partial differential equations, , which satisfies where is an embedding parameter, and are nonzero auxiliary functions for and , , and, clearly, we have where are convergence control parameters. The selection of functions might be different type of polynomial and so on. It is very important to choose the functions, since the convergence of the solution very much depends on these functions. The auxiliary function provides us with a simple way to adjust and control the convergence. It also increases the accuracy of the results and effectiveness of the method [17, 18]. The convergence control parameters adjust and control the convergence, which provides optimal accuracy and effectiveness of the method. The presence of ensures the fast convergence.(c)Expand and in Taylor’s series about , to get approximate solutions as(i)It has been observed that the convergence of the series (8) and (9) depends upon the convergence control parameters.(ii)If it converges at , one has(d)Equate the coefficients of like powers of after substituting (8) and (9) in (5) and (6), respectively; one can get zeroth-order system given by (12), first-order and second-order system given by (13)-(14), respectively, and the higher order system if needed: and so on.(e)The above system of nonlinear equations, that is, zeroth-order, first-order, and higher order systems (if needed), can be easily solved. Put these solutions of different order problems in (10) and (11); one can obtain the approximate solutions and , respectively.(f)Determine the convergence control parameters, , by using one of the methods given in [8–11]. Using auxiliary constants in (10) and (11), one can get the approximate solutions and , respectively.
3. Application
Example 1. Let us consider with , , , , and two-dimensional Brusselator system (1) can be written in the following form:subject to the initial conditions:The exact solutions of (15) is found to be [19]OHAM formulation presented in Section 2 generates the series of problems, which can be written asThe above zeroth-order, first-order, second-order, and third-order problems are given in (18) which can easily be solved. We getSubstituting (19) in (10) and (11), respectively, we obtainAuxiliary constants shown in (20) can be found by using step (f) in Section 2.
4. Results and Discussions
The formulation of OHAM for two-dimensional nonlinear Brusselator system with , , and is presented in Section 2 and the demonstration of the formulation is presented in Example 1.
The third-order approximate solutions of , are given in (19). These solutions depend upon the optimal convergence control parameters which are given in Table 1. The simplicity and accuracy of the presented method are illustrated by computing , . Tables 2 and 3 at different gird points show the comparisons of absolute error of and between OHAM and the exact solutions, respectively. In this study, only up to third-order solutions are considered. From Tables 2 and 3, it is clear that OHAM achieves accurate solutions at only third-order term of approximation without any spatial discretization. Thus the third-order approximate solutions of Brusselator reaction-diffusion system converge.
Figures 1, 2, and 3 show the convergence of first-order, second-order, and third-order approximate solutions of obtained by OHAM, respectively. Figures 4, 5, and 6 show the convergence of first-order, second-order, and third-order approximate solutions of obtained by OHAM, respectively. Figure 7 shows the error convergence with the order of approximation of and . It is clear from the figures that the behaviour of approximate solutions is highly same as exact solution. It can also be clear from Figures 3 and 6 that third-order approximate solution converges and there is no need to compute extra terms when OHAM is used. It is observed that as we move along the domain we get consistent accuracy. It can also be observed that approximate solutions by formulation are in excellent agreement with the exact solutions. Thus the series solutions for fractional equations converge. Results indicate the performance of the method nonlinear system of partial differential equations in precisely approximating the solution at less computational cost.
5. Conclusion
In this work, we employed a new powerful semianalytic technique, optimal homotopy asymptotic method to solve Brusselator reaction-diffusion system. The objective of this work is to illustrate the usefulness of the technique. This technique is simple in applicability, as it does not need discretization like numerical methods. Furthermore, this technique delivers an appropriate way to control the convergence by optimally shaping the convergence control parameters. Additionally, this method converges rapidly at lower order of approximations. Therefore, OHAM exhibits its concealed supremacy and is latent for the solution of nonlinear problems in real life applications. One substantial objective of this effort is the investigation of convergence and practicality of the method. By using this method, we acquire a new effective recurrent relation to solve nonlinear Brusselator system. The results demonstrate that the method is a prevailing mathematical tool for solving systems of nonlinear partial differential equations having extensive applications in science and engineering.
Competing Interests
The author declares that there is no conflict of interests regarding the publication of this paper.