Regularity theory and high order numerical methods for the (1d)-Fractional Laplacian

Abstract

This paper presents regularity results and associated high-order numerical methods for one-dimensional Fractional-Laplacian boundary-value problems. On the basis of a factorization of solutions as a product of a certain edge-singular weight $\omega$ times a "regular" unknown, a characterization of the regularity of solutions is obtained in terms of the smoothness of the corresponding right-hand sides. In particular, for right-hand sides which are analytic in a Bernstein Ellipse, analyticity in the same Bernstein Ellipse is obtained for the "regular" unknown. Moreover, a sharp Sobolev regularity result is presented which completely characterizes the co-domain of the Fractional-Laplacian operator in terms of certain weighted Sobolev spaces introduced in (Babuška and Guo, SIAM J. Numer. Anal. 2002). The present theoretical treatment relies on a full eigendecomposition for a certain weighted integral operator in terms of the Gegenbauer polynomial basis. The proposed Gegenbauer-based Nyström numerical method for the Fractional-Laplacian Dirichlet problem, further, is significantly more accurate and efficient than other algorithms considered previously. The sharp error estimates presented in this paper indicate that the proposed algorithm is spectrally accurate, with convergence rates that only depend on the smoothness of the right-hand side. In particular, convergence is exponentially fast (resp. faster than any power of the mesh-size) for analytic (resp. in nitely smooth) right-hand sides. The properties of the algorithm are illustrated with a variety of numerical results.

References added and typos corrected. To appear in Mathematics of Computation

References in corpus (3)

• From the long jump random walk to the fractional Laplacian
• Interior regularity of solutions of non-local equations in Sobolev and Nikol'skii spaces
• Gevrey regularity for integro-differential operators

Cited by in corpus (16)

• Numerical methods for nonlocal and fractional models
• Numerical Methods for Fractional Diffusion
• A short FE implementation for a 2d homogeneous Dirichlet problem of a Fractional Laplacian
• $hp$-Finite Elements for Fractional Diffusion
• Fast implicit difference schemes for time-space fractional diffusion equations with the integral fractional Laplacian
• Finite element approximations of the nonhomogeneous fractional Dirichlet problem
• Structure Probing Neural Network Deflation
• Numerical algorithm for the space-time fractional Fokker-Planck system with two internal states
• Exponential Convergence of hp FEM for the Integral Fractional Laplacian in Polygons
• A simple solver for the fractional Laplacian in multiple dimensions
• Spectral Method for the Fractional Laplacian in 2D and 3D
• Numerical approximations for fractional elliptic equations via the method of semigroups
• Evaluation of resonances: adaptivity and AAA rational approximation of randomly scalarized boundary integral resolvents
• Unbiased `walk-on-spheres' Monte Carlo methods for the fractional Laplacian
• On spectral Petrov-Galerkin method for solving fractional initial value problems in weighted Sobolev space
• Local discontinuous Galerkin method for the fractional diffusion equation with integral fractional Laplacian