An Efficient Method for Systems of Variable Coefficient Coupled Burgers' Equation with Time-Fractional Derivative

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.


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][2][3][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 1 ( ), 2 ( ), 1 ( ), 2 ( ), 1 ( ), and 2 ( ) 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.

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 1 [ , ] is given by It has the following properties: (i) exists for any ∈ [ , ], , , ≥ 0, and > −1.
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.
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 > 0 is defined as

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: ( , ) = ( ) + ( ) + ( , ) . (10) For solving system (7) by NHPM, we construct the following homotopy: where ∈ [0, 1] is an embedding or homotopy parameter, Now, let us write (12) in the following form: ( , ) = ,0 + ( ( ) + ( ) + ( , )) . (13) Applying the inverse operator, , which is the Riemann-Liouville fractional integral of order ≥ 0, on both sides of (13), we have Suppose that the initial approximation of (10) has the form where , ( ), = 0, 1, 2, . . . are unknown coefficients and ( ), = 0, 1, 2, . . . 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 ,1 ( , ) = 0. Therefore, the approximate solution may be obtained as

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

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