Numerical simulation of the motion of neutrally buoyant particles in the plane Poiseuille flow of a Newtonian fluid

In this note we discuss the generalization of a Lagrange multiplier based fictitious domain method to the simulation of migration of neutrally buoyant particles in plane Poiseuille flow of a Newtonian fluid. The migration away from the center of the channel is believed to be an effect of the curvature of the velocity profile. We found this effect is not weakened by the presence of many particles, but instead by the collisions among the particles. Experiments and simulations show that the particles concentrate in the central region where the shear rate is low.