Project Details
Projekt Print View

Higher order numerical methods for acoustic scattering problems with locally perturbed periodic structures

Subject Area Mathematics
Term from 2019 to 2024
Project identifier Deutsche Forschungsgemeinschaft (DFG) - Project number 433126998
 
Final Report Year 2024

Final Report Abstract

Time-harmonic wave scattering problems in (locally perturbed) bi-periodic structures are important topics in application areas such as optics or photonic crystals. From the mathematical point of view, the problem is modeled by the Helmholtz equation in a half space above a bounded bi-periodic surface (assume it is periodic in x1 and x2 directions) in a three-dimensional space, with a boundary condition on the bi-periodic surface and a proper radiation condition in the x3 direction which guarantees the well-posedness of the problem. The numerical simulation for such a problem is particularly interesting but also very challenging due to the unbounded domain. To solve the problem by the finite element method, we need to approximate the scattering problem by one defined in a bounded domain. With a proper truncation technique, such as the transparent boundary condition (TBC) or the perfectly matched layers (PML), we first reduce the unbounded domain into a bi-periodic waveguide which is bounded in the x3 direction. Then applying the Floquet-Bloch transform, the new problem is rewritten into an equivalent family of quasi-periodic problems, which are suitable to be solved by the finite element method. The TBC results in a complex singular integral in a square due to the Rayleigh anomalies. To approximate the integral numerically, we developed a high-order tailor-made quadrature method. The PML, on the other hand, is no longer exact but gets rid of singularities from the TBC. But since it is an approximation, we also work on the error estimations to guarantee the convergence of this method. Finally, the method is also extended to locally perturbed periodic surfaces based on a domain transformation. An efficient numerical method is also developed based on the Schur complement. We have also presented numerous numerical examples to illustrate our theoretical results.

Publications

 
 

Additional Information

Textvergrößerung und Kontrastanpassung