On the application of quaternion‐based approaches in discrete element methods