Adaptation of an existing filament model as a two-wall sliding filament model, using spring forces to model passive links and active motor protein links as seen in biological axonemal architecture of flagellum and cilia.
Original repository (accessed 1 October 2019): https://github.com/ekeaveny/filaments/tree/master/matlab-2d-rpy
This original code demonstrates the use of the method described here in simulating a single flexible filament falling under gravity in an infinite domain.
It uses the 'EJBb' version of Broyden's method (Algorithm 2 in the paper) with a reduced 'robot arm' system of nonlinear equations. For the hydrodynamic solver, it uses the RPY tensors.
In this version, external forces such as gravity and internal forces such as motor protein links have been adjusted to replicate sliding filament models seen in the literature.