Modeling and direct simulation of velocity fluctuations and particle-velocity correlations in sedimentation

In this paper we present direct numerical simulations of monodisperse and polydisperse suspensions of non-Brownian particles sedimenting (it low Reynolds number We describe a scheme to generate ergodic ensembles of random particidate Systems and a numerical procedure for computing interactions among spherical particles based oil Elvald summation technique for hydrodynamic moility tensors. From the generation process truly random both monodisperse and multimodal size distributions of particles were obtained for dilute and moderate densities based oil a minimum energy criterion. Concerned with computations of the Ewald sums our numerical procedure drastically reduces the CPU simulation time providing results of the hindered settling function in good agreement with available experimental data and asymptotic results for ordered and random periodic arrays of particles. We show new computer simulations with no flux boundary perpendicular to gravity and periodic boundary conditions in horizontal direction. The simulations reproduce the experimental correlation-time and anisotropy of the velocity fluctuations, but have the magnitude of these fluctuations increasing proportional to the size of the system.