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

TL;DR

This paper presents an advanced immersed boundary method for fully resolving dense suspensions of buoyant particles, extending previous techniques to handle high volume fractions and comparing results with experimental bubbly flow data.

Contribution

It introduces an improved virtual force stabilization technique for simulating high volume fraction buoyant particle suspensions with accurate results.

Findings

The model constant $C_v$ increases with volume fraction.

Simulations of single and multiple rising particles match experimental observations.

A new drag correlation for buoyant particle suspensions is proposed.

Abstract

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.