#### Abstract

A new homotopy perturbation method (NHPM) is applied to system of variable coefficient coupled Burgers' equation with time-fractional derivative. The fractional derivatives are described in the Caputo fractional derivative sense. The concept of new algorithm is introduced briefly, and NHPM is examined for two systems of nonlinear Burgers' equation. In this approach, the solution is considered as a power series expansion that converges rapidly to the nonlinear problem. The new approximate analytical procedure depends on two iteratives. The modified algorithm provides approximate solutions in the form of convergent series with easily computable components. Results indicate that the introduced method is promising for solving other types of systems of nonlinear fractional-order partial differential equations.

#### 1. Introduction

In recent years, the differential equations of fractional order have been the focus of many studies due to their frequent appearance in various applications in fluid mechanics, medical sciences, biological research, as well as various chemical, biochemical, and physical fields, viscoelasticity, biology, physics, and engineering. Consequently, considerable attention has been given to the solutions of fractional differential equations and integral equations of physical interest [1–4]. Various powerful methods have been presented so far such as homotopy perturbation method [5, 6], variational iteration method [7], differential transform method [8], homotopy analysis method [9], and Adomian decomposition method [10, 11] for solving different kinds of fractional partial differential equations. In this paper, we construct the solution of a system of variable coefficient coupled Burgers’ equation with time-fractional derivative by extending the idea of [12, 13]. A new version of homotopy perturbation method is proposed, which we called it NHPM and then applied it to the nonlinear systems of variable coefficient coupled Burgers’ equation with time-fractional derivative that can be written as the following basic form: subject to the initial condition where the subscripts , , , , , and are arbitrary smooth functions of .

The paper is organized as follows. In Section 2, we begin with an introduction to some necessary definitions of fractional calculus theory. In Section 3, we illustrated a basic idea of the new method. In Section 4, the uses of the new method for solving nonlinear variable coefficient coupled Burgers’ equation are presented. Two examples are solved by the proposed method in this section. Conclusion will appear in Section 5.

#### 2. Fractional Calculus

We give some basic definitions and properties of the fractional calculus theory used in this work. Some of these are Riemann-Liouville, Grunwald-Letnikov, Caputo, and generalized functions approach. The most commonly used definitions are the Riemann-Liouville and Caputo derivatives.

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

The Riemann-Liouville fractional derivative is mostly used by mathematicians, but this approach is not suitable for physical problems of the real world since it requires the definition of fractional order initial conditions which have no physically meaningful explanation yet. Caputo introduced an alternative definition, which has the advantage of defining integer-order initial conditions for fractional order differential equations.

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

Lemma 3. *If , , and , then
**
The Caputo fractional derivative is considered here because it allows traditional initial and boundary conditions to be included in the formulation of the problem. In this paper, we have considered some systems of linear and nonlinear FPDEs, where fractional derivatives are taken in Caputo sense as follows.*

*Definition 4. *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

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, we construct the following homotopy: where is an embedding or homotopy parameter, , and are the initial approximation of solution of the problem in (10).

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 (12) in the following form: Applying the inverse operator, , which is the Riemann-Liouville fractional integral of order , on both sides of (13), we have Suppose that the initial approximation of (10) has the form where , are unknown coefficients and , are specific functions on the problem. By substituting (12) and (15) into (14), 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. Examples

In this section, to illustrate the method and to show the ability of the method, two examples are presented.

*Example 1. *Consider the following variable coefficient coupled Burgers’ equation:
subject to the initial condition
The exact solutions of (19) for the special case are and .

To obtain the solution of (19) by NHPM, we construct the following homotopy:
Applying the inverse operator of on both sides of the above equation, we obtain
For solving system (22), by new homotopy perturbation method, we use the Taylor series of
The solution of (19) has the following form:
Substituting (23) and (24) in (22) and equating the coefficients of like powers of , we get the following set of equations:
Assuming , , , , and and solving the above equation for and lead to the result
Vanishing and lets the coefficients , have the following values:
Therefore, we obtain the solutions of (19) as
If we put in (28) or solve (19) with , we obtain the exact solution

*Example 2. *Consider the following variable coefficient coupled Burgers’ equation:
subject to the initial condition
The exact solution for is and .

To obtain the solution of (30) by NHPM, we construct the following homotopy:
Applying the inverse operator of on both sides of the above equation, we obtain
For solving system (33), by new homotopy perturbation method, we use the Taylor series of
The solution of (30) has the following form:
Substituting (34) and (35) in (33) and equating the coefficients of like powers of , we get the following set of equations:
Assuming , , , , and and solving the above equation for and lead to the result
Vanishing and lets the coefficients , have the following values:
Therefore, we obtain the solutions of (30) as
If we put in (39) or solve (30) with , we obtain the exact solution

#### 5. Concluding Remarks

In this paper, we have used a new homotopy perturbation method for solving a system of two nonlinear time-fractional partial differential equations. The NHPM for solving system of variable coefficient coupled Burgers’ equation with time-fractional derivative is based on two-component procedure and polynomial initial condition. The Computations finally lead to a set of nonlinear equations with one unspecified value in each equation. This set can be readily solved using Maple, and putting these values into the first approximate solution yields the analytical approximate solution. The present study has confirmed that NHPM offers significant advantages in terms of its straightforward applicability, computational efficiency, and accuracy. Thus, we conclude that the new method can be considered as an efficient method for solving linear and nonlinear problems.

#### Conflict of Interests

The authors of the paper do not have a direct financial relation that might lead to a “conflict of interests” for any of the authors.