We present recent advances in the development of a Finite Element Time Domain (Discontinuous Galerkin) solver for computational nanophotonics, focusing on metallic nano structures irradiated by laser pulses. Metallic nano structures of sizes between 2 nm and 25 nm are well known to show spatial dispersion which is modeled by a nonlocal dispersion model for the electron gas [1, 2]. We briefly review and classify the dispersion model used in this work and then step to the numerical modeling of the resulting linear hydrodynamic equations when coupled to the 3D Maxwell's equations in time-domain. The considered 3D finite element time domain method combines the advantage of high order polynomial solutions on very flexible unstructured tetrahedral meshes and a broad band frequency solution due to short pulses [3]. We assess the performance of the numerical method on multiple setups. Among them, we considered a spherical dimer system. Such systems are well known for high field enhancements in the gap and are well suited for single molecule detection [4]. Our numerical simulations investigate the influence of spatial dispersion on the scattering cross-section of spherical dimers with a diameter of 20 nm and 2 nm gap size. Figure 1 compares the scattering cross-section for an incident polarization of the electric field along the dimer.