- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 728736, 8 pages

http://dx.doi.org/10.1155/2014/728736

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

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

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

Academic Editor: Abdon Atangana

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., [1–7] 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 [15–19], finite difference method [20–22], and wavelet methods [23–25].

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 [26–32]). 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 ,

According to Legendre-Gauss quadratures

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*

*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*

*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*

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

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

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

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

*6. Conclusion*

*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*

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

*Acknowledgment*

*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*

*References*

- D. Baleanu, K. Diethelm, E. Scalas, and J. J. Trujillo,
*Fractional Calculus*, vol. 3 of*Series on Complexity, Nonlinearity and Chaos*, World Scientific, 2012. View at Publisher · View at Google Scholar · View at 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - 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. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Google Scholar · View at MathSciNet - 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. - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar - 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. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - 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. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Google Scholar - 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. View at Publisher · View at Google Scholar · View at MathSciNet - C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang,
*Spectral Methods in Fluid Dynamics*, Springer, New York, NY, USA, 1988. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Fornberg,
*A Practical Guide to Pseudospectral Methods*, vol. 1, Cambridge University Press, Cambridge, Mass, USA, 1996. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Peyret,
*Spectral Methods for Incompressible Viscous Flow*, vol. 148, Springer, New York, NY, USA, 2002. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at MathSciNet - L. N. Trefethen,
*Spectral Methods in MATLAB*, vol. 10, SIAM, Philadelphia, Pa, USA, 2000. View at Publisher · View at Google Scholar · View at 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 · View at Google Scholar · View at MathSciNet - A. Rababah, “Jacobi-Bernstein basis transformation,”
*Computational Methods in Applied Mathematics*, vol. 4, no. 2, pp. 206–214, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet

*
*