Variational Iteration Method for a Fractional-Order Brusselator System

and Applied Analysis 3


Introduction
In recent years, fractional differential equations (FDEs) have been the focus of many studies due to their appearance in various fields such as physics, chemistry, and engineering [1][2][3].On the other hand, much attention has been paid to the solutions of fractional differential equations.Since most fractional differential equations do not have exact analytic solutions and approximate and numerical techniques, therefore, they are used extensively.Recently, the Adomian decomposition method, homotopy perturbation method, homotopy analysis method, and differential transform method have been used for solving a wide range of problems [4][5][6][7][8][9][10].
Another powerful analytical method, called the variational iteration method (VIM), was first introduced in [11].This technique has successfully been applied to many situations: for example, see [12][13][14][15][16][17].Reference [18] was the first where the variational iteration method was applied to fractional differential equations.Odibat and Momani [19] implemented the variational iteration method to solve partial differential equations of fractional order.
In this paper, we introduce a new application of the variational iteration method to provide approximate solutions of the fractional-order Brusselator system in the following form: subject to the initial conditions with  > 0,  > 0, 0 <   ≤ 1 ( = 1, 2), and  1 ,  2 are constants.
is used to represent the Caputo-type fractional derivative of order   .
The Riemann-Liouville definition of the fractional integration [2] is given by Abstract and Applied Analysis For our purpose in this paper, we adopt Caputo's fractional derivative [2]: where  is a positive integer and Γ(⋅) is the Gamma function.
In particular, 0 <   < 1, and we have The fractional-order Brusselator system has been considered by several authors recently [20][21][22].Gafiychuk and Datsko investigated its stability [20].Wang and Li proved by numerical method that the solutions of the fractionalorder Brusselator system have a limit cycle [22].We used the variational iteration method to investigate the approximate solutions of the fractional-order Brusselator system.

Variational Iteration Method
The principles of the variational iteration method and its applicability for various kinds of differential equations are given in [23,24].In [18], it was shown that the variational iteration method is also valid for fractional differential equations.In this section, following the discussion presented in [18], we extend the application of the variational iteration method to solve the fractional Brusselator equation: According to the variational iteration method, we can construct the correction functional for (6) as where   , ( = 1, 2) are the general Lagrange multiplier, which can be identified optimally via variational theory [25,26].
To identify approximately Lagrange multiplier, some approximations must be made.The correction functional equation ( 7) can be approximately expressed as follows: where x and ỹ are considered as restricted variations, in which  x = ỹ  = 0. To find the optimal  1 and  2 , we proceed as follows: The stationary conditions can be obtained as follows: We substitute   () = −1, ( = 1, 2) into the functional equation (11) to obtain the following iteration formula: The initial approximations  0 () and  0 () can be freely chosen if they satisfy the initial conditions of the problem.Finally, we approximate the solutions () = lim  → 1   () and () = lim  → 1   () by the th terms   () and   ().

Illustrative Examples
For purposes of illustration of (VIM) for solving Brusselator equation, we present two examples.
Example 1.Consider the following fractional-order Brusselator system: with the initial conditions: According to the variational iteration method and (11), the iteration formula for ( 12) is given by By using the above variational iteration formula, if we start with the initial approximations  0 () = 1 and  0 () = 1, we can obtain directly the other components as and so on; in the same way the rest of the components of the iteration formula can be obtained.Figures 1(a) and 1(b) show comparison between the approximate solutions (() ≅  5 ()), ((() ≅  5 ()) of ( 12) obtained using VIM for the special case  1 =  2 = 0.98 and the numerical solutions for the special case  1 =  2 = 1, respectively.Figures 2(a) and 2(b) show the approximate solutions (() ≅  5 ()), ((() ≅  5 ()) of ( 12) using VIM for the special case  1 =  2 = 1 and the numerical solutions, respectively.

Conclusions
The variational iteration method is a powerful method which is able to handle linear/nonlinear fractional differential equations.The method has been applied to fractional-order Brusselator system in order to find its approximate solutions.
The results show that the applied method is suitable and inexpensive for obtaining the approximate solutions.