A unified Petrov-Galerkin spectral method for fractional PDEs

Existing numerical methods for fractional PDEs suffer from low accuracy and inefficiency in dealing with three-dimensional problems or with long-time integrations. We develop a unified and spectrally accurate Petrov–Galerkin (PG) spectral method for a weak formulation of the general linear Fractional Partial Differential Equations (FPDEs) of the form 0Dt2τu+∑j=1dcj[ajDxj2μju]+γu=f, where 2τ2τ, μj∈(0,1)μj∈(0,1), in a (1+d1+d)-dimensional space–time domain subject to Dirichlet initial and boundary conditions. We perform the stability analysis (in 1-D) and the corresponding convergence study of the scheme (in multi-D). The unified PG spectral method applies to the entire family of linear hyperbolichyperbolic-, parabolicparabolic- and ellipticelliptic-like equations. We develop the PG method based on a new spectral theory for fractional Sturm–Liouville problems (FSLPs), recently introduced in Zayernouri and Karniadakis (2013). Specifically, we employ the eigenfunctions of the FSLP of first kind (FSLP-I), called Jacobi poly-fractonomials, as temporal/spatial bases. Next, we construct a different space for test functions from poly-fractonomial eigenfunctions of the FSLP of second kind (FSLP-II). Besides the high-order spatial accuracy of the PG method, we demonstrate its efficiency and spectral accuracy in time-integration schemes for solving time-dependent FPDEs as well, rather than employing algebraically accurate traditional methods, especially when 2τ=12τ=1. Finally, we formulate a general fast linear solver based on the eigenpairs of the corresponding temporal and spatial mass matrices with respect to the stiffness matrices, which reduces the computational cost drastically. We demonstrate that this framework can reduce to hyperbolic FPDEs such as time- and space-fractional advection (TSFA), parabolic FPDEs such as time- and space-fractional diffusion (TSFD) model, and elliptic FPDEs such as fractional Helmholtz/Poisson equations with the same ease and cost. Several numerical tests confirm the efficiency and spectral convergence of the unified PG spectral method for the aforementioned families of FPDEs. Moreover, we demonstrate the computational efficiency of the new approach in higher-dimensions e.g., (1+3), (1+5) and (1+9)-dimensional problems.

[1]  J. P. Roop.  Variational Solution of the Fractional Advection Dispersion Equation , 2004 .

[2]  George E. Karniadakis,et al.  Fractional Spectral Collocation Method , 2014, SIAM J. Sci. Comput..

[3]  J. M. Sanz-Serna,et al.  A numerical method for a partial integro-differential equation , 1988 .

[4]  Cécile Piret,et al.  A radial basis functions method for fractional diffusion equations , 2013, J. Comput. Phys..

[5]  I. Podlubny.  Fractional differential equations , 1998 .

[6]  Giuseppina Autuori,et al.  Elliptic problems involving the fractional Laplacian in RN , 2013 .

[7]  R. Magin.  Fractional Calculus in Bioengineering , 2006 .

[8]  J. Guermond,et al.  Theory and practice of finite elements , 2004 .

[9]  G. Karniadakis,et al.  Spectral/hp Element Methods for CFD , 1999 .

[10]  Ali H Bhrawy,et al.  A shifted Legendre spectral method for fractional-order multi-point boundary value problems , 2012 .

[11]  George E. Karniadakis,et al.  Fractional spectral collocation methods for linear and nonlinear variable order FPDEs , 2015, J. Comput. Phys..

[12]  Jan S. Hesthaven,et al.  Stable multi-domain spectral penalty methods for fractional partial differential equations , 2014, J. Comput. Phys..

[13]  J. P. Roop.  Computational aspects of FEM approximation of fractional advection dispersion equations on bounded domains in R 2 , 2006 .

[14]  B. Henry,et al.  The accuracy and stability of an implicit solution method for the fractional diffusion equation , 2005 .

[15]  C. Lubich.  Discretized fractional calculus , 1986 .

[16]  Jakob Keller,et al.  Propagation of simple non-linear waves in gas filled tubes with friction , 1981 .

[17]  George E. Karniadakis,et al.  Exponentially accurate spectral and spectral element methods for fractional ODEs , 2014, J. Comput. Phys..

[18]  E. A. Rawashdeh,et al.  Numerical solution of fractional integro-differential equations by collocation method , 2006, Appl. Math. Comput..

[19]  G. Fix,et al.  Least squares finite-element solution of a fractional order two-point boundary value problem , 2004 .

[20]  Eid H. Doha,et al.  A Chebyshev spectral method based on operational matrix for initial and boundary value problems of fractional order , 2011, Comput. Math. Appl..

[21]  O. C. Zienkiewicz,et al.  The Finite Element Method: Its Basis and Fundamentals , 2005 .

[22]  F. Mainardi.  Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models , 2010 .

[23]  Luise Blank,et al.  Numerical Treatment of Differential Equations of Fractional Order , 1996 .

[24]  W. Chester,et al.  Resonant oscillations in closed tubes , 1964, Journal of Fluid Mechanics.

[25]  Dumitru Baleanu,et al.  Two Efficient Generalized Laguerre Spectral Algorithms for Fractional Initial Value Problems , 2013 .

[26]  Nobumasa Sugimoto.  Burgers equation with a fractional derivative; hereditary effects on nonlinear acoustic waves , 1991, Journal of Fluid Mechanics.

[27]  J. Klafter,et al.  The random walk's guide to anomalous diffusion: a fractional dynamics approach , 2000 .

[28]  George E. Karniadakis,et al.  Fractional Sturm-Liouville eigen-problems: Theory and numerical approximation , 2013, J. Comput. Phys..

[29]  Igor M. Sokolov,et al.  Physics of Fractal Operators , 2003 .

[30]  H. Kreiss,et al.  Time-Dependent Problems and Difference Methods , 1996 .

[31]  A. H. Bhrawy,et al.  A New Legendre Spectral Galerkin and Pseudo-Spectral Approximations for Fractional Initial Value Problems , 2013 .

[32]  D. Benson,et al.  Application of a fractional advection‐dispersion equation , 2000 .

[33]  Zhi‐zhong Sun,et al.  A fully discrete difference scheme for a diffusion-wave system , 2006 .

[34]  S. Wearne,et al.  Fractional Reaction-Diffusion , 2000 .

[35]  C. Lubich,et al.  On the Stability of Linear Multistep Methods for Volterra Convolution Equations , 1983 .

[36]  Jan S. Hesthaven,et al.  Spectral Methods for Time-Dependent Problems: Contents , 2007 .

[37]  T. Kakutani,et al.  ‘Generalized Burgers' equation’ for nonlinear viscoelastic waves , 1985 .

[38]  M. Khader.  On the numerical solutions for the fractional diffusion equation , 2011 .

[39]  X. Li,et al.  Existence and Uniqueness of the Weak Solution of the Space-Time Fractional Diffusion Equation and a Spectral Method Approximation , 2010 .

[40]  J. Bouchaud,et al.  Anomalous diffusion in disordered media: Statistical mechanisms, models and physical applications , 1990 .

[41]  Junying Cao,et al.  A high order schema for the numerical solution of the fractional ordinary differential equations , 2013, J. Comput. Phys..

[42]  George E. Karniadakis,et al.  Discontinuous Spectral Element Methods for Time- and Space-Fractional Advection Equations , 2014, SIAM J. Sci. Comput..

[43]  K. Miller,et al.  An Introduction to the Fractional Calculus and Fractional Differential Equations , 1993 .

[44]  Chuanju Xu,et al.  Finite difference/spectral approximations for the time-fractional diffusion equation , 2007, J. Comput. Phys..

[45]  MOHSEN ZAYERNOURI,et al.  Spectral and Discontinuous Spectral Element Methods for Fractional Delay Equations , 2014, SIAM J. Sci. Comput..

[46]  Saeid Abbasbandy,et al.  An Adaptive Pseudospectral Method for Fractional Order Boundary Value Problems , 2012 .

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

[48]  Mathias Kluge.  Parallel Scientific Computing In C And Mpi A Seamless Approach To Parallel Algorithms And Their Implementation , 2016 .

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