An Effective Finite Element Method with Shifted Fractional Powers Bases for Fractional Boundary Value Problems

In this paper, an effective finite element method with shifted fractional powers bases is developed for fractional convection diffusion equations involving a Riemann–Liouville derivative of order α∈(3/2,2)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \in (3/2,2)$$\end{document}. A Petrov-Galerkin variational formulation is constructed on the domain H~α-1(Ω)×H~1(Ω)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tilde{H}^{\alpha -1}(\Omega )\times \tilde{H}^{1}(\Omega )$$\end{document}, based on which the finite element approximation scheme is developed by employing shifted fractional power functions and continuous piecewise polynomials of degree up to m(m∈N+)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m~(m\in \mathbb {N}^+)$$\end{document} for trial and test finite element spaces, respectively. The approximation property of trial finite element space and inf\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\inf $$\end{document}-sup\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sup $$\end{document} condition for discrete variational form are derived, which enables us to derive the error estimates in L2(Ω)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L^2(\Omega )$$\end{document} and Hα-1(Ω)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H^{\alpha -1}(\Omega )$$\end{document} norms. Numerical examples are included to verify the theoretical findings and demonstrate an actual convergence rate of order α-1+m\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha -1+m$$\end{document}, where m equals to 1 or 2.