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