An Efficient Space–Time Method for Time Fractional Diffusion Equation

A space–time Petrov–Galerkin spectral method for time fractional diffusion equations is developed in this paper. The Petrov–Galerkin method is used to simplify the computation of stiffness matrix but leads to full non-symmetric mass matrix. However, the matrix decomposition method based on eigen-decomposition is numerically unstable for non-symmetric linear systems. A QZ decomposition is adopted instead of eigen-decomposition. The QZ decomposition has essentially the same computational complexity as the eigen-decomposition but is numerically stable. Moreover, the enriched Petrov–Galerkin method is developed to resolve the weak singularity at the initial time. We also carry out the error analysis for the proposed methods and present ample numerical results to validate the accuracy and robustness of our numerical schemes.

[1]  Bangti Jin,et al.  A Petrov-Galerkin Finite Element Method for Fractional Convection-Diffusion Equations , 2015, SIAM J. Numer. Anal..

[2]  B. Guo,et al.  Spectral Methods and Their Applications , 1998 .

[3]  I. Podlubny.  Fractional differential equations : an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications , 1999 .

[4]  Bangti Jin,et al.  Two Fully Discrete Schemes for Fractional Diffusion and Diffusion-Wave Equations with Nonsmooth Data , 2016, SIAM J. Sci. Comput..

[5]  Weihua Deng,et al.  Polynomial spectral collocation method for space fractional advection–diffusion equation , 2012, 1212.3410.

[6]  Jan S. Hesthaven,et al.  A multi-domain spectral method for time-fractional differential equations , 2015, J. Comput. Phys..

[7]  Jie Shen,et al.  Fourierization of the Legendre--Galerkin method and a new space--time spectral method , 2007 .

[8]  Bangti Jin,et al.  Space-Time Petrov–Galerkin FEM for Fractional Diffusion Problems , 2017, Comput. Methods Appl. Math..

[9]  Jie Shen,et al.  Enriched Spectral Methods and Applications to Problems with Weakly Singular Solutions , 2018, Journal of Scientific Computing.

[10]  William Rundell,et al.  Variational formulation of problems involving fractional order differential operators , 2013, Math. Comput..

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

[12]  Dale B. Haidvogel,et al.  The Accurate Solution of Poisson's Equation by Expansion in Chebyshev Polynomials , 1979 .

[13]  Jie Shen,et al.  Spectral and High-Order Methods with Applications , 2006 .

[14]  Jie Shen,et al.  Spectral Methods: Algorithms, Analysis and Applications , 2011 .

[15]  Chuanju Xu,et al.  A fractional spectral method with applications to some singular problems , 2016, Adv. Comput. Math..

[16]  Li-Lian Wang Shuhuang Xiang Wenjie Liu.  A New Spectral Method Using Nonstandard Singular Basis Functions for Time-Fractional Diferential Equations , 2019 .

[17]  T. A. Zang,et al.  Spectral Methods: Fundamentals in Single Domains , 2010 .

[18]  Zhongqiang Zhang,et al.  Fast difference schemes for solving high-dimensional time-fractional subdiffusion equations , 2016, J. Comput. Phys..

[19]  Masahiro Yamamoto,et al.  Initial value/boundary value problems for fractional diffusion-wave equations and applications to some inverse problems , 2011 .

[20]  Fawang Liu,et al.  The Use of Finite Difference/Element Approaches for Solving the Time-Fractional Subdiffusion Equation , 2013, SIAM J. Sci. Comput..

[21]  Changtao Sheng and Jie Shen,et al.  A Space-Time Petrov-Galerkin Spectral Method for Time Fractional Diffusion Equation , 2018, Numerical Mathematics: Theory, Methods and Applications.

[22]  M. Stynes,et al.  An analysis of the Grünwald–Letnikov scheme for initial-value problems with weakly singular solutions , 2019, Applied Numerical Mathematics.

[23]  Jie Shen.  Generalized Jacobi Spectral-Galerkin Method for Nonlinear Volterra Integral Equations with Weakly Singular Kernels , 2015 .

[24]  Chi-Wang Shu,et al.  Foreword by the Editor-in-Chief , 2019, Communication on Applied Mathematics and Computation.

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

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

[27]  Xiaoping Xie,et al.  Analysis of a Time-Stepping Scheme for Time Fractional Diffusion Problems with Nonsmooth Data , 2018, SIAM J. Numer. Anal..

[28]  Zhimin Zhang,et al.  Optimal Fractional Integration Preconditioning and Error Analysis of Fractional Collocation Method Using Nodal Generalized Jacobi Functions , 2016, SIAM J. Numer. Anal..

[29]  Jie Shen,et al.  Generalized Jacobi functions and their applications to fractional differential equations , 2014, Math. Comput..

[30]  Zhongqiang Zhang,et al.  A Generalized Spectral Collocation Method with Tunable Accuracy for Variable-Order Fractional Differential Equations , 2015, SIAM J. Sci. Comput..

[31]  K. Diethelm.  The Analysis of Fractional Differential Equations: An Application-Oriented Exposition Using Differential Operators of Caputo Type , 2010 .

[32]  Chuanju Xu,et al.  Müntz Spectral Methods for the Time-Fractional Diffusion Equation , 2017, Comput. Methods Appl. Math..

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

[34]  Zhongqiang Zhang,et al.  Optimal Error Estimates of Spectral Petrov-Galerkin and Collocation Methods for Initial Value Problems of Fractional Differential Equations , 2015, SIAM J. Numer. Anal..

[35]  Masahiro Yamamoto,et al.  Time-fractional diffusion equation in the fractional Sobolev spaces , 2015 .

[36]  Mark Ainsworth,et al.  A unified Petrov-Galerkin spectral method for fractional PDEs , 2015 .

[37]  Bangti Jin,et al.  Numerical Analysis of Nonlinear Subdiffusion Equations , 2017, SIAM J. Numer. Anal..

[38]  Weihua Deng,et al.  Time Discretization of a Tempered Fractional Feynman-Kac Equation with Measure Data , 2018, SIAM J. Numer. Anal..

[39]  John R. Rice,et al.  Direct solution of partial difference equations by tensor product methods , 1964 .

[40]  Bangti Jin,et al.  Error Analysis of a Finite Element Method for the Space-Fractional Parabolic Equation , 2014, SIAM J. Numer. Anal..

[41]  Zhi-Zhong Sun,et al.  Alternating direction implicit schemes for the two-dimensional fractional sub-diffusion equation , 2011, J. Comput. Phys..