Journal
JOURNAL OF SCIENTIFIC COMPUTING
Volume 81, Issue 2, Pages 1088-1110Publisher
SPRINGER/PLENUM PUBLISHERS
DOI: 10.1007/s10915-019-01052-8
Keywords
Fractional derivative; Spectral method; QZ decomposotion; Generalized Jacobi functions; Error analysis
Categories
Funding
- NSF [DMS-1620262, DMS-1720442]
- AFOSR [FA9550-16-1-0102]
Ask authors/readers for more resources
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.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available