Abstract: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. The key idea is a novel 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 local $l_1$-type optimization procedure that automatically guarantees the stability and, therefore, the convergence of the meshfree discretization for linear elliptic equations. 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 (Comput Methods Appl Mech Eng 198(3-4):592-601, 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. At the end, we present numerical results for the performance of our method in both 2d and 3d, examining a range of ellipticity constants including the near-degenerate regime.
What problem does this paper attempt to address?
This paper attempts to solve the numerical solution problem of linear elliptic partial differential equations (PDEs) in non - divergence form on point cloud data. Specifically, the article focuses on how to design a monotonic mesh - free finite - difference method to handle such equations and ensure the stability and convergence of the numerical solution.
### Problem Background
The linear elliptic equation in non - divergence form can be expressed as:
\[
\begin{cases}
- Lu(x) := -\sum_{i,j = 1}^d a_{ij}(x)\partial_{ij}u(x)=f(x), & x\in\Omega\\
u(x) = g(x), & x\in\partial\Omega
\end{cases}
\]
where $\Omega\subset\mathbb{R}^d$ is a bounded region, $A(x)=(a_{ij}(x))_{d\times d}$ is a symmetric positive - definite matrix, satisfying the uniform elliptic condition:
\[
\lambda|\xi|^2\leq\xi^T A(x)\xi\leq\Lambda|\xi|^2,\quad\forall\xi\in\mathbb{R}^d,\forall x\in\Omega
\]
Here, $\lambda$ and $\Lambda$ are positive constants, and $\rho:=\frac{\lambda}{\Lambda}\leq1$.
### Main Challenges
1. **Requirement for Wide Templates**: In order to construct consistent and monotonic finite - difference schemes, wide templates are usually required, especially when the ellipticity constant is small.
2. **Non - Smooth Coefficient Matrix**: When the coefficient matrix $A(x)$ is not smooth or continuous, it cannot be transformed into divergence form, which makes the traditional variational methods no longer applicable.
### Solution
The author proposes a mesh - free finite - difference method based on the non - local relaxation method to solve the problem through the following steps:
1. **Non - Local Integral Approximation**: Define a family of non - local operators $\{L_\delta\}_{\delta > 0}$ to approximate the original elliptic operator $L$, where $\delta$ is a non - local parameter.
2. **Elliptical Neighborhood Search**: Define an elliptical neighborhood around each interior point and find the minimum positive template through a local $l_1$ - like optimization method.
3. **Theoretical Guarantee**: Prove that there exist consistent and positive templates under the given point - cloud geometric conditions and give an estimate of the elliptical search area.
### Main Contributions
- **Improved Template Width Estimation**: Compared with existing works, this paper significantly reduces the required template width in the near - degenerate case ($\rho\ll1$).
- **Guidance for Practical Applications**: Provide numerical algorithms and implementation guidance, especially suitable for cases with small ellipticity constants.
- **Two - Dimensional and Three - Dimensional Verification**: Show the numerical results of this method in two - dimensional and three - dimensional cases, verifying its effectiveness and practicality.
Through these improvements, this method not only improves the accuracy and efficiency of the numerical solution but also provides more stringent theoretical guarantees.