A geometric optics ansatz-based plane wave method for two dimensional Helmholtz equations with variable wave numbers

Qiya Hu,Zezhong Wang

DOI: https://doi.org/10.48550/arXiv.2012.14250

2021-07-21

Abstract:In this paper we develop a plane wave type method for discretization of homogeneous Helmholtz equations with variable wave numbers. In the proposed method, local basis functions (on each element) are constructed by the geometric optics ansatz such that they approximately satisfy a homogeneous Helmholtz equation without boundary condition. More precisely, each basis function is expressed as the product of an exponential plane wave function and a polynomial function, where the phase function in the exponential function approximately satisfies the eikonal equation and the polynomial factor is recursively determined by transport equations associated with the considered Helmholtz equation. We prove that the resulting plane wave spaces have high order $h$-approximations as the standard plane wave spaces (which are available only to the case with constant wave number). We apply the proposed plane wave spaces to the discretization of nonhomogeneous Helmholtz equations with variable wave numbers and establish the corresponding error estimates of their finite element solutions. We report some numerical results to illustrate the efficiency of the proposed method.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop a new plane - wave method for discretizing the homogeneous Helmholtz equation with variable wave numbers. Specifically, the author aims to overcome the limitations of existing plane - wave methods in dealing with inhomogeneous Helmholtz equations and general variable wave numbers, especially the pollution effect at high wave numbers. The following is a detailed description of the problem:

1. **Background and Challenges**:
- The Helmholtz equation is a fundamental model for acoustic wave propagation, and high - precision solution in the case of large wave numbers is an important research topic.
- The standard finite - element method has the so - called "pollution effect" in the case of large wave numbers, that is, when the wave number increases, the accuracy of the approximate solution will be affected, unless the mesh size \(h\) satisfies \(h\omega = O(\omega^{-\delta})\), where \(\delta>0\).

2. **Limitations of Existing Methods**:
- The standard plane - wave method is suitable for cases of constant or piecewise - constant wave numbers, but not for general variable wave numbers.
- High - order finite - element methods can reduce the pollution effect, but not completely eliminate it.
- The standard plane - wave method is also not weakly pollution - free for inhomogeneous Helmholtz equations.

3. **Objectives of This Paper**:
- Propose a plane - wave method based on geometric - optics approximation to construct local basis functions that can approximately satisfy the homogeneous Helmholtz equation without boundary conditions on each element.
- Ensure that the new method can provide high - order \(h\)-approximation properties, similar to the standard plane - wave space (only applicable to constant wave numbers).
- Apply the new method to the discretization of inhomogeneous Helmholtz equations, and establish corresponding error estimates to prove its weak pollution - free property.

4. **Specific Methods**:
- Use the geometric - optics approximation (geometric optics ansatz) to represent the basis functions as the product of exponential plane - wave functions and polynomial functions.
- Among them, the phase function approximately satisfies the eikonal equation in the exponential function, and the polynomial factor is recursively determined by the transport equation derived from the Helmholtz equation.
- The basis functions constructed by this method can approximately satisfy the homogeneous Helmholtz equation on each element.

5. **Expected Results**:
- Prove that the proposed plane - wave space has high - order \(h\)-approximation properties.
- Apply it to the discretization of inhomogeneous Helmholtz equations and give error estimates, indicating that the method is weakly pollution - free.
- Verify the effectiveness of the proposed method through numerical experiments.

In summary, the main objective of this paper is to develop a new plane - wave method to effectively solve the discretization problem of the Helmholtz equation with variable wave numbers, especially to reduce the pollution effect at high wave numbers.