TY - JOUR

T1 - A simple finite element method for boundary value problems with a Riemann–Liouville derivative

AU - Jin, Bangti

AU - Lazarov, Raytcho

AU - Lu, Xiliang

AU - Zhou, Zhi

N1 - KAUST Repository Item: Exported on 2020-10-01
Acknowledged KAUST grant number(s): KUS-C1-016-04
Acknowledgements: The authors are grateful to the anonymous referees for their insightful comments, which have led to improved presentation of the paper. The research of R. Lazarov was supported in parts by National Science Foundation Grant DMS-1016525 and also by Award No. KUS-C1-016-04, made by King Abdullah University of Science and Technology (KAUST). X. Lu is supported by National Natural Science Foundation of China Nos. 91230108 and 11471253. Z. Zhou was partially supported by National Science Foundation Grant DMS-1016525.
This publication acknowledges KAUST support, but has no KAUST affiliated authors.

PY - 2016/2

Y1 - 2016/2

N2 - © 2015 Elsevier B.V. All rights reserved. We consider a boundary value problem involving a Riemann-Liouville fractional derivative of order α∈(3/2,2) on the unit interval (0,1). The standard Galerkin finite element approximation converges slowly due to the presence of singularity term xα-$^{1}$ in the solution representation. In this work, we develop a simple technique, by transforming it into a second-order two-point boundary value problem with nonlocal low order terms, whose solution can reconstruct directly the solution to the original problem. The stability of the variational formulation, and the optimal regularity pickup of the solution are analyzed. A novel Galerkin finite element method with piecewise linear or quadratic finite elements is developed, and $^{L2}$(D) error estimates are provided. The approach is then applied to the corresponding fractional Sturm-Liouville problem, and error estimates of the eigenvalue approximations are given. Extensive numerical results fully confirm our theoretical study.

AB - © 2015 Elsevier B.V. All rights reserved. We consider a boundary value problem involving a Riemann-Liouville fractional derivative of order α∈(3/2,2) on the unit interval (0,1). The standard Galerkin finite element approximation converges slowly due to the presence of singularity term xα-$^{1}$ in the solution representation. In this work, we develop a simple technique, by transforming it into a second-order two-point boundary value problem with nonlocal low order terms, whose solution can reconstruct directly the solution to the original problem. The stability of the variational formulation, and the optimal regularity pickup of the solution are analyzed. A novel Galerkin finite element method with piecewise linear or quadratic finite elements is developed, and $^{L2}$(D) error estimates are provided. The approach is then applied to the corresponding fractional Sturm-Liouville problem, and error estimates of the eigenvalue approximations are given. Extensive numerical results fully confirm our theoretical study.

UR - http://hdl.handle.net/10754/597407

UR - https://linkinghub.elsevier.com/retrieve/pii/S0377042715001454

UR - http://www.scopus.com/inward/record.url?scp=84941260935&partnerID=8YFLogxK

U2 - 10.1016/j.cam.2015.02.058

DO - 10.1016/j.cam.2015.02.058

M3 - Article

SN - 0377-0427

VL - 293

SP - 94

EP - 111

JO - Journal of Computational and Applied Mathematics

JF - Journal of Computational and Applied Mathematics

ER -