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

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.

Keywords

Fractional Laplacian / Meshfree / Finite difference (FD) / Arbitrary domain / Overlay grid / 65N06 / 35R11

References

[1] Acosta G, Bersetche FM, Borthagaray JP. A short FE implementation for a 2D homogeneous Dirichlet problem of a fractional Laplacian. Comput. Math. Appl.. 2017, 74: 784-816

[2] Acosta G, Borthagaray JP. A fractional Laplace equation: regularity of solutions and finite element approximations. SIAM J. Numer. Anal.. 2017, 55: 472-495

[3] Ainsworth M, Glusa C. Aspects of an adaptive finite element method for the fractional Laplacian: a priori and a posteriori error estimates, efficient implementation and multigrid solver. Comput. Methods Appl. Mech. Eng.. 2017, 327: 4-35

[4] Ainsworth M, Glusa C. Towards an efficient finite element method for the integral fractional Laplacian on polygonal domains. Contemporary Computational Mathematics—a Celebration of the 80th Birthday of Ian Sloan. 2018, Cham, Springer: 1757

[5] Antil H, Brown T, Khatri R, Onwunta A, Verma D, Warma M. Chapter 3 — Optimal control, numerics, and applications of fractional PDEs. Handbook of Numerical Analysis. 2022, Amsterdam, Elsevier: 8711423

[6] Antil H, Dondl P, Striet L. Approximation of integral fractional Laplacian and fractional PDEs via sinc-basis. SIAM J. Sci. Comput.. 2021, 43: A2897-A2922

[7] Bonito A, Lei W, Pasciak JE. Numerical approximation of the integral fractional Laplacian. Numer. Math.. 2019, 142: 235-278

[8] Burkardt J, Wu Y, Zhang Y. A unified meshfree pseudospectral method for solving both classical and fractional PDEs. SIAM J. Sci. Comput.. 2021, 43: A1389-A1411

[9] Chew LP. Constrained Delaunay triangulations. Algorithmica. 1989, 4: 97-108 Computational geometry (Waterloo, ON) (1987)

[10] Du N, Sun H-W, Wang H. A preconditioned fast finite difference scheme for space-fractional diffusion equations in convex domains. Comput. Appl. Math.. 2019, 38: 14

[11] Du Q, Ju L, Lu J. A discontinuous Galerkin method for one-dimensional time-dependent nonlocal diffusion problems. Math. Comput.. 2019, 88: 123-147

[12] Duo S, van Wyk HW, Zhang Y. A novel and accurate finite difference method for the fractional Laplacian and the fractional Poisson problem. J. Comput. Phys.. 2018, 355: 233-252

[13] Dyda B, Kuznetsov A, Kwaśnicki M. Fractional Laplace operator and Meijer G-function. Constr. Approx.. 2017, 45: 427-448

[14] Faustmann M, Karkulik M, Melenk JM. Local convergence of the FEM for the integral fractional Laplacian. SIAM J. Numer. Anal.. 2022, 60: 1055-1082

[15] Hao Z, Zhang Z, Du R. Fractional centered difference scheme for high-dimensional integral fractional Laplacian. J. Comput. Phys.. 2021, 424 109851

[16] Huang W, Russell RD. Adaptive Moving Mesh Methods. Applied Mathematical Sciences Series. 2011, New York, Springer 174

[17] Huang, W., Shen, J.: A grid-overlay finite difference method for the fractional Laplacian on arbitrary bounded domains. SIAM J. Sci. Comput. (to appear). http://arxiv.org/abs/arXiv:2307.14437 (2023)

[18] Huang Y, Oberman A. Numerical methods for the fractional Laplacian: a finite difference-quadrature approach. SIAM J. Numer. Anal.. 2014, 52: 3056-3084

[19] Huang, Y., Oberman, A.: Finite difference methods for fractional Laplacians. arXiv:1611.00164 (2016)

[20] Li H, Liu R, Wang L-L. Efficient Hermite spectral-Galerkin methods for nonlocal diffusion equations in unbounded domains. Numer. Math. Theory Methods Appl.. 2022, 15: 1009-1040

[21] Lischke A, Pang G, Gulian M, Song F, Glusa C, Zheng X, Mao Z, Cai W, Meerschaert MM, Ainsworth M, Karniadakis GE. What is the fractional Laplacian? A comparative review with new results. J. Comput. Phys.. 2020, 404 109009

[22] Liu GR. Meshfree Methods: Moving Beyond the Finite Element Method. 20102Boca Raton, CRC Press

[23] Minden V, Ying L. A simple solver for the fractional Laplacian in multiple dimensions. SIAM J. Sci. Comput.. 2020, 42: A878-A900

[24] Ortigueira MD. Riesz potential operators and inverses via fractional centred derivatives. Int. J. Math. Math. Sci.. 2006, 2006: 48391

[25] Ortigueira MD. Fractional central differences and derivatives. J. Vib. Control. 2008, 14: 1255-1266

[26] Pang G, Chen W, Fu Z. Space-fractional advection-dispersion equations by the Kansa method. J. Comput. Phys.. 2015, 293: 280-296

[27] Pang H-K, Sun H-W. Multigrid method for fractional diffusion equations. J. Comput. Phys.. 2012, 231: 693-703

[28] Shewchuk JR. General-dimensional constrained Delaunay and constrained regular triangulations. I. Combinatorial properties. Discrete Comput. Geom.. 2008, 39: 580-637

[29] Somasekhar M, Vivek S, Malagi KS, Ramesh V, Deshpande SM. Adaptive cloud refinement (ACR)-adaptation in meshless framework. Commun. Comput. Phys.. 2012, 11: 1372-1385

[30] Song F, Xu C, Karniadakis GE. Computing fractional Laplacians on complex-geometry domains: algorithms and simulations. SIAM J. Sci. Comput.. 2017, 39: A1320-A1344

[31] Suchde P, Jacquemin T, Davydov O. Point cloud generation for meshfree methods: an overview. Arch. Comput. Methods Eng.. 2023, 30: 889-915

[32] Sun J, Nie D, Deng W. Algorithm implementation and numerical analysis for the two-dimensional tempered fractional Laplacian. BIT. 2021, 61: 1421-1452

[33] Tian X, Du Q. Analysis and comparison of different approximations to nonlocal diffusion and linear peridynamic equations. SIAM J. Numer. Anal.. 2013, 51: 3458-3482

[34] Trobec R, Kosec G. Parallel Scientific Computing: Theory, Algorithms, and Applications of Mesh Based and Meshless Methods. 2015, Cham, Springer

[35] Wang H, Basu TS. A fast finite difference method for two-dimensional space-fractional diffusion equations. SIAM J. Sci. Comput.. 2012, 34: A2444-A2458

Funding

National Natural Science Foundation of China(12101509)

University of Kansas General Research(FY23)

Simons Foundation(MP-TSM-00002397)

RIGHTS & PERMISSIONS

Shanghai University