An -robust finite element method for a multi-term time-fractional diffusion problem
Chaobao Huang,Martin Stynes,Hu Chen
DOI: https://doi.org/10.1016/j.cam.2020.113334
IF: 2.872
2021-01-01
Journal of Computational and Applied Mathematics
Abstract:A time-fractional initial-boundary problem is considered on a bounded spatial domain Omega subset of R-d, where d is an element of {1, 2, 3} and Omega is convex or smooth. The differential equation is Sigma(l)(i=1) q(i)D(t)(alpha i) u(x, t) - Delta u(x, t) = f, where each D-t(alpha i) is a Caputo derivative with 0 < alpha(l) < ... < alpha(1) < 1 and the q(i) are positive constants. A new error analysis is given for the numerical method (L1 scheme in time, finite elements in space) of Huang and Stynes (2020). The error bounds that are derived here remain valid (i.e., are "alpha-robust") if alpha(1) -> 1(-), unlike the bounds in Huang and Stynes (2020), which blow up as alpha(1) -> 1(-). Thus these new error bounds are a more natural and desirable result, since the unknown solution is well behaved as alpha(1) -> 1(-). (C) 2020 Elsevier B.V. All rights reserved.