An integral model based on slender body theory, with applications to curved rigid fibers

Abstract

We propose a novel integral model describing the motion of both flexible and rigid slender fibers in viscous flow and develop a numerical method for simulating dynamics of curved rigid fibers. The model is derived from nonlocal slender body theory (SBT), which approximates flow near the fiber using singular solutions of the Stokes equations integrated along the fiber centerline. In contrast to other models based on (singular) SBT, our model yields a smooth integral kernel which incorporates the (possibly varying) fiber radius naturally. The integral operator is provably negative definite in a nonphysical idealized geometry, as expected from the partial differential equation theory. This is numerically verified in physically relevant geometries. We discuss the convergence and stability of a numerical method for solving the integral equation. The accuracy of the model and method is verified against known models for ellipsoids. Finally, we develop an algorithm for computing dynamics of rigid fibers with complex geometries in the case where the fiber density is much greater than that of the fluid, for example, in turbulent gas-fiber suspensions.