Analysis of stability and convergence for L-type formulas combined with a spatial finite element method for solving subdiffusion problems. ETNA - Electronic Transactions on Numerical Analysis
Fiche du document
10 juin 2022
- handle: 10670/1.qdmsmc
- issn: 1068-9613
- doi: 10.1553/etna_vol55s568
info:eu-repo/semantics/openAccess
Sujets proches
Methods of analysis Analysis and chemistry Analytical methods Chemical analysis Analysis methods Analysis and examinationCiter ce document
Gundolf Haase et al., « Analysis of stability and convergence for L-type formulas combined with a spatial finite element method for solving subdiffusion problems. ETNA - Electronic Transactions on Numerical Analysis », Elektronisches Publikationsportal der Österreichischen Akademie der Wissenschafte, ID : 10.1553/etna_vol55s568
Métriques
Partage / Export
Résumé
A time-fractional diffusion equation with the Caputo fractional derivative of order α∈(0,1) is considered on a bounded polygonal domain. Some numerical methods are presented based on the finite element method (FEM) in space on a quasi-uniform mesh and L-type discretizations (i.e., L1, L1-2, and L1-2-3 formulas) to approximate the Caputo derivative. Stability and convergence of the L1-2-3 FEM as well as L1-2 FEM are proved rigorously. The lack of positivity of the coefficients of these formulas is the main difficulty in the analysis of the proposed methods. This has hampered the analysis of methods using finite elements mixed with L1-2 and L1-2-3 discretizations. Our proofs are based on the concept of a special kind of discrete Grönwall's inequality and the energy method. Numerical examples confirm the theoretical analysis.