Communications on Applied Mathematics and Computation ›› 2025, Vol. 7 ›› Issue (2): 589-605.doi: 10.1007/s42967-024-00368-z

Previous Articles     Next Articles

Meshfree Finite Difference Solution of Homogeneous Dirichlet Problems of the Fractional Laplacian

Jinye Shen1, Bowen Shi1, Weizhang Huang2   

  1. 1 School of Mathematics, Southwestern University of Finance and Economics, Chengdu 611130, Sichuan, China;
    2 Department of Mathematics, University of Kansas, Jayhawk Boulevard, Lawrence, KS 66045, USA
  • Received:2023-08-17 Revised:2024-01-05 Accepted:2024-01-06 Online:2025-06-20 Published:2025-04-21
  • Supported by:
    W.H. was supported in part by the University of Kansas General Research Fund FY23 and the Simons Foundation through Grant MP-TSM-00002397. J.S. was supported in part by the National Natural Science Foundation of China through Grant 12101509.

Abstract: A so-called grid-overlay finite difference method (GoFD) was proposed recently for the numerical solution of homogeneous Dirichlet boundary value problems (BVPs) of the fractional Laplacian on arbitrary bounded domains. It was shown to have advantages of both finite difference (FD) and finite element methods, including their efficient implementation through the fast Fourier transform (FFT) and the ability to work for complex domains and with mesh adaptation. The purpose of this work is to study GoFD in a meshfree setting, a key to which is to construct the data transfer matrix from a given point cloud to a uniform grid. Two approaches are proposed, one based on the moving least squares fitting and the other based on the Delaunay triangulation and piecewise linear interpolation. Numerical results obtained for examples with convex and concave domains and various types of point clouds are presented. They show that both approaches lead to comparable results. Moreover, the resulting meshfree GoFD converges in a similar order as GoFD with unstructured meshes and finite element approximation as the number of points in the cloud increases. Furthermore, numerical results show that the method is robust to random perturbations in the location of the points.

Key words: Fractional Laplacian, Meshfree, Finite difference (FD), Arbitrary domain, Overlay grid

CLC Number: