Minimal positive stencils in meshfree finite difference methods for linear elliptic equations in non-divergence form

by   Qihao Ye, et al.

We design a monotone meshfree finite difference method for linear elliptic equations in the non-divergence form on point clouds via a nonlocal relaxation method. Nonlocal approximations of linear elliptic equations are first introduced to which a meshfree finite difference method applies. Minimal positive stencils are obtained through a local l_1-type optimization procedure that automatically guarantees the stability and, therefore, the convergence of the meshfree discretization for linear elliptic equations. The key to the success of the method relies on the existence of 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 study improves the known theoretical results on the existence of positive stencils for linear elliptic equations when the ellipticity constant becomes small. Numerical algorithms and practical guidance are provided with an eye on the case of small ellipticity constant. We present numerical results in 2d and 3d at the end.


page 1

page 2

page 3

page 4


Convergent Finite Difference Methods for Fully Nonlinear Elliptic Equations in Three Dimensions

We introduce a generalized finite difference method for solving a large ...

Numerical analysis of optimal control problems governed by fourth-order linear elliptic equations using the Hessian discretisation method

This paper focusses on the optimal control problems governed by fourth-o...

Scheduled Relaxation Jacobi schemes for non-elliptic partial differential equations

The Scheduled Relaxation Jacobi (SRJ) method is a linear solver algorith...

Error Bounds for a Least Squares Meshless Finite Difference Method on Closed Manifolds

We present an error bound for a least squares version of the kernel base...

MGM: A meshfree geometric multilevel method for systems arising from elliptic equations on point cloud surfaces

We develop a new meshfree geometric multilevel (MGM) method for solving ...

On the Reduction in Accuracy of Finite Difference Schemes on Manifolds without Boundary

We investigate error bounds for numerical solutions of divergence struct...

Please sign up or login with your details

Forgot password? Click here to reset