A Novel Effective Approach for Solving Fractional Nonlinear PDEs

The present work introduces an effective modification of homotopy perturbation method for the solution of nonlinear time-fractional biological population model and a system of three nonlinear time-fractional partial differential equations. In this approach, the solution is considered a series expansion that converges to the nonlinear problem. The new approximate analytical procedure depends only on two iteratives. The analytical approximations to the solution are reliable and confirm the ability of the new homotopy perturbation method as an easy device for computing the solution of nonlinear equations.


Introduction
In recent years, fractional differential equations have played an important role in different research areas such as mechanics, electricity, biology, economics, notably control theory, and signal and image processing [1][2][3][4][5]. Many phenomena can be described very successfully by models using mathematical tools from fractional calculus, that is, the theory of derivatives and integrals of fractional order. The most important advantage of using fractional differential equations in these and other applications is their nonlocal property. It is well known that the integer order differential operator is a local operator, but the fractional-order differential operator is nonlocal. This means that the next state of a system depends not only upon its current state but also upon all of its historical states. In fact, this is the main reason why differential operators of fractional order provide an excellent instrument for description of memory and hereditary properties of various mathematical, physical, and engineering processes. The reader is asked to refer to [6][7][8][9] in order to know more details about the fractional DEs, including their history and kinds, their existence and uniqueness of solutions, and their applications and methods of solutions.
During the last few years, the numerical methods and exact solution methods have been proposed to solve fractional differential equations, for example, the Adomian decomposition method [10], the homotopy perturbation method [11,12], the variational iteration method [13,14], the differential transform method [15,16], the / method [17,18], the first integral method [19], and the exp-function method [20]. Aminikhah and Hemmatnezhad [21] and Aminikhah and Biazar [22] proposed a new homotopy perturbation method to obtain the exact and numerical solutions of ordinary differential equation.
In this paper, we will implement the new homotopy perturbation method to obtain the approximate solution for the following time-fractional derivative nonlinear partial differential equation.
(i) The fractional-order biological population equation as a nonlinear model is as follows [23,24]: 2 International Scholarly Research Notices real numbers, is given initial condition, and denotes the differential operator in the sense of Caputo.
(ii) The system of three nonlinear time-fractional partial equations which have been widely discussed in the literature [25,26] is as follows: with initial conditions ( , , 0) = ℎ 1 ( , ) , where denotes the differential operator in the sense of Caputo.
The rest of this paper is organized as follows. In Section 2, the properties of fractional derivative are described. In Section 3, basic ideas of new homotopy perturbation method for solving fractional partial differential equations are presented. In Section 4, the application of new homotopy perturbation method for solving the nonlinear fractional biological population equation and the system of three nonlinear timefractional partial equations are presented. Finally, we give a brief conclusion in the last section.

Background on Fractional Derivatives
The history of fractional calculus is more than three centuries old, yet only in the past 20 years has the field received much attention and interest. The reader may refer to [8,9]. There are several definitions for fractional differential equations. These definitions include Grunwald-Letnikov, Riemann-Liouville, Caputo, Weyl, Marchaud, and Riesz fractional derivatives; among those Riemann-Liouville and Caputo fractional derivatives are the most popular. The differential equations defined in terms of Riemann-Liouville derivatives require fractional initial conditions whereas the differential equations defined in terms of Caputo derivatives require boundary conditions involving integer order derivatives, which have clear physical interpretations. For this reason, Caputo fractional derivatives are popular among scientists and engineers. Now, we give some basic definitions and properties of the fractional calculus theory used in this work.

Remark 3.
For to be the smallest integer that exceeds , the Caputo time-fractional derivative operator of > 0 is defined as

Analysis of New Homotopy Perturbation Method
In this section, we construct the solution of system of partial differential equations with time-fractional derivative by extending the idea of [21,22]. Let us consider the system of nonlinear fractional differential equations: with the following initial conditions: where are the operators, are known functions, and are sought functions. Assume that operators can be written as ( ) = ( )+ ( ), where are the linear operators and are the nonlinear operators. Hence, (7) can be rewritten as follows: International Scholarly Research Notices 3 For solving system (7) by NHPM [22] we construct the following homotopy: where ∈ [0, 1] is an embedding or homotopy parameter, ( , x; ) : Ω × [0, 1] → , and ,0 are the initial approximation of solution of the problem in (9). Clearly Now, let us write (11) in the following form: Applying the inverse operator which is the Riemann-Liouville fractional integral of order ≥ 0, on both sides of (12), we have Suppose that the initial approximation of (9) has the form where , (x), = 0, 1, 2, . . ., are unknown coefficients and ( ), = 0, 1, 2, . . ., are specific functions on the problem. By substituting (11) and (14) into (13), we get Equating the coefficients of like powers of , we get the following set of equations: Now, we solve these equations in such a way that ,1 (x, ) = 0. Therefore, the approximate solution may be obtained as

Applications
In this section, we present some examples with analytical solution to show the efficiency of methods described in the previous section for solving (1) and (2).

Conclusion
In the present work, we proposed a new homotopy perturbation method to solve nonlinear time-fractional biological population model and a system of three nonlinear timefractional partial differential equations. The new method for solving fractional-order differential equations is based on two component procedures and polynomial initial condition.
The new homotopy perturbation method is very simple in application and is less computational and more accurate in comparison with other mentioned methods. Unlike the Adomian decomposition method, the new homotopy perturbation method is free from the need to use Adomian polynomials. In this method, we do not need the Lagrange multiplier, correction functional, stationary conditions, and calculating integrals, which eliminate the complications that exist in the variational iteration method. The obtained results show that these approaches can solve the problem effectively and can be applied to many nonlinear differential equations.
International Scholarly Research Notices 9 The computations associated with the examples discussed above were performed by MAPLE.