E ect of flexibility on the growth of concentration fluctuations in a suspension of sedimenting fibers: Particle simulations

Three-dimensional numerical simulations are performed to study the stability of a sedimenting suspension of weakly flexible fibers. It is well known that a suspension of rigid rods sedimenting under gravity at low Reynolds number is unstable to concentration fluctuations owing to hydrodynamic interactions. Flexible fibers, however, reorient while settling and even weak flexibility can alter their collective dynamics. In our recent work [Manikantan et al., “The instability of a sedimenting suspension of weakly flexible fibres,” J. Fluid Mech. 756, 935–964 (2014)], we developed a mean-field theory to predict the linear stability of such a system. Here, we verify these predictions using accurate and e cient particle simulations based on a slender-body model. We also demonstrate the mechanisms by which flexibility-induced reorientation alters suspension microstructure, and through it, its stability. Specifically, we first show that the anisotropy of the base state in the case of a suspension of flexible fibers has a destabilizing e↵ect compared to a suspension of rigid rods. Second, a conflicting e↵ect of flexibility is also shown to suppress particle clustering and slow down the growth of the instability. The relative magnitude of filament flexibility and rotational Brownian motion dictates which e↵ect dominates, and our simulations qualitatively follow theoretically predicted trends. The mechanism for either e↵ects is tied to the flexibilityinduced reorientation of particles, which we illustrate using velocity and orientation statistics from our simulations. Finally, we also show that, in the case of an initially homogeneous and isotropic suspension, flexibility always acts to suppress the growth of the instability. C 2016 AIP Publishing LLC. [http://dx.doi.org/10.1063/1.4938493]