Abstract
A new spectral method is constructed for the linear and semilinear subdiffusion equations with possibly discontinuous rough initial data. The new method effectively combines several computational techniques, including the contour integral representation of the solutions, the quadrature approximation of contour integrals, the exponential integrator using the de la Vallée Poussin means of the source function, and a decomposition of the time interval geometrically refined towards the singularity of the solution and the source function. Rigorous error analysis shows that the proposed method has spectral convergence for the linear and semilinear subdiffusion equations with bounded measurable initial data and possibly singular source functions under the natural regularity of the solutions. © 2023 Society for Industrial and Applied Mathematics Publications. All rights reserved.
Original language | English |
---|---|
Pages (from-to) | 2305-2326 |
Journal | SIAM Journal on Numerical Analysis |
Volume | 61 |
Issue number | 5 |
Online published | 9 Oct 2023 |
DOIs | |
Publication status | Published - 2023 |
Research Keywords
- contour integral
- convolution quadrature
- exponential integrator
- geometric decomposition
- quadrature approximation
- semilinear subdiffusion equation
- singularity
- spectral method
- VP means
Publisher's Copyright Statement
- COPYRIGHT TERMS OF DEPOSITED FINAL PUBLISHED VERSION FILE: © 2023 Society for Industrial and Applied Mathematics.