3D Simulation of Active thin Structures in a Viscous Fluid and Application to Mucociliary Transport

¹ Université de Bordeaux, Institut de Mathématiques de Bordeaux (CNRS-UMR 5251), 351 cours de la Libération, 33405 Talence, France
² Université Paris Cité, Laboratoire MAP5 (CNRS-UMR 8145), 45 rue des Saints-Pères, 75270 Paris cedex 06, France
³ Université Paris-Saclay, Laboratoire de Mathématiques d’Orsay (CNRS-UMR 8628), Bâtiment 307, rue Michel Magat, Faculté des Sciences d’Orsay, 91405 Orsay cedex, France

^* Corresponding author: sebastien.martin@u-paris.fr

Received: 27 October 2023
Accepted: 17 May 2024

Abstract

We propose a hierarchy of mathematical models for the numerical simulation of active thin structures in a viscous fluid and its application to mucociliary clearance. Our aim is to simulate large forests of cilia and analyze the collective dynamics arising in the flow, as well as their impact on the efficiency of the mucus transport. In a 3D model we describe the cilia individually and study their joint actions on the fluid. The model is built upon a 3D Stokes problem with singular source terms that represent the action of the 1D cilia on the fluid, including the influence of the background flow (making the problem nonlocal). Surface tension between the periciliary layer and the mucus is taken into account. From the 3D model we also derive a 1D space averaged model, describing the dynamics of the mean velocity of the mucus that is propelled by the cilia, hence allowing lower computational costs and still providing useful characterization of the efficiency of the transport. Mathematical properties of the models (existence and uniqueness of solutions in suitable functional spaces) are analyzed. Numerical simulations highlight the influence of critical parameters on the efficiency of the mucociliary transport in the case of dense forests of cilia.