Nonlocal hydrodynamic Drude model overcomes the limitations of the conventional Drude model in describing light-matter interactions in nanoplasmonic structures. Here we present a weak formulation based rigorous numerical method for it, which avoids spurious resonances as in the case of the curl-free approximation setting. The simulated results agree with Mie results, and the method is capable of handling arbitrarily shaped scatterers. © OSA 2014.