Two-step perfectly matched layer for arbitrary-order pseudo-spectral analytical time-domain methods

Olga Shapoval 1 Jean-Luc Vay 1 Henri Vincenti 2, 3
3 PHI - Physique à Haute Intensité
IRAMIS - Institut Rayonnement Matière de Saclay, LIDyl - Laboratoire Interactions, Dynamique et Lasers (ex SPAM)
Abstract : Numerical simulation of an electrodynamic system in empty space requires the implementation of open boundary conditions (BC) to terminate the solution of Maxwell’s equations on the boundaries of the computational domain. The Perfectly Matched Layer (PML) has become the method of choice for open BC with wave equations, as it is straightforward and relatively easy to implement, and offers very efficient and user-adjustable absorption rates. PMLs are most often employed with the Finite-Difference Time-Domain (FDTD) algorithm, which in its most common implementation offers second-order accuracy in space and time on Cartesian grids. Yet, simulations (including some class of electromagnetic Particle-In-Cell simulations) that require higher precision may resort to higher-order Maxwell solvers employing extended finite-difference stencils, or even to pseudo-spectral Maxwell solvers, for which a general, versatile and efficient formulation of the PML has been missing so far. In this paper, we propose a novel “two-step” formulation of the PML that is simple, very versatile and can be used as is with any Maxwell solver. In particular, it is applicable to a large class of Maxwell solvers including the arbitrary-order Pseudo-Spectral Analytical Time-Domain (PSATD) solver, which offers arbitrarily low numerical dispersion when increasing solver order and becomes dispersion-free at infinite order. Analysis and numerical simulations demonstrate that the new formulation is as efficient as the standard PML formulation, both for the FDTD and the PSATD implementations.
Document type :
Journal articles
Complete list of metadatas

https://hal-cea.archives-ouvertes.fr/cea-01934810
Contributor : Caroline Lebe <>
Submitted on : Monday, November 26, 2018 - 11:38:32 AM
Last modification on : Thursday, February 7, 2019 - 5:14:51 PM

Identifiers

Citation

Olga Shapoval, Jean-Luc Vay, Henri Vincenti. Two-step perfectly matched layer for arbitrary-order pseudo-spectral analytical time-domain methods. Computer Physics Communications, Elsevier, 2019, 235, pp.102 - 110. ⟨10.1016/j.cpc.2018.09.015⟩. ⟨cea-01934810⟩

Share

Metrics

Record views

55