Abstract and Applied Analysis

Volume 2014 (2014), Article ID 636191, 10 pages

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

## A New Legendre Collocation Method for Solving a Two-Dimensional Fractional Diffusion Equation

^{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 11 March 2014; Accepted 5 April 2014; Published 4 May 2014

Academic Editor: Dumitru Baleanu

Copyright © 2014 A. H. Bhrawy. 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 spectral shifted Legendre Gauss-Lobatto collocation (SL-GL-C) method is developed and analyzed to solve a class of two-dimensional initial-boundary fractional diffusion equations with variable coefficients. The method depends basically on the fact that an expansion in a series of shifted Legendre polynomials , for the function and its space-fractional derivatives occurring in the partial fractional differential equation (PFDE), is assumed; the expansion coefficients are then determined by reducing the PFDE with its boundary and initial conditions to a system of ordinary differential equations (SODEs) for these coefficients. This system may be solved numerically by using the fourth-order implicit Runge-Kutta (IRK) method. This method, in contrast to common finite-difference and finite-element methods, has the exponential rate of convergence for the two spatial discretizations. Numerical examples are presented in the form of tables and graphs to make comparisons with the results obtained by other methods and with the exact solutions more easier.

#### 1. Introduction

In recent years there has been a high level of interest of employing spectral methods for numerically solving many types of integral and differential equations, due to their ease of applying them for finite and infinite domains [1–8]. The speed of convergence is one of the great advantages of spectral method. Besides, spectral methods have exponential rates of convergence; they also have high level of accuracy. From the overview of spectral approximation to differential equations, the spectral methods have been divided to four types, namely, collocation [9–11], tau [12, 13], Galerkin [14, 15], and Petrov Galerkin [16, 17] methods. The main idea of all versions of spectral methods is to express the spectral solution of the problem as a finite sum of certain basis functions (combination of orthogonal polynomials or functions) and then to choose the coefficients in order to minimize the difference between the exact and numerical solutions as possible. The spectral collocation method is a specific type of spectral methods, which is more applicable and widely used to solve almost types of differential equations [18–21].

Fractional differential equations [22–31] model many phenomena in several fields such as biology, viscoelasticity, finance, fluid mechanics, physics, chemistry, and engineering. Most fractional differential equations do not have exact solutions, so approximation and numerical techniques must be used. Here, we focus on the application of SL-GL-C scheme for the numerical solution of the time dependent fractional diffusion equation in two dimensions. This problem models phenomena exhibiting [32–34] anomalous diffusion that cannot be modeled accurately by the classical two-dimensional diffusion equations with variable coefficients. For the previous numerical methods used for time dependent fractional diffusion equation, see, for example, the papers in [35–40].

The solution of this equation is approximated as which can be expressed as a finite expansion of shifted Legendre polynomials for the space variables and then evaluate the spatial fractional partial derivatives of the finite expansion of at the shifted Legendre-Gauss-Lobatto (SL-GL) quadrature points. Substituting these approximations in the underlined PFDE provides a SODEs in time. This system may be solved by IRK scheme. This scheme is one of the suitable schemes for solving SODEs of first order. Since the IRK scheme has excellent stability properties, then it allows us to reduce computational costs. Indeed, this is the first work for employing SL-GL-C scheme to solve two-dimensional fractional diffusion equations.

This paper is organized as follows. We present few relevant properties of fractional integration in the Riemann-Liouville and shifted Legendre polynomials in the coming section. Section 3 presents the numerical solution of the two-dimensional fractional diffusion equations with initial-boundary conditions. Section 4 is devoted to solve two test examples. Finally, some concluding remarks are given in the last section.

#### 2. Preliminaries and Notation

##### 2.1. The Fractional Integration in the Riemann-Liouville Sense

There are several definitions of the fractional integration of order and not necessarily equivalent to each other; see [41, 42]. The most used definition is due to Riemann-Liouville, which is defined as One of the basic properties of the operator is The Riemann-Liouville fractional derivative of order will be denoted by . The next relation defines the Riemann-Liouville fractional derivative of order as where , and is the smallest integer greater than .

Lemma 1. *If , then
*

*2.2. Properties of Shifted Legendre Polynomials*

*The well-known Legendre polynomials are defined on the interval . Some properties about the standard Legendre polynomials will be recalled in this section. The Legendre polynomials satisfy the following Rodrigue's formula:
We recall also that is a polynomial of degree , and therefore (the th derivative of ) is
where
The Legendre polynomials satisfy the following relations:
and the orthogonality relation
where , , and is the Kronecker delta function. The Legendre-Gauss-Lobatto quadrature is used to evaluate the previous integral accurately. For any , we get
We introduce the following discrete inner product:
where () and () are used as usual nodes and the corresponding Christoffel numbers in the interval , respectively. Let the shifted Legendre polynomials be denoted by ; . Then can be obtained with the aid of the following recurrence formula:
The analytic form of the shifted Legendre polynomials of degree is given by
The Riemann-Liouville fractional integration of shifted Legendre polynomials of degree is given by
where . The orthogonality condition is
where and .*

*The function , square integrable in , may be expressed in terms of shifted Legendre polynomials as
where the coefficients are given by
In practice, only the first -terms shifted Legendre polynomials are considered. Hence can be expressed in the form
*

*3. Shifted Legendre Spectral Collocation Method*

*In this section, we derive the SL-GL-C method for solving two-dimensional fractional diffusion equation of the form
which is subject to the initial condition
and boundary conditions
*

*Denote by , and , the nodes and Christoffel numbers of the standard (shifted) Legendre-Gauss-Lobatto quadrature on the intervals , respectively. Then one can clearly deduce that
and if denotes the set of all polynomials of degree at most , then it follows for any that
*

*Here, the SL-GL-C method will be employed to transform the previous two-dimensional fractional diffusion equation into SODEs. Let us expand the dependent variables in the form
We observe from (15) that the coefficients can be approximated by
The expression (24) can be rewritten as
where . In what follows, the approximation of the fractional spatial partial derivative with respect to for can be computed as
*

*Also, the fractional spatial partial derivative of fractional order with respect to of the approximate solution is
The previous fractional derivatives can be computed at the SL-GL interpolation nodes as
where
In the proposed SL-GL-C method the residual of (19) is set to zero at of SL-GL points. Moreover, the values of , , , and are given by
Therefore, adapting (26)–(31) enables one to write (19)–(21) as a SODEs:
which is subject to the initial conditions:
where
Finally, we can arrange (32) and (33) to their matrix formulation, which can be solved by IRK scheme:
where
*

*The Runge-Kutta method can be expressed as one of powerful numerical integrations tools used for initial value SODEs of first order. The IRK method presents a subclass of the well-known family of Runge-Kutta methods and has many applications in the numerical solution of systems of ordinary differential equations.*

*4. Numerical Results*

*In this section, we present two numerical examples to show the accuracy, robustness, and applicability of the proposed method. We compare the results obtained by our method with those obtained in [43].*

*The difference between the measured value of the approximate solution and the exact solution (absolute error) is given by
where , , , , and are the space vectors, time, exact, and numerical solutions, respectively. Moreover, the maximum absolute error (MAEs) is given by
Also, we define the norm infinity as
*

*Example 1. *Consider the two-dimensional problem [43] is

*
with the diffusion coefficients
and consider the forcing function is
with the initial condition
and Dirichlet boundary conditions
The exact solution of this two-dimensional fractional diffusion equation is given by
*

*From the first look on Table 1, we see high accurate results based on MAEs at different choices of and , compared with the results given in [43]. Moreover, results listed in Table 2 confirm that the proposed method is very accurate.*

*In Figure 1, we display the numerical solution of problem (40) computed by the present method at , while we plot the space graph of absolute error of problem (40) computed by the present method at in Figure 2. As shown in Figures 3 and 4, we see the agreement of numerical and exact solutions curves.*

*Example 2. *Consider the two-dimensional fractional diffusion equation is

*
with the diffusion coefficients
and consider the forcing function is
with the initial condition
and the Dirichlet boundary conditions
The exact solution of this two-dimensional fractional diffusion equation is given by
*

*More accurate results for maximum absolute and absolute errors of (46) have been obtained in Tables 3 and 4, respectively, with values of parameters listed in their captions, the approximate solution at and is plotted in Figure 5. Moreover, we plot the absolute error for this problem at and in Figure 6. Finally, the - and -directions curves of absolute error at and are displayed in Figures 7 and 8.*

*5. Concluding Remarks and Future Work*

*The SL-GL-C method was investigated successfully in spatial discretizations to get accurate approximate solution of the two-dimensional fractional diffusion equations with variable coefficients. This problem was transformed into SODEs in time variable which greatly simplifies the problem. The IRK scheme was then applied to the resulting system. From the numerical experiments, the obtained results demonstrated the effectiveness and the high accuracy of SL-GL-C method for solving the mentioned problem.*

*The technique can be extended to more sophisticated problems in two- and three-dimensional spaces. In principle, this method may be extended to related problems in mathematical physics. It is possible to use other orthogonal polynomials, say, Chebyshev polynomials or Jacobi polynomials, to solve the mentioned problem in this paper. Furthermore, the proposed spectral method might be developed by considering the Legendre pseudospectral approximation in both temporal and spatial discretizations. We should note that, as a numerical method, we are restricted to solving problem over a finite domain. Also, the pseudospectral approximation might be employed based on generalized Laguerre or modified generalized Laguerre polynomials [44] to solve similar problems in a semi-infinite spatial intervals.*

*Conflict of Interests*

*The author declares that there is no conflict of interests regarding the publication of this paper.*

*Acknowledgments*

*This paper was funded by the Deanship of Scientific Research DSR, King Abdulaziz University, Jeddah, Grant no. 130-53-D1435. The author, therefore, acknowledges with thanks DSR for the technical and financial support.*

*References*

- C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang,
*Spectral Methods: Fundamentals in Single Domains*, Springer, New York, NY, USA, 2006. - H. Schamel and K. Elsaesser, “The application of the spectral method to nonlinear wave propagation,”
*Journal of Computational Physics*, vol. 22, no. 4, pp. 501–516, 1976. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - W. M. Abd-Elhameed, E. H. Doha, and Y. H. Youssri, “Efficient spectral-Petrov-Galerkin methods for third- and fifth-order differential equations using general parameters generalized Jacobi polynomials,”
*Quaestiones Mathematicae*, vol. 36, pp. 15–38, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - E. H. Doha, A. H. Bhrawy, M. A. Abdelkawy, and R. M. Hafez, “A Jacobi collocation approximation for nonlinear coupled viscous Burgers' equation,”
*Central European Journal of Physics*, vol. 12, pp. 111–122, 2014. View at Publisher · View at Google Scholar - H. Adibi and A. M. Rismani, “On using a modified Legendre-spectral method for solving singular IVPs of Lane-Emden type,”
*Computers and Mathematics with Applications*, vol. 60, no. 7, pp. 2126–2130, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - 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 Schrodinger equations,”
*Journal of Computational Physics*, vol. 261, pp. 244–255, 2014. View at Publisher · View at Google Scholar - E. H. Doha, A. H. Bhrawy, D. Baleanu, and R. M. Hafez, “A new Jacobi rational-Gauss collocation method for numerical solution of generalized Pantograph equations,”
*Applied Numerical Mathematics*, vol. 77, pp. 43–54, 2014. View at Publisher · View at Google Scholar - F. M. Mahfouz, “Numerical simulation of free convection within an eccentric annulus filled with micropolar fluid using spectral method,”
*Applied Mathematics and Computation*, vol. 219, pp. 5397–5409, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J. Ma, B.-W. Li, and J. R. Howell, “Thermal radiation heat transfer in one- and two-dimensional enclosures using the spectral collocation method with full spectrum k-distribution model,”
*International Journal of Heat and Mass Transfer*, vol. 71, pp. 35–43, 2014. View at Publisher · View at Google Scholar - X. Ma and C. Huang, “Spectral collocation method for linear fractional integro-differential equations,”
*Applied Mathematical Modelling*, vol. 38, pp. 1434–1448, 2014. View at Publisher · View at Google Scholar - W. M. Abd-Elhameed, E. H. Doha, and Y. H. Youssri, “New wavelets collocation method for solving second-order multipoint boundary value problems using Chebyshev polynomials of third and fourth,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 542839, 9 pages, 2013. View at Publisher · View at Google Scholar - S. R. Lau and R. H. Price, “Sparse spectral-tau method for the three-dimensional helically reduced wave equation on two-center domains,”
*Journal of Computational Physics*, vol. 231, pp. 7695–7714, 2012. View at Publisher · View at Google Scholar - F. Ghoreishi and S. Yazdani, “An extension of the spectral Tau method for numerical solution of multi-order fractional differential equations with convergence analysis,”
*Computers and Mathematics with Applications*, vol. 61, no. 1, pp. 30–43, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - 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 and Mathematics with Applications*, vol. 64, pp. 558–571, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - T. Boaca and I. Boaca, “Spectral galerkin method in the study of mass transfer in laminar and turbulent flows,”
*Computer Aided Chemical Engineering*, vol. 24, pp. 99–104, 2007. View at Publisher · View at Google Scholar · View at Scopus - E. H. Doha, A. H. Bhrawy, and R. M. Hafez, “A Jacobi-Jacobi dual-Petrov-Galerkin method for third- and fifth-order differential equations,”
*Mathematical and Computer Modelling*, vol. 53, no. 9-10, pp. 1820–1832, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - W. M. Abd-Elhameed, E. H. Doha, and M. A. Bassuony, “Two Legendre-Dual-Petrov-Galerkin algorithms for solving the integrated forms of high odd-order boundary value problems,”
*The Scientific World Journal*, vol. 2013, Article ID 309264, 11 pages, 2013. View at Publisher · View at Google Scholar - L. Wang, Y. Ma, and Z. Meng, “Haar wavelet method for solving fractional partial differential equations numerically,”
*Applied Mathematics and Computation*, vol. 227, pp. 66–76, 2014. View at Publisher · View at Google Scholar - A. Golbabai and M. Javidi, “A numerical solution for non-classical parabolic problem based on Chebyshev spectral collocation method,”
*Applied Mathematics and Computation*, vol. 190, no. 1, pp. 179–185, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - 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 - E. H. Doha, A. H. Bhrawy, and R. M. Hafez, “On shifted Jacobi spectral method for high-order multi-point boundary value problems,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 17, no. 10, pp. 3802–3810, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - R. Garrappa and M. Popolizio, “On the use of matrix functions for fractional partial differential equations,”
*Mathematics and Computers in Simulation*, vol. 81, no. 5, pp. 1045–1056, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - A. A. Pedas and E. Tamme, “Numerical solution of nonlinear fractional differential equations by spline collocation methods,”
*Journal of Computational and Applied Mathematics*, vol. 255, pp. 216–230, 2014. View at Publisher · View at Google Scholar - F. Gao, X. Lee, F. Fei, H. Tong, Y. Deng, and H. Zhao, “Identification time-delayed fractional order chaos with functional extrema model via differential evolution,”
*Expert Systems With Applications*, vol. 41, pp. 1601–1608, 2014. View at Publisher · View at Google Scholar - Y. Zhao, D. F. Cheng, and X. J. Yang, “Approximation solutions for local fractional Schrödinger equation in the one-dimensional Cantorian system,”
*Advances in Mathematical Physics*, vol. 2013, Article ID 291386, 5 pages, 2013. View at Publisher · View at Google Scholar - 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 - J. W. Kirchner, X. Feng, and C. Neal, “Frail chemistry and its implications for contaminant transport in catchments,”
*Nature*, vol. 403, no. 6769, pp. 524–526, 2000. View at Publisher · View at Google Scholar · View at Scopus - A. Pedas and E. Tamme, “Piecewise polynomial collocation for linear boundary value problems of fractional differential equations,”
*Journal of Computational and Applied Mathematics*, vol. 236, no. 13, pp. 3349–3359, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - A. Ahmadian, M. Suleiman, S. Salahshour, and D. Baleanu, “A Jacobi operational matrix for solving a fuzzy linear fractional differential equation,”
*Advances in Difference Equations*, vol. 2013, article 104, 29 pages, 2013. View at Publisher · View at Google Scholar - 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 - R. Hilfer,
*Applications of Fractional Calculus in Physics*, Word Scientific, Singapore, 2000. - E. Kotomin and V. Kuzovkov,
*Modern Aspects of Diffusion-Controlled Reactions: Cooperative Phenomena in Bimolecular Processes*, Comprehensive Chemical Kinetics, Elsevier, 1996. - D. Baleanu, K. Diethelm, E. Scalas, and J. J. Trujillo,
*Fractional Calculus Models and Numerical Methods*, Series on Complexity, Nonlinearity and Chaos, World Scientific, Singapore, 2012. - R. S. Cantrell and C. Cosner,
*Spatial Ecology via Reaction Diffusion Equations*, Wiley, 2004. - H. Wang and N. Du, “Fast solution methods for space-fractional diffusion equations,”
*Journal of Computational and Applied Mathematics*, vol. 255, pp. 376–383, 2014. View at Publisher · View at Google Scholar - A. H. Bhrawy and D. Baleanu, “A spectral Legendre-Gauss-Lobatto collocation method for a space-fractional advection diffusion equations with variable coefficients,”
*Reports on Mathematical Physics*, vol. 72, pp. 219–233, 2013. View at Publisher · View at Google Scholar - E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “Numerical approximations for fractional diffusion equations via a Chebyshev spectral-tau method,”
*Central European Journal of Physics*, vol. 11, pp. 1494–1503, 2013. View at Publisher · View at Google Scholar - F. Liu, P. Zhuang, I. Turner, K. Burrage, and V. Anh, “A new fractional finite volume method for solving the fractional diffusion equation,”
*Applied Mathematical Modelling*, 2013. View at Publisher · View at Google Scholar - S. B. Yuste and J. Quintana-Murillo, “A finite difference method with non-uniform timesteps for fractional diffusion equations,”
*Computer Physics Communications*, vol. 183, pp. 2594–2600, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - K. Wang and H. Wang, “A fast characteristic finite difference method for fractional advection-diffusion equations,”
*Advances in Water Resources*, vol. 34, no. 7, pp. 810–816, 2011. View at Publisher · View at Google Scholar · View at Scopus - K. Miller and B. Ross,
*An Introduction to the Fractional Calaulus and Fractional Differential Equations*, John Wiley & Sons, New York, NY, USA, 1993. - I. Podluny,
*Fractional Differential Equations*, Academic Press, San Diego, Calif, USA, 1999. - C. Tadjeran and M. M. Meerschaert, “A second-order accurate numerical method for the two-dimensional fractional diffusion equation,”
*Journal of Computational Physics*, vol. 220, no. 2, pp. 813–823, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - 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

*
*