Inertial migration in dilute and semi-dilute suspensions of rigid particles in laminar square duct flow

Abstract

We study the inertial migration of finite-size neutrally buoyant spherical particles in dilute and semi-dilute suspensions in laminar square duct flow. We perform several direct numerical simulations using an immersed boundary method to investigate the effects of the bulk Reynolds number $Re_b$, particle Reynolds number $Re_p$ and duct to particle size ratio $h/a$ at different solid volume fractions $\phi$, from very dilute conditions to $20\%$. We show that the bulk Reynolds number $Re_b$ is the key parameter in inertial migration of particles in dilute suspensions. At low solid volume fraction ($\phi=0.4\%$) and low bulk Reynolds number $(Re_b= 144)$, particles accumulate at the center of the duct walls. As $Re_b$ is increased, the focusing position moves progressively towards the corners of the duct. At higher volume fractions, $\phi=5, 10$ and $20\%$, and in wider ducts with $Re_b=550$, particles are found to migrate away from the duct core towards the walls. In particular, for $\phi=5$ and $10\%$, particles accumulate preferentially at the corners. At the highest volume fraction considered, $\phi=20\%$, particles sample all the volume of the duct, with a lower concentration at the duct core. The presence of particles induces secondary cross-stream motions in the duct cross-section, for all $\phi$. The intensity of these secondary flows depends strongly on particle rotation rate, on the maximum concentration of particles in focusing positions, and on the solid volume fraction. We find that the secondary flow intensity increases with the volume fraction up to $\phi=5\%$. However, beyond $\phi=5\%$ excluded volume effects lead to a strong reduction of cross-stream velocities.