We demonstrate that simulating plasmonic nanostructures by means of curved elements (CEs) significantly increases the accuracy and computation speed not only in the linear but also in the nonlinear regime. We implemented CEs within the discontinuous Galerkin (DG) method and, as an example of a nonlinear effect, investigated second-harmonic generation (SHG) at a silver nanoparticle. The second-harmonic response of the material is simulated by an extended Lorentz model (ELM). In the linear regime the CEs are ? 9 times faster than ordinary elements for the same accuracy, provide a much better convergence and show fewer unphysical field artifacts. For DG-SHG calculations CEs are almost indispensable to obtain physically reasonable results at all. Additionally, their boundary approximation has to be of the same order as their polynomial degree to achieve artifact-free field distributions. In return, the use of such CEs with the DG method pays off more than evidently, since the additional computation time is only 1%.