# Efficient Implementation of Elastohydrodynamics via Integral Operators

## TL;DR

This paper introduces an efficient integral operator-based method for simulating the complex dynamics and interactions of flexible filaments in fluid flows, addressing computational challenges and enabling new biological and physical insights.

## Contribution

The manuscript develops a computationally efficient integral equation framework using regularized stokeslets for simulating multiple interacting filaments, including active and passive cases.

## Key findings

- Filament interactions can promote buckling instabilities.
- The method accurately models filament dynamics in various flow conditions.
- Open-source MATLAB implementation facilitates further research.

## Abstract

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.