Numerical Methods for Quasi-Periodic Incident Fields Scattered by Locally Perturbed Periodic Surfaces
For researchers in computational wave scattering, this work extends numerical methods to a practically important class of incident fields (plane waves) that were previously excluded due to non-decay at infinity.
This paper develops Bloch transform-based finite element methods for scattering of quasi-periodic incident fields (including plane waves) by locally perturbed periodic surfaces, achieving high-order accuracy for non-decaying incident fields. Numerical experiments validate both standard and high-order methods.
Waves scattering from unbounded structures are always complicated problems for numerical simulations. For the case that the non-periodic incident field scattered by (locally perturbed) periodic surfaces, with the help of the Bloch transform, the problem could be solved by some finite element methods, if the incident fields decay at certain rate at the infinity. For faster decaying incident fields, a high order numerical method is also available. However, in these cases, the plain waves, which belong to a very important family of incident fields but do not decay at the infinity, are not included. In this paper, we aim to develop the Bloch transform based standard finite method for this certain case, and then establish the high order method afterwards. Numerical experiments have been carried out for both the standard and high order numerical methods. Based on the algorithms for incident plain waves, we could also extend the numerical methods to more generalized cases when only the not so efficient standard method is available.