Abstract:In this paper, we present a Localized Orthogonal Decomposition (LOD) in Petrov-Galerkin formulation for a two-scale Helmholtz-type problem. The two-scale problem is, for instance, motivated from the homogenization of the Helmholtz equation with high contrast, studied together with a corresponding multiscale method in (Ohlberger, Verfürth. A new Heterogeneous Multiscale Method for the Helmholtz equation with high contrast, <a class="link-https" data-arxiv-id="1605.03400" href="https://arxiv.org/abs/1605.03400">arXiv:1605.03400</a>, 2016). There, an unavoidable resolution condition on the mesh sizes in terms of the wave number has been observed, which is known as "pollution effect" in the finite element literature. Following ideas of (Gallistl, Peterseim. Comput. Methods Appl. Mech. Engrg. 295:1-17, 2015), we use standard finite element functions for the trial space, whereas the test functions are enriched by solutions of subscale problems (solved on a finer grid) on local patches. Provided that the oversampling parameter $m$, which indicates the size of the patches, is coupled logarithmically to the wave number, we obtain a quasi-optimal method under a reasonable resolution of a few degrees of freedom per wave length, thus overcoming the pollution effect. In the two-scale setting, the main challenges for the LOD lie in the coupling of the function spaces and in the periodic boundary conditions.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the "pollution effect" encountered in the numerical solution of high - frequency Helmholtz problems. Specifically, when using the standard finite - element method to solve two - scale Helmholtz - type problems with high contrast, due to the oscillatory nature of the solution, a mesh size much smaller than the wavelength is required to ensure the stability and convergence of the numerical scheme. This usually leads to a resolution condition \(k\alpha H = O(1)\) instead of \(kH = O(1)\), that is, more degrees of freedom are required per wavelength.
To overcome this challenge, the author introduced a Localized Orthogonal Decomposition (LOD) method, specifically for two - scale Helmholtz - type problems. This method is implemented through the Petrov - Galerkin formulation and combines the ideas of multi - scale methods. Specifically, this method uses standard finite - element functions as the trial space, while the test functions are enriched by solving sub - scale problems on a fine grid. As long as the oversampling parameter \(m\) is logarithmically related to the wave number \(k\), this method can provide quasi - optimal results under reasonable degrees of freedom per wavelength conditions, thereby overcoming the pollution effect.
### Main contributions:
1. **Theoretical analysis**: The author conducted a rigorous theoretical analysis of LOD in a two - scale setting, especially in dealing with coupled function spaces and sesquilinear forms.
2. **Oversampling parameter selection**: By selecting an appropriate oversampling parameter \(m\approx\log(k)\), the stability and quasi - optimality of the method are ensured.
3. **Periodic boundary conditions**: The periodic boundary conditions on the unit cube are considered, which are crucial for defining the interpolation operator and oversampling patches.
### Key formulas:
- Relationship between wave number \(k\) and mesh sizes \(H_c, h_c\): \(k(H_c + h_c)\lesssim1\)
- Selection of oversampling parameter \(m\): \(m\approx\log(k)\)
Through these improvements, the LOD method can effectively solve the pollution effect in high - frequency Helmholtz problems while maintaining a relatively low computational complexity.