Abstract
A coupled system of hydrodynamic and volume integral equations is solved for analyzing electromagnetic wave interactions with non-local dispersion effects on nanostructures. The proposed scheme discretizes the scatterer into a mesh of tetrahedral elements and expands (unknown) electric flux and hydrodynamic current using Schaubert-Wilton-Glisson basis functions defined on this mesh. Inserting these expansions into the coupled equations and applying Galerkin testing yield a matrix system. An iterative scheme is used to solve this matrix system for unknown expansion coefficients. Numerical results show additional resonance peaks in the scattering cross section spectrum of nanospheres, which can be explained by non-local dispersion effects accounted for by the proposed method.