An accurate and efficient space-time Galerkin spectral method for the subdiffusion equation  

作　　者：Wei Zeng Chuanju Xu 

机构地区：[1]School of Mathematical Sciences and Fujian Provincial Key Laboratory of Mathematical Modeling and High Performance Scientific Computing,Xiamen University,Xiamen 361005,China

出　　处：《Science China Mathematics》2024年第10期2387-2408,共22页中国科学（数学）（英文版）

基　　金：supported by National Natural Science Foundation of China (Grant No. 11971408)。

摘　　要：In this paper, we design and analyze a space-time spectral method for the subdiffusion equation.Here, we are facing two difficulties. The first is that the solutions of this equation are usually singular near the initial time. Consequently, traditional high-order numerical methods in time are inefficient. The second obstacle is that the resulting system of the space-time spectral approach is usually large and time-consuming to solve. We aim at overcoming the first difficulty by proposing a novel approach in time, which is based on variable transformation techniques. Suitable ψ-fractional Sobolev spaces and a new variational framework are introduced to establish the well-posedness of the associated variational problem. This allows us to construct our space-time spectral method using a combination of temporal generalized Jacobi polynomials(GJPs) and spatial Legendre polynomials. For the second difficulty, we propose a fast algorithm to effectively solve the resulting linear system. The fast algorithm makes use of a matrix diagonalization in space and QZ decomposition in time. Our analysis and numerical experiments show that the proposed method is exponentially convergent with respect to the polynomial degrees in both space and time directions, even though the exact solution has very limited regularity.

关 键 词：subdiffusion equations variable transformation Ψ-Sobolev spaces WELL-POSEDNESS space-time Galerkin spectral method error estimate fast algorithm 

分 类 号：O241.82[理学—计算数学]