International Journal of Mathematics and Mathematical Sciences

Volume 2013 (2013), Article ID 421378, 8 pages

http://dx.doi.org/10.1155/2013/421378

## Using Homo-Separation of Variables for Solving Systems of Nonlinear Fractional Partial Differential Equations

Division of Informatics, Logistics and Management, School of Technology and Health (STH), Royal Institute of Technology (KTH), 100 44 Stockholm, Sweden

Received 18 March 2013; Accepted 24 May 2013

Academic Editor: Irena Lasiecka

Copyright © 2013 Abdolamir Karbalaie et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### 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) [1–3]. 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.

#### References

- R. L. Magin,
*Fractional Calculus in Bioengineering,*, Begell House, West Redding, Conn, USA, 2006. - R. L. Magin and M. Ovadia, “Modeling the cardiac tissue electrode interface using fractional calculus,”
*JVC/Journal of Vibration and Control*, vol. 14, no. 9-10, pp. 1431–1442, 2008. View at Publisher · View at Google Scholar · View at Scopus - J. Keener and J. Sneyd,
*Mathematical Physiology*, Springer, New York, NY, USA, 2004. - A. A. M. Arafa, S. Z. Rida, and M. Khalil, “Fractional modeling dynamics of HIV and CD4+ T-cells during primary infection,”
*Nonlinear Biomedical Physics*, vol. 6, no. 1, article 1, 2012. View at Publisher · View at Google Scholar · View at Scopus - I. S. Jesus, J. A. Tenreiro MacHado, and J. Boaventure Cunha, “Fractional electrical impedances in botanical elements,”
*Journal of Vibration and Control*, vol. 14, no. 9-10, pp. 1389–1402, 2008. View at Publisher · View at Google Scholar · View at Scopus - K. S. Cole, “Sketches of general and comparative demography,” in
*Cold Spring Harbor Symposia on Quantitative Biology*, vol. 1, pp. 107–116, Cold Spring Harbor Laboratory Press, 1933. View at Google Scholar - V. D. Djordjevi, J. Jari, B. Fabry, J. J. Fredberg, and D. Stamenovi, “Fractional derivatives embody essential features of cell rheological behavior,”
*Annals of Biomedical Engineering*, vol. 31, no. 6, pp. 692–699, 2003. View at Google Scholar - C. Coussot, S. Kalyanam, R. Yapp, and M. Insana, “Fractional derivative models for ultrasonic characterization of polymer and breast tissue viscoelasticity,”
*IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control*, vol. 56, no. 4, pp. 715–725, 2009. View at Publisher · View at Google Scholar · View at Scopus - F. Mainardi, “Fractional diffusive waves in viscoelastic solids,” in
*Nonlinear Waves in Solids*, pp. 93–97, Fairfield, 1995. View at Google Scholar - V. G. Gupta and B. Sharma, “Application of Sumudu transform in reaction-diffusion systems and nonlinear waves,”
*Applied Mathematical Sciences*, vol. 4, no. 9–12, pp. 435–446, 2010. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - I. Podlubny,
*Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications*, vol. 198 of*Mathematics in Science and Engineering*, Academic Press, San Diego, Calif, USA, 1999. View at MathSciNet - S. Momani and Z. Odibat, “Numerical approach to differential equations of fractional order,”
*Journal of Computational and Applied Mathematics*, vol. 207, no. 1, pp. 96–110, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - O. Abdulaziz, I. Hashim, and S. Momani, “Solving systems of fractional differential equations by homotopy-perturbation method,”
*Physics Letters A*, vol. 372, no. 4, pp. 451–459, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Ganjiani and H. Ganjiani, “Solution of coupled system of nonlinear differential equations using homotopy analysis method,”
*Nonlinear Dynamics of Nonlinear Dynamics and Chaos in Engineering Systems*, vol. 56, no. 1-2, pp. 159–167, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Jafari, M. Nazari, D. Baleanu, and C. M. Khalique, “A new approach for solving a system of fractional partial differential
equations,”
*Computers & Mathematics with Applications*, 2012. View at Publisher · View at Google Scholar - S. Kumar, A. Yildirim, and L. Wei, “A fractional model of the diffusion equation and its analytical solution using Laplace transform,”
*Scientia Iranica*, vol. 19, no. 4, pp. 1117–1123, 2012. View at Publisher · View at Google Scholar - D. Kumar, J. Singh, and S. Rathore, “Sumudu decomposition method for nonlinear equations,”
*International Mathematical Forum*, vol. 7, no. 9–12, pp. 515–521, 2012. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Singh, D. Kumar, and Sushila, “Homotopy perturbation sumudu transform method for nonlinear equations,”
*Advances in Theoretical and Applied Mechanics*, vol. 4, no. 4, pp. 165–175, 2011. View at Google Scholar - A. Karbalaie, M. M. Montazeri, and H. H. Muhammed, “New approach to find the exact solution of fractional partial differential equation,”
*WSEAS Transactions on Mathematics*, vol. 11, no. 10, pp. 908–917, 2012. View at Google Scholar - S. G. Samko, A. A. Kilbas, and O. I. Marichev,
*Fractional Integrals and Derivatives: Theory and Applications*, Gordon and Breach Science Publishers, Yverdon, Switzerland, 1993. View at MathSciNet - A. A. M. Arafa, S. Z. Rida, and H. Mohamed, “Homotopy analysis method for solving biological population model,”
*Communications in Theoretical Physics*, vol. 56, no. 5, pp. 797–800, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Peng and K. Li, “A note on property of the Mittag-Leffler function,”
*Journal of Mathematical Analysis and Applications*, vol. 370, no. 2, pp. 635–638, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J.-F. Cheng and Y.-M. Chu, “Solution to the linear fractional differential equation using Adomian decomposition method,”
*Mathematical Problems in Engineering*, vol. 2011, Article ID 587068, 14 pages, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z. M. Odibat, “Analytic study on linear systems of fractional differential equations,”
*Computers & Mathematics with Applications*, vol. 59, no. 3, pp. 1171–1183, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. He, “Application of homotopy perturbation method to nonlinear wave equations,”
*Chaos, Solitons & Fractals*, vol. 26, no. 3, pp. 695–700, 2005. View at Publisher · View at Google Scholar · View at Scopus - J. H. He, “Homotopy perturbation method for bifurcation of nonlinear problems,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 6, pp. 207–208, 2005. View at Google Scholar