I want to simulate bead collision in fluid, the bead size is similar to the fluid flow characteristic length, e.g. 1/5 of flow depth. I found that most of existing CFD-DEM coupling code or software can only deal particles which are much smaller that mesh size. And I do believe the wake and other flow features around bead should be resolved in my problem. I am looking for a more direct and simple method than linking a cfd code with dem code using immersed boundary method, anyone get a clue?