A continuum active structure model for the interaction of cilia with a viscous fluid

Abstract

Abstract This paper presents a model for a thin active structure interacting with a viscous fluid, as well as a discretization and numerical simulations of the arising fluid‐structure interaction problem. The developed model allows to reproduce the behavior of cilia or flagella immersed in a viscous flow. In the context of linear or nonlinear elasticity, the model is based upon the definition of a suitable internal Piola‐Kirchoff tensor mimicking the action of the internal dyneins that induce the motility of the structure. In the subsequent fluid‐structure interaction problem, two difficulties arise and are discussed: on the one hand the internal activity of the structure leads to more restrictive well‐posedness conditions and, on the other hand, the coupling conditions between the fluid and the structure require a specific numerical treatment. A weak formulation of the time‐discretized problem is derived in functional spaces that include the coupling conditions, but for numerical purposes, an equivalent formulation using Lagrange multipliers is introduced in order to get rid of the constraints in the functional spaces. This new formulation allows for the use of standard (fluid and structure) solvers, up to an iterative procedure. Numerical simulations are presented, including the beating of one or two cilia in 2d, discussing the competition between the magnitude of the internal activity and the viscosity of the surrounding fluid.