High-Order Methods for High-Frequency Scattering Applications

Abstract

The problem of simulating the scattering of (acoustic, electromagnetic, elastic) waves has provided a particularly sustained, demanding and motivating challenge for the development of efficient and accurate numerical methods since the advent of computers. The classical issues present in most other applications, such as those related to the environmental and/or geometrical intricacies of the media in which quantities of interest are defined, are augmented in the context of wave propagation by the intrinsic complexities (i.e. oscillations) of the quantities themselves. Still, very efficient methodologies have been devised, particularly in the last twenty years, to simulate wave processes in rather complex settings. These techniques can be based, for instance, on finite elements (see e.g. [34, 35, 46] and the references therein), finite differences [40, 49] or boundary integral equations [4, 8, 11, 16, 25], and they can, today, effectively address these problems, with a high degree of accuracy, in domains that can span tens or perhaps even a few hundred wavelengths. The very nature of these classical approaches, however, limits their applicability at higher frequencies since the numerical resolution of field oscillations translates in a commensurately higher number of degrees of freedom and this, in turn, can easily lead to impractical computational times. In this chapter we review some recently proposed methodologies [5, 13–15, 29] that can overcome these limitations while retaining the mathematical rigor of classical numerical procedures.