A novel coarse space applying to the weighted Schwarz method for Helmholtz equations

Qiya Hu,Ziyi Li

2024-02-10

Abstract:In this paper we are concerned with restricted additive Schwarz with local impedance transformation conditions for a family of Helmholtz problems in two dimensions. These problems are discretized by the finite element method with conforming nodal finite elements. We design and analyze a new adaptive coarse space for this kind of restricted additive Schwarz method. This coarse space is spanned by some eigenvalue functions of local generalized eigenvalue problems, which are defined by weighted positive semi-definite bilinear forms on subspaces consisting of local discrete Helmholtz-harmonic functions from impedance boundary data. We proved that a two-level hybrid Schwarz preconditioner with the proposed coarse space possesses uniformly convergence independent of the mesh size, the subdomain size and the wave numbers under suitable assumptions. We also introduce an economic coarse space to avoid solving generalized eigenvalue problems. Numerical experiments confirm the theoretical results.

Numerical Analysis

What problem does this paper attempt to address?

### What problems does this paper attempt to solve?

This paper aims to solve several challenges in the numerical solution of the Helmholtz equation, especially for Helmholtz problems with large wave numbers (wave number). Specifically, the authors focus on how to improve the existing Weighted Schwarz Method and design a new coarse space to improve the solution efficiency and stability.

#### Main problems:
1. **Wave Number Pollution**: For the Helmholtz equation, when the wave number \(\kappa\) increases, in order to ensure the accuracy of the solution, the mesh size \(h\) must be small enough, which leads to a sharp increase in the size of the linear system and the system is highly indefinite and ill - conditioned. Therefore, the traditional finite - element discretization method is very expensive and difficult when dealing with large - wave - number problems.

2. **Limitations of the Standard Domain Decomposition Methods**: The standard Domain Decomposition Methods are inefficient when dealing with the Helmholtz equation with large wave numbers, unless a very fine coarse mesh (coarse mesh) is selected, that is, requiring \(d\kappa = O(1)\), which is impractical in practical applications.

3. **Insufficient Convergence Analysis of Existing Methods**: Although there have been some studies on the Weighted Schwarz Method, there is a lack of strict convergence analysis of the two - level weighted Schwarz method combined with an adaptive coarse space.

### Solutions:
The authors propose a new adaptive coarse space, which is spanned by the eigenfunctions of local generalized eigenvalue problems. These eigenvalue problems are defined on the subspaces composed of local discrete Helmholtz - harmonic functions defined by impedance boundary data.

#### Key innovation points:
- **New Weighted Positive Definite Bilinear Form**: Used to define local generalized eigenvalue problems to ensure the effectiveness of the new coarse space.
- **Theoretical Proof**: It is proved that the two - level hybrid WASI method using the proposed coarse space has uniform convergence independent of the mesh size \(h\), the sub - domain size \(d\) and the wave number \(\kappa\).
- **Economical Coarse Space**: A simplified version of the coarse space is introduced to avoid solving generalized eigenvalue problems, thus improving the computational efficiency.

### Conclusion:
Through the above methods, the authors have successfully solved the problems of wave number pollution and inefficiency of the domain decomposition method encountered by traditional methods when dealing with large - wave - number Helmholtz problems. Numerical experiments have verified the effectiveness and superiority of the proposed method.

---

If you have more specific questions or need further explanations, please feel free to let me know!