A fictitious domain approach for the simulation of dense suspensions

  • Authors:

  • Stany Gallier;Elisabeth Lemaire;Laurent Lobry;François Peters

  • Affiliations:

  • SAFRAN-Herakles, Le Bouchet Research Center, 91710 Vert le Petit, France and CNRS, University of Nice LPMC, UMR, 6622, Parc Valrose, 06108 Nice, France;CNRS, University of Nice LPMC, UMR, 6622, Parc Valrose, 06108 Nice, France;CNRS, University of Nice LPMC, UMR, 6622, Parc Valrose, 06108 Nice, France;CNRS, University of Nice LPMC, UMR, 6622, Parc Valrose, 06108 Nice, France

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2014

Quantified Score

Hi-index 31.45

Visualization

Abstract

Low Reynolds number concentrated suspensions do exhibit an intricate physics which can be partly unraveled by the use of numerical simulation. To this end, a Lagrange multiplier-free fictitious domain approach is described in this work. Unlike some methods recently proposed, the present approach is fully Eulerian and therefore does not need any transfer between the Eulerian background grid and some Lagrangian nodes attached to particles. Lubrication forces between particles play an important role in the suspension rheology and have been properly accounted for in the model. A robust and effective lubrication scheme is outlined which consists in transposing the classical approach used in Stokesian Dynamics to our present direct numerical simulation. This lubrication model has also been adapted to account for solid boundaries such as walls. Contact forces between particles are modeled using a classical Discrete Element Method (DEM), a widely used method in granular matter physics. Comprehensive validations are presented on various one-particle, two-particle or three-particle configurations in a linear shear flow as well as some O(10^3) and O(10^4) particle simulations.