A Novel Effective Approach for Solving Fractional Nonlinear PDEs

Abstract

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.

1. 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–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–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]: where denotes population density, represents the population supply due to births and deaths, , , , and are 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 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 time-fractional partial equations are presented. Finally, we give a brief conclusion in the last section.

2. 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.

Definition 1. The Riemann-Liouville fractional integral operator of order on the usual Lebesgue space is given by It has the following properties:(i) exists for any ,(ii),(iii),(iv),(v),

where , , and .

Definition 2. The Caputo definition of fractal derivative operator is given by where ,  , , and .

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

3. 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: For solving system (7) by NHPM [22] we construct the following homotopy: where is an embedding or homotopy parameter, , and are the initial approximation of solution of the problem in (9).

Clearly, the homotopy equations and are equivalent to the equations and , respectively. Thus, a monotonous change of parameter from zero to one corresponds to a continuous change of the trivial problem to the original problem. Next, we assume that the solution of equation can be written as a power series in embedding parameter , as follows: Now, let us write (11) in the following form: Applying the inverse operator which is the Riemann-Liouville fractional integral of order , on both sides of (12), we have

Suppose that the initial approximation of (9) has the form where , are unknown coefficients and , 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 . Therefore, the approximate solution may be obtained as

4. 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).

Example 1. Let us consider (1) with , , and , then we have
The exact solution of (18), for the special case , is .
To obtain the solution of (18) by new homotopy perturbation method, we construct the following homotopy: Applying the inverse operator of on both sides of (19), we obtain
The solution of (21) has the following form: Substituting (21) into (20) and equating the coefficients of like powers of , we get the following set of equations:
Assuming , , and and solving (22) for lead to the following result: Vanishing lets the coefficients , have the following values: Therefore, we obtain the solutions of (18) as If , we have which is an exact solution.

Example 2. Consider (1) with , and , then we have subject to the initial condition The exact solution of (27), for the special case , is To obtain the solution of (27) by the new homotopy perturbation method, we construct the following homotopy: Applying the inverse operator of on both sides of (30), we obtain The solution of (31) has the following form: Substituting (32) into (33) and equating the coefficients of like powers of , we get the following set of equations: Assuming , , , and solving (33) for lead to the following result: Vanishing lets the coefficients , have the following values: Therefore, we obtain the solution of (27) as If we put in (36) or solve (30) with , we obtain the exact solution

Example 3. Let us consider (2) with , , and . Therefore, we have subject to the initial conditions The exact solution of (38), for the special case , is To obtain the solution of (38) by new homotopy perturbation method, we construct the following homotopy: Applying the inverse operator of on both sides of the above equation, we obtain The solution of (42) has the following form: Substituting (43) into (42) and equating the coefficients of like powers of , we get the following set of equations:
Assuming and solving the above equation for and and lead to the following result: Vanishing , , and lets the coefficients , have the following values: Therefore, we obtain the solutions of (38) as
If we solve (38) with , we obtain the exact solution

5. 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 time-fractional 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. The computations associated with the examples discussed above were performed by MAPLE.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors are very grateful to anonymous referees for their careful reading and valuable comments which led to the improvement of this paper.