A New Grünwald-Letnikov Derivative Derived from a Second-Order Scheme

Abstract

A novel derivation of a second-order accurate Grünwald-Letnikov-type approximation to the fractional derivative of a function is presented. This scheme is shown to be second-order accurate under certain modifications to account for poor accuracy in approximating the asymptotic behavior near the lower limit of differentiation. Some example functions are chosen and numerical results are presented to illustrate the efficacy of this new method over some other popular choices for discretizing fractional derivatives.

B. A. Jacobs 1, 2

Academic Editor:Tiecheng Xia

1, School of Computer Science and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag 3, Wits 2050, South Africa
2, DST-NRF Centre of Excellence in Mathematical and Statistical Sciences (CoE-MaSS), University of the Witwatersrand, Johannesburg, Private Bag 3, Wits 2050, South Africa

Received 25 May 2015; Accepted 16 August 2015; 26 August 2015

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.

1. Introduction

Fractional differential equations and fractional order derivatives have become increasingly popular in recent years. The use of fractional order derivatives has become abundant in modeling techniques as well as computational methods for the numerical solution of these models. Zeng et al. [1] derive a generalized fractional diffusion equation by addressing deficiencies in previously proposed models for anomalous diffusion. An application of modelling plumes moving at a fractional rate was investigated by Benson et al. [2] which showed a fractional advection-dispersion model to be appropriate as the mean squared displacement followed [figure omitted; refer to PDF] . He [3] presented a more exact model for seepage flow by relaxing the continuity assumption and went on to solve the arising equation using the variational iteration method (VIM). Metzler et al. [4] provide a review of anomalous diffusion equations presented as well as the formulation of a general diffusion case overcoming previous inconsistencies. Moving past anomalous diffusion Gafiychuk et al. [5, 6] studied a fractional reaction-diffusion system modelling an activator-inhibitor dynamic. The application of fractional order models is various and with this increased complexity requires both analytic and numerical methods to find useful solutions.

Meerschaert and Tadjeran [7] construct an implicit scheme for two-sided space-fractional partial differential equations based on two first-order approximations to the left- and right-shifted Grünwald-Letnikov fractional derivative. Scherer et al. [8] also derive an implicit scheme for the fractional heat equation in the Caputo sense, making use of the standard Grünwald-Letnikov scheme as well as the relationship between the Riemann-Liouville derivative and the Caputo derivative. There has been work conducted on comparing numerical methods with analytical methods; El-Sayed et al. [9], for example, compare a modified numerical method, wherein the authors apply a transform to homogenize the initial conditions, with the Adomian decomposition method. Another approach taken by Karatay et al. [10] derives an unconditionally stable implicit scheme based on a first-order approximation to the Caputo derivative. An interesting study is presented by Su et al. [11] which combines the first-order shifted Grünwald-Letnikov scheme with characteristic methods and a departure from the Riemann-Liouville and Caputo derivatives leads Chen et al. [12] to derive a conditionally stable explicit scheme and unconditionally stable implicit scheme for the Caputo-Riesz fractional diffusion equation that are superlinearly convergent.

Tadjeran and Meerschaert [13] implement a second-order accurate Crank-Nicolson-type finite difference scheme for the space-fractional diffusion equation, making use of the right-shifted Grünwald-Letnikov fractional derivative. Diethelm et al. [14] propose a predictor-corrector-type approach for the numerical solution of both linear and nonlinear fractional differential equations including systems of equations. Mohebbi et al. [15] study a high-order accurate numerical method for obtaining the solution to the subdiffusion equation with a nonlinear source term. Cui [16] goes on to solve the two-dimensional time-fractional diffusion equation using a high-order compact alternating direction implicit method. Meerschaert et al. [17] solve the two-dimensional space-fractional dispersion equation using a right-shifted Grünwald-Letnikov discretization since the order of the fractional derivative is above 1 but less than 2. Lynch et al. [18] implement a finite difference approach for the space-fractional dispersion problem but implement discretization switching to maximize accuracy based on the order of derivative. Momani and Odibat [19] depart from the finite-difference approach and instead compare two analytic methods for the solution of linear fractional differential equations, namely, the variational iteration method and the Adomian decomposition method. These authors go on to explore the application of these methods to nonlinear partial differential equations of fractional order in [20].

The present work derives a difference approximation in a novel way that produces a form similar to the known result presented by Oldham and Spanier in [21]. This text cites an error in the approximation of the fractional derivative of a constant and linear function of order [figure omitted; refer to PDF] that is proportional to the number of terms in the series approximation. A novel technique is coupled with this approximation to attain order [figure omitted; refer to PDF] accuracy. Although the result obtained is known, to the best of the author's knowledge this derivation is novel.

Another important technique used in the numerical application of fractional differential equations is the "short-memory" principle [22]. This technique truncates the approximation sum to avoid high computational costs. It is shown in Section 4 that the present method is highly accurate even with very few terms in the approximation sum.

The following section presents some preliminary results used in the remainder of the paper. Section 3 goes on to derive the new method of discretizing a fractional derivative with results on the performance of this technique reported in Section 4. We make some concluding remarks and discuss opportunities for further work in Section 5.

2. Preliminaries

This section presents some preliminary results that are put to use in subsequent sections. However, in the interest of brevity, we omit some useful properties but direct the interested reader to the work by Li and Deng [23]. Often for numerical simulation of fractional order problems the Grünwald-Letnikov derivative is used as a means of discretizing a fractional order derivative. The Grünwald-Letnikov discretization has been used for numerical schemes for fractional partial differential equations [13, 24, 25] and is defined as follows.

Definition 1.

The Grünwald-Letnikov derivative of order [figure omitted; refer to PDF] of a function [figure omitted; refer to PDF] is [figure omitted; refer to PDF] where [figure omitted; refer to PDF] and [figure omitted; refer to PDF] .

However the limit definition is only practically useful for finite-difference implementations; hence the following definition is used.

Definition 2.

The Grünwald-Letnikov derivative of order [figure omitted; refer to PDF] of a function [figure omitted; refer to PDF] is [22] [figure omitted; refer to PDF] where [figure omitted; refer to PDF] .

Definition 2 is used as a means of calculating the exact fractional derivative of a function. The exact solution obtained here is then used as a measure for accuracy when compared with the discrete approximations. It is prudent to mention that if [figure omitted; refer to PDF] is sufficiently smooth, [figure omitted; refer to PDF] , then the Grünwald-Letnikov derivative is equivalent to the popular Riemann-Liouville defined below.

Definition 3.

The Riemann-Liouville derivative of fractional order [figure omitted; refer to PDF] of a function [figure omitted; refer to PDF] is [22] [figure omitted; refer to PDF]

The right-shifted Grünwald-Letnikov approximation introduced by Meerschaert and Tadjeran [25], as used in [13], leads to stable numerical schemes and is defined below.

Definition 4.

The right-shifted Grünwald-Letnikov approximation of a function [figure omitted; refer to PDF] for [figure omitted; refer to PDF] is [25] [figure omitted; refer to PDF]

Oldham and Spanier obtain a second-order Grünwald-Letnikov approximation simply by noting that the Riemann sum present in the derivative would converge faster if a simple modification was made.

Definition 5.

The Grünwald-Letnikov derivative of order [figure omitted; refer to PDF] of a function [figure omitted; refer to PDF] with a central difference modification is [21] [figure omitted; refer to PDF] where [figure omitted; refer to PDF] .

The primary difference between Definition 5 and the present work is the derivation and the choice of [figure omitted; refer to PDF] , the upper limit of the Riemann sum. Although these discrepancies are slight we also address the issue of the asymptotic behavior of the fractional derivative as [figure omitted; refer to PDF] approaches zero, a difficult behavior for a difference scheme to capture.

3. Method

In much the same way the Grünwald-Letnikov derivative is generalized from the first-order backward difference scheme, as described by Podlubny in [22], we derive a new approximate fractional derivative in the form of a difference scheme which is derived from the second-order central difference scheme. Yuste and Acedo [24] state that the standard Grünwald-Letnikov derivative is accurate to [figure omitted; refer to PDF] , while we derive a scheme based on the central-difference scheme; we hope to achieve an approximation that is [figure omitted; refer to PDF] . Podlubny [26] shows, by example, that the Grünwald-Letnikov derivative of a function [figure omitted; refer to PDF] or [figure omitted; refer to PDF] is accurate to [figure omitted; refer to PDF] and hence conclude that the Grünwald-Letnikov derivative of any function that can be written as a power series will be accurate to [figure omitted; refer to PDF] . The central difference approximation to the first-order derivative is given by [figure omitted; refer to PDF] or [figure omitted; refer to PDF] Similarly the second-order derivative can be approximated by [figure omitted; refer to PDF] Following this trend we may write [figure omitted; refer to PDF] Generalizing the order of the derivative to be a positive real value [figure omitted; refer to PDF] and substituting into (10) we obtain a general formula: [figure omitted; refer to PDF] where the binomial coefficient is calculated by the use of the Gamma function that generalizes the factorial operator: [figure omitted; refer to PDF] It is well known [26, 27] that the fractional differentiation operator is a nonlocal operator. This effect is seen in the discretization being affected by every point in the function and not just a neighborhood of the point [figure omitted; refer to PDF] . This insight guides our choice of [figure omitted; refer to PDF] as the upper limit of the summation. We choose [figure omitted; refer to PDF] so that the last term in the series corresponds with the terminating point of the function [figure omitted; refer to PDF] , and in terms of fractional differential equations and fractional partial differential equations this would correspond with a boundary or initial condition: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] represents the ceiling function. Due to the singularity of a fractional derivative [figure omitted; refer to PDF] when [figure omitted; refer to PDF] it is convenient to select [figure omitted; refer to PDF] such that [figure omitted; refer to PDF] to alleviate the asymptotic behavior. However we do not want to restrict ourselves only to functions that obey this property. Hence we define [figure omitted; refer to PDF] so that [figure omitted; refer to PDF] and then [figure omitted; refer to PDF] since the fractional derivative is a linear operator. We may now exploit Property 2.2 in [23] which yields [figure omitted; refer to PDF] which then leads to [figure omitted; refer to PDF] Finally (15) becomes [figure omitted; refer to PDF] While correcting for the asymptotic behavior at [figure omitted; refer to PDF] produces sufficiently accurate results, we may go one step further and substitute the exact solution of the linear component of [figure omitted; refer to PDF] . This is done to ensure that the accuracy of the scheme is at least [figure omitted; refer to PDF] . Let [figure omitted; refer to PDF] then [figure omitted; refer to PDF]

3.1. Accuracy

As described by Oldham and Spanier in [21] the error in the above scheme for functions [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] is given by [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] , respectively. Since we make the exact substitutions for the linear components we are able to increase the accuracy of the scheme to [figure omitted; refer to PDF] for any analytic function.

Theorem 6.

The error of [figure omitted; refer to PDF] , denoted by [figure omitted; refer to PDF] is at least [figure omitted; refer to PDF] accurate for the function [figure omitted; refer to PDF] for all [figure omitted; refer to PDF] given that the exact derivative of order [figure omitted; refer to PDF] of [figure omitted; refer to PDF] is [figure omitted; refer to PDF]

Proof.

Let [figure omitted; refer to PDF] . Then (see [21], Chapter 8) (while Oldham and Spanier cite the relative error it is easy to show that the more restrictive, absolute error is also [figure omitted; refer to PDF] . More information can be found in [22], Chapter 7) [figure omitted; refer to PDF] Assume that the order of error is at least [figure omitted; refer to PDF] for [figure omitted; refer to PDF] ; that is, [figure omitted; refer to PDF] Now let [figure omitted; refer to PDF] ; then by the product rule for fractional derivatives (see [21], Chapter 5) we have [figure omitted; refer to PDF] Then using the recursive property of the Gamma function [figure omitted; refer to PDF] we have [figure omitted; refer to PDF]

Using Theorem 6 we can deduce that the accuracy of [figure omitted; refer to PDF] for any analytic function [figure omitted; refer to PDF] is [figure omitted; refer to PDF] provided that the linear and constant components of that function are dealt with accordingly.

4. Results

This section presents some numerical results for different functions. Each example shows the absolute error over the domain of the function, the consistency of each derivative as [figure omitted; refer to PDF] tends to zero, the absolute error of the truncated series, and the absolute error over a range of [figure omitted; refer to PDF] . We compare the aforementioned criteria for four different methods: the standard Grünwald-Letnikov approximation [figure omitted; refer to PDF] , the second-order Grünwald-Letnikov approximation [figure omitted; refer to PDF] , the second-order Grünwald-Letnikov approximation with constant correction described in Section 3 [figure omitted; refer to PDF] , and finally the right-shifted Grünwald-Letnikov approximation [figure omitted; refer to PDF] .

4.1. Example 1

Let [figure omitted; refer to PDF] and let the order of derivative be [figure omitted; refer to PDF] , with time step [figure omitted; refer to PDF] . Since this choice of [figure omitted; refer to PDF] satisfies [figure omitted; refer to PDF] the approximation does not suffer from poor accuracy near [figure omitted; refer to PDF] and hence the second-order Grünwald-Letnikov [figure omitted; refer to PDF] and the second-order Grünwald-Letnikov with constant correction [figure omitted; refer to PDF] approximations are similar.

4.2. Example 2

Let [figure omitted; refer to PDF] the order of derivative [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] . As mentioned by Oldham and Spanier [21] this method is robust for all [figure omitted; refer to PDF] ; this example shows highly accurate results for [figure omitted; refer to PDF] .

4.3. Example 3

Let [figure omitted; refer to PDF] the order of derivative [figure omitted; refer to PDF] , and again [figure omitted; refer to PDF] . This function suffers from the asymptotic behavior of the fractional derivative near [figure omitted; refer to PDF] and hence the constant correction method enjoys a vastly superior accuracy over the other methods.

4.4. Example 4

Let [figure omitted; refer to PDF] the order of derivative [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] .

4.5. Example 5

Let [figure omitted; refer to PDF] the order of derivative [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] .

The above results indicate that the present method enjoys a high accuracy. The nature of the fractional derivative near the lower limit of differentiation results in a high gradient. Difference schemes with fixed interval steps are unable to accurately capture very steep gradients, similar to those at the asymptote of the fractional derivative. The use of the constant correction allows one to avoid the asymptotic behavior near [figure omitted; refer to PDF] and the provable accuracy.

Figures 1, 5, 9, 13, and 17 show the increasing accuracy of the approximations as [figure omitted; refer to PDF] deviates from 0, as well as the robustness of our method to the asymptotic behavior. The consistency of the methods is evident in Figures 2, 6, 10, 14, and 18. These figures indicate that as [figure omitted; refer to PDF] tends to 0 the absolute error of the approximations also tends to 0. As mentioned in Section 1 the "short-memory" principle allows for the truncation of the approximating series in order to alleviate computational cost; our present method is highly accurate even for small values of [figure omitted; refer to PDF] as is evident in Figures 3, 7, 11, 15, and 19. The final set of figures, Figures 4, 8, 12, 16, and 20, indicate the robustness of our method to varying values of [figure omitted; refer to PDF] . It is important to note here the log scale of the graph.

Figure 1: Absolute error over time domain for four different methods.

[figure omitted; refer to PDF]

Figure 2: Absolute error over varying values of [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

Figure 3: Absolute error for increasing [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

Figure 4: Absolute error over varying values of [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

Figure 5: Absolute error over time domain for four different methods.

[figure omitted; refer to PDF]

Figure 6: Absolute error over varying values of [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

Figure 7: Absolute error for increasing [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

Figure 8: Absolute error over varying values of [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

Figure 9: Absolute error over time domain for four different methods.

[figure omitted; refer to PDF]

Figure 10: Absolute error over varying values of [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

Figure 11: Absolute error for increasing [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

Figure 12: Absolute error over varying values of [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

Figure 13: Absolute error over time domain for four different methods.

[figure omitted; refer to PDF]

Figure 14: Absolute error over varying values of [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

Figure 15: Absolute error for increasing [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

Figure 16: Absolute error over varying values of [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

Figure 17: Absolute error over time domain for four different methods.

[figure omitted; refer to PDF]

Figure 18: Absolute error over varying values of [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

Figure 19: Absolute error for increasing [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

Figure 20: Absolute error over varying values of [figure omitted; refer to PDF] for four different methods.

[figure omitted; refer to PDF]

5. Conclusion

This work has shown a novel derivation of a second-order accurate scheme for approximating the fractional order derivative of a function. To the best of the author's knowledge performing a linear transformation on the function to avoid asymptotic behavior near the origin is also novel. Section 4 provides several examples with varying parameters and behaviors to illustrate the efficacy of the present method.

Using a higher-order accurate approximation leads to the potential for higher-order accurate numerical schemes for fractional differential equations and fractional partial differential equations. Tadjeran and Meerschaert [13] implement the right-shifted Grünwald-Letnikov derivative for fractional derivatives of order [figure omitted; refer to PDF] which lead to an unconditionally stable numerical scheme. According to Oldham and Spanier [21] the present method is robust for all values of alpha . Potential extensions to this work include constructing numerical schemes centered on this approximation and the fractional order derivative for nonlinear fractional differential equations as well as time-marching schemes for time or space fractional partial differential equations.

The accuracy of the method has been proven to be at least of [figure omitted; refer to PDF] , which suggests that the accuracy should be of order [figure omitted; refer to PDF] for [figure omitted; refer to PDF] . However in these examples the error is observed to be in the region of [figure omitted; refer to PDF] . Nevertheless the mathematical proof of accuracy as well as the compelling results presented substantiates the efforts toward this field.

Acknowledgments

The author would like to acknowledge and thank Professor Charis Harley for her helpful suggestions and useful insight. The DST-NRF Centre of Excellence in Mathematical and Statistical Sciences (CoE-MaSS) is acknowledged for funding under Grant no. 94005. Thanks are due to the reviewers for their comments which helped to elevate the level of this work.

Conflict of Interests

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

References

[1] Q. Zeng, H. Li, D. Liu, "Anomalous fractional diffusion equation for transport phenomena," Communications in Nonlinear Science and Numerical Simulation , vol. 4, no. 2, pp. 99-104, 1999.

[2] D. A. Benson, S. W. Wheatcraft, M. M. Meerschaert, "Application of a fractional advection-dispersion equation," Water Resources Research , vol. 36, no. 6, pp. 1403-1412, 2000.

[3] 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. 1-2, pp. 57-68, 1998.

[4] R. Metzler, W. G. Glöckle, T. F. Nonnenmacher, "Fractional model equation for anomalous diffusion," Physica A: Statistical Mechanics and its Applications , vol. 211, no. 1, pp. 13-24, 1994.

[5] V. Gafiychuk, B. Datsko, V. Meleshko, "Mathematical modeling of time fractional reaction-diffusion systems," Journal of Computational and Applied Mathematics , vol. 220, no. 1-2, pp. 215-225, 2008.

[6] B. Y. Datsko, V. V. Gafiychuk, "Mathematical modeling of fractional reaction-diffusion systems with different order time derivatives," Journal of Mathematical Sciences , vol. 165, no. 3, pp. 392-402, 2010.

[7] M. M. Meerschaert, C. Tadjeran, "Finite difference approximations for two-sided space-fractional partial differential equations," Applied Numerical Mathematics , vol. 56, no. 1, pp. 80-90, 2006.

[8] R. Scherer, S. L. Kalla, L. Boyadjiev, B. Al-Saqabi, "Numerical treatment of fractional heat equations," Applied Numerical Mathematics , vol. 58, no. 8, pp. 1212-1223, 2008.

[9] A. M. El-Sayed, I. L. El-Kalla, E. A. Ziada, "Analytical and numerical solutions of multi-term nonlinear fractional orders differential equations," Applied Numerical Mathematics , vol. 60, no. 8, pp. 788-797, 2010.

[10] I. Karatay, R. S. Bayramoglu, A. Sahin, "Implicit difference approximation for the time fractional heat equation with the nonlocal condition," Applied Numerical Mathematics , vol. 61, no. 12, pp. 1281-1288, 2011.

[11] L. Su, W. Wang, H. Wang, "A characteristic difference method for the transient fractional convection-diffusion equations," Applied Numerical Mathematics , vol. 61, no. 8, pp. 946-960, 2011.

[12] M. Chen, W. Deng, Y. Wu, "Superlinearly convergent algorithms for the two-dimensional space-time caputo-riesz fractional diffusion equation," Applied Numerical Mathematics , vol. 70, pp. 22-41, 2013.

[13] C. Tadjeran, 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.

[14] K. Diethelm, N. J. Ford, A. D. Freed, "A predictor-corrector approach for the numerical solution of fractional differential equations," Nonlinear Dynamics , vol. 29, no. 1-4, pp. 3-22, 2002.

[15] A. Mohebbi, M. Abbaszadeh, M. Dehghan, "A high-order and unconditionally stable scheme for the modified anomalous fractional sub-diffusion equation with a nonlinear source term," Journal of Computational Physics , vol. 240, pp. 36-48, 2013.

[16] M. Cui, "Compact alternating direction implicit method for two-dimensional time fractional diffusion equation," Journal of Computational Physics , vol. 231, no. 6, pp. 2621-2633, 2012.

[17] M. M. Meerschaert, H.-P. Scheffler, C. Tadjeran, "Finite difference methods for two-dimensional fractional dispersion equation," Journal of Computational Physics , vol. 211, no. 1, pp. 249-261, 2006.

[18] V. E. Lynch, B. A. Carreras, D. del-Castillo-Negrete, K. M. Ferreira-Mejias, H. R. Hicks, "Numerical methods for the solution of partial differential equations of fractional order," Journal of Computational Physics , vol. 192, no. 2, pp. 406-421, 2003.

[19] S. Momani, Z. Odibat, "Numerical comparison of methods for solving linear differential equations of fractional order," Chaos, Solitons & Fractals , vol. 31, no. 5, pp. 1248-1255, 2007.

[20] Z. Odibat, S. Momani, "Numerical methods for nonlinear partial differential equations of fractional order," Applied Mathematical Modelling , vol. 32, no. 1, pp. 28-39, 2008.

[21] K. B. Oldham, J. Spanier The Fractional Calculus , Academic Press, New York, NY, USA, 1974.

[22] I. Podlubny Fractional Dfferential Equations , vol. 198, of Mathematics in Science and Engineering, Academic Press, 1999.

[23] C. Li, W. Deng, "Remarks on fractional derivatives," Applied Mathematics and Computation , vol. 187, no. 2, pp. 777-784, 2007.

[24] S. B. Yuste, 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.

[25] M. M. Meerschaert, 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.

[26] I. Podlubny, "Geometric and physical interpretation of fractional integration and fractional differentiation," Fractional Calculus & Applied Analysis , vol. 5, no. 4, pp. 367-386, 2002.

[27] A. Compte, R. Metzler, "The generalized Cattaneo equation for the description of anomalous transport processes," Journal of Physics A: Mathematical and General , vol. 30, no. 21, pp. 7277-7289, 1997.