Fully resolved simulation of dense suspensions of freely evolving buoyant particles using an improved immersed boundary method.

Fully resolved simulation of flows with buoyant particles is a challenging problem since buoyant particles are lighter than the surrounding fluid, and as a result, the two phases are strongly coupled together. In this work, the virtual force stabilization technique introduced by Schwarz et al. [Schwarz, S., Kempe, T., & Fr\"ohlich, J. (2015). A temporal discretization scheme to compute the motion of light particles in viscous flows by an immersed boundary method. J. Comput. Phys., 281, 591-613] is extended to simulate buoyant particle suspensions with high volume fractions (up to $40 \%$). It is concluded that the dimensionless numerical model constant $C_v$ in the virtual force technique should increase with volume fraction. The behavior of a single rising particle, two in-line rising particles, and buoyant particle suspensions are studied. In each case, results are compared with experimental works on bubbly flows to highlight the differences and similarities between buoyant particles and bubbles. Finally, the drag coefficient is extracted from simulations of buoyant particle suspensions at different volume fractions and based on that a drag correlation is presented.