Numerical Simulation of Complex Multi-Fluid Flows using a Combined Immersed Boundary and Volume of Fluid Approach

In this paper a simulation model is presented for the Direct Numerical Simulation (DNS) of complex multi-fluid flows in which simultaneously (moving) deformable (drops or bubbles) and non-deformable (moving) elements (particles) are present, possibly with the additional presence of free surfaces. Our model combines the VOF model developed by van Sint Annaland et al. (2005) and the Immersed Boundary (IB) model developed by van der Hoef et al. (2006). The Volume of Fluid (VOF) part features i) an interface reconstruction technique based on piecewise linear interface representation ii) a three-dimensional version of the CSF model of Brackbill et al. (1992). The Immersed Boundary (IB) part incorporates both particle-fluid and particle-particle interaction via a Direct Forcing Method (DFM) and a hard sphere Discrete Particle (DP) approach. In our model a fixed (Eulerian) grid is utilized to solve the Navier-Stokes equations for the entire computational domain. The no-slip condition at the surface of the moving particles is enforced via a momentum source term which only acts in the vicinity of the particle surface. For the enforcement of the no-slip condition Lagrangian force points are used which are distributed evenly over the surface of the particle. Dissipative particle-particle and/or particle-wall collisions are accounted via a hard sphere DP approach (Hoomans et al., 1996) using a three-parameter particle-particle interaction model accounting for normal and tangential restitution and tangential friction. The capabilities of the hybrid VOF-IB model are demonstrated with a number of examples in which complex topological changes in the interface are encountered.