Analytical and Numerical Methods for Solving Partial Differential Equations and Integral Equations Arising in Physical Models
View this Special IssueResearch Article  Open Access
Numerical Solution of the Fractional Partial Differential Equations by the TwoDimensional FractionalOrder Legendre Functions
Abstract
A numerical method is presented to obtain the approximate solutions of the fractional partial differential equations (FPDEs). The basic idea of this method is to achieve the approximate solutions in a generalized expansion form of twodimensional fractionalorder Legendre functions (2DFLFs). The operational matrices of integration and derivative for 2DFLFs are first derived. Then, by these matrices, a system of algebraic equations is obtained from FPDEs. Hence, by solving this system, the unknown 2DFLFs coefficients can be computed. Three examples are discussed to demonstrate the validity and applicability of the proposed method.
1. Introduction
Fractional partial differential equations play a significant role in modeling physical and engineering processes. Therefore, there is an urgent need to develop efficient and fast convergent methods for FPDEs. Recently, several different techniques, including Adomian’s decomposition method (ADM) [1, 2], homotopy perturbation method (HPM) [3–5], variational iteration method (VIM) [6–8], spectral methods [9–13], orthogonal polynomials method [14–17], and wavelets method [18–21] have been presented and applied to solve FPDEs.
The method based on the orthogonal functions is a wonderful and powerful tool for solving the FDEs and has enjoyed many successes in this realm. The operational matrix of fractional integration has been determined for some types of orthogonal polynomials, such as Chebyshev polynomials [16], Legendre polynomials [22], Laguerre polynomials [23–25], and Jacobi polynomials [26]. Moreover, the operational matrix of fractional derivative for Chebyshev polynomials [9] and Legendre polynomials [9, 14] also has been derived. However, since these polynomials using integer power series to approximate fractional ones, it cannot accurately represent properties of fractional calculus. Recently, Rida and Yousef [27] presented a fractional extension of the classical Legendre polynomials by replacing the integer order derivative in Rodrigues formula with fractional order derivatives. The defect is that the complexity of these functions made them unsuitable for solving FDEs. Subsequently, Kazem et al. [28] presented the orthogonal fractional order Legendre functions based on shifted Legendre polynomials to find the numerical solution of FDEs and drew a conclusion that their method is accurate, effective, and easy to implement.
Benefiting from their “exponentialconvergence” property when smooth solutions are involved, spectral methods have been widely and effectively used for the numerical solution of partial differential equations. The basic idea of spectral methods is to expand a function into sets of smooth global functions, called the trial functions. Because of their special properties, the orthogonal polynomials are usually chosen to be trial functions. Spectral methods can obtain very accurate approximations for a smooth solution while only need a few degrees of freedom. Recently, Chebyshev spectral method [9], Legendre spectral method [10], and adaptive pseudospectral method [11] were proposed for solving fractional boundary value problems. Moreover, generalized Laguerre spectral algorithms and Legendre spectral Galerkin method were developed by Baleanu et al. [12] and Bhrawy and Alghamdi [13] for fractional initial value problems, respectively.
Motivated and inspired by the ongoing research in orthogonal polynomials methods and spectral methods, we construct twodimensional fractionalorder Legendre functions and derive the operational matrices of integration and derivative for the solution of FPDEs. To the best of the authors’ knowledge, such approach has not been employed for solving FPDEs.
The rest of the paper is organized as follows. In Section 2, we introduce some mathematical preliminaries of the fractional calculus theory and fractionalorder Legendre functions. In Section 3, a basis of 2DFLFs is defined and some properties are given. Section 4 is devoted to the operational matrices of fractional derivative and integration for 2DFLFs. Some numerical examples are presented in Section 5. Finally, we conclude the paper with some remarks.
2. Preliminaries and Notations
2.1. Fractional Calculus Theory
Some necessary definitions and Lemma of the fractional calculus theory [29, 30] are listed here for our subsequent development.
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 and only if , .
Definition 2. RiemannLiouville fractional integral operator () of order , of a function , is defined as where is the wellknown Gamma function. Some properties of the operator can be found, for example, in [29, 30].
Definition 3. The fractional derivative of in the Caputo sense is defined as where , denotes a continuous (but not necessarily differentiable) function.
Lemma 4. Let , , , , . Then
2.2. FractionalOrder Legendre Functions
In this section, we introduce the fractionalorder Legendre functions which were first proposed by Kazem et al. [28]. The normalized eigenfunctions problem for FLFs is which is a singular SturmLiouville problem. The fractionalorder Legendre polynomials, denoted by , are defined on the interval and can be determined with the aid of following recurrence formulae: and the analytic form of of degree is given by where and . The orthogonality condition is where is the weight function and is the Kronecker delta. For more details, please see [28].
3. 2DFLFs
In this section, the definitions and theorems of 2DFLFs are given by Liu’s method described in [31].
3.1. Definitions and Properties of the 2DFLFs
Definition 5. Let be the fractional Legendre polynomials on ; we call the twodimensional fractional Legendre polynomials on .
Theorem 6. The basis is orthogonal on with the weight function .
Proof. Let or
Theorem 7. Consider
3.2. 2DFLFs Expansion
Definition 8. A function of two independent variables which is integrable in square can be expanded as where
Theorem 9. If the series converges uniformly to on the square , then we have
Proof. By multiplying on both sides of (10), where and are fixed and integrating termwise with regard to and on , then Finally one can get (11).
If the infinite series in (10) is truncated, then it can be written as where and are given by where , , and .
According to the definition of FLFs, one can find that fractional Legendre polynomials are identical to Legendre polynomials shifted to when using the transform , . Therefore, in a similar method described in [31], we can easily get the convergence and stability theorems of proposed method.
Lemma 10. If the function is a continuous function on and the series converges uniformly to , then is the 2DFLFs expansion of .
Proof (by contradiction). Let Then there is at least one coefficient such that . However,
Lemma 11. If two continuous functions defined on have the identical 2DFLFs expansions, then these two functions are identical.
Proof. Suppose that and can be expanded by 2DFLFs as follows: By subtracting the above two equations with each other, one has Then Lemma 11 can be proved.
Theorem 12. If the 2DFLFs expansion of a continuous function converges uniformly, then the 2DFLFs expansion converges to the function .
Proof. Theorem 12 can be proved by Theorems 7 and 9.
Theorem 13. If the sum of the absolute values of the 2DFLFs coefficients of a continuous function forms a convergent series, then the 2DFLFs expansion is absolutely uniformly convergent, and converges to the function .
Proof. Consider Then converges uniformly to the function .
Theorem 14. If a continuous function , defined on , has bounded mixed partial derivative , then the 2DFLFs expansion of the function converges uniformly to the function.
Proof. Let be a function defined on such that where is a positive constant and By employing the transform and , one can obtain Consequently, in a similar method described in [31], Theorem 14 can be proved.
4. Operational Matrices of 2DFLFs
4.1. Integration Operational Matrices of 2DFLFs
Lemma 15. The RiemannLiouville fractional integration of order of the 2DFLFs can be obtained in the form of
Proof. Consider
Lemma 16. Let ; then one has
Proof. Using previous Lemma 15 and (6), one can have
Theorem 17. Let be the 2DFLFs vector defined in (16); then one has where is the operational matrix of RiemannLiouville fractional integration of order , and has the form as follows: in which is matrix and the elements are defined as follows: and is identity matrix.
Proof. Using (29) and orthogonality property of FLFs, one can get
where and are two matrices defined as
Now by substituting above equations in (32), Theorem 12 can be proved.
In a similar way as previous, one can obtain the operational matrix of RiemannLiouville fractional integration with respect to variable .
Theorem 18. Let be the 2DFLFs vector defined in (16); one has where is the operational matrix of RiemannLiouville fractional integration of order , and has the form as follows: in which is matrix and the elements are defined as follows:
4.2. Derivative Operational Matrices of 2DFLFs
Lemma 19. The FLFs Caputo fractional derivative of can be obtained in the form of where when and in other case .
Proof. Consider
Lemma 20. Let , ; then one has
Proof. Using previous Lemma 19 and (6), one can have
Theorem 21. Let be the 2DFLFs vector defined in (16); one has where is the operational matrix of Caputo fractional derivative of order , and has the form as follows: in which is matrix and the elements are defined as follows: and is a identity matrix.
Proof. Using (41) and the orthogonality property of FLFs, one can have
where and are two matrices defined as
Now by substituting above equations in (44), Theorem 21 can be proved.
In a similar way as above, one can get Caputo fractional derivative of order with respect to variable .
Theorem 22. Let be the 2DFLFs vector defined in (16); one can have where is the operational matrix of Caputo fractional derivative of order , and has the form as follows: in which is matrix and the elements are defined as follows:
5. Applications and Results
Consider the following FPDEs: where and are linear operator and nonlinear operator; respectively. and are the Caputo fractional derivatives of order and , respectively; is a known analytic function.
By employing operator on both sides of (49) and then using the Lemma 4, one can have
We first express unknown function and derivative term as
Now for the nonlinear part, by employing the nonlinear term approximation method described in [32] and then by using transform , , one can get the 2DFLFs expansion of nonlinear term as where is coefficient matrix of nonlinear term which must be computed and its order is .
For the linear part, we have where is a matrix of order .
After substituting (51)–(53) into (50), one can obtain According to the Wu’s [33] technology for determining the initial iteration value, the initial iteration value is chosen as . The coefficient matrix can be computed by using the MATLAB function fsolve( ) or the method described in [34].
Now, the present method is applied to solve the linear and nonlinear FPDEs, and their results are compared with the solution of other methods. The accuracy of our approach is estimated by the following error functions:
Example 23. Consider the onedimensional linear inhomogeneous fractional Burger’s equation [35]:
with the initial condition and the exact solution being .
By employing 2DFLFs method, one can get
where . Then we can get .
Figures 1(a) and 1(b) show the numerical results for with , and with , , respectively. It should be found that the accuracy of 2DFLFs method is very high while only a small number of 2DFLFs are needed.
(a)
(b)
(c)
(d)
Example 24. Consider nonlinear fractional KleinGordon equation [36, 37]:
subject to the initial conditions
and . The exact solution of (58) is .
By employing 2DFLFs method with and , one can have
The numerical results of Example 24 for different values of and are shown in Figure 2. In addition, and errors are presented in Table 1. From Table 1, one can conclude that the solutions of 2DFLFs method are in good agreement with the exact results. Compared with homotopy analysis method (HAM) [36] and homotopy perturbation method (HPM) [37], 2DFLFs method can get high accuracy solution while only need a few terms of 2DFLFs.

(a)
(b)
(c)
(d)
Example 25. Consider the nonlinear timefractional advection partial differential equation [37–39]
subject to the initial condition
Figure 3 gives the approximation solutions of (61) for with , and with , . Moreover, Table 2 shows the approximate solutions for (61) obtained for different values of using the fractional variational iteration method (FVIM) [39] and 2DFLFs method. The values of are the only case for which we know the exact solution . It should be noted that only the fourthorder term of the FVIM was used in evaluating the approximate solutions for Table 2. From Table 2, it clearly appears that 2DFLFs method is more accurate than FVIM and the obtained results are in good agreement with exact solution.

(a)
(b)
Example 26. We finally consider the linear timefractional wave equation:
subject to the initial conditions
Table 3 gives a comparison of the approximate solutions at different values of using the FVIM [39] and 2DFLFs method. Figure 4 shows the numerical solutions of 2DFLFs method for (63) at different values of with , . The values of are the only case for which we know the exact solution . As previous, only the fourthorder term of the FVIM was used in evaluating the numerical solutions for Table 3. In the case of , it can be found that absolute error of 2DFLFs is not bigger than which is very small compared with that obtained by FVIM.

(a)
(b)
(c)
(d)
6. Conclusion
We define a basis of 2DFLFs and derived its operational matrices of fractional derivative and integration, which are used to approximate the numerical solution of FPDEs. Compared with other numerical methods, 2DFLFs method can accurately represent properties of fractional calculus. Moreover, only a small number of 2DFLFs are needed to obtain a satisfactory result. The obtained results demonstrate the validity and applicability of proposed method for solving the FPFEs.
Acknowledgments
This work is supported by the National Natural Science Foundation of China (Grant no. 11272352). The authors are grateful to the anonymous referees for their comments which substantially improved the quality of this paper.
References
 Z. Odibat and S. Momani, “Numerical methods for nonlinear partial differential equations of fractional order,” Applied Mathematical Modelling, vol. 32, no. 1, pp. 28–39, 2008. View at: Publisher Site  Google Scholar
 S. Momani and Z. Odibat, “Analytical approach to linear fractional partial differential equations arising in fluid mechanics,” Physics Letters A, vol. 355, no. 45, pp. 271–279, 2006. View at: Publisher Site  Google Scholar
 Z. Odibat and S. Momani, “Modified homotopy perturbation method: application to quadratic Riccati differential equation of fractional order,” Chaos, Solitons & Fractals, vol. 36, no. 1, pp. 167–174, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. H. Hosseinnia, A. Ranjbar, and S. Momani, “Using an enhanced homotopy perturbation method in fractional differential equations via deforming the linear part,” Computers & Mathematics with Applications, vol. 56, no. 12, pp. 3138–3149, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 O. Abdulaziz, I. Hashim, and S. Momani, “Solving systems of fractional differential equations by homotopyperturbation method,” Physics Letters A, vol. 372, no. 4, pp. 451–459, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. H. He, “Approximate analytical solution for seepage flow with fractional derivatives in porous media,” Computer Methods in Applied Mechanics and Engineering, vol. 167, no. 12, pp. 57–68, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. H. He, “Variational iteration method—a kind of nonlinear analytical technique: some examples,” International Journal of NonLinear Mechanics, vol. 34, no. 4, pp. 699–708, 1999. View at: Google Scholar
 J. H. He, “Variational iteration method—some recent results and new interpretations,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 3–17, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. H. Doha, A. H. Bhrawy, and S. S. EzzEldien, “A Chebyshev spectral method based on operational matrix for initial and boundary value problems of fractional order,” Computers & Mathematics with Applications, vol. 62, no. 5, pp. 2364–2373, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. H. Bhrawy and M. M. AlShomrani, “A shifted Legendre spectral method for fractionalorder multipoint boundary value problems,” Advances in Difference Equations, vol. 2012, article 8, 19 pages, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 M. Maleki, I. Hashim, M. Tavassoli Kajani, and S. Abbasbandy, “An adaptive pseudospectral method for fractional order boundary value problems,” Abstract and Applied Analysis, vol. 2012, Article ID 381708, 19 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. Baleanu, A. H. Bhrawy, and T. M. Taha, “Two efficient generalized Laguerre spectral algorithms for fractional initial value problems,” Abstract and Applied Analysis, vol. 2013, Article ID 546502, 10 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 A. H. Bhrawy and M. A. Alghamdi, “A new legendre spectral galerkin and pseudospectral approximations for fractional initial value problems,” Abstract and Applied Analysis, vol. 2013, Article ID 306746, 10 pages, 2013. View at: Publisher Site  Google Scholar
 A. Saadatmandi and M. Dehghan, “A new operational matrix for solving fractionalorder differential equations,” Computers & Mathematics with Applications, vol. 59, no. 3, pp. 1326–1336, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. H. Bhrawy, A. S. Alofi, and S. S. EzzEldien, “A quadrature tau method for fractional differential equations with variable coefficients,” Applied Mathematics Letters, vol. 24, no. 12, pp. 2146–2152, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. H. Bhrawy and A. S. Alofi, “The operational matrix of fractional integration for shifted Chebyshev polynomials,” Applied Mathematics Letters, vol. 26, no. 1, pp. 25–31, 2013. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Yüzbasi, “Numerical solution of the BagleyTorvik equation by the Bessel collocation method,” Mathematical Methods in the Applied Sciences, vol. 36, no. 3, pp. 300–312, 2013. View at: Publisher Site  Google Scholar
 J. L. Wu, “A wavelet operational method for solving fractional partial differential equations numerically,” Applied Mathematics and Computation, vol. 214, no. 1, pp. 31–40, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Saadatmandi and M. Dehghan, “A new operational matrix for solving fractionalorder differential equations,” Computers & Mathematics with Applications, vol. 59, no. 3, pp. 1326–1336, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Jafari, S. A. Yousefi, M. A. Firoozjaee, S. Momani, and C. M. Khalique, “Application of Legendre wavelets for solving fractional differential equations,” Computers & Mathematics with Applications, vol. 62, no. 3, pp. 1038–1045, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Li, “Solving a nonlinear fractional differential equation using Chebyshev wavelets,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 9, pp. 2284–2292, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Ahmadian, M. Suleiman, and S. Salahshour, “An operational matrix based on Legendre polynomials for solving fuzzy fractionalorder differential equations,” Abstract and Applied Analysis, vol. 2013, Article ID 505903, 29 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 A. H. Bhrawy, D. Baleanu, L. M. Assas, and J. A. T. Machado, “On a generalized Laguerre operational matrix of fractional integration,” Mathematical Problems in Engineering, Article ID 569286, 7 pages, 2013. View at: Google Scholar  MathSciNet
 A. H. Bhrawy and T. M. Taha, “An operational matrix of fractional integration of the Laguerre polynomials and its application on a semiinfinite interval,” Mathematical Sciences, vol. 6, article 41, 7 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. H. Bhrawy, M. M. Alghamdi, and T. M. Taha, “A new modified generalized Laguerre operational matrix of fractional integration for solving fractional differential equations on the half line,” Advances in Difference Equations, p. 2012:179, 12, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 S. Kazem, “An integral operational matrix based on Jacobi polynomials for solving fractionalorder differential equations,” Applied Mathematical Modelling, vol. 37, no. 3, pp. 1126–1136, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 S. Z. Rida and A. M. Yousef, “On the fractional order Rodrigues formula for the Legendre polynomials,” Advances and Applications in Mathematical Sciences, vol. 10, no. 5, pp. 509–517, 2011. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 S. Kazem, S. Abbasbandy, and S. Kumar, “Fractionalorder Legendre functions for solving fractionalorder differential equations,” Applied Mathematical Modelling, vol. 37, no. 7, pp. 5498–5510, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 I. Podlubny, Fractional Differential Equations, vol. 198, Academic Press, San Diego, Calif, USA, 1999. View at: MathSciNet
 I. Podlubny, “Geometric and physical interpretation of fractional integration and fractional differentiation,” Fractional Calculus & Applied Analysis, vol. 5, no. 4, pp. 367–386, 2002. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 N. Liu and E. B. Lin, “Legendre wavelet method for numerical solutions of partial differential equations,” Numerical Methods for Partial Differential Equations, vol. 26, no. 1, pp. 81–94, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. Yin, J. Song, X. Cao, and F. Lu, “Couple of the variational iteration method and Legendre wavelets for nonlinear partial differential equations,” Journal of Applied Mathematics, vol. 2013, Article ID 157956, 11 pages, 2013. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. C. Wu, “Challenge in the variational iteration method—a new approach to identification of the Lagrange multipliers,” Journal of King Saud University, vol. 25, pp. 175–178, 2013. View at: Publisher Site  Google Scholar
 F. Yin, J. Song, and F. Lu, “A coupled method of Laplace transform and Legendre wavelets for nonlinear KleinGordon equations,” Mathematical Methods in the Applied Sciences, 2013. View at: Publisher Site  Google Scholar
 N. Imran and S. T. MohyudDin, “Decomposition method for fractional partial differential equation (PDEs) using Laplace transformation,” International Journal of Physical Sciences, vol. 8, no. 16, pp. 684–688, 2013. View at: Publisher Site  Google Scholar
 K. A. Gepreel and M. S. Mohamed, “Analytical approximate solution for nonlinear spacetime fractional Klein Gordon equation,” Chinese Physics B, vol. 22, no. 1, Article ID 010201, 2013. View at: Publisher Site  Google Scholar
 A. M. A. ElSayed, A. Elsaid, and D. Hammad, “A reliable treatment of homotopy perturbation method for solving the nonlinear KleinGordon equation of arbitrary (fractional) orders,” Journal of Applied Mathematics, vol. 2012, Article ID 581481, 13 pages, 2012. View at: Publisher Site  Google Scholar
 S. Momani and Z. Odibat, “A novel method for nonlinear fractional partial differential equations: combination of DTM and generalized Taylor’s formula,” Journal of Computational and Applied Mathematics, vol. 220, no. 12, pp. 85–95, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Song, F. Yin, X. Cao, and F. Lu, “Fractional variational iteration method versus Adomian’s decomposition method in some fractional partial differential equations,” Journal of Applied Mathematics, vol. 2013, Article ID 392567, 10 pages, 2013. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2013 Fukang Yin 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.