Numerical simulations of granular dynamics: I. Hard-sphere discrete element method and tests

► We present a new numerical method for simulating granular dynamics. ► The technique has application to small-body surface particle motions. ► Our approach is fast, scalable, and adaptable to different problems. ► Full derivations of all collision equations are provided. ► Two tests are presented showing correct functionality of the method.