Andrea Adriani,Rosita Luisa Sormani,Cristina Tablino-Possio,Rolf Krause,Stefano Serra-Capizzano
Abstract:The current study investigates the asymptotic spectral properties of a finite difference approximation of nonlocal Helmholtz equations with a Caputo fractional Laplacian and a variable coefficient wave number $\mu$, as it occurs when considering a wave propagation in complex media, characterized by nonlocal interactions and spatially varying wave speeds. More specifically, by using tools from Toeplitz and generalized locally Toeplitz theory, the present research delves into the spectral analysis of nonpreconditioned and preconditioned matrix-sequences. We report numerical evidences supporting the theoretical findings. Finally, open problems and potential extensions in various directions are presented and briefly discussed.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: to study the asymptotic spectral characteristics of the non - local Helmholtz equation with the fractional - order Laplace operator and the variable - coefficient wavenumber \(\mu(x, y)\), and propose effective preconditioning methods. Specifically, the author focuses on the non - local interactions and spatially varying wave speeds that occur when describing wave propagation problems in complex numerical domains.
### Specific problem description
1. **Non - local Helmholtz equation**:
The core of the research is the two - dimensional form of the non - local Helmholtz equation, whose mathematical expression is:
\[
(-\Delta)^{\alpha/2} u(x, y) + \mu(x, y) u(x, y)=v(x, y), \quad (x, y) \in \Omega \subset \mathbb{R}^2, \quad \alpha \in (1, 2)
\]
Here, \((-\Delta)^{\alpha/2}\) is the fractional - order Laplace operator, defined as:
\[
(-\Delta)^{\alpha/2} u(x, y)=c_\alpha \text{P.V.} \int_{\mathbb{R}^2} \frac{u(x, y)-u(\tilde{x}, \tilde{y})}{[(x - \tilde{x})^2+(y - \tilde{y})^2]^{(2+\alpha)/2}} d\tilde{x} d\tilde{y}
\]
Here, \(c_\alpha\) is a constant related to \(\alpha\), \(\text{P.V.}\) represents the Cauchy principal value integral, and \(\Gamma\) is the Gamma function.
2. **Variable - coefficient wavenumber \(\mu(x, y)\)**:
The wavenumber \(\mu(x, y)\) in the equation is a complex - valued function and may be different at different positions. This variable coefficient makes the problem more complex because the traditional Toeplitz structure and real symmetry no longer apply.
3. **Asymptotic spectral analysis**:
The main objective of the paper is to study the asymptotic spectral properties of the coefficient matrix sequence after the discretization of the above equation. By introducing the generalized local Toeplitz (GLT) theoretical tool, the author analyzes the singular value and eigenvalue distributions of the matrix sequences before and after preconditioning.
4. **Preconditioning methods**:
In order to improve the solution efficiency, the paper proposes several preconditioning methods, especially the preconditioning techniques based on τ - algebra. These methods aim to improve the condition number of the linear system, thereby accelerating the iterative solution process.
### Main contributions
- **Theoretical analysis**: Provide a complete spectral characteristic analysis of the matrix sequence when \(\mu(x, y)\) is non - constant and complex - valued.
- **Numerical verification**: Verify the validity of the theoretical results through numerical experiments.
- **Potential expansion**: Discuss possible future research directions and open problems, such as further generalization to more complex geometric domains or higher - dimensional problems.
In summary, this paper aims to solve the efficient solution problem of the non - local Helmholtz equation with the fractional - order Laplace operator and the variable - coefficient wavenumber through in - depth theoretical analysis and numerical verification.