Efficient implementation of elastohydrodynamics via integral operators

The dynamics of geometrically non-linear 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 non-trivial 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 non-local 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 non-local 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 non-dimensional 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 sperm-like swimmers for different governing parameters. The MATLAB(R) implementation is made available as an open-source library, enabling flexible extension for alternate discretizations and different surrounding flows.

[1]  H. Gadêlha On the optimal shape of magnetic swimmers , 2013 .

[2]  G. G. Stokes "J." , 1890, The New Yale Book of Quotations.

[3]  D. Bates,et al.  Direct detection and measurement of wall shear stress using a filamentous bio-nanoparticle , 2015, Nano Research.

[4]  Michael Shelley,et al.  Simulating the dynamics and interactions of flexible fibers in Stokes flows , 2004 .

[5]  C. Brokaw,et al.  Bend propagation by a sliding filament model for flagella. , 1971, The Journal of experimental biology.

[6]  E A Gaffney,et al.  Bend propagation in the flagella of migrating human sperm, and its modulation by viscosity. , 2009, Cell motility and the cytoskeleton.

[7]  Jackson Kirkman-Brown,et al.  Human spermatozoa migration in microchannels reveals boundary-following navigation , 2012, Proceedings of the National Academy of Sciences.

[8]  David B. Stein,et al.  Coarse graining the dynamics of immersed and driven fiber assemblies , 2019, Physical Review Fluids.

[9]  G. J. HANCOCKf,et al.  THE PROPULSION OF SEA-URCHIN SPERMATOZOA , 2005 .

[10]  R. Adhikari,et al.  Autonomous motility of active filaments due to spontaneous flow-symmetry breaking. , 2012, Physical review letters.

[11]  C. Brokaw,et al.  Computer simulation of flagellar movement. I. Demonstration of stable bend propagation and bend initiation by the sliding filament model. , 1972, Biophysical journal.

[12]  S. Chu,et al.  Effect of hydrodynamic interactions on DNA dynamics in extensional flow: Simulation and single molecule experiment , 2004 .

[13]  Franck Plouraboué,et al.  A general formulation of Bead Models applied to flexible fibers and active filaments at low Reynolds number , 2015, J. Comput. Phys..

[14]  Anke Lindner,et al.  Morphological transitions of elastic filaments in shear flow , 2018, Proceedings of the National Academy of Sciences.

[15]  T. Ishikawa,et al.  Nodal cilia-driven flow: Development of a computational model of the nodal cilia axoneme. , 2017, Journal of biomechanics.

[16]  L. Fauci,et al.  Dynamics of a macroscopic elastic fibre in a polymeric cellular flow , 2017, Journal of Fluid Mechanics.

[17]  E. Gaffney,et al.  The counterbend phenomenon in flagellar axonemes and cross-linked filament bundles , 2013, Proceedings of the National Academy of Sciences.

[18]  J. Higdon,et al.  A spectral boundary element approach to three-dimensional Stokes flow , 1995, Journal of Fluid Mechanics.

[19]  H. Gadêlha,et al.  The counterbend dynamics of cross-linked filament bundles and flagella , 2017, Journal of The Royal Society Interface.

[20]  C. Lindemann A model of flagellar and ciliary functioning which uses the forces transverse to the axoneme as the regulator of dynein activation. , 1994, Cell motility and the cytoskeleton.

[21]  Ricardo Cortez,et al.  Regularized Stokeslet segments , 2018, J. Comput. Phys..

[22]  Ricardo Cortez,et al.  Enhanced flagellar swimming through a compliant viscoelastic network in Stokes flow , 2016, Journal of Fluid Mechanics.

[23]  D. Woolley,et al.  A study of helical and planar waves on sea urchin sperm flagella, with a theory of how they are generated. , 2001, The Journal of experimental biology.

[24]  S. F. Schoeller,et al.  Methods for suspensions of passive and active filaments , 2019, J. Comput. Phys..

[25]  H. Gadêlha,et al.  The asymptotic coarse-graining formulation of slender-rods, bio-filaments and flagella , 2017, Journal of The Royal Society Interface.

[26]  D. Smith,et al.  Modelling the fluid mechanics of cilia and flagella in reproduction and development , 2012, The European physical journal. E, Soft matter.

[27]  J. Gower Generalized procrustes analysis , 1975 .

[28]  W. Marsden I and J , 2012 .

[29]  D J Smith,et al.  Model-based image analysis of a tethered Brownian fibre for shear stress sensing , 2017, Journal of The Royal Society Interface.

[30]  David J. Smith,et al.  Meshfree and efficient modelling of swimming cells , 2018, 1801.04142.

[31]  J J Blum,et al.  Bend propagation in flagella. I. Derivation of equations of motion and their simulation. , 1978, Biophysical journal.

[32]  Lawrence F. Shampine,et al.  The MATLAB ODE Suite , 1997, SIAM J. Sci. Comput..

[33]  L. E. Becker,et al.  Instability of elastic filaments in shear flow yields first-normal-stress differences. , 2001, Physical review letters.

[34]  E. Gaffney,et al.  Nonlinear instability in flagellar dynamics: a novel modulation mechanism in sperm migration? , 2010, Journal of The Royal Society Interface.

[35]  L. Fauci,et al.  The method of regularized Stokeslets in three dimensions : Analysis, validation, and application to helical swimming , 2005 .

[36]  Thomas D. Montenegro-Johnson,et al.  Microtransformers: Controlled microscale navigation with flexible robots , 2018, Physical Review Fluids.

[37]  D. Smith,et al.  A boundary element regularized Stokeslet method applied to cilia- and flagella-driven flow , 2009, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[38]  global sci A Treecode Algorithm for 3D Stokeslets and Stresslets , 2019 .

[39]  Y. Young Hydrodynamic interactions between two semiflexible inextensible filaments in Stokes flow. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[40]  D. J. Smith,et al.  Spermatozoa scattering by a microchannel feature: an elastohydrodynamic model , 2014, Royal Society Open Science.

[42]  J. Casademunt,et al.  Nonlinear amplitude dynamics in flagellar beating , 2016, Royal Society Open Science.

[43]  L. Fauci,et al.  Hydrodynamic interactions of sheets vs filaments: Synchronization, attraction, and alignment , 2015 .

[44]  M. Hicks,et al.  Detection of pathogenic bacteria using a homogeneous immunoassay based on shear alignment of virus particles and linear dichroism. , 2012, Analytical chemistry.

[45]  M Cosentino Lagomarsino,et al.  Hydrodynamic induced deformation and orientation of a microscopic elastic filament. , 2005, Physical review letters.

[46]  Svetlana Tlupova,et al.  A Treecode Algorithm for 3D Stokeslets and Stresslets , 2018, Advances in Applied Mathematics and Mechanics.

[47]  M. Shelley,et al.  Bistability in the synchronization of actuated microfilaments , 2017, Journal of Fluid Mechanics.

[48]  J. Kirkman-Brown,et al.  CASA: tracking the past and plotting the future. , 2018, Reproduction, fertility, and development.

[49]  Franck Plouraboué,et al.  Large-scale simulation of steady and time-dependent active suspensions with the force-coupling method , 2015, J. Comput. Phys..

[50]  Ricardo Cortez,et al.  The Method of Regularized Stokeslets , 2001, SIAM J. Sci. Comput..

[51]  F. Jülicher,et al.  Physical limits of flow sensing in the left-right organizer , 2017, eLife.