skip to main content

Title: Slender body theories for rotating filaments
Slender fibres are ubiquitous in biology, physics and engineering, with prominent examples including bacterial flagella and cytoskeletal fibres. In this setting, slender body theories (SBTs), which give the resistance on the fibre asymptotically in its slenderness $\epsilon$ , are useful tools for both analysis and computations. However, a difficulty arises when accounting for twist and cross-sectional rotation: because the angular velocity of a filament can vary depending on the order of magnitude of the applied torque, asymptotic theories must give accurate results for rotational dynamics over a range of angular velocities. In this paper, we first survey the challenges in applying existing SBTs, which are based on either singularity or full boundary integral representations, to rotating filaments, showing in particular that they fail to consistently treat rotation–translation coupling in curved filaments. We then provide an alternative approach which approximates the three-dimensional dynamics via a one-dimensional line integral of Rotne–Prager–Yamakawa regularized singularities. While unable to accurately resolve the flow field near the filament, this approach gives a grand mobility with symmetric rotation–translation and translation–rotation coupling, making it applicable to a broad range of angular velocities. To restore fidelity to the three-dimensional filament geometry, we use our regularized singularity model to inform a simple empirical equation which relates the mean force and torque along the filament centreline to the translational and rotational velocity of the cross-section. The single unknown coefficient in the model is estimated numerically from three-dimensional boundary integral calculations on a rotating, curved filament.  more » « less
Award ID(s):
2052515 1646339
Author(s) / Creator(s):
Date Published:
Journal Name:
Journal of Fluid Mechanics
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This paper presents a theory to obtain the force per unit length acting on a slender filament with a non-circular cross-section moving in a fluid at low Reynolds number. Using a regular perturbation of the inner solution, we show that the force per unit length has $O(1/\ln (2A))+O(\unicode[STIX]{x1D6FC}/\ln ^{2}(2A))$ contributions driven by the relative motion of the particle and the local fluid velocity and an $O(\unicode[STIX]{x1D6FC}/(\ln (2A)A))$ contribution driven by the gradient in the imposed fluid velocity. Here, the aspect ratio ( $A=l/a_{0}$ ) is defined as the ratio of the particle size ( $l$ ) to the cross-sectional dimension ( $a_{0}$ ) and $\unicode[STIX]{x1D6FC}$ is the amplitude of the non-circular perturbation. Using thought experiments, we show that two-lobed and three-lobed cross-sections affect the response to relative motion and velocity gradients, respectively. A two-dimensional Stokes flow calculation is used to extend the perturbation analysis to cross-sections that deviate significantly from a circle (i.e. $\unicode[STIX]{x1D6FC}\sim O(1)$ ). We demonstrate the ability of our method to accurately compute the resistance to translation and rotation of a slender triaxial ellipsoid. Furthermore, we illustrate novel dynamics of straight rods in a simple shear flow that translate and rotate quasi-periodically if they have two-lobed cross-section, and rotate chaotically and translate diffusively if they have a combination of two- and three-lobed cross-sections. Finally, we show the remarkable ability of our theory to accurately predict the motion of rings, retaining great accuracy for moderate aspect ratios ( ${\sim}10$ ) and cross-sections that deviate significantly from a circle, thereby making our theory a computationally inexpensive alternative to other Stokes flow solvers. 
    more » « less
  2. Dynamic organization of the cytoskeletal filaments and rod-like proteins in the cell membrane and other biological interfaces occurs in many cellular processes. Previous modeling studies have considered the dynamics of a single rod on fluid planar membranes. We extend these studies to the more physiologically relevant case of a single filament moving in a spherical membrane. Specifically, we use a slender-body formulation to compute the translational and rotational resistance of a single filament of length L moving in a membrane of radius R and 2D viscosity ηm, and surrounded on its interior and exterior with Newtonian fluids of viscosities η− and η+. We first discuss the case where the filament's curvature is at its minimum κ=1/R. We show that the boundedness of spherical geometry gives rise to flow confinement effects that increase in strength with increasing the ratio of filament's length to membrane radius L/R. These confinement flows only result in a mild increase in filament's resistance along its axis, ξ∥, and its rotational resistance, ξΩ. As a result, our predictions of ξ∥ and ξΩ can be quantitatively mapped to the results on a planar membrane. In contrast, we find that the drag in perpendicular direction, ξ⊥, increases superlinearly with the filament's length, when L/R>1 and ultimately ξ⊥→∞ as L/R→π. Next, we consider the effect of the filament's curvature, κ, on its parallel motion, while fixing the membrane's radius. We show that the flow around the filament becomes increasingly more asymmetric with increasing its curvature. These flow asymmetries induce a net torque on the filament, coupling its parallel and rotational dynamics. This coupling becomes stronger with increasing L/R and κ. 
    more » « less
  3. Semiflexible slender filaments are ubiquitous in nature and cell biology, including in the cytoskeleton, where reorganization of actin filaments allows the cell to move and divide. Most methods for simulating semiflexible inextensible fibers/polymers are based on discrete (bead-link or blob-link) models, which become prohibitively expensive in the slender limit when hydrodynamics is accounted for. In this paper, we develop a novel coarse-grained approach for simulating fluctuating slender filaments with hydrodynamic interactions. Our approach is tailored to relatively stiff fibers whose persistence length is comparable to or larger than their length and is based on three major contributions. First, we discretize the filament centerline using a coarse non-uniform Chebyshev grid, on which we formulate a discrete constrained Gibbs–Boltzmann (GB) equilibrium distribution and overdamped Langevin equation for the evolution of unit-length tangent vectors. Second, we define the hydrodynamic mobility at each point on the filament as an integral of the Rotne–Prager–Yamakawa kernel along the centerline and apply a spectrally accurate “slender-body” quadrature to accurately resolve the hydrodynamics. Third, we propose a novel midpoint temporal integrator, which can correctly capture the Ito drift terms that arise in the overdamped Langevin equation. For two separate examples, we verify that the equilibrium distribution for the Chebyshev grid is a good approximation of the blob-link one and that our temporal integrator for overdamped Langevin dynamics samples the equilibrium GB distribution for sufficiently small time step sizes. We also study the dynamics of relaxation of an initially straight filament and find that as few as 12 Chebyshev nodes provide a good approximation to the dynamics while allowing a time step size two orders of magnitude larger than a resolved blob-link simulation. We conclude by applying our approach to a suspension of cross-linked semiflexible fibers (neglecting hydrodynamic interactions between fibers), where we study how semiflexible fluctuations affect bundling dynamics. We find that semiflexible filaments bundle faster than rigid filaments even when the persistence length is large, but show that semiflexible bending fluctuations only further accelerate agglomeration when the persistence length and fiber length are of the same order. 
    more » « less
  4. Abstract

    The distortion of the charge cloud around a uniformly charged, dielectric, rigid sphere that translates and rotates in an unbounded binary, symmetric electrolyte at zero Reynolds number is examined. The zeta potential of the particle ζ is assumed small relative to the thermal voltage scale. It is assumed that the equilibrium structure of the cloud is slightly distorted, which requires that the Péclet numbers characterizing distortion due to particle translation,, and rotation,, are small compared to unity. Here,ais radius of the particle;Dis the ionic diffusion coefficient;and, whereUandΩare the rectilinear and angular velocities of the particle, respectively. Perturbation expansions for smallandare employed to calculate the nonequilibrium structure of the cloud, whence the force and torque on the particle are determined. In particular, we predict that the sphere experiences a force orthogonal to its directions of translation and rotation. This “lift” force arises from the nonlinear distortion of the cloud under the combined actions of particle translation and rotation. The lift force is given by. Here, ε is the permittivity of the electrolyte;is the Debye length; andis a negative function that decreases in magnitude with increasing. The lift force implies that an unconstrained particle would follow a curved path; an electrokinetic analog of the inertial Magnus effect. Finally, the implication of the lift force on cross‐streamline migration of an electrophoretic particle in shear flow is discussed.

    more » « less
  5. Abstract

    We explore the effects of rapid rotation on the properties of neutrino-heated winds from proto-neutron stars (PNS) formed in core-collapse supernovae or neutron-star mergers by means of three-dimensional general-relativistic hydrodynamical simulations with M0 neutrino transport. We focus on conditions characteristic of a few seconds into the PNS cooling evolution when the neutrino luminosities obeyLνe+Lν¯e7×1051erg s−1, and over which most of the wind mass loss will occur. After an initial transient phase, all of our models reach approximately steady-state outflow solutions with positive energies and sonic surfaces captured on the computational grid. Our nonrotating and slower rotating models (angular velocity relative to Keplerian Ω/ΩK≲ 0.4; spin periodP≳ 2 ms) generate approximately spherically symmetric outflows with properties in good agreement with previous PNS wind studies. By contrast, our most rapidly spinning PNS solutions (Ω/ΩK≳ 0.75;P≈ 1 ms) generate outflows focused in the rotational equatorial plane with much higher mass-loss rates (by over an order of magnitude), lower velocities, lower entropy, and lower asymptotic electron fractions, than otherwise similar nonrotating wind solutions. Although such rapidly spinning PNS are likely rare in nature, their atypical nucleosynthetic composition and outsized mass yields could render them important contributors of light neutron-rich nuclei compared to more common slowly rotating PNS birth. Our calculations pave the way to including the combined effects of rotation and a dynamically important large-scale magnetic field on the wind properties within a three-dimensional GRMHD framework.

    more » « less