Abstract:A grid-overlay finite difference method is proposed for the numerical approximation of the fractional Laplacian on arbitrary bounded domains. The method uses an unstructured simplicial mesh and an overlay uniform grid for the underlying domain and constructs the approximation based on a uniform-grid finite difference approximation and a data transfer from the unstructured mesh to the uniform grid. The method takes full advantages of both uniform-grid finite difference approximation in efficient matrix-vector multiplication via the fast Fourier transform and unstructured meshes for complex geometries and mesh adaptation. It is shown that its stiffness matrix is similar to a symmetric and positive definite matrix and thus invertible if the data transfer has full column rank and positive column sums. Piecewise linear interpolation is studied as a special example for the data transfer. It is proved that the full column rank and positive column sums of linear interpolation is guaranteed if the spacing of the uniform grid is smaller than or equal to a positive bound proportional to the minimum element height of the unstructured mesh. Moreover, a sparse preconditioner is proposed for the iterative solution of the resulting linear system for the homogeneous Dirichlet problem of the fractional Laplacian. Numerical examples demonstrate that the new method has similar convergence behavior as existing finite difference and finite element methods and that the sparse preconditioning is effective. Furthermore, the new method can readily be incorporated with existing mesh adaptation strategies. Numerical results obtained by combining with the so-called MMPDE moving mesh method are also presented.

Arbitrary order of convergence for Riesz fractional derivative via central difference method

Efficient Fourth-Order Weights in Kernel-Type Methods without Increasing the Stencil Size with an Application in a Time-Dependent Fractional PDE Problem

A refined convergence estimate for a fourth order finite difference numerical scheme to the Cahn-Hilliard equation

On the order of accuracy for finite difference approximations of partial differential equations using stencil composition

Sharp error estimates for spatial-temporal finite difference approximations to fractional sub-diffusion equation without regularity assumption on the exact solution

Finite Difference and Chebyshev Collocation for Time-Fractional and Riesz Space Distributed-Order Advection–Diffusion Equation with Time-Delay

Efficient Numerical Algorithms for Three-Dimensional Fractional Partial Differential Equations

A grid-overlay finite difference method for the fractional Laplacian on arbitrary bounded domains

A fast method for variable-order Caputo fractional derivative with applications to time-fractional diffusion equations

Fast time-stepping discontinuous Galerkin method for the subdiffusion equation

Improved Stencil Selection for Meshless Finite Difference Methods in 3D

-analysis of H3N3-2 -based difference method for fractional hyperbolic equations

Fast second-order implicit difference schemes for time distributed-order and Riesz space fractional diffusion-wave equations

Two fast and unconditionally stable finite difference methods for Riesz fractional diffusion equations with variable coefficients

High Order Finite Difference Discretization for Composite Grid Hierarchy and Its Applications

High-Order Central Difference Scheme for Caputo Fractional Derivative

A Corrected L1 Method for a Time-Fractional Subdiffusion Equation

An optimal frequency-domain finite-difference operator with a flexible stencil and its application in discontinuous-grid modeling

Enhancing RBF-FD Efficiency for Highly Non-Uniform Node Distributions via Adaptivity

Convergence analysis of a fast second‐order time‐stepping numerical method for two‐dimensional nonlinear time–space fractional Schrödinger equation

Efficient numerical method for multi-term time-fractional diffusion equations with Caputo-Fabrizio derivatives