A finite volume scheme with preconditioned Lanczos method for two-dimensional space-fractional reaction–diffusion equations

Abstract Fractional differential equations have been increasingly used as a powerful tool to model the non-locality and spatial heterogeneity inherent in many real-world problems. However, a constant challenge faced by researchers in this area is the high computational expense of obtaining numerical solutions of these fractional models, owing to the non-local nature of fractional derivatives. In this paper, we introduce a finite volume scheme with preconditioned Lanczos method as an attractive and high-efficiency approach for solving two-dimensional space-fractional reaction–diffusion equations. The computational heart of this approach is the efficient computation of a matrix-function-vector product f ( A ) b , where A is the matrix representation of the Laplacian obtained from the finite volume method and is non-symmetric. A key aspect of our proposed approach is that the popular Lanczos method for symmetric matrices is applied to this non-symmetric problem, after a suitable transformation. Furthermore, the convergence of the Lanczos method is greatly improved by incorporating a preconditioner. Our approach is show-cased by solving the fractional Fisher equation including a validation of the solution and an analysis of the behaviour of the model.

[1]  Manuel Duarte Ortigueira,et al.  Fractional Calculus for Scientists and Engineers , 2011, Lecture Notes in Electrical Engineering.

[2]  Fawang Liu,et al.  Novel numerical methods for time-space fractional reaction diffusion equations in two dimensions , 2011 .

[3]  K. Burrage,et al.  Restarted GMRES preconditioned by deflation , 1996 .

[4]  I. Sokolov,et al.  Anomalous transport : foundations and applications , 2008 .

[5]  J. A. Tenreiro Machado,et al.  New Trends in Nanotechnology and Fractional Calculus Applications , 2010 .

[6]  J. Crawford,et al.  A mass balance based numerical method for the fractional advection‐dispersion equation: Theory and application , 2005 .

[7]  Qianqian Yang,et al.  Efficient solution of two-sided nonlinear space-fractional diffusion equations using fast Poisson preconditioners , 2013, J. Comput. Phys..

[8]  Gene H. Golub,et al.  Adaptively Preconditioned GMRES Algorithms , 1998, SIAM J. Sci. Comput..

[9]  Qianqian Yang,et al.  A banded preconditioner for the two-sided, nonlinear space-fractional diffusion equation , 2013, Comput. Math. Appl..

[10]  H. V. D. Vorst,et al.  An iterative solution method for solving f ( A ) x = b , using Krylov subspace information obtained for the symmetric positive definite matrix A , 1987 .

[11]  Fawang Liu,et al.  Computationally efficient methods for solving time-variable-order time-space fractional reaction-diffusion equation , 2012 .

[12]  M. Meerschaert,et al.  Stochastic Models for Fractional Calculus , 2011 .

[13]  Fawang Liu,et al.  Stability and convergence of a finite volume method for the space fractional advection-dispersion equation , 2014, J. Comput. Appl. Math..

[14]  YuanTong Gu,et al.  AN ADVANCED MESHLESS METHOD FOR TIME FRACTIONAL DIFFUSION EQUATION , 2011 .

[15]  Tao Lin,et al.  On the Accuracy of the Finite Volume Element Method Based on Piecewise Linear Polynomials , 2001, SIAM J. Numer. Anal..

[16]  Hans Engler,et al.  On the Speed of Spread for Fractional Reaction-Diffusion Equations , 2009, 0908.0024.

[17]  O. Agrawal,et al.  Advances in Fractional Calculus , 2007 .

[18]  Fawang Liu,et al.  Numerical solution of the space fractional Fokker-Planck equation , 2004 .

[19]  V. Ervin,et al.  Variational formulation for the stationary fractional advection dispersion equation , 2006 .

[20]  A. Luo,et al.  Fractional Dynamics and Control , 2011 .

[21]  M. T. Cicero.  FRACTIONAL CALCULUS AND WAVES IN LINEAR VISCOELASTICITY , 2012 .

[22]  Fawang Liu,et al.  An implicit RBF meshless approach for time fractional diffusion equations , 2011 .

[23]  I. Turner,et al.  Numerical methods for fractional partial differential equations with Riesz space fractional derivatives , 2010 .

[24]  I. Turner,et al.  Numerical Approximation of a Fractional-In-Space Diffusion Equation, I , 2005 .

[25]  K. Burrage,et al.  Fourier spectral methods for fractional-in-space reaction-diffusion equations , 2014 .

[26]  Y. Saad,et al.  Numerical Methods for Large Eigenvalue Problems , 2011 .

[27]  Ralf Metzler,et al.  Fractional dynamics : recent advances , 2011 .

[28]  Nicholas Hale,et al.  An Efficient Implicit FEM Scheme for Fractional-in-Space Reaction-Diffusion Equations , 2012, SIAM J. Sci. Comput..

[29]  Willem Hundsdorfer,et al.  Stability of implicit-explicit linear multistep methods , 1997 .

[30]  M. Meerschaert,et al.  Finite difference approximations for fractional advection-dispersion flow equations , 2004 .

[31]  Fawang Liu,et al.  Novel Numerical Methods for Solving the Time-Space Fractional Diffusion Equation in Two Dimensions , 2011, SIAM J. Sci. Comput..

[32]  Vo Anh,et al.  A Numerical Solution Using an Adaptively Preconditioned Lanczos Method for a Class of Linear Systems Related with the Fractional Poisson Equation , 2008 .

[33]  Xianjuan Li,et al.  A Space-Time Spectral Method for the Time Fractional Diffusion Equation , 2009, SIAM J. Numer. Anal..

[34]  Fawang Liu,et al.  Novel numerical methods for solving the time-space fractional diffusion equation in 2D , 2011 .

[35]  M. Meerschaert,et al.  Finite difference approximations for two-sided space-fractional partial differential equations , 2006 .