Numerical Simulation of The Settling of Polydisperse Suspensions of Spheres

. The extension of Kynch's kinematical theory of ideal suspensions to polydisperse suspensions of spheres leads to a nonlinear system of conservation laws for the volumetric concentration of each species. In this work, we consider particle species different in sizes and densities, including the buoyant case. We show that modern shock-capturing numerical schemes for the solution of systems of conservation laws can be employed as an efficient tool for the simulation of the settling and separation of polydisperse suspensions. This is demonstrated by comparison with published experimental and theoretical results and by simulating some hypothetical configurations. Particular attention is focused on the emergence of rarefaction waves. 1. Introduction The one-dimensional sedimentation of ideal suspensions of monosized spheres under the influence of gravity can be described quite well by Kynch's sedimentation theory [1, 2, 3]. The principal assumption of this theory is that the local settling..