Abstract:We study the numerical evaluation of the integral fractional Laplacian and its application in solving fractional diffusion equations. We derive a pseudo-spectral formula for the integral fractional Laplacian operator based on fractional order-dependent, generalized multi-quadratic radial basis functions (RBFs) to address efficient computation of the hyper-singular integral. We apply the proposed formula to solving fractional diffusion equations and design a simple, easy-to-implement and nearly integration-free meshless method. We discuss the convergence of the novel meshless method through equivalent Galerkin formulations. We carry out numerical experiments to demonstrate the accuracy and efficiency of the proposed approach compared to the existing method using Gaussian RBFs.
What problem does this paper attempt to address?
This paper attempts to address the challenges of the fractional Laplacian in numerical calculations, especially in dealing with singular integrals and non - locality. Specifically, the main contributions of the article include:
1. **Proposing a new mesh - free method**: By using generalized multi - quadratic radial basis functions (GMQ RBFs), this method avoids the complex calculations of singular integrals and series truncation for the fractional Laplacian.
2. **Improving computational efficiency**: By using the generalized multi - quadratic function as a radial basis function, the integral form of the fractional Laplacian can be solved explicitly and simply, thus significantly improving the computational efficiency, especially in large - scale numerical evaluations.
3. **Proving the convergence of the method**: The article discusses the convergence of the proposed method and proves that the method has spectral convergence through an equivalent Galerkin formula.
4. **Applying to the fractional - order diffusion equation**: The proposed formula is applied to solve the fractional - order diffusion equation, and a simple, easy - to - implement, and almost integration - free mesh - free method is designed.
5. **Numerical experiment verification**: Numerical experiments are used to show the advantages of this method in terms of accuracy and efficiency compared to existing methods (such as the method using Gaussian radial basis functions).
### Main problems
The main problem solved in the paper is how to calculate the fractional Laplacian efficiently and accurately, especially for the fractional - order diffusion equation on general bounded domains. Traditional numerical methods face problems of high computational complexity and difficult handling of singular integrals when dealing with the fractional Laplacian, while the method proposed in this paper effectively solves these problems by introducing generalized multi - quadratic radial basis functions.
### Formula representation
The key formulas involved in the article include:
- Definition of the fractional Laplacian:
\[
(-\Delta)^{\alpha/2}u(x) := c_{d,\alpha}\int_{\mathbb{R}^d}\frac{u(x) - u(y)}{|x - y|^{d+\alpha}}\,dy,
\]
where \(c_{d,\alpha}=\frac{2^\alpha\Gamma(\frac{\alpha + d}{2})}{\pi^{d/2}|\Gamma(-\alpha/2)|}\).
- Fractional Laplacian of the generalized multi - quadratic function:
\[
(-\Delta)^{\alpha/2}\Phi_{d,\alpha}(x)=\mu_{d,\alpha}\Phi_{d,-\alpha}(x),
\]
where \(\Phi_{d,\alpha}(x)=(1 + |x|^2)^{\alpha/2 - d/2}\), \(\mu_{d,\alpha}=\frac{2^\alpha\Gamma(d/2+\alpha/2)}{\Gamma(d/2-\alpha/2)}\).
The introduction of these formulas makes the calculation of the fractional Laplacian more concise and efficient.