Abstract:We generate Gaussian radial function based higher order compact RBF-FD formulas for some differential operators. Analytical expressions for weights associated to first and second derivative formulas (up to order 10) and 2D-Laplacian formulas (up to order 6) are derived. Then these weights are used to obtain analytical expression for local truncation errors. The weights are obtained by symbolic computation of a linear system in Mathematica. Often such linear systems are not directly amenable to symbolic computation. We make use of symmetry of formula stencil along with Taylor series expansions for performing the computation. In the flat limit, the formulas converge to their respective order polynomial based compact FD formulas. We validate the formulas with standard test functions and demonstrate improvement in approximation accuracy with respect to corresponding order multiquadric based compact RBF-FD formulas and compact FD schemes. We also compute optimal value of shape parameter for each formula.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop higher - order compact radial basis function finite - difference (RBF - FD) formulas, especially those based on the Gaussian radial basis function (Gaussian RBF). Specifically, the author aims to: 1. **Derive higher - order RBF - FD formulas**: Generate higher - order compact RBF - FD formulas based on the Gaussian radial basis function for some differential operators, including first - and second - order derivative formulas (up to 10th order) and 2D - Laplacian formulas (up to 6th order). The weights of these formulas are obtained through symbolic computation, and Taylor series expansion is used to handle the non - constant terms of the linear system. 2. **Analyze local truncation error**: By substituting the derived weights into the local truncation error formula and using Taylor series expansion, obtain the expressions of the local truncation error of these formulas. This helps to evaluate the accuracy and convergence of the formulas. 3. **Verify the convergence of the formulas in the limit case**: Study whether these formulas converge to the corresponding polynomial - based compact finite - difference (compact FD) formulas when the shape parameter \(\epsilon\) approaches zero. The author verifies this by comparing with the classical results of Collatz [2] and Lele [3]. 4. **Compare the accuracy of different methods**: Verify the effectiveness of these new formulas through standard test functions, and compare them with the multi - quadratic (Multiquadric, MQ) - based RBF - FD formulas and the classical compact finite - difference schemes, showing their approximation accuracy under different shape parameter values. 5. **Calculate the optimal shape parameter**: For each formula, calculate the optimal shape parameter value through optimization techniques to ensure the best approximation effect under a specific test function. ### Key contributions - Proposed higher - order compact RBF - FD formulas based on the Gaussian radial basis function. - Obtained the exact weights and local truncation error expressions of these formulas through symbolic computation and Taylor series expansion. - Verified that these formulas converge to the corresponding polynomial - based compact finite - difference formulas when \(\epsilon \to 0\). - Demonstrated the superiority of these formulas over existing methods through numerical experiments and provided a method for selecting the optimal shape parameter. ### Formulas involved - **First - order derivative formula**: \[ u'(x_0) \approx \sum_{i = - n}^{n} \alpha_i u(x_0 + ih) + \sum_{k = - m}^{m} \beta_k u'(x_0 + kh) \] - **Second - order derivative formula**: \[ u''(x_0) \approx \sum_{i = - n}^{n} \alpha_i u(x_0 + ih) + \sum_{k = - m}^{m} \beta_k u''(x_0 + kh) \] - **2D - Laplacian formula**: \[ \Delta u(x_0, y_0) \approx \sum_{(i, j) \in S} \alpha_{i, j} u(x_0 + ih, y_0 + jh) + \sum_{(i, j) \in \mu} \beta_{i, j} \Delta u(x_0 + ih, y_0 + jh) \] The weights \(\alpha\) and \(\beta\) in these formulas are obtained by solving the linear system and converge to the weights of the corresponding polynomial - based compact finite - difference formulas when \(\epsilon \to 0\).

Comparative study on higher order compact RBF-FD formulas with Gaussian and Multiquadric radial functions

Multiquadric based RBF-HFD approximation formulas and convergence properties

Stable and Accurate Least Squares Radial Basis Function Approximations on Bounded Domains

Efficient Fourth-Order Weights in Kernel-Type Methods without Increasing the Stencil Size with an Application in a Time-Dependent Fractional PDE Problem

On the Fractional Derivatives of Radial Basis Functions

On a generalized Gaussian radial basis function: Analysis and applications

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

High precision compact numerical approximation in exponential form for the system of 2D quasilinear elliptic BVPs on a discrete irrational region

Efficient Scheme for the Economic Heston–Hull–White Problem Using Novel RBF-FD Coefficients Derived from Multiquadric Function Integrals

Application of Randomized Quadrature Formulas to the Finite Element Method for Elliptic Equations

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

On the numerical rank of radial basis function kernels in high dimension

Generalized Karagiannidis–Lioumpas Approximations and Bounds to the Gaussian Q-Function With Optimized Coefficients

Fractional-order dependent Radial basis functions meshless methods for the integral fractional Laplacian

Explicit radial basis function Runge-Kutta methods

Gaussian Variational Schemes on Bounded and Unbounded Domains

A superconvergence result in the RBF-FD method

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

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

On Newton–Cotes Formula-Type Inequalities for Multiplicative Generalized Convex Functions via Riemann–Liouville Fractional Integrals with Applications to Quadrature Formulas and Computational Analysis

Hyperviscosity-Based Stabilization for Radial Basis Function-Finite Difference (RBF-FD) Discretizations of Advection-Diffusion Equations