TY - JOUR
T1 - Meshfree and efficient modelling of swimming cells
AU - Gallagher, Meurig Thomas
AU - Smith, David
PY - 2018/5/31
Y1 - 2018/5/31
N2 - Locomotion in Stokes flow is an intensively-studied problem because it describes important bi- ological phenomena such as the motility of many species’ sperm, bacteria, algae and protozoa. Numerical computations can be challenging, particularly in three dimensions, due to the presence of moving boundaries and complex geometries; methods which combine ease-of-implementation and computational efficiency are therefore needed. A recently-proposed method to discretise the regularised Stokeslet boundary integral equation without the need for a connected ‘mesh’ is ap- plied to the inertialess locomotion problem in Stokes flow. The mathematical formulation and key aspects of the computational implementation in Matlab® /GNU Octave are described, followed by numerical experiments with biflagellate algae and multiple uniflagellate sperm swimming between no-slip surfaces, for which both swimming trajectories and flow fields are calculated. These compu- tational experiments required minutes of time on modest hardware; an extensible implementation is provided in a github repository. The nearest neighbour discretisation dramatically improves convergence and robustness, a key challenge in extending the regularised Stokeslet method to complicated, three dimensional, biological fluid problems.
AB - Locomotion in Stokes flow is an intensively-studied problem because it describes important bi- ological phenomena such as the motility of many species’ sperm, bacteria, algae and protozoa. Numerical computations can be challenging, particularly in three dimensions, due to the presence of moving boundaries and complex geometries; methods which combine ease-of-implementation and computational efficiency are therefore needed. A recently-proposed method to discretise the regularised Stokeslet boundary integral equation without the need for a connected ‘mesh’ is ap- plied to the inertialess locomotion problem in Stokes flow. The mathematical formulation and key aspects of the computational implementation in Matlab® /GNU Octave are described, followed by numerical experiments with biflagellate algae and multiple uniflagellate sperm swimming between no-slip surfaces, for which both swimming trajectories and flow fields are calculated. These compu- tational experiments required minutes of time on modest hardware; an extensible implementation is provided in a github repository. The nearest neighbour discretisation dramatically improves convergence and robustness, a key challenge in extending the regularised Stokeslet method to complicated, three dimensional, biological fluid problems.
KW - cell locomotion
KW - biological fluid dynamics
KW - collective behaviour
KW - flows in microchannels
KW - swiming
U2 - 10.1103/PhysRevFluids.3.053101
DO - 10.1103/PhysRevFluids.3.053101
M3 - Article
SN - 2469-990X
VL - 3
JO - Physical Review Fluids
JF - Physical Review Fluids
IS - 5
M1 - 053101
ER -