On the rectification of oscillatory flows by flexible leaflets in a confined geometry

Énoncé du problème
Les interactions fluide-structure (IFS) dans les canaux compliants et les géométries confinées sont fondamentales pour le transport biologique et la conception microfluidique. Bien que les non-linéarités locales induites par des structures flexibles uniques soient bien caractérisées, prédire le comportement collectif de multiples structures en interaction sous des écoulements variables dans le temps reste un défi important. La littérature existante se concentre principalement sur les régimes stationnaires ou sur des structures uniques, négligeant souvent la richesse des dynamiques temporelles qui émergent de l'interaction de plusieurs éléments compliants. Plus précisément, les mécanismes régissant la manière dont des réseaux de lamelles flexibles redressent les écoulements oscillatoires dans des environnements à faible nombre de Reynolds — brisant la symétrie de renversement du temps inhérente au régime de Stokes pour générer un transport unidirectionnel net — ne sont pas pleinement compris.

Méthodologie
Les auteurs étudient ce problème à l'aide d'un modèle bidimensionnel minimal d'un canal contenant un lit de lamelles asymétriques ancrées à la paroi inférieure. Le système est piloté par l'oscillation verticale de la paroi supérieure, simulant un écoulement de type « squeeze flow ». L'étude emploie une approche numérique entièrement couplée combinant :

1. La méthode de Boltzmann sur réseau (LBM) : Utilisée pour résoudre les équations de Navier-Stokes incompressibles dans la limite de faible nombre de Mach. Pour gérer la nature dépendante du temps du problème sans restrictions de pas de temps prohibitives, un schéma de double pas de temps (DTS) accéléré par une approche multigrid est utilisé.
2. La méthode de la frontière immergée (IBM) : Utilisée pour intégrer la géométrie flexible des lamelles dans le domaine fluide. Les lamelles sont modélisées comme des plaques rigides pivotant autour d'un axe, attachées à un ressort de rotation, répondant aux torques hydrodynamiques.
3. Un cadre analytique : Un modèle de limite continue est développé pour isoler l'influence de l'élasticité et de la densité des lamelles. Ce modèle utilise l'équilibre du couple local et la conservation du moment angulaire global pour prédire le redressement de l'écoulement dans la limite stationnaire, ce qui est ensuite étendu via une approximation quasi-stationnaire au régime oscillatoire.

Contributions clés et résultats
L'étude caractérise la génération d'un transport de fluide net à travers l'interaction entre l'élasticité, la densité des lamelles et la géométrie du canal. Les principales conclusions incluent :

• Rôle du nombre elasto-visqueux ($\eta$) : Le comportement du système est régi par le nombre elasto-visqueux sans dimension, défini comme le rapport entre les forces visqueuses hydrodynamiques et les forces élastiques de rappel des lamelles. Le flux net présente une dépendance non monotone vis-à-vis de $\eta$ :
  • Lamelles rigides ($\eta \ll 1$) : Les déformations sont négligeables, ce qui entraîne un redressement de l'écoulement minimal.
  • Flexibilité optimale : Une valeur intermédiaire de $\eta$ maximise le transport. Ici, les lamelles se déforment suffisamment pour créer une asymétrie géométrique entre les phases de contraction et de relaxation sans perdre leur intégrité structurelle.
  • Lamelles excessivement souples ($\eta \gg 1$) : Les lamelles se déforment excessivement (dépassant souvent 90 degrés), réduisant l'asymétrie géométrique et diminuant le transport net.
• Interactions collectives et densité : Les densités de lamelles élevées maximisent les interactions collectives et le transport net. À mesure que la densité augmente, le réseau discret de lamelles approche une limite continue où le profil de la frontière effective devient lisse. Le débit net augmente avec la densité jusqu'à saturation, un comportement capturé par une fonction d'échelle phénoménologique $Q(\eta, \phi) \sim \phi/(\phi + K(\eta))$.
• Régimes stationnaires vs dynamiques : Dans la limite stationnaire, le flux net est indépendant du nombre de Reynolds et de la fréquence d'oscillation, dépendant uniquement de $\eta$ et de la densité. Cependant, dans le régime pleinement oscillatoire, le système présente un déphasage entre le mouvement de la paroi et la réponse des lamelles. Ce retard dynamique est contrôlé par un rapport de l'échelle de temps sans dimension $T_r$, représentant le rapport entre le temps de relaxation visqueuse de la lamelle et la période d'oscillation de la paroi. Lorsque $T_r$ est significatif, des configurations transitoires apparaissent, ce qui diminue le transport redressé net, un phénomène non capturé par les modèles purement quasi-statiques.
• Dépendance géométrique : Bien que le modèle quasi-stationnaire prédise avec précision l'écoulement dans la limite de haute densité, il ne parvient pas à capturer la dépendance du flux net vis-à-vis de la longueur du canal ($L_x$) observée dans les simulations, soulignant l'importance des effets dynamiques dans les canaux courts.

Signification
Cet article établit un cadre fondamental pour analyser comment des structures élancées collectives interagissent dynamiquement au sein d'environnements visqueux. En développant un modèle de continuum analytique validé par des simulations IFS entièrement couplées, les auteurs fournissent un outil prédictif pour comprendre le redressement de l'écoulement piloté par des réseaux flexibles. Ce travail démontre que les interactions collectives peuvent briser les symétries du régime de Stokes pour générer un transport unidirectionnel, offrant des perspectives pertinentes pour les systèmes de transport de fluides biologiques et la conception de redresseurs microfluidiques passifs. L'identification d'un nombre elasto-visqueux optimal et la caractérisation des déphasages dynamiques fournissent des paramètres de conception spécifiques pour contrôler l'écoulement dans les géométries compliantes.