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.

An Optimal Adaptive Grid Method Based on L1 Scheme for a Nonlinear Caputo Fractional Differential Equation

A Compact Difference Scheme on Graded Meshes for the Nonlinear Fractional Integro-differential Equation with Non-smooth Solutions

α -robust analysis of fast and novel two-grid FEM with nonuniform L1 scheme for semilinear time-fractional variable coefficient diffusion equations

A Corrected L1 Method for a Time-Fractional Subdiffusion Equation

A fully discrete local discontinuous Galerkin method based on generalized numerical fluxes to variable-order time-fractional reaction-diffusion problem with the Caputo fractional derivative

The Stability and Convergence Analysis of Finite Difference Methods for the Fractional Neutron Diffusion Equation

Mathematical Analysis and a Second-Order Compact Scheme for Nonlinear Caputo–Hadamard Fractional Sub-diffusion Equations

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

Finite difference methods for nonlinear fractional differential equation with ψ -Caputo derivative

Two-grid methods for nonlinear time fractional diffusion equations by L 1 -Galerkin FEM

Numerical Methods for Caputo–Hadamard Fractional Differential Equations with Graded and Non-Uniform Meshes

A High-Order Accurate Numerical Scheme for the Caputo Derivative with Applications to Fractional Diffusion Problems

Fully Finite Element Adaptive Algebraic Multigrid Method for Time-Space Caputo-Riesz Fractional Diffusion Equations

A Robust and higher order numerical technique for a time-fractional equation with nonlocal condition

Two‐grid finite element method on grade meshes for time‐fractional nonlinear Schrödinger equation

An analysis of the L1 Scheme for the subdiffusion equation with nonsmooth data

An efficient numerical method for a time-fractional diffusion equation

A high-order adaptive numerical algorithm for fractional diffusion wave equation on non-uniform meshes

Mathematical analysis and efficient finite element approximation for variable‐order time‐fractional reaction–diffusion equation with nonsingular kernel

Efficient Numerical Schemes for Fractional Sub-Diffusion Equation with the Spatially Variable Coefficient

A new robust compact difference scheme on graded meshes for the time-fractional nonlinear Kuramoto–Sivashinsky equation