Abstract
© 2014 Society for Industrial and Applied Mathematics We consider an initial boundary value problem for a one-dimensional fractional-order parabolic equation with a space fractional derivative of Riemann-Liouville type and order α ∈ (1, 2). We study a spatial semidiscrete scheme using the standard Galerkin finite element method with piecewise linear finite elements, as well as fully discrete schemes based on the backward Euler method and the Crank-Nicolson method. Error estimates in the L2(D)- and Hα/2 (D)-norm are derived for the semidiscrete scheme and in the L2(D)-norm for the fully discrete schemes. These estimates cover both smooth and nonsmooth initial data and are expressed directly in terms of the smoothness of the initial data. Extensive numerical results are presented to illustrate the theoretical results.
| Original language | English (US) |
|---|---|
| Pages (from-to) | 2272-2294 |
| Number of pages | 23 |
| Journal | SIAM Journal on Numerical Analysis |
| Volume | 52 |
| Issue number | 5 |
| DOIs | |
| State | Published - Jan 2014 |
| Externally published | Yes |
Bibliographical note
KAUST Repository Item: Exported on 2020-10-01Acknowledged KAUST grant number(s): KUS-C1-016-04
Acknowledgements: This work was supported by award KUS-C1-016-04, made by King Abdullah University of Science and Technology (KAUST).This author's work was supported by NSF grant DMS-1319052.The second author's work was supported in part by NSF grant DMS-1016525, and the third author's work was supported by NSF grant DMS-1216551.
This publication acknowledges KAUST support, but has no KAUST affiliated authors.