average normal, shear and penetration
int rather than std::size_t
pure virtual function
relocate calcTimeStep before findContact
calculate time steps
synchronize particle and fluid time step
test cvs
use drag coefficient Cd
print and apply penalization and pressure gradient forces to a particle
Further corrections in current configuration