The hydrodynamics of a twisting, bending, inextensible fiber in Stokes flow

by   Ondrej Maxian, et al.

In swimming microorganisms and the cell cytoskeleton, inextensible fibers resist bending and twisting, and interact with the surrounding fluid to cause or resist large-scale fluid motion. In this paper, we develop a novel numerical method for the simulation of cylindrical fibers by extending our previous work on inextensible bending fibers [Maxian et al., Phys. Rev. Fluids 6 (1), 014102] to fibers with twist elasticity. In our "Euler" model, twist is a scalar function that measures the deviation of the fiber cross section relative to a twist-free frame, the fiber exerts only torque parallel to the centerline on the fluid, and the perpendicular components of the rotational fluid velocity are discarded in favor of the translational velocity. In the first part of this paper, we justify this model by comparing it to another commonly-used "Kirchhoff" formulation where the fiber exerts both perpendicular and parallel torque on the fluid, and the perpendicular angular fluid velocity is required to be consistent with the translational fluid velocity. We then develop a spectral numerical method for the hydrodynamics of the Euler model. We define hydrodynamic mobility operators using integrals of the Rotne-Prager-Yamakawa tensor, and evaluate these integrals through a novel slender-body quadrature, which requires on the order of 10 points along the fiber to obtain several digits of accuracy. We demonstrate that this choice of mobility removes the unphysical negative eigenvalues in the translation-translation mobility associated with asymptotic slender body theories, and ensures strong convergence of the fiber velocity and weak convergence of the fiber constraint forces. We pair the spatial discretization with a semi-implicit temporal integrator to confirm the negligible contribution of twist elasticity to the relaxation dynamics of a bent fiber and study the instability of a twirling fiber.



There are no comments yet.


page 1

page 2

page 3

page 4


An integral-based spectral method for inextensible slender fibers in Stokes flow

Every animal cell is filled with a cytoskeleton, a dynamic gel made of i...

A single-layer based numerical method for the slender body boundary value problem

Fluid flows containing dilute or dense suspensions of thin fibers are wi...

Continuous data assimilation applied to a velocity-vorticity formulation of the 2D Navier-Stokes equations

We study a continuous data assimilation (CDA) algorithm for a velocity-v...

An immersed boundary method with subgrid resolution and improved numerical stability applied to slender bodies in Stokes flow

The immersed boundary method is a numerical and mathematical formulation...

The Fictitious Domain Method Based on Navier Slip Boundary Condition for Simulation of Flow-Particle Interaction

In this article, we develop a least–squares/fictitious domain method for...

A-ULMPM: An Arbitrary Updated Lagrangian Material Point Method for Efficient Simulation of Solids and Fluids

We present an arbitrary updated Lagrangian Material Point Method (A-ULMP...

A novel formulation for the evolution of relativistic rotating stars

We present a new formulation to construct numerically equilibrium config...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.