Direct simulation of rigid particles in a viscoelastic fluid

Astrid Decoene (1) , Sébastien Martin (2) , Bertrand Maury (1, 3)

1 LM-Orsay - Laboratoire de Mathématiques d'Orsay
2 MAP5 - UMR 8145 - Mathématiques Appliquées Paris 5
3 DMA - Département de Mathématiques et Applications - ENS Paris

In this paper, we present a numerical method which performs the direct simulation of 2D viscoelastic suspensions. Objectives. Interactions between viscoelastic fluids of Oldroyd type, including inertial effects, and a large number of rigid spheres or ellipsoids are addressed. Method. The numerical method is built upon a fictitious domain approach which consists in defining the fluid-structure problem over the whole domain (fluid+solid). In the Newtonian case, the resulting variational formulation, although classical, uses non-standard functional spaces which include the rigidity constraints. From the numerical point of view, these constraints can be addressed by penalty methods, leading to the possible use of standard finite elements solvers with fixed structured meshes. We show how this approach can be adapted to viscoelastic fluids of Oldroyd type. Results. Two types of configurations are investigated in dilute and dense regimes: sedimentation of spherical particles and shear flows with ellipsoidal particles.

Mots clés: characteristics method; dense suspensions; viscoelastic fluid; Oldroyd model; fluid-structure interaction; fictitious domain; penalty method; characteris-; tics method

Domaines: Equations aux dérivées partielles [math.AP]; Analyse numérique [math.NA]; Mécanique des fluides [physics.class-ph]

Astrid Decoene, Sébastien Martin, Bertrand Maury. Direct simulation of rigid particles in a viscoelastic fluid. 2018. ⟨hal-01456089v2⟩