Fluid-particle simulations with FreeFem++. (English) Zbl 1354.76109

Lefebvre, Aline

ESAIM, Proc. 18, 120-132 (2007).

Summary: We present here a method to simulate the motion of rigid particles in a Newtonian fluid. The rigid motion is enforced by penalizing the strain tensor on the rigid domain and the time discretization is performed using the method of characteristics. This leads to a generalized Stokes variational formulation on the whole domain which can easily be implemented from any Finite Element Stokes/Navier-Stokes solver. In order to ensure robustness, we describe a strategy to take collisions into account. To validate the method, simulations implemented with FreeFem++ are presented.

MSC:

76M10	Finite element methods applied to problems in fluid mechanics
76T20	Suspensions
76F65	Direct numerical and large eddy simulation of turbulence
76D05	Navier-Stokes equations for incompressible viscous fluids
65M60	Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs

Keywords:

particulate flows; fluid-particle flows; direct numerical simulation; penalty method; Navier-Stokes equations; finite-element method; inelastic contacts

Software:

FreeFem++