- About this Journal ·
- Abstracting and Indexing ·
- 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

Abstract and Applied Analysis

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

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

## Numerical Solution of a Kind of Fractional Parabolic Equations via Two Difference Schemes

^{1}Institute for Groundwater Studies, Faculty of Natural and Agricultural Sciences, University of the Free State, Bloemfontein 9300, South Africa^{2}Department of Chemical and Materials Engineering, Faculty of Engineering, King Abdulaziz University, P.O. Box 80204, Jeddah 21589, Saudi Arabia^{3}Department of Mathematics and Computer Sciences, Cankaya University, Faculty of Art and Sciences, Balgat, 06530 Ankara, Turkey^{4}Institute of Space Sciences, P.O. Box MG-23, Magurele, 76900 Bucharest, Romania

Received 5 August 2013; Accepted 23 August 2013

Academic Editor: Soheil Salahshour

Copyright © 2013 Abdon Atangana and Dumitru Baleanu. 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 kind of parabolic equation was extended to the concept of fractional calculus. The resulting equation is, however, difficult to handle analytically. Therefore, we presented the numerical solution via the explicit and the implicit schemes. We presented together the stability and convergence of this time-fractional parabolic equation with two difference schemes. The explicit and the implicit schemes in this case are stable under some conditions.

#### 1. Introduction

A parabolic partial differential equation is a type of second-order partial differential equations (PDEs), describing a wide family of problems in science including heat diffusion and ocean acoustic propagation, in physical or mathematical systems with a time variable, which behave essentially like heat diffusing through a solid [1–4]. This mathematical model is a simplified description of physical reality expressed in mathematical terms. Thus, the investigation of the exact or approximate solution helps us to understand the means of these mathematical models. In most cases, it is difficult, or infeasible, to find the analytical solution or good numerical solution of the problems. Numerical solutions or approximate analytical solutions become necessary. Numerical methods typically yield approximate solutions to the governing equation through the discretization of space and time and can relax the rigid idealized conditions of analytical models or lumped-parameter models. They can, therefore, be more realistic and flexible for simulating field conditions. Within the discredited problem domain, the variable internal properties, boundaries, and stresses of the system are approximated. One of the most important aspects of this numerical method is the study of the stability and convergence of the numerical method [5, 6].

The purpose of this work is to study the stability and the convergence of the numerical scheme of the parabolic equation of the following form: where is sufficiently given smooth function and . Here, is a sufficiently large positive constant.

#### 2. Useful Tools for the Fractional Calculus

*Definition 1 (see [7–17]). *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 space if , .

*Definition 2 (see [7–16]). *The Riemann-Liouville fractional integral operator of order , of a function , , is defined as
Properties of the operator can be found in [7–15]; we only mention the following:

*Definition 3. *The Caputo fractional-order derivative is given as follows [7–10]:

*Definition 4. *The Riemann-Liouville fractional-order derivative is given as follows [8–16]:

*Definition 5. *The Jumarie fractional-order derivative is given as follows [16]:

Lemma 6. *If , , and , , then
*

*Definition 7 (partial derivatives of fractional order [7, 8, 11, 18]). *Assume now that is a function of variables , , ??also of class on . We define partial derivative of order for in respect to the the function as follows:
where is the usual partial derivative of integer order .

#### 3. Examination of the Numerical Solution via Difference Schemes

This section is devoted to the discussion underpinning the numerical simulation of the solution above (1) via the explicit scheme [19–23] and the implicit scheme [23–30]. However, before we present the numerical schemes, we must assume that (1) has a unique and sufficiently smooth solution [23]. In addition, to present the numerical schemes, we let , , , , , and ; is the space step size, and and are grid points. We will start with the implicit scheme.

##### 3.1. Implicit Scheme for the Main Problem

It is important to recall that the finite difference approximation for the second-order spatial derivative is known as follows [23]: The discretization of the Caputo-type time-fractional-order derivative can be presented as follows: Now, substituting (9) and (10) into (1), we obtained the following expression: For ease, let Then, (11) can be rephrased as follows: It is important to inform that if , then the term of the sum of the right-hand side automatically vanished. Then, (13) can be divided as follows: where .

The above equation can be written in matrix form as follows:

###### 3.1.1. Stability of the Implicit Difference Scheme

In this subsection, we present the stability analysis of the implicit difference scheme for solving the time-fractional parabolic equation (1). To accomplish this, we let , with being the approximate solution of the main problem at the point , ; ; in addition, is the transpose of the matrix .

To analyze the stability, we exploit the Fourier method [22], and the expression of can be defined as follows: Then, the function can be expressed in Fourier series as follows: It was proven in [22, 29, 30] that Let us now examine the stability of the implicit scheme of the main problem We next assume that with being the real space wave number and . Then, by replacing (21) into (20), we obtain the following:

Lemma 8. *Assume that verified (14) and that for all , ; then, the following inequality is satisfied:
*

*Proof. *To prove this lemma, we make use of the recursive method on the natural number from (3) and (15), we have that, for all and for ,
And for , we have
Observe that then, the following is obtained:
Applying the absolute value on both sides of (20), we have
Now, let us assume that (18) is true for all ; then,
Now, applying the absolute value on both sides of (28), and making further use of the inequality theorem, we arrived at the following:
Making use of the induction hypothesis and factorizing , we obtain
But we have that
Therefore,
Thus,

Theorem 9. *The implicit difference scheme for the time-fractional parabolic equation (1) is stable providing that, for all??.*

*Proof. *From (19) and Lemma 8, we obtain
and this proves that the implicit difference scheme for the time-fractional parabolic equation (1) is stable.

*Remark 10. *It is observed that, from Theorem 9, the statement of stability of implicit difference scheme for time-fractional parabolic equation (1) depends on the evolution of the function . It follows that the stability condition can change in time advancement and space position.

###### 3.1.2. Convergence Analysis of the Implicit Difference Scheme

Assuming that is the exact solution of (1) at the point (), then, by defining, , and is the transpose of the matrix ; here, is neglected because of being equal to zero. Therefore, we have the following relation for the implicit difference scheme for the time-fractional parabolic equation (1) and from (13): Here, is the truncate or the remainder term of the approximation and has the following expression: It follows from (3), (4), and (5) that, Thus, from (36) and (37), we have the following: where , , , and are constant. The interested reader can find the error analysis of the chosen fractional derivative (the Caputo fractional derivative) in [29, 30].

Lemma 11. *One has that is true for all .*

*Proof. *Again we employ the induction method to achieve this, so that, for ,
Now, assuming that, for all , , then,
Making use of the induction hypothesis, we obtain the following:
This completes the proof.

Theorem 12. *The implicit difference scheme of the time-fractional parabolic equation (1) is convergent, and there exists a constant such that,
*

##### 3.2. Explicit Difference Scheme for Time-Fractional Parabolic Equation

It is important to recall that the finite-difference approximation for the second-order spatial derivative is known as follows: The discretization of the Caputo-type time-fractional-order derivative can be presented as follows: Now, substituting (44) and (45) into (1), we obtained the following expression: For simplicity, let Then, (37) can be rewritten as follows where and .

###### 3.2.1. Stability of the Explicit Difference Scheme of the Time-Fractional Parabolic Equation

Following the discussion presented earlier for the analysis of the implicit scheme, we obtain the following roundoff error equation from explicit scheme equation (1): where Again, we suppose that in (39) can be given in exponential form as follows: Then, replacing this form into (49), we arrive at the following expression:

Lemma 13. *Suppose that is the solution of (40) for all , , and , then, the following inequality holds:
*

*Proof. *To prove this, we make use of the recurrence technique on the natural number , and we will examine first the case for . Therefore, if , we have that
Since and , then, we have
We next assume that, for all , ; then
Now, making use of the triangular inequality, we arrive at the following inequality:
Now, making use of the induction hypothesis, we arrive at the following:
And this produces, since ,
This completes the proof.

Theorem 14. *The explicit difference scheme of time-fractional parabolic equation (1) is stable under the condition that, for all , and .*

*Proof. *It is straightforward from (11) and Lemma 13, that,
and this proves that the explicit difference scheme for the time-fractional parabolic equation (3) is stable.

###### 3.2.2. Convergence Analysis of the Explicit Scheme

Theorem 15. *The explicit scheme is convergent if, for all (; ), and and there exists a constant such that
*

Using the same method in implicit difference, the above theorem can be proven.

#### 4. Conclusions

Partial differential equations are sometimes very difficult to be solved analytically. It is, therefore, sometimes easy to solve them via numerical techniques. In this paper, we solve the time-fractional parabolic equation (1) via the implicit and explicit difference schemes. We study the stability and the convergence of implicit and explicit difference schemes.

#### References

- A. Ashyralyev, F. Dal, and Z. Pinar, “On the numerical solution of fractional hyperbolic partial differential equations,”
*Mathematical Problems in Engineering*, vol. 2009, Article ID 730465, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Ashyralyev, “A note on fractional derivatives and fractional powers of operators,”
*Journal of Mathematical Analysis and Applications*, vol. 357, no. 1, pp. 232–236, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - I. Podlubny and A. M. A. El-Sayed,
*On Two Definitions of Fractional Calculus*, Slovak Academy of Science, Institute of Experimental Physics, 1996. - A. Yakar and M. E. Koksal, “Existence results for solutions of nonlinear fractional differential equations,”
*Abstract and Applied Analysis*, vol. 2012, Article ID 267108, 12 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Atangana and J. F. Botha, “Generalized groundwater flow equation using the concept of variable order derivative,”
*Boundary Value Problems*, vol. 2013, article 53, 2013. View at Publisher · View at Google Scholar - A. Atangana and S. C. Oukouomi Noutchie, “Stability and convergence of a tme-fractional variable order Hantush equation for a deformable aquifer,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 691060, 8 pages, 2013. View at Publisher · View at Google Scholar - D. Baleanu, K. Diethelm, E. Scalas, and J. J. Trujillo,
*Fractional Calculus Models and Numerical Methods Series on Complexity, Nonlinearity and Chaos*, World Scientific, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo,
*Theory and Applications of Fractional Differential Equations*, vol. 204 of*North-Holland Mathematics Studies*, Elsevier Science B.V., Amsterdam, The Netherlands, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. B. Oldham and J. Spanier,
*The Fractional Calculus*, vol. 111 of*Mathematics in Science and Engineering*, Academic Press, New York, NY, USA, 1974. View at Zentralblatt MATH · View at MathSciNet - B. Ross, Ed.,
*Fractional Calculus and Its Applications*, vol. 457 of*Lecture Notes in Mathematics*, Springer, Berlin, Germany, 1975. View at MathSciNet - Y. Luchko and R. Gorenflo, “The initial value problem for some fractional differential equations with the Caputo derivative,”
*Preprint Series*A08–98, Freic Universitat Berlin, 1998, Fachbreich Mathematik and Informatik. View at Google Scholar - A. Atangana and E. Alabaraoye, “Solving system of fractional partial differential equations arisen in the model of HIV infection of CD4
^{+}cells and attractor one-dimensional Keller-Segel equation,”*Advances in Difference Equations*, vol. 2013, article 94, 2013. View at Publisher · View at Google Scholar - I. Podlubny,
*Fractional Differential Equations*, vol. 198 of*Mathematics in Science and Engineering*, Academic Press, San Diego, Calif, USA, 1999. View at Zentralblatt MATH · View at MathSciNet - A. Atangana and A. Secer, “A Note on fractional order derivatives and table of fractional derivatives of some special functions,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 279681, 8 pages, 2013. View at Publisher · View at Google Scholar - K. S. Miller and B. Ross,
*An Introduction to the Fractional Calculus and Fractional Differential Equations*, John Wiley & Sons, New York, NY, USA, 1993. View at MathSciNet - A. Atangana and A. Kilicman, “Analytical solutions of the space-time fractional derivative of advection dispersion equation,”
*Mathematical Problems in Engineering*, vol. 2013, Article ID 853127, 9 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - X.-J. Yang, D. Baleanu, and J. A. Tenreiro Machado, “Systems of navier-stokes equations on cantor sets,”
*Mathematical Problems in Engineering*, vol. 2013, Article ID 769724, 8 pages, 2013. View at Publisher · View at Google Scholar - F. Mainardi, “Fractional calculus: some basic problems in continuum and statistical mechanics,” in
*Fractals and Fractional Calculus in Continuum Mechanics*, vol. 378 of*CISM Courses and Lectures*, pp. 291–348, Springer, Vienna, Austria, 1997. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. M. Meerschaert and C. Tadjeran, “Finite difference approximations for fractional advection-dispersion flow equations,”
*Journal of Computational and Applied Mathematics*, vol. 172, no. 1, pp. 65–77, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - Y. Zhang, “A finite difference method for fractional partial differential equation,”
*Applied Mathematics and Computation*, vol. 215, no. 2, pp. 524–529, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - C. Tadjeran, M. M. Meerschaert, and H. P. Scheffler, “A second-order accurate numerical approximation for the fractional diffusion equation,”
*Journal of Computational Physics*, vol. 213, no. 1, pp. 205–213, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - C. M. Chen, F. Liu, I. Turner, and V. Anh, “A Fourier method for the fractional diffusion equation describing sub-diffusion,”
*Journal of Computational Physics*, vol. 227, no. 2, pp. 886–897, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - S. B. Yuste and L. Acedo, “An explicit finite difference method and a new von Neumann-type stability analysis for fractional diffusion equations,”
*SIAM Journal on Numerical Analysis*, vol. 42, no. 5, pp. 1862–1874, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - I. Podlubny, A. Chechkin, T. Skovranek, Y. Chen, and B. M. Vinagre Jara, “Matrix approach to discrete fractional calculus II: partial fractional differential equations,”
*Journal of Computational Physics*, vol. 228, no. 8, pp. 3137–3153, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - E. Hanert, “On the numerical solution of space-time fractional diffusion models,”
*Computers and Fluids*, vol. 46, no. 1, pp. 33–39, 2011. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - P. Zhuang, F. Liu, V. Anh, and I. Turner, “Numerical methods for the variable-order fractional advection-diffusion equation with a nonlinear source term,”
*SIAM Journal on Numerical Analysis*, vol. 47, no. 3, pp. 1760–1781, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Lin, F. Liu, V. Anh, and I. Turner, “Stability and convergence of a new explicit finite-difference approximation for the variable-order nonlinear fractional diffusion equation,”
*Applied Mathematics and Computation*, vol. 212, no. 2, pp. 435–445, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - J. Crank and P. Nicolson, “A practical method for numerical evaluation of solutions of partial differential equations of the heat conduction type,”
*Proceedings of the Cambridge Philosophical Society*, vol. 43, no. 1, pp. 50–67, 1947. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. Diethelm, N. J. Ford, and A. D. Freed, “Detailed error analysis for a fractional Adams method,”
*Numerical Algorithms*, vol. 36, no. 1, pp. 31–52, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - C. P. Li and C. X. Tao, “On the fractional Adams method,”
*Computers and Mathematics with Applications*, vol. 58, no. 8, pp. 1573–1588, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus