Efficient implementation of elastohydrodynamics via integral operators

A. L. Hall-McNair, T. D. Montenegro-Johnson, H. Gadêlha, D. J. Smith, and M. T. Gallagher
Phys. Rev. Fluids 4, 113101 – Published 12 November 2019
PDFHTMLExport Citation

Abstract

The dynamics of geometrically nonlinear flexible filaments play an important role in a host of biological processes, from flagella-driven cell transport to the polymeric structure of complex fluids. Such problems have historically been computationally expensive due to numerical stiffness associated with the inextensibility constraint, as well as the often nontrivial boundary conditions on the governing high-order PDEs. Formulating the problem for the evolving shape of a filament via an integral equation in the tangent angle has recently been found to greatly alleviate this numerical stiffness. The contribution of the present manuscript is to enable the simulation of nonlocal interactions of multiple filaments in a computationally efficient manner using the method of regularized stokeslets within this framework. The proposed method is benchmarked against a nonlocal bead and link model, and recent code utilizing a local drag velocity law. Systems of multiple filaments (1) in a background fluid flow, (2) under a constant body force, and (3) undergoing active self-motility are modeled efficiently. Buckling instabilities are analyzed by examining the evolving filament curvature, as well as by coarse graining the body frame tangent angles using a Chebyshev approximation for various choices of the relevant nondimensional parameters. From these experiments, insight is gained into how filament-filament interactions can promote buckling, and further reveal the complex fluid dynamics resulting from arrays of these interacting fibers. By examining active moment-driven filaments, we investigate the speed of worm- and spermlike swimmers for different governing parameters. The matlab® implementation is made available as an open-source library, enabling flexible extension for alternate discretizations and different surrounding flows.

  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
9 More
  • Received 8 March 2019

DOI:https://doi.org/10.1103/PhysRevFluids.4.113101

©2019 American Physical Society

Physics Subject Headings (PhySH)

Fluid DynamicsPhysics of Living Systems

Authors & Affiliations

A. L. Hall-McNair1, T. D. Montenegro-Johnson1, H. Gadêlha2, D. J. Smith1, and M. T. Gallagher1,*

  • 1School of Mathematics, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom
  • 2Faculty of Engineering, University of Bristol, Bristol BS8 1UB, United Kingdom

  • *Corresponding author: m.t.gallagher@bham.ac.uk

Article Text (Subscription Required)

Click to Expand

Supplemental Material (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 4, Iss. 11 — November 2019

Reuse & Permissions
Access Options
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×
×

Images

×

Sign up to receive regular email alerts from Physical Review Fluids

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×