An Eulerian-based immersed boundary method for particle suspensions with implicit lubrication model. (30th March 2022)