Deformable bodies in a 3-dimensional viscous flow: Vorticity-Stream vector formulation

Résumé Technique : Corps Déformables en Écoulement Visqueux 3D via la Formulation Vecteur Vorticité-Flux

Énoncé du Problème
La simulation de l'interaction entre des écoulements visqueux incompressibles tridimensionnels (3D) et des obstacles déformables ou élastiques (tels que des vésicules, des gouttelettes et des cellules) présente des défis importants. Les méthodes existantes, incluant les méthodes de Boltzmann sur réseau (LBM), les techniques d'intégrale de bord et les approches basées sur les particules comme la dynamique des particules dissipatives lissées (SDPD), souffrent souvent d'une complexité computationnelle élevée, de mises en œuvre numériques complexes ou de difficultés à gérer les mécaniques d'interface élaborées. Bien que les simulations 2D utilisant des formulations vorticité-fonction de courant couplées à des modèles de champ de phase aient prouvé leur robustesse, un cadre comparable, léger et fondé sur la théorie classique pour les interfaces déformables 3D régies par des énergies de type Canham–Hérlif hers (flexion de membrane) ou Cahn–Hilliard (tension superficielle) manquait à l'appel. Plus précisément, il existe une lacune dans la littérature concernant les formulations basées sur la vorticité pour les interfaces 3D.

Méthodologie
Les auteurs introduisent une nouvelle formulation numérique basée sur l'approche vecteur vorticité-flux, spécifiquement adaptée aux écoulements à faible nombre de Reynolds. La méthode généralise la fonction de courant 2D en un potentiel vecteur 3D ($\vec{\psi}$), où le champ de vitesse est défini par $\vec{v} = \nabla \times \vec{\psi}$.

Le cœur de la méthodologie consiste à résoudre un système couplé d'équations pour trois observables primaires :

1. Potentiel Chimique de l'Interface ($\mu$) : Dérivé de la dérivée fonctionnelle d'une énergie libre choisie ($F$) par rapport au paramètre d'ordre du champ de phase ($\phi$). Les auteurs démontrent le cadre en utilisant à la fois l'énergie de Canham–Hérlif hers (pour la flexion de la membrane) et l'énergie de Cahn–Hilliard (pour la tension superficielle).
2. Vorticité du Fluide ($\vec{\omega}$) : Définie par $\vec{\omega} = \nabla \times \vec{v}$.
3. Vecteur Flux ($\vec{\psi}$) : Le potentiel vecteur satisfaisant la condition de gauge $\nabla \cdot \vec{\psi} = 0$.

Les équations directrices se réduisent à un ensemble de problèmes de Poisson couplés :

• L'évolution du champ de phase $\phi$ est régie par une équation d'advection-diffusion impliquant le potentiel chimique et le champ de vitesse.
• L'équation de la vorticité est dérivée du rotationnel de l'équation de Stokes (ou de Navier-Stokes pour les extensions inertielles), résultant en une équation de Poisson où le terme source dépend du produit vectoriel du gradient de phase et du gradient de potentiel chimique : $\nabla^2 \vec{\omega} = \frac{1}{\eta} (\nabla \phi) \times (\nabla \mu)$.
• Le vecteur flux est obtenu en résolvant $\nabla^2 \vec{\psi} = -\vec{\omega}$.

Implémentation Numérique
Le système est résolu à l'aide d'une méthode de différences finies sur une grille cubique fixe et non ajustée (unfitted). L'interface est traitée comme une région diffusive via l'approche de champ de phase, éliminant ainsi le besoin de déformation complexe de maillage ou de techniques de frontière immergée. L'évolution temporelle utilise une intégration d'Euler explicite pour le champ de phase et un solveur itératif pour les équations de Poisson, utilisant un schéma ségrégué (séquentiel et décalé). Le code est implémenté en Python et est disponible sous forme de dépôt open-source.

Contributions Clés

• Simplification des Solveurs 3D : La formulation transforme le solveur de fluide en un ensemble de problèmes de Poisson couplés, évitant le calcul direct de la pression et garantissant une incompressibilité stricte par construction.
• Cadre Unifié pour les Énergies d'Interface : La méthode intègre de manière transparente différents modèles d'énergie libre (Canham–Hérlif hers pour la flexion et Cahn–Hilliard pour la tension superficielle) au sein de la même structure numérique, permettant une personnalisation aisée de la mécanique des membranes.
• Observables Diagnostiques : Les auteurs proposent de nouveaux observables basés sur les fluctuations quadratiques des champs de vorticité et de vecteur flux par rapport à l'écoulement non perturbé. Ces métriques se sont révélées très sensibles aux déformations de forme et aux transitions morphologiques.

Résultats
La méthode a été validée par des simulations de vésicules et de gouttelettes uniques dans des écoulements de Poiseuille et de Couette newtoniens :

• Écoulement de Poiseuille : Les simulations ont récupéré avec succès les formes axisymétriques canoniques des globules rouges (RBC), incluant les morphologies en « parachute » et en « slipper », conformément à la littérature. Les observables diagnostiques ont suivi efficacement la dynamique de relaxation, montrant des pics distincts correspondant aux changements de forme (par exemple, la disparition d'une fossette) et se stabilisant une fois que l'objet a atteint un état stationnaire.
• Écoulement de Couette : Le cadre a capturé la migration latérale, un phénomène de portance non inertiel où les vésicules migrent perpendiculairement à la direction de l'écoulement. Les simulations ont démontré que la position d'équilibre et la dynamique de migration dépendent systématiquement du contraste de viscosité ($\zeta$) entre les fluides interne et externe. Des contrastes de viscosité plus élevés ont entraîné un dérive plus prononcée vers la paroi mobile.
• Analyse de Contrainte : La méthode a calculé avec succès les distributions de contraintes de cisaillement dans le plan, révélant que des taux de cisaillement plus élevés induisent une plus grande hétérogénéité de contrainte et peuvent mener à des instabilités morphologiques.

Signification et Revendications
L'article affirme que cette formulation offre une alternative légère, flexible et efficace sur le plan computationnel aux solveurs existants mésoscopiques ou à haute complexité pour l'hydrodynamique 3D à faible nombre de Reynolds. En s'appuyant sur les principes classiques de la mécanique des fluides plutôt que sur la théorie cinétique (comme dans la LBM), la méthode évite le besoin de preuves de cohérence avec la mécanique des milieux continus.

Les auteurs soulignent que bien que le présent travail se concentre sur la dynamique d'un corps unique dans des fluides newtoniens, le cadre est conçu pour être extensible. Il peut naturellement accommoder des forces de corps, des effets inertiels, des milieux viscoélastiques et des interactions multi-corps. La capacité à capturer la dynamique de cisaillement dans le plan en 3D — une caractéristique souvent perdue dans les approximations 2D — positionne cette méthode comme un outil précieux pour l'étude de la biomécanique cellulaire et de la physique de la matière molle. La nature open-source du code est destinée à favoriser l'étalonnage et le développement futur dans ce domaine émergent.