Approximate Solutions for Nonlinear Initial Value Problems Using the Modified Variational Iteration Method

We have used themodified variational iterationmethod MVIM to find the approximate solutions for some nonlinear initial value problems in the mathematical physics, via the Burgers-Fisher equation, the Kuramoto-Sivashinsky equation, the coupled Schrodinger-KdV equations, and the long-short wave resonance equations together with initial conditions. The results of these problems reveal that the modified variational iteration method is very powerful, effective, convenient, and quite accurate to systems of nonlinear equations. It is predicted that this method can be found widely applicable in engineering and physics.


Introduction
Nonlinear partial differential equations are known to describe a wide variety of phenomena not only in physics, where applications extend over magnetofluid dynamics, water surface gravity waves, electromagnetic radiation reactions, and ion acoustic waves in plasma, but also in biology, chemistry, and several other fields. It is one of the important tasks in the study of the nonlinear partial differential equations to seek exact and explicit solutions. In the past several decades both mathematicians and physicists have made many attempts in this direction. Various methods for obtaining exact solutions to nonlinear partial differential equations have been proposed. Among these methods are the Bäcklund transformation method 1, 2 , the Hirota's bilinear method 3 , the inverse scattering transform method 4 , extended tanh method 5-7 , the Adomian-Pade approximation 8-10 , the variational method 11-14 , the variational iteration method 15, 16 , the various Lindstedt-Poincare methods 17-20 , the Adomian decomposition method 8, 21, 22 , the F-expansion method 2 Journal of Applied Mathematics 23, 24 , the Exp-function method 25-27 and others 28-35 . Zayed et al. 36 investigated the travelling wave solutions for nonlinear initial value problems using the homotopy perturbation method. The modified variational iteration method is the couples of the variational iteration method with the homotopy pertirbation method. Recently Akbarzade and Langari 37 and Mei and Zhang 38 had used the modified variational iteration method for some nonlinear partial differential equations.
The main objective of the present paper is to use the modified variational iteration method MVIM for constructing the traveling wave solutions of the following nonlinear partial differential equations in mathematical physics: i the nonlinear Burgers-Fisher equation 39 : ii the nonlinear Kuramoto-Sivashinsky equation 40 : u t auu x bu xx cu xxxx 0, 1.2 iii the nonlinear coupled Schrodinger KdV equations 41 : iv the nonlinear long-short wave resonance equations 42 : together with initial conditions, where a, b, c, and β are arbitrary constants while i √ −1. It is interesting to point out that 1.1 includes the convection term uu x and the dissipation term u xx . Equation 1.2 describes the fluctuations of the position of a flame front, the motion of a fluid going down a vertical wall, or a spatially uniform oscillating chemical reaction in a homogeneous medium. Equation 1.3 describe various processes in dusty plasma such as Langmuri, dust-acoustic wave and electromagnetic waves, while in 1.4 u is the envelope of the short wave and is a complex function, and v is the amplitude of the long wave which is a real function.

Basic Idea of He's Homotopy Perturbation Method
We illustrate the following nonlinear differential equation [43][44][45][46][47][48][49][50][51][52][53][54] : Journal of Applied Mathematics 3 with the boundary conditions: where A is a general differential operator, B is a boundary operator, f r is an analytic function, and Γ is the boundary of the domain Ω. Generally speaking, the operator A can be divided into two parts L and N, where L is linear but N is nonlinear. Therefore, 2.1 can be rewritten in the following form: By the homotopy technique, we construct a homotopy V r, p : Ω× 0, 1 → R which satisfies where p ∈ 0, 1 is an embedding parameter and u 0 is an initial approximation of 2.1 which satisfies the boundary conditions 2.2 . Obviously, from 2.4 and 2.5 , we have The changing process of p from zero to unity is just that of V r, p from u 0 r to u r . In topology, this is called the deformation but L V − L u 0 and A V − f r are called the homotopies. According to the homotopy perturbation method, we can first use the embedding parameter "p as a small parameter and assume that 2.4 or 2.5 can be written as a power series in "p as follows: Letting p → 1 in 2.7 , the approximate solution of 2.3 takes the following form: The combination of the perturbation method and the homotopy method is called the homotopy perturbation method which has eliminated the limitations of the traditional perturbation methods. On the other hand, this technique can have full advantage of the traditional perturbation techniques.
which is variational iteration algorithm I, and λ is a general Lagrange multipliers. The variational iteration method can be identified optimally via variational theory 6, 7 . The second term on the right-hand side in 3.2 is called the corrections, the subscript n denotes the nth order approximation, and u n is restricted variations. We can assume that the aforementioned correctional functionals are stationary i.e., δu n 1 0 , and then the Lagrange multipliers can be identified. Now we can start with the given initial approximation and by the previous iteration formulas we can obtain the approximate solutions. He 55 has used the fractional iteration method to obtain the approximate solutions for nonlinear fractional differential equations.

The Modified Variational Iteration Method
To convey the basic idea of the variational homotopy perturbation method 2, 3 , we consider the following general differential equation: where L is a linear differential operator, N is a nonlinear operator and f x, t is an inhomogeneous term. According to the variational iteration method 4-13 , we can construct a correct functional as follows: where λ is a Lagrange multipliers, which can be identified optimally via variational theory 6, 7 . The subscripts n denote the nth approximation, and u n is considered as a restricted variation. That is, δ u n 0 is called a correct functional. Now, we apply the homotopy perturbation method to 4.2 : Journal of Applied Mathematics 5 which is the variational iteration algarithm II and is formulated by the modified variational iteration method. The embedding parameter p ∈ 0, 1 can be considered as an expanding parameter 14-19 .
The homotopy perturbation method uses the homotopy parameter p as an expanding parameter 14-19 to obtain 4.4 If p → 1, then 4.4 becomes the approximate solution of the following form: A comparison of like powers of p gives solutions of various orders. The application of the Adomain polynomial is too complex so that we consider the variational iteration method and He's polynomial to calculate the approximate solutions see, e.g., 56-60 .

Applications
In this section, we construct the approximate solutions for some nonlinear evolution equations in the mathematical physics, namely, the Burgers-Fisher equation 1.1 , the Kuramoto-Sivashinsky equation 1.2 , the coupled Schrodinger-KdV equations 1.3 , and the long-short wave resonance equations 1.4 together with initial conditions by using the the modified variational iteration method. Applications of this method to similar equations can be found in 61-67 .

Approximate Solution of Burgers-Fisher Equation with Initial Conditions Using Modified Variational Iteration Method
In this subsection, we use the MVIM to find the solution u x, t of an initial value problem consisting of the nonlinear Burgers-Fisher equation 1.1 and the following initial condition 39 : This initial condition follows by setting t 0 in the following exact solution of 1.1 : This exact solution has been derived by Wazwaz 39 using the tanh-coth method. To this end, we construct the modified variational iteration method for the nonlinear Burgers-Fisher equation 1.1 which satisfies

5.3
Comparing the different coefficient of like power of p, we have and so on. Consequently after some reduction with help of Maple or Mathematica, we get:

5.5
In this manner the other components can be obtained. Substituting from 5.5 into 4.5 , we obtain the following approximate solution of the initial value problem 1.1 and 5.1 : 5.6  Note that if we expand the exact solution 5.2 in Taylor series near t 0, we obtain the approximate solution 5.6 . To demonstrate the convergence of the variational homotopy perturbation method, the results of the numerical example are presented and only few terms are required to obtain accurate solutions. The accuracy of the modified variational iteration method for the nonlinear Burgers-Fisher equation is controllable and absolute errors are very small with the present choice of x and t. These results are listed in Table 1. Both the exact solution 5.2 and the approximate solution 5.6 obtained for the first three approximations are plotted in Figure 1. There are no visible differences in diagrams. It is also evident that when more terms for the modified variational iteration method are computed, the numerical results get much more closer to the corresponding exact solution with the initial condition 5.1 .

Approximate Solution of the Nonlinear Kuramoto-Sivashinsky Equation with Initial Conditions Using MVIM
In this subsection, we use the MVIM to find the solution u x, t of an initial-value problem consisting of the nonlinear Kuramoto-Sivashinsky equation 1.2 with the following initial condition 40 : where λ −b/19c, b/c < 0, and a, b, and c are constants. This initial condition follows by setting t 0 in the following exact solution of 1.2 : This exact solution has been derived by Wazwaz Comparing the different coefficient of like power of p, we have

5.10
Journal of Applied Mathematics 9 and so on. Consequently after some reduction with help of Maple or Mathematica, we get

5.11
and so on. Substituting from 5.11 -5.12 into 4.5 , we obtain the approximate solution of the initial value problem 1.2 :

5.12
which is in agreement with the exact solution 5.8 using Taylor series expansion near t 0. The comparison between the exact solution 5.8 and the approximate solution 5.12 is shown in Table 2

Approximate Solutions for the Nonlinear-Coupled Schrodinger-KdV Equations with Initial Conditions Using MVIM
In this subsection, we find the solutions u x, t and v x, t satisfying the nonlinear coupled Schrodinger-KdV equations 1.3 with the following initial conditions 41 : where k, α, and c are arbitrary constants and α 2k 2 c/2. These initial conditions follow by setting t 0 in the following exact solutions of 1.3 : 14 v x, t −2k 2 sech 2 ik x − ct .

5.15
These exact solutions have been derived by Zhang 41 using a direct algebraic approach. Let us now apply the MVIM to the initial value problem 1.3 and 5.13 :

5.16
Journal of Applied Mathematics

11
Comparing the different coefficient of like power of p, we have p 0 : u x, 0 −cα 2αk tanh ikx ,

5.17
and so on. Consequently after some reduction with help of Maple or Mathematica, we get x In this manner the other components can be obtained. Substituting 5.18 into 4.5 , we obtain the approximate solutions of the initial value problem 1.3 and 5.13 :

5.18
which are in the closed form of the exact solutions 5.14 and 5.15 using Taylor series expansion near t 0.
The comparison between the exact solutions 5.14 , 5.15 and the approximate solutions 5.19 , 5.20 respectively, are shown in Table 3

Approximate Solution of the Nonlinear Long and Short Wave Resonance Equations with Initial Conditions Using MVIM
In this subsection, we find the solutions u x, t and v x, t satisfying the nonlinear long-short wave resonance equations 1.4 with the following initial conditions 42 : where β, k, C, ξ 0 and ξ 1 are arbitrary constants. These initial conditions follow by setting t 0 in the following exact solutions of 1.4 : where w C − k 2 is constant. These exact solutions have been derived by Shang 42 using the extended hyperbolic function method, which describes the resonance interaction between the long wave and the short wave. Let us now apply the MVIM to the initial value problem 1.4 and 5.21 :

5.26
In this manner the other components can be obtained. Consequently, we obtain the following approximate solutions of the initial value problem 1.4 and 5.  Table 4 and Figures 5 and 6. It seems that the errors are very small if t 0.5, k 0.1, c 0.2, ξ 0 2, and ξ 1 3.

Conclusions
In this paper, the modified variational iteration method was applied for finding the approximate solutions for some nonlinear evolution equations in mathematical physics via the nonlinear Burgers-Fisher equation, nonlinear Kuramoto-Sivashinsky equation, nonlinear coupled Schrodinger KdV equations, and nonlinear long-short wave resonance equations with well-known initial conditions. It seems to us that the modified variational iteration method presents a rapid convergence solutions. It can be concluded that this method is very powerful and efficient technique in finding approximate solutions for wide classes of nonlinear problems.