Exponentially accurate spectral and spectral element methods for fractional ODEs

• Published in: Journal of computational physics Vol. 257; pp. 460 - 480
• Main Authors: Zayernouri, Mohsen, Karniadakis, George Em
• Format: Journal Article
• Language: English
• Published: Elsevier Inc 01.01.2014
• Subjects: Asymptotic properties; Basis functions; Dirichlet problem; Discontinuous spectral element method; Distributed memory; Eigenfunctions; Exponential convergence; Jacobi polyfractonomials; Mathematical analysis; Mathematical models; Petrov–Galerkin spectral methods; Spectral methods; Texts; Exponential convergence; Jacobi polyfractonomials; Petrov–Galerkin spectral methods; Discontinuous spectral element method
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2013.09.039

Current discretizations of fractional differential equations lead to numerical solutions of low order of accuracy. Here, we present different methods for fractional ODEs that lead to exponentially fast decay of the error. First, we develop a Petrov–Galerkin (PG) spectral method for Fractional Initial-Value Problems (FIVPs) of the form Dtν0u(t)=f(t) and Fractional Final-Value Problems (FFVPs) DTνtu(t)=g(t), where ν∈(0,1), subject to Dirichlet initial/final conditions. These schemes are developed based on a new spectral theory for fractional Sturm–Liouville problems (FSLPs), which has been recently developed in [1]. Specifically, we obtain solutions to FIVPs and FFVPs in terms of the new fractional (non-polynomial) basis functions, called Jacobi polyfractonomials, which are the eigenfunctions of the FSLP of first kind (FSLP-I). Correspondingly, we employ another space of test functions as the span of polyfractonomial eigenfunctions of the FSLP of second kind (FSLP-II). Subsequently, we develop a Discontinuous Spectral Method (DSM) of Petrov–Galerkin sense for the aforementioned FIVPs and FFVPs, where the basis functions do not satisfy the initial/final conditions. Finally, we extend the DSM scheme to a Discontinuous Spectral Element Method (DSEM) for efficient longer time-integration and adaptive refinement. In these discontinuous schemes, we employ the asymptotic eigensolutions to FSLP-I & -II, which are of Jacobi polynomial forms, as basis and test functions. Our numerical tests confirm the exponential/algebraic convergence, respectively, in p- and h-refinements, for various test cases with integer- and fractional-order solutions.