Dear Dr. Mathias,
I think the current hard sphere collision model will result in spheres overlapping if the suspension is not dilute or multiple contacts occur. Would you mind to share your code used for elastic collision (the implementation of parameter gamma) showing in your paper?
Thanks,
Junwei Guo