Abstract:We propose a novel method to compute a finite difference stencil for Riesz derivative for artibitrary speed of convergence. This method is based on applying a pre-filter to the Grünwald-Letnikov type central difference stencil. The filter is obtained by solving for the inverse of a symmetric Vandemonde matrix and exploiting the relationship between the Taylor's series coefficients and fast Fourier transform. The filter costs O\left(N^{2}\right) operations to evaluate for O\left(h^{N}\right) of convergence, where h is the sampling distance. The higher convergence speed should more than offset the overhead with the requirement of the number of nodal points for a desired error tolerance significantly reduced. The benefit of progressive generation of the stencil coefficients for adaptive grid size for dynamic problems with the Grünwald-Letnikov type difference scheme is also kept because of the application of filtering. The higher convergence rate is verified through numerical experiments.
What problem does this paper attempt to address?
### What problem does this paper attempt to solve?
This paper aims to solve the problem of arbitrary - order convergence of the Riesz fractional - order derivative. Specifically, the author proposes a novel method to calculate the finite - difference template of the Riesz derivative, enabling its convergence rate to reach an arbitrary order. Traditional methods can usually only achieve lower - order convergence (such as second - order or fourth - order), while this paper can significantly improve the convergence rate on the basis of the central difference method by introducing a pre - filter.
#### Key points:
1. **Riesz derivative**: The Riesz derivative is the symmetric version of the fractional - order Riemann - Liouville derivative, which can describe non - local effects and has important applications in many physical problems.
2. **Grüwald - Letnikov - type central difference scheme**: This is a commonly used discretization method, but traditional methods can only reach a lower convergence order.
3. **Pre - filter**: The author proposes to improve the Grüwald - Letnikov - type central difference scheme by applying a pre - filter, thereby achieving an arbitrary - order convergence rate.
4. **Filter design**: The filter is designed based on solving the inverse of a symmetric Vandemonde matrix and using the relationship between Taylor series coefficients and the fast Fourier transform. The computational complexity of the filter is \(O(N^2)\), where \(N\) is the desired convergence order.
#### Core idea of the method:
- By introducing a pre - filter, on the basis of the central difference scheme, low - order error terms can be eliminated, thereby improving the convergence rate.
- The design of the filter ensures that in the frequency domain, the response of the difference operator is closer to the response of the Riesz derivative in the ideal case.
#### Numerical experiment verification:
- The author verifies the effectiveness of the proposed method through numerical experiments, indicating that this method can indeed achieve the expected high - order convergence rate.
- The experimental results also show that for some specific types of functions, even if the smoothness requirements are not fully met, this method can still exhibit good convergence performance.
In general, the main contribution of this paper is to provide an efficient and flexible method that can achieve arbitrary - order convergence of the Riesz derivative in the finite - difference framework, thereby providing a powerful tool for solving complex physical problems involving fractional - order calculus.