MODELING AND SIMULATION OF RIGID PARTICULATE FLOW IN AN INCOMPRESSIBLE NEWTONIAN FLUID

Abstract

We present a penalization method for the simulation of the interaction of an incompressible Newtonian fluid flow with rigid bodies. This method is based on a variational formulation throughout the fluid/solid domain, with constraints on the unknown and on the test functions. The rigid motion of the particle is enforced by penalizing the strain tensor on the rigid domain for canceling the deformation rate in the volume occupied by the particle. A simulations , implemented with FreeFem++, of a rigid particles in a incompressible Newtonian fluid flow are performed. The method is first validated by simulating the Sheared particle and then aplicated for simulating the sedimentation of one single particle and two identical particles in an incompressible fluid.