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 timedomain method, along with a low-storage Runge-Kutta time scheme of order four. The semidiscrete 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.