On Spectral Petrov–Galerkin Method for Solving Time-Fractional Diffusion Equations in Weighted Sobolev Space
- 1College of Sciences, Nanjing Agricultural University, Nanjing 210095, China
- 2School of Mathematics, Southeast University, Nanjing 210096, China
Abstract
Time-fractional diffusion equations with a Caputo derivative of order α∈(0,1) provide effective models for anomalous diffusion processes with memory effects, while the weakly singular fractional kernel and rough source terms generally lead to limited temporal regularity. In this paper, we develop a fully discrete scheme that combines a second-order central difference approximation in space with a spectral Petrov–Galerkin method in time. The spatial semidiscretization leads to fractional initial value problems with reaction parameters given by the eigenvalues of the spatial discretization matrix. Unlike the standard reaction-free case, the presence of these reaction parameters complicates the weighted regularity analysis. An ultraweak formulation and a bootstrapping argument are employed to establish the temporal regularity of the semidiscrete solution in weighted Sobolev spaces for source terms with regularity index r≥−α, including rough inputs. The dependence of the resulting regularity estimates on the reaction coefficients is explicitly characterized for the subsequent fully discrete error analysis. Under suitable spatial regularity of the source term, the fully discrete error is bounded by C(h2+N−m), where m=α+min{r,2α+1−ϵ} and ϵ>0 is an arbitrarily small positive constant. No coupling condition between the spatial mesh size h and the temporal polynomial degree N is required. A fast iterative solver is also developed using fast polynomial transforms and matrix–vector products. Numerical experiments support the theoretical convergence rates. For the tested scalar FIVP, the fast solver requires less computational time than the direct solver at the larger temporal polynomial degrees.
Related
- Highly Accurate Galerkin Approach for Modified Anomalous Time-Fractional Sub-Diffusion Equations Based on Fibonacci Coefficient Polynomials
- Galerkin finite element method for time fractional diffusion-wave equations with constant time delay
- Preconditioning for the time-dependent Stokes equations: stage-weighted Radau IIA and weighted discontinuous Galerkin time-stepping
- A weighted polygonal discontinuous Galerkin method for hierarchically coupled reaction-diffusion systems with cubic interactions
- A stable and efficient Galerkin spectral method with asymmetrically-weighted Hermite functions for the Vlasov-Poisson system
BibTeX
@article{li2026spectral,
title = {{On Spectral Petrov–Galerkin Method for Solving Time-Fractional Diffusion Equations in Weighted Sobolev Space}},
author = {Shengyue Li and Wanrong Cao and Zhaopeng Hao},
journal = {Fractal and Fractional},
year = {2026},
month = oct,
doi = {10.3390/fractalfract10110743},
url = {https://www.mdpi.com/2504-3110/10/11/743},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - On Spectral Petrov–Galerkin Method for Solving Time-Fractional Diffusion Equations in Weighted Sobolev Space
AU - Shengyue Li
AU - Wanrong Cao
AU - Zhaopeng Hao
PY - 2026
DA - 2026/10/25
JO - Fractal and Fractional
DO - 10.3390/fractalfract10110743
UR - https://www.mdpi.com/2504-3110/10/11/743
AB - Time-fractional diffusion equations with a Caputo derivative of order α∈(0,1) provide effective models for anomalous diffusion processes with memory effects, while the weakly singular fractional kernel and rough source terms generally lead to limited temporal regularity. In this paper, we develop a fully discrete scheme that combines a second-order central difference approximation in space with a spectral Petrov–Galerkin method in time. The spatial semidiscretization leads to fractional initial value problems with reaction parameters given by the eigenvalues of the spatial discretization matrix. Unlike the standard reaction-free case, the presence of these reaction parameters complicates the weighted regularity analysis. An ultraweak formulation and a bootstrapping argument are employed to establish the temporal regularity of the semidiscrete solution in weighted Sobolev spaces for source terms with regularity index r≥−α, including rough inputs. The dependence of the resulting regularity estimates on the reaction coefficients is explicitly characterized for the subsequent fully discrete error analysis. Under suitable spatial regularity of the source term, the fully discrete error is bounded by C(h2+N−m), where m=α+min{r,2α+1−ϵ} and ϵ>0 is an arbitrarily small positive constant. No coupling condition between the spatial mesh size h and the temporal polynomial degree N is required. A fast iterative solver is also developed using fast polynomial transforms and matrix–vector products. Numerical experiments support the theoretical convergence rates. For the tested scalar FIVP, the fast solver requires less computational time than the direct solver at the larger temporal polynomial degrees.
SN - 2504-3110
ER -