A numerical study of the settling of non-spherical particles in quiescent water

Abstract

Settling of non-spherical particles is poorly understood with previous studies having focused mainly on spherical particles. Here, a series of particle-resolved direct numerical simulations are conducted using FLOW-3D (commercial computational fluid dynamics software) for spheres and five regular, non-spherical shapes of sediment particles, i.e., prolate spheroid, oblate spheroid, cylinder, disk, and cube. The Galileo number varies from 0.248 to 360, and the particle Reynolds number Rep ranges from 0.002 77 to 562. The results show that a non-spherical particle may experience larger drag and, consequently, attain a lower terminal velocity than an equivalent sphere. If Rep is sufficiently small, the terminal velocity is less affected by particle shape as characterized by the particle aspect ratio. For relatively large Rep, the shape effect (represented by the Corey shape factor) becomes more significant. Empirical correlations are derived for the dimensionless characteristic time t95∗ and displacement s95∗ of particle settling, which show that t95∗ remains constant in the Stokes regime (Rep < 1) and decreases with increasing Rep in the intermediate regime (1 ≤ Rep < 103), whereas s95∗ increases progressively with increasing Rep over the simulated range. It is also found that in the Stokes regime, particle orientation remains essentially unchanged during settling, and so the terminal velocity is governed by the initial orientation. In the intermediate regime, a particle provisionally settling at an unstable orientation self-readjusts to a stable equilibrium state, such that the effect of initial orientation on the terminal velocity is negligible. Moreover, an unstable initial orientation can enhance the vertical displacement and may promote vortex shedding.