`Mathematical Problems in EngineeringVolume 2014 (2014), Article ID 728736, 8 pageshttp://dx.doi.org/10.1155/2014/728736`
Research Article

## An Accurate Spectral Galerkin Method for Solving Multiterm Fractional Differential Equations

1Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
2Department of Mathematics, Faculty of Science, Beni-Suef University, Beni-Suef 62511, Egypt

Received 4 May 2014; Accepted 20 May 2014; Published 12 June 2014

Copyright © 2014 A. H. Bhrawy and A. S. Alofi. 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

This paper reports a new formula expressing the Caputo fractional derivatives for any order of shifted generalized Jacobi polynomials of any degree in terms of shifted generalized Jacobi polynomials themselves. A direct solution technique is presented for solving multiterm fractional differential equations (FDEs) subject to nonhomogeneous initial conditions using spectral shifted generalized Jacobi Galerkin method. The homogeneous initial conditions are satisfied exactly by using a class of shifted generalized Jacobi polynomials as a polynomial basis of the truncated expansion for the approximate solution. The approximation of the spatial Caputo fractional order derivatives is expanded in terms of a class of shifted generalized Jacobi polynomials with , and is the polynomial degree. Several numerical examples with comparisons with the exact solutions are given to confirm the reliability of the proposed method for multiterm FDEs.

#### 1. Introduction

Fractional differential equations (FDEs), as generalizations of classical integer order differential equations, are increasingly used to model several real phenomena emerging in engineering and science fields. Owing to the increasing applications, there has been important interest in developing analytical and numerical methods for the solution of fractional differential equations (see e.g., [17] and the references therein). These methods include variational iteration method [8, 9], Adomian decomposition method [10, 11], generalized differential transform method [12], Laplace decomposition method [13], homotopy analysis method [14], spectral method [1519], finite difference method [2022], and wavelet methods [2325].

Spectral method is one of the principal methods of discretization for the numerical solution of most types of differential equations. The three most widely used spectral versions are the Galerkin, Tau, and collocation methods (see, for instance [2632]). Recently, spectral method is a class of important tools for obtaining the numerical solutions of fractional differential equations. They have excellent error properties and they offer exponential rates of convergence for smooth problems. In the present paper we intend to extend the application of Galerkin method based on generalized Jacobi polynomials form solving linear problems to solve multiterm FDEs. To the best of our knowledge, there are not so many results on using this technique to solve such problems arising in mathematical physics. This partially motivated our interest in such a method.

Spectral Galerkin method for the numerical solution of fractional differential equations is characterized by expanding the solution by a truncated series of the trial functions. The unknown coefficients of this expansion will be determined by minimizing the error between the exact and numerical solutions in appropriate weighted space. This method provides exponential rates of convergence. An explicit expression for the derivatives of an infinitely differentiable function of any degree and for any fractional order in terms of the function itself is needed. Doha et al. [16] have obtained such a relation in the case of the basis functions of expansion that are shifted Jacobi polynomials. Another formula for shifted Legendre coefficients is obtained by Bhrawy et al. [17]. Moreover, in [33] the authors expressed explicitly the Caputo fractional derivatives of generalized Laguerre polynomials of any degree in terms of the generalized Laguerre polynomials themselves to solve fractional initial value problems on the half line.

An explicit expression for any Caputo fractional order derivative of the shifted generalized Jacobi polynomials of any degree in terms of the shifted generalized Jacobi polynomials themselves is the first goal of this paper. The fundamental goal of this paper is to develop a direct solution technique based on shifted generalized Jacobi-Galerkin method (SGJG) for solving multiterm FDEs with homogeneous and nonhomogeneous initial conditions. Finally, we present some numerical results exhibiting the accuracy and efficiency of our numerical algorithm.

The next section of this paper is for fractional preliminaries. Section 3 is devoted to proving a formula that expresses the Caputo fractional order derivative of the shifted generalized Jacobi polynomials. In Section 4, we construct and develop algorithms for solving linear FDEs by using shifted generalized Jacobi Galerkin spectral method. In Section 5, several examples are presented. Finally, some concluding remarks are given in the last section.

#### 2. Preliminaries and Notations

In this section, we present some basic knowledge of fractional calculus, orthogonal shifted Jacobi polynomials, and generalized Jacobi polynomials these are most relevant to spectral approximations.

##### 2.1. The Fractional Derivative in the Caputo Sense

In this section, we state the definition and preliminaries of fractional calculus.

Definition 1. For to be the smallest integer that exceeds , Caputo’s fractional derivative operator of order is defined as where For the Caputo derivative we have
Similar to the integer-order differentiation, the Caputo’s fractional differentiation is a linear operation; that is, where and are constants.

##### 2.2. Classical Jacobi Polynomials

The Jacobi polynomials with the real parameters are a sequence of polynomials , satisfying the orthogonality relation where It is convenient to standardize the Jacobi polynomials so that where . In this form the polynomials may be generated using the standard recurrence relation of Jacobi polynomials starting from and , or obtained from Rodrigue’s formula where .

The shifted Jacobi polynomials [34, 35] are orthogonal polynomials on with respect to the weight function . Note that the shifted Jacobi polynomials satisfy the orthogonality relation where

If we denote by , and , to the nodes and Christoffel numbers of the standard (shifted) Legendre-Gauss-Lobatto quadratures on the intervals , respectively, then one can easily show that and if denotes the set of all polynomials of degree at most , then it follows that for any ,

We define the discrete inner product and norm as follows:

Obviously,

##### 2.3. Generalized Jacobi Polynomials

Recently, Guo et al. [36] presented and developed the generalized Jacobi approximation, in which the parameters and considered in the generalized Jacobi polynomials might be any real numbers. In this section, we give some properties of such polynomials. Let and . We denote the set of integers by . For any , the generalized Jacobi polynomials are defined by (see [36, 37]) For our present purposes it is convenient to use the shifted Jacobi polynomials ; let and . We define the shifted GJPs and separate them into four cases as follows.

Case 1.

Case 2.

Case 3.

Case 4. where , .

Lemma 2. Each of GJPs and forms a complete orthogonal system in , respectively. And the square of the norm of each of the four GJPs cases is defined as where is the square of the norm of the classical shifted Jacobi polynomials and

Proof. Firstly, Secondly, Thirdly, And lastly,

#### 3. The Fractional Derivatives of

The main objective of this section is to prove the following theorem for the fractional derivatives of the shifted generalized Jacobi polynomials. The analytic form of the shifted generalized Jacobi polynomials of degree is given by A function , square integrable in , can be expressed in terms of shifted generalized Jacobi polynomials as where the coefficients are given by

Lemma 3. Let be a shifted generalized Jacobi polynomial of degree ; then

Proof. This lemma can be easily proved by making use of relations (3)-(4) with relation (26).

Theorem 4. The fractional derivative of order in the Caputo sense for the shifted generalized Jacobi polynomials is given by where

Proof. The analytic form of the shifted generalized Jacobi polynomials of degree is given by (26). Using (3)-(4) and (26), we have Now, approximating by terms of shifted generalized Jacobi series, we have where is given from (28) with , and this immediately gives
Employing (32)–(34), we get where is given as in (30), and this proves the theorem.

#### 4. Shifted Generalized Jacobi Galerkin Method for FDEs

In this section, we are interested in employing the SGJG method to solve the linear multiterm FDE subject to the homogeneous initial conditions where and are constants, denotes the Caputo fractional derivative of order for , and is a given source function. Let us first introduce some basic notation that will be used in the upcoming sections. We set where denotes th-order differentiation of with respect to . Then the shifted generalized Jacobi-Galerkin approximation to (36) is to find such that where and is the inner product in the weighted space . The norm in will be denoted by . Let Then we can write (39) as follows: Let us denote

By virtue of (31) and making use of the orthogonality relation of shifted generalized Jacobi polynomials (21), and after some rather lengthy calculation, we get

Thereby, we can write (41) in the following matrix system form

##### 4.1. Treatment of the Nonhomogeneous Initial Conditions

In the following we can always modify the right-hand side to take care of the nonhomogeneous initial conditions. Let us consider for instance the one-dimensional fractional differential equation (36) subject to the nonhomogeneous initial conditions: In such a case we proceed as follows.

Setting where the transformation (46) turns the nonhomogeneous initial conditions (45) into the homogeneous initial conditions Hence it suffices to solve the following modified multiterm fractional differential equation: subject to the homogeneous initial conditions (48), where is given by (46), and

#### 5. Illustrative Examples

Several test examples are carried out in this section. The results obtained by the present methods reveal that the present method is very effective and convenient for linear FDEs.

Example 1. Consider the linear FDE equation with homogeneous initial conditions whose exact solution is given by .

Table 1 lists the maximum absolute errors, using the shifted generalized Jacobi Galerkin (SGJG) method with various choices of and. Accuracy and stability of the SGJG method for all choices of are achieved in this table.

Table 1: Maximum absolute errors with various choices of and , for Example 1.

Example 2. Consider the equation whose exact solution is given by .

In Table 2, we present the maximum absolute errors, using SGJG method with various choices of and.

Table 2: Maximum pointwise errors using SGJG method for .

Example 3. Consider the equation whose exact solution is given by .

In Table 3, we present the maximum absolute errors, using SGJG method with various choices of and.

Table 3: Maximum pointwise errors using SGJG method for .

Example 4. Consider the equation whose exact solution is given by .

In Table 4, we exhibit maximum pointwise error using SGJG method with two choices of the shifted generalized Jacobi parameters and . We observe from this table that the suggested algorithm provides accurate and stable numerical results. This numerical experiment demonstrates the utility of the method.

Table 4: Maximum pointwise errors using SGJG method for .

#### 6. Conclusion

We have derived a new formula expressing explicitly the Caputo fractional derivatives for any fractional-order of shifted generalized Jacobi polynomials of any degree in terms of shifted generalized Jacobi polynomials themselves. We have derived a Galerkin method, involving a specified class of the shifted generalized Jacobi polynomials, which permits us to numerically solve an important class of FDEs. Indeed, in Section 5, we demonstrated that for all parameter shifted generalized Jacobi considered, the method results in rather small errors with relatively few modes are considered. Since the method is rather robust, it is likely that it may be applied to other types of FDEs. For instance, one- and two-dimensional time-dependent FDEs

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

This paper was funded by the Deanship of Scientific Research DSR, King Abdulaziz University, Jeddah. The authors, therefore, acknowledge with thanks DSR technical and financial support.

#### References

1. D. Baleanu, K. Diethelm, E. Scalas, and J. J. Trujillo, Fractional Calculus, vol. 3 of Series on Complexity, Nonlinearity and Chaos, World Scientific, 2012.
2. 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.
3. A. H. Bhrawy and M. A. Alghamdi, “A shifted Jacobi-Gauss-Lobatto collocation method for solving nonlinear fractional Langevin equation involving two fractional orders in different intervals,” Boundary Value Problems, vol. 2012, article 62, 2012.
4. Y. Yang, Y. Chen, and Y. Huang, “Convergence analysis of the Jacobi spectral-collocation method for fractional integro-differential equations,” Acta Mathematica Scientia B, vol. 34, no. 3, pp. 673–690, 2014.
5. A. Atangana and N. Bildik, “The use of fractional order derivative to predict the groundwater flow,” Mathematical Problems in Engineering, vol. 2013, Article ID 543026, 9 pages, 2013.
6. A. Atangana and S. B. Belhaouari, “Solving partial differential equation with space- and time-fractional derivatives via homotopy decomposition method,” Mathematical Problems in Engineering, vol. 2013, Article ID 318590, 9 pages, 2013.
7. X.-J. Yang, D. Baleanu, and J.-H. He, “Transport equations in fractal porous media within fractional complex transform method,” Proceedings of the Romanian Academy A: Mathematics, Physics, Technical Sciences, Information Science, vol. 14, no. 4, pp. 287–292, 2013.
8. A. M. Yang, J. Li, H. M. Srivastava, G. Xie, and X. J. Yang, “The local fractional Laplace variational iteration method for solving linear partial differential equations with local fractional derivatives,” Discrete Dynamics in Nature and Society. In press.
9. F. Yin, J. Song, and X. Cao, “A general iteration formula of VIM for fractional heat- and wave-like equations,” Journal of Applied Mathematics, vol. 2013, Article ID 428079, 9 pages, 2013.
10. 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.
11. L. Song and W. Wang, “A new improved Adomian decomposition method and its application to fractional differential equations,” Applied Mathematical Modelling, vol. 37, no. 3, pp. 1590–1598, 2013.
12. J. Liu and G. Hou, “Numerical solutions of the space- and time-fractional coupled Burgers equations by generalized differential transform method,” Applied Mathematics and Computation, vol. 217, no. 16, pp. 7001–7008, 2011.
13. S. Kumar, D. Kumar, S. Abbasbandy, and M. M. Rashidi, “Analytical solution of fractional Navier-Stokes equation by using modified Laplace decomposition method,” Ain Shams Engineering Journal, vol. 5, no. 2, pp. 569–574, 2014.
14. S. Kumar, J. Singh, D. Kumar, and S. Kapoor, “New homotopy analysis transform algorithm to solve volterra integral equation,” Ain Shams Engineering Journal, vol. 5, pp. 243–246, 2014.
15. A. H. Bhrawy and M. M. Al-Shomrani, “A shifted Legendre spectral method for fractional-order multi-point boundary value problems,” Advances in Difference Equations, vol. 2012, article 8, 2012.
16. E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “A new Jacobi operational matrix: an application for solving fractional differential equations,” Applied Mathematical Modelling, vol. 36, no. 10, pp. 4931–4943, 2012.
17. A. H. Bhrawy, A. S. Alofi, and S. S. Ezz-Eldien, “A quadrature tau method for fractional differential equations with variable coefficients,” Applied Mathematics Letters, vol. 24, no. 12, pp. 2146–2152, 2011.
18. M. A. Mohamed and M. Sh. Torky, “Approximate solution of fractional nonlinear partial differential equations by the Legendre multiwavelet Galerkin method,” Journal of Applied Mathematics, vol. 2014, Article ID 192519, 12 pages, 2014.
19. B. A. Jacobs and C. Harley, “Two hybrid methods for solving two-dimensional linear time-fractional partial differential equations,” Abstract and Applied Analysis, vol. 2014, Article ID 757204, 10 pages, 2014.
20. M. Y. Ongun, D. Arslan, and R. Garrappa, “Nonstandard finite difference schemes for a fractional-order Brusselator system,” Advances in Difference Equations, vol. 2013, article 102, 2013.
21. L. Su, W. Wang, and Z. Yang, “Finite difference approximations for the fractional advection-diffusion equation,” Physics Letters A: General, Atomic and Solid State Physics, vol. 373, no. 48, pp. 4405–4408, 2009.
22. L. Su and P. Cheng, “A weighted average finite difference method for the fractional convection-diffusion equation,” Advances in Mathematical Physics, vol. 2013, Article ID 129404, 5 pages, 2013.
23. M. H. Heydari, M. R. Hooshmandasl, C. Cattani, and M. Li, “Legendre wavelets method for solving fractional population growth model in a closed system,” Mathematical Problems in Engineering, vol. 2013, Article ID 161030, 8 pages, 2013.
24. M. Mahalakshmi and G. Hariharan, “An efficient wavelet based approximation method to steady state reaction diffusion model arising in mathematical chemistry,” The Journal of Membrane Biology, vol. 247, pp. 263–271, 2013.
25. M. H. Heydari, M. R. Hooshmandasl, F. M. Maalek Ghaini, and M. Li, “Chebyshev wavelets method for solution of nonlinear fractional integrodifferential equations in a large interval,” Advances in Mathematical Physics, vol. 2013, Article ID 482083, 12 pages, 2013.
26. C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods in Fluid Dynamics, Springer, New York, NY, USA, 1988.
27. A. H. Bhrawy, E. Tohidi, and F. Soleymani, “A new Bernoulli matrix method for solving high-order linear and nonlinear Fredholm integro-differential equations with piecewise intervals,” Applied Mathematics and Computation, vol. 219, no. 2, pp. 482–497, 2012.
28. B. Fornberg, A Practical Guide to Pseudospectral Methods, vol. 1, Cambridge University Press, Cambridge, Mass, USA, 1996.
29. E. H. Doha and A. H. Bhrawy, “An efficient direct solver for multidimensional elliptic Robin boundary value problems using a Legendre spectral-Galerkin method,” Computers & Mathematics with Applications, vol. 64, no. 4, pp. 558–571, 2012.
30. R. Peyret, Spectral Methods for Incompressible Viscous Flow, vol. 148, Springer, New York, NY, USA, 2002.
31. A. H. Bhrawy, “A Jacobi-Gauss-Lobatto collocation method for solving generalized Fitzhugh-Nagumo equation with time-dependent coefficients,” Applied Mathematics and Computation, vol. 222, pp. 255–264, 2013.
32. L. N. Trefethen, Spectral Methods in MATLAB, vol. 10, SIAM, Philadelphia, Pa, USA, 2000.
33. 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.
34. A. Rababah, “Jacobi-Bernstein basis transformation,” Computational Methods in Applied Mathematics, vol. 4, no. 2, pp. 206–214, 2004.
35. E. H. Doha, A. H. Bhrawy, M. A. Abdelkawy, and R. A. Van Gorder, “Jacobi-Gauss-Lobatto collocation method for the numerical solution of $1+1$ nonlinear Schrödinger equations,” Journal of Computational Physics, vol. 261, pp. 244–255, 2014.
36. B.-Y. Guo, J. Shen, and L.-L. Wang, “Optimal spectral-Galerkin methods using generalized Jacobi polynomials,” Journal of Scientific Computing, vol. 27, no. 1–3, pp. 305–322, 2006.
37. B.-Y. Guo and Y.-G. Yi, “Generalized Jacobi rational spectral method and its applications,” Journal of Scientific Computing, vol. 43, no. 2, pp. 201–238, 2010.