Simulation of semidilute suspensions of non-Brownian fibers in shear flow

Particle-level simulations are performed to study semidilute suspensions of monodispersed non-Brownian fibers in shear flow, with a Newtonian fluid medium. The incompressible three-dimensional Navier-Stokes equations are used to describe the motion of the medium, while fibers are modeled as chains of fiber segments, interacting with the fluid through viscous drag forces. The two-way coupling between the solids and the fluid phase is taken into account by enforcing momentum conservation. The model includes long-range and short-range hydrodynamic fiber-fiber interactions, as well as mechanical interactions. The simulations rendered the time-dependent fiber orientation distribution, whose time average was found to agree with experimental data in the literature. The viscosity and first normal stress difference was calculated from the orientation distribution using the slender body theory of Batchelor [J. Fluid Mech. 46, 813 (1971)], with corrections for the finite fiber aspect ratios. The viscosity was also obtained from direct computation of the shear stresses of the suspension for comparison. These two types of predictions compared well in the semidilute regime. At higher concentrations, however, a discrepancy was seen, most likely due to mechanical interactions, which are only accounted for in the direct computation method. The simulated viscosity determined directly from shear stresses was in fair agreement with experimental data found in the literature. The first normal stress difference was found to be proportional to the square of the volume concentration of fibers in the semidilute regime. As concentrations approached the concentrated regime, the first normal stress difference became proportional to volume concentration. It was also found that the coefficient of friction has a strong influence on the tendency for flocculation as well as the apparent viscosity of the suspension in the semidilute regime.