Simulation of three-dimensional nanoscale light interaction with spatially dispersive metals using a high order curvilinear DGTD method

In this work, we present and study a flexible and accurate numerical solver in the context of three-dimensional computational nanophotonics. More precisely, we focus on the propagation of electromagnetic waves through metallic media described by a non-local dispersive model. For this model, we propose a discretization based on a high-order Discontinuous Galerkin time-domain method, along with a low-storage Runge-Kutta time scheme of order four. The semi-discrete stability of the scheme is analyzed for classical numerical fluxes, i.e. centered and upwind. Furthermore, the numerical treatment is enriched with an enhanced approximation of the geometry based on isoparametric curvilinear meshes. We finally assess our approach on several test cases, from academic to more physical ones.