Abstract:In this paper, we present how high-order accurate solutions to elliptic partial differential equations can be achieved in arbitrary spatial domains using radial basis function-generated finite differences (RBF-FD) on unfitted node sets (i.e., not adjusted to the domain boundary). In this novel method, we only collocate on nodes interior to the domain boundary and enforce boundary conditions as constraints by means of Lagrange multipliers. This combination enables full geometric flexibility near boundaries without compromising the high-order accuracy of the RBF-FD method. The high-order accuracy and robustness of two formulations of this approach are illustrated by numerical experiments.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to use the Radial Basis Function - Generated Finite Difference method (RBF - FD) combined with the Lagrange multiplier method on a spatial domain of arbitrary shape to achieve high - precision numerical solutions of elliptic partial differential equations (PDEs) while maintaining the meshfree characteristic. Specifically, the paper proposes two new RBF - FD methods based on Lagrange multipliers (RBF - FD - LM1 and RBF - FD - LM2), and these methods can still maintain high - order accuracy without adjusting the node set to fit the boundary. ### Main Problems and Solutions 1. **Problem Background**: - Elliptic partial differential equations play a crucial role in mathematical models that describe physical systems and processes. - Traditional numerical methods usually require the node set to be aligned with the boundary, which limits the geometric flexibility and is more complicated when dealing with moving boundaries. - Mesh - free methods (such as RBF - FD) provide flexibility, but there are challenges in maintaining high - order accuracy. 2. **Proposed Methods**: - **RBF - FD - LM1**: It is achieved by minimizing the sum of squared residuals on internal nodes and using the Lagrange multiplier method with the boundary conditions as constraints. - **RBF - FD - LM2**: High - order accuracy is achieved by adding boundary constraints to the discrete Laplace operator, similar to the methods in optimization theory. 3. **Advantages of the Methods**: - **Geometric Flexibility**: There is no need to adjust the node set to fit the boundary, which simplifies the generation of the node set. - **High - Order Accuracy**: Combined with the Lagrange multiplier method, it can handle complex boundary conditions without sacrificing high - order accuracy. - **Computational Efficiency**: It is applicable to static and dynamic boundary problems and simplifies the handling of moving boundaries. ### Numerical Experiments The paper verifies the effectiveness and robustness of the proposed methods through several test problems: 1. **Test Problem 1**: Poisson's equation on the unit disk with Dirichlet boundary conditions. 2. **Test Problem 2**: Poisson's equation on a butterfly - shaped domain with Dirichlet and Neumann boundary conditions. 3. **Test Problem 3**: Poisson's equation on the unit sphere with Dirichlet boundary conditions. 4. **Test Problem 4**: Poisson's equation on the unit sphere with Dirichlet and Neumann boundary conditions. ### Conclusion The paper shows how to combine the Lagrange multiplier method and RBF - FD to achieve high - precision numerical solutions of elliptic partial differential equations without adjusting the node set to fit the boundary. The two proposed RBF - FD - LM methods show high - order accuracy in both two - dimensional and three - dimensional problems and have good robustness when dealing with complex boundary conditions. Future research can further explore the application of these methods in time - dependent problems.

High-order numerical method for solving elliptic partial differential equations on unfitted node sets

An unfitted RBF-FD method in a least-squares setting for elliptic PDEs on complex geometries

Adaptive residual refinement in an RBF finite difference scheme for 2D time-dependent problems

Highly localized RBF Lagrange functions for finite difference methods on spheres

Solve Partial Differential Equations by Two or More Radial Basis Functions

A RBF partition of unity collocation method based on finite difference for initial-boundary value problems

Radial boundary elements method, a new approach on using radial basis functions to solve partial differential equations, efficiently

Enhancing RBF-FD Efficiency for Highly Non-Uniform Node Distributions via Adaptivity

A novel Hermite RBF-based differential quadrature method for solving two-dimensional variable-order time fractional advection-diffusion equation with Neumann boundary condition

On solving elliptic boundary value problems using a meshless method with radial polynomials

A novel local Hermite radial basis function‐based differential quadrature method for solving two‐dimensional variable‐order time fractional advection–diffusion equation with Neumann boundary conditions

An RBF-FD polynomial method based on polyharmonic splines for the Navier-Stokes equations: Comparisons on different node layouts

A fully implicit method using nodal radial basis functions to solve the linear advection equation

High order unfitted finite element discretizations for explicit boundary representations

Solving multiscale elliptic problems by sparse radial basis function neural networks

Solve Partial Differential Equations by Meshless Subdomains Method Combined with RBFs

Adaptive radial basis function–generated finite differences method for contact problems

Numerical investigation of high-dimensional option pricing PDEs by utilizing a hybrid radial basis function - finite difference procedure

A stabilized radial basis-finite difference (RBF-FD) method with hybrid kernels

Randomized Radial Basis Function Neural Network for Solving Multiscale Elliptic Equations

ON SPECIAL FINITE DIFFERENCE APPROXIMATIONS FOR SOLVING SECOND ORDER DIFFERENTIAL EQUATIONS