Optimizing perfectly matched layers in discrete contexts

Axel Modave, Eric Delhez and Christophe Geuzaine
august, 2014
Publication type:
Paper in peer-reviewed journals
International Journal for Numerical Methods in Engineering, vol. 99 (6), pp. 410 - 437
Perfectly matched layers (PMLs) are widely used for the numerical simulation of wave-like problems defined on large or infinite spatial domains. However, for both time-dependent and time-harmonic cases, their performance critically depends on the so-called absorption function. This paper deals with the choice of this function when classical numerical methods are used (based on finite differences, finite volumes, continuous finite elements and discontinuous finite elements). After reviewing the properties of the PMLs at the continuous level, we analyze how they are altered by the different spatial discretizations. In the light of these results, different shapes of absorption function are optimized and compared by means of both one-dimensional and two-dimensional representative time-dependent cases. This study highlights the advantages of the so-called shifted hyperbolic function, which is efficient in all cases and does not require the tuning of a free parameter, by contrast with the widely used polynomial functions.
    author={Axel Modave and Eric Delhez and Christophe Geuzaine },
    title={Optimizing perfectly matched layers in discrete contexts },
    doi={10.1002/nme.4690 },
    journal={International Journal for Numerical Methods in Engineering },
    year={2014 },
    volume={99 (6) },