simulation of charge motion in ball mills part 1