Energy and accuracy issues in numerical simulations of a non-linear impact model