Monotone meshfree methods for linear elliptic equations in non-divergence form via nonlocal relaxation
We design a monotone meshfree finite difference method for linear elliptic PDEs in non-divergence form on point clouds via a nonlocal relaxation method. The key idea is a combination of a nonlocal integral relaxation of the PDE problem with a robust meshfree discretization on point clouds. Minimal positive stencils are obtained through a linear optimization procedure that automatically guarantees the stability and, therefore, the convergence of the meshfree discretization. A major theoretical contribution is the existence of consistent and positive stencils for a given point cloud geometry. We provide sufficient conditions for the existence of positive stencils by finding neighbors within an ellipse (2d) or ellipsoid (3d) surrounding each interior point, generalizing the study for Poisson’s equation by Seibold in 2008. It is well-known that wide stencils are in general needed for constructing consistent and monotone finite difference schemes for linear elliptic equations. Our result represents a significant improvement in the stencil width estimate for positive-type finite difference methods for linear elliptic equations in the near-degenerate regime (when the ellipticity constant becomes small), compared to previously known works in this area. Numerical algorithms and practical guidance are provided with an eye on the case of small ellipticity constant. Numerical results will be presented in both 2d and 3d, examining a range of ellipticity constants including the near-degenerate regime.
Tuesday, February 14, 2023
11:00AM AP&M 2402 and Zoom ID 994 0149 1091
Center for Computational Mathematics9500 Gilman Dr. #0112La Jolla, CA 92093-0112Tel: (858)534-9056