Abstract

A new method proposed and coined by the authors as the homo-separation of variables method is utilized to solve systems of linear and nonlinear fractional partial differential equations (FPDEs). The new method is a combination of two well-established mathematical methods, namely, the homotopy perturbation method (HPM) and the separation of variables method. When compared to existing analytical and numerical methods, the method resulting from our approach shows that it is capable of simplifying the target problem at hand and reducing the computational load that is required to solve it, considerably. The efficiency and usefulness of this new general-purpose method is verified by several examples, where different systems of linear and nonlinear FPDEs are solved.

1. Introduction

In the recent years, it has turned out that many phenomena in various technical and scientific fields can be described very successfully by using fractional calculus. In particular, fractional calculus can be employed to solve many problems within the biomedical research field and get better results. Such a practical application of fractional order models is to use these models to improve the behavior and efficiency of bioelectrodes. The importance of this application is based on the fact that bioelectrodes are usually needed to be used for all types of biopotential recording and signal measurement purposes such as electroencephalography (EEG), electrocardiography (ECG), and electromyography (EMG) [13]. Another promising biomedical application field is proposed by Arafa et al. [4] where a fractional-order model of HIV-1 infection of CD4 T cells is introduced. Other examples of applications of fractional calculus in life sciences and technology can be found in botanics [5], biology [6], rheology [7], and elastography [8].

Numerous analytical methods have been presented in the literature to solve FPDEs, such as the fractional Greens function method [9], the Fourier transform method [2], the Sumudu transform method [10], the Laplace transform method, and the Mellin transform method [11]. Some numerical methods have also widely been used to solve systems of FPDEs, such as the variational iteration method [12], the Adomian decomposition method [2], the homotopy perturbation method [13] and the homotopy analysis method [14]. Some of these methods use specific transformations and others give the solution as a series which converges to the exact solution.

In addition, some numerical methods use a combination of utilizing specific transformations and obtaining series which converge to the exact solutions. An example of such a method is the iterative Laplace transform method which is a combination of the Laplace transform method and an iterative method [15]. Another such a combination is the homotopy perturbation transformation method, which is constructed by combining two powerful methods, namely, the Laplace transform method and the homotopy perturbation method [16]. A third example is the Sumudu decomposition method, which is a combination of the Sumudu transform method and Adomian decomposition method [17]. A fourth such an approach is combining the Sumudu transformation method with the homotopy perturbation method, which gives a new method called the homotopy perturbation Sumudu transform method [18].

Recently, the homotopy perturbation method and the Adomian decomposition method are frequently used for solving nonlinear FPDE problems.

According to the combinational methods mentioned above, it is possible to notice that HPM has a strong potential to be combined with another method to produce a more efficient approach. The main reason is that HPM is an efficient method for solving PDEs and ODEs (ordinary differential equations) with integer or fractional order. Recently, Karbalaie et al. [19] found the exact solution of one-dimensional FPDEs by using truncated versions of modified HPM. Therefore, by getting inspiration of the ideas, methods, and tools of previous works, the novel approach called the homo-separation of variables method is developed and utilized to find the exact solutions of systems of FPDEs. The methods presented in Karbalaie et al. [19] are extended and developed further to achieve methods that can solve n-dimensional systems of FPDEs. This new approach is constructed by a smart combination of HPM and the separation of variables method. By using this method, the system of FPDEs to be solved is changed into a system of FODEs. Consequently, the original problem is simplified considerably which makes it more straightforward and easier to solve.

The structure of the current paper is as follows. For easy reader-friendly reference, some basic definitions and properties of fractional differential equations are presented in Section 2. In Section 3, the homo-separation of variables method is described. After that, in Section 4, three examples are presented to show the simplicity and efficiency (at the same time) of the homo-separation of variables method. Finally, in Section 5, relevant conclusions are drawn.

2. Basic Definitions and Properties

In this section, to make it easier to introduce the new approach of this research work, a number of good-to-know concepts and properties within the field of fractional differential equations are defined and presented briefly.

Definition 1. A real function , is said to be in the space , if there exists a real number, such that , where , and it is said to be in the space if .

Definition 2. The left sided Riemann-Liouville fractional integral of order , of a function , is defined as where , and is the gamma function.
Properties of the operator (for ) can be found in [20] and are defined as follows:(1), (2), (3).

Definition 3. The left sided Caputo fractional derivative of , in Caputo sense is defined by [11] as follows:
Note that, according to [21], (1) and (2) become

Definition 4. The single-parameter and the two-parameter variants of the Mittag-leffler function are denoted by and , respectively. These two variants are relevant to be considered here because of their connection with fractional calculus and are defined as
The th derivatives of these two variants are
Some special cases and properties of the Mittag-Leffler functions are the following: (1), (2), (3), (4), (5), (6), (7).
Further properties of the Mittag-Leffler functions can be found in [21]. These functions are generalizations of the exponential function; therefore, most linear differential equations of fractional order (FPDEs) have solutions that are expressed in terms of these functions.

Important Note. According to [22], cannot satisfy the following semigroup property: for any noninteger . The above property is not true unless in the special cases where or . Unfortunately, this misunderstanding can be found in a number of publications, which can be misleading and confusing.

Definition 5. Consider the following relaxation-oscillation equation, which is discussed in [23]: subject to the initial condition where are real constants and . The solution of (7) is obtained according to [23] as follows: where is the Mittag-Leffler function. It is easy to realize that if , then the solution of (7) becomes as follows: which will be used to solve the systems of FPDEs in the examples discussed in this paper.

Definition 6. According to [24], consider the following linear system of fractional differential equations with its initial condition: where the matrix and is the Caputo fractional derivative of order, where . The general analytical solution for (11) is given by where are arbitrary constants, is the Mittag-Leffler function, and and are the eigenvalues and the corresponding eigenvectors of the characteristic equation , where and is an identity matrix.

3. The Homo-Separation of Variables Method

This approach is achieved by a novel combination of HPM and separation of variables. In this section, the algorithm of this new method is presented briefly by considering the following system of FPDEs: subject to the initial conditions where , and for , the terms represent linear operators, the terms are nonlinear operators, while the terms represent known analytical functions. In addition, is the Caputo fractional derivative of order , where for .

According to the homotopy perturbation technique (HPM), we can construct the following homotopies: where , , and the homotopy parameter, denoted by , is considered as small . In addition, for, is an initial approximation of the solution of (13) which satisfies the initial condition in (14).

We can assume that the solution of (15) can be expressed as a power series in , as given below in (16).

By substituting (16) into (15) and equating the terms with identical powers of , a set of equations is obtained as follows: where the terms and are obtained from (17)–(19) by equating the coefficients of the nonlinear operators and with the identical powers of , that is, the terms containing .

In case the -parameter is considered as small, the best approximate solution of (13) can be readily obtained as follows: where .

If in (20) there exists some , where and , then the exact solution can be written in the following form:

For simplicity, we assume that in (21), which means that the exact solution in (13) is simplified as follows:

Now when solving (17), we obtain the following result:

By using (22) and (23) we have

The core of the new method proposed in this research work is to formulate the initial approximation of (13) in the form of separation of variables, as follows:

The task now is to find the terms and to obtain the exact solution of the system of FPDEs in (13). Since (25) satisfies the initial conditions in (14), we get

By substituting (25) into (18), we get where and

Obviously, our goal is finally achieved and the system of FPDEs is changed into a system of FODEs. Consequently, the problem at hand is simplified considerably as it is well known that solving a system of FODEs is generally more straightforward and much easier than solving the corresponding system of FPDEs.

Finally, the target unknowns and can be obtained by utilizing and solving the system of FODEs and the initial conditions in (26).

4. Applications

In this section, we illustrate the applicability, simplicity, and efficiency of the homo-separation of variables method for solving systems of linear and nonlinear FPDEs.

Example 7. Consider the following system of linear FPDEs: where , subject to the initial conditions
To solve (29) and (30) by using the proposed homo-separation of variables method, we choose the initial approximations for (29) as follows:
By substituting the two initial conditions given by (30) into these two initial approximations, we get the following initial approximations:
By considering and following the approach which resulted in (27) and substituting the results presented by (32) into the system of linear FPDEs in (29), we get
By solving the latter system of linear FPDEs presented in (33), we obtain the following system of linear FODEs: together with the initial conditions
By solving (34) and (35) by considering and applying Definition 6 which is discussed in Section 2, we obtain the following results:
By substituting (36) into (32) and applying Definition 4, in particular properties , and which are discussed in Section 2, we obtain the exact solution for the system of linear FPDEs in (29), as follows:
If we now put in (37) or solve (29) and (30) for , we obtain the following exact solution for the corresponding system of linear PDEs:

Example 8. Consider the following system of inhomogeneous FPDEs: where and , subject to the initial conditions
To solve (39) and (40) by using the new homo-separation of variables method, we choose the initial approximations for (39) as follows:
By substituting the two initial conditions given by (40) into these two initial approximations, we get the following initial approximations:
By considering and following the approach which resulted in (27) and substituting the results presented by (42) into the system of inhomogeneous FPDEs that we would like to solve, which is presented by (39), we get
By solving the latter system of FPDEs presented in (43), we obtain the following system of FODEs: together with the initial conditions
By solving (44) and (45) by considering and applying Definition 5 which is discussed in Section 2, we obtain the following results:
By substituting (46) into (42), we obtain the exact solution for the system of inhomogeneous FPDEs presented in (39), as follows:
Thus, this means that . If we now put in (47), or solve (39) and (40) for , we obtain the following exact solution for the corresponding system of inhomogeneous PDEs:

Example 9. Consider the following system of nonlinear FPDEs: where and ,  subject to the initial conditions
To solve (49) and (50) by using the new homo-separation of variables method, we choose the initial approximations for (49) as follows:
By substituting the initial conditions given by (50) into these initial approximations, we get the following initial approximations:
By considering and following the approach which resulted in (27) and substituting the results presented by (52) into the system of nonlinear FPDEs that we would like to solve, which is presented by (49), we get
By solving the latter system of FPDEs presented in (53), we obtain the following system of FODEs: together with the initial conditions
By solving (54) and (55) by considering and applying Definition 5 which is discussed in Section 2 of this paper, we obtain the following results:
By substituting (56) into (52), we obtain the exact solution for the system of nonlinear FPDEs in (49), as follows:
If we now put in (57), or solve (49) and (50) for , we obtain the following exact solution for the corresponding system of nonlinear PDEs:

5. Conclusions

In this paper, a novel approach was introduced and utilized to solve linear and nonlinear systems of FPDEs. The new technique was coined by the authors of this paper as the homo-separation of variables method. In this research work, it was demonstrated through different examples how this new method can be used for solving various systems of FPDEs.

When compared with the existing published methods, it is easy to notice that the new method has many advantages. It is straightforward, easy to understand, and fast, requiring much less computations to perform a limited number of steps of the simple procedure that can be applied to find the exact solution of a wide range of types of FPDE equations’ systems. Furthermore, there is no need for using linearization or restrictive assumptions when employing this new method.

The basic principle of the current method is to perform an elegant combination (which is simple and smart) of two existing techniques. These two techniques are the popular homotopy perturbation method (HPM) that was originally proposed by He [25, 26] and another popular approach called separation of variables. Finally, the resulting homo-separation of variables method, which is analytical, can be used to solve various systems of PDEs with integer and fractional order, for example, with respect to time.

Acknowledgments

The authors would like to thank Maryam Shabani and Mehran Shafaiee for their help and support as well as the inspiration the authors got through friendly scientific discussions.