A hybrid Volume of Fluid Phase-Field method for Direct Numerical Simulations of soluble surfactant-laden interfacial flows

Ce papier présente une méthode hybride Volume-of-Fluid Phase-Field avec raffinement adaptatif de maillage pour les simulations numériques directes d'écoulements chargés en tensioactifs solubles, qui capture avec précision le couplage entre le transport en volume et à l'interface afin de démontrer comment les contraintes de Marangoni modifient significativement la dynamique de montée de bulles dans des géométries tridimensionnelles.

L'essentiel

La vue d'ensemble : Une nouvelle façon de simuler l'eau « savonneuse »

Imaginez que vous observez une bulle monter à travers un verre d'eau. Si l'eau est pure, la bulle monte droit, en oscillant un peu mais rapidement. Mais si vous ajoutez du savon (des tensioactifs) à l'eau, la bulle se comporte différemment. Elle peut ralentir, osciller davantage, ou même changer de trajectoire.

Cela se produit parce que les molécules de savon adorent se coller à la surface de la bulle. Lorsque la bulle se déplace, ces molécules sont poussées, créant une « tension » inégale sur la peau de la bulle. Cette tension inégale agit comme des mains invisibles poussant la bulle dans différentes directions, modifiant ainsi son mouvement.

Le Problème :
Simuler cela sur un ordinateur est incroyablement difficile. C'est comme essayer de filmer une bulle de savon avec un appareil photo ayant deux paramètres contradictoires :

1. L'appareil photo net : Il doit voir le bord de la bulle comme une ligne aussi fine qu'un rasoir (pour calculer la pression et la forme).
2. L'appareil photo flou : Il doit voir les molécules de savon s'étaler de manière lisse sur ce bord (pour calculer le mouvement du savon).

La plupart des méthodes informatiques vous obligent à choisir un seul paramètre, rendant la simulation soit physiquement inexacte, soit impossible à exécuter pour des formes 3D complexes.

La Solution :
Les auteurs de ce document ont créé une méthode hybride. Imaginez une simulation en « écran divisé » qui utilise le meilleur des deux mondes simultanément :

• Le bord net (Volume-of-Fluid) : Ils utilisent une méthode qui maintient le bord de la bulle net et conserve parfaitement la quantité de liquide (comme un contour haute définition).
• Le savon lisse (Phase-Field) : Ils utilisent une seconde couche « floue » qui agit comme une autoroute lisse pour le déplacement des molécules de savon. Cela permet au savon de se déplacer naturellement entre l'eau et la surface de la bulle sans rester coincé ni se perdre.

Comment cela fonctionne : L'analogie du « contrôleur de trafic »

Pour rendre cela fonctionnel, les auteurs ont créé un système de trafic numérique pour les molécules de savon :

1. L'autoroute (L'interface) : La surface de la bulle est une autoroute bondée. Les molécules de savon sont les voitures.
2. Les bretelles d'entrée et de sortie (Adsorption/Désorption) :
   • Adsorption : Les molécules de savon provenant de l'eau (le volume) veulent sauter sur l'autoroute (la surface de la bulle).
   • Désorption : Les molécules de savon se fatiguent et sautent de l'autoroute pour retourner dans l'eau.
   • La nouvelle méthode calcule exactement combien de voitures montent ou descendent à chaque instant, garantissant qu'aucune voiture ne disparaît ni n'apparaît de nulle part.
3. L'embouteillage (Contrainte de Marangoni) : Lorsque trop de voitures de savon s'accumulent à un endroit précis sur la bulle, cet endroit devient « collant » (tension élevée). La peau de la bulle tente de s'éloigner de ce point collant, créant une force qui ralentit la bulle ou la fait osciller. La simulation capture parfaitement ce jeu de tiraillement.

Ce qu'ils ont testé (Les « tests de conduite »)

Avant de laisser leur nouvelle voiture rouler sur l'autoroute, ils l'ont fait passer par une école de conduite avec trois tests spécifiques :

1. Le test d'étirement (Sphère en expansion) : Ils ont gonflé une bulle recouverte de savon. Ils ont vérifié si le savon se répartissait uniformément à mesure que la bulle grossissait. La simulation correspondait parfaitement aux mathématiques.
2. Le test de rotation (Bulle en rotation) : Ils ont fait tourner une bulle avec du savon dessus. Ils ont vérifié si le savon se déplaçait correctement autour du cercle sans fuite. Encore une fois, la simulation était parfaite.
3. Le test d'échange (Mur plat) : Ils ont observé le déplacement du savon de l'eau vers un mur plat et inversement. Ils ont testé trois scénarios :
   • Uniquement montée : Le savon colle-t-il ? Oui.
   • Uniquement descente : Le savon part-il ? Oui.
   • Les deux : Trouve-t-il un équilibre ? Oui.

L'événement principal : La bulle montante

Enfin, ils ont laissé leur nouvelle méthode simuler une bulle montant dans un réservoir d'eau en 3D.

• La bulle propre : Elle montait relativement vite et droit.
• La bulle de savon « insoluble » : Le savon était coincé à la surface et ne pouvait pas partir. Cela créait un fort « embouteillage » à l'arrière de la bulle, la ralentissant considérablement.
• La bulle de savon « soluble » (Le vrai cas) : C'est là que la nouvelle méthode brille. Le savon pouvait monter et descendre de la bulle.
  • Si le savon descendait facilement (désorption élevée), la bulle se comportait presque comme une bulle propre.
  • Si le savon montait facilement (adsorption élevée), la bulle se comportait comme la version au savon coincé.
  • Au milieu, la bulle montrait une danse complexe : elle ralentissait, changeait de trajectoire, et laissait une « traînée » de savon dans l'eau derrière elle en montant.

Pourquoi cela compte (Selon le document)

Les auteurs affirment que cette méthode est robuste, évolutive et précise.

• Robuste : Elle ne plante pas lorsque la bulle prend des formes étranges ou se brise.
• Évolutive : Elle peut s'exécuter sur des supercalculateurs pour gérer efficacement d'énormes simulations 3D complexes.
• Précise : Elle prédit correctement la vitesse de montée des bulles et leurs oscillations, en accord avec la physique réelle.

En bref : Ils ont construit un nouveau moteur numérique capable enfin de simuler le comportement des bulles de savon dans l'espace 3D, gérant la danse délicate entre la forme de la bulle et le mouvement des molécules de savon entrant et sortant de sa peau, le tout sans perdre en précision ni faire planter l'ordinateur.

Résumé technique

Résumé technique : Une méthode hybride Volume de Fluide et Champ de Phase pour les écoulements interfaciaux chargés de tensioactifs solubles

Énoncé du problème
La simulation numérique précise des écoulements interfaciaux impliquant des tensioactifs solubles présente un défi redoutable en raison de la nature multi-échelle et multiphysique du problème. Contrairement aux tensioactifs insolubles, les espèces solubles sont régies par un système couplé impliquant la convection-diffusion en volume, le transport interfacial et les processus d'échange cinétique (adsorption/désorption). Ces interactions génèrent des champs de tension de surface variant spatialement et temporellement, conduisant à des contraintes de Marangoni qui modifient considérablement la topologie, la stabilité et le mouvement de l'écoulement. Les méthodes numériques existantes présentent des limitations distinctes :

• Méthodes intégrales/éléments de frontière : Limitées aux régimes d'écoulement de Stokes et peinent à gérer les changements de topologie (rupture/coalescence).
• Suivi de front (Front-Tracking) : Haute précision mais mise en œuvre complexe en 3D, nécessitant un traitement topologique explicite et confrontée à des problèmes d'équilibrage de charge en calcul parallèle.
• Méthodes Level-Set : Gèrent naturellement les changements de topologie mais souffrent souvent de problèmes de conservation de masse nécessitant une réinitialisation et des corrections a posteriori.
• Méthodes Champ de Phase : Gèrent naturellement les changements de topologie mais reposent souvent sur des hypothèses d'équilibre thermodynamique qui peuvent être violées par une forte convection, entraînant une non-conservation de la masse.
• Méthodes Volume de Fluide (VoF) : Excellentes pour la conservation de masse et les interfaces nettes, mais traditionnellement en difficulté pour définir les opérateurs différentiels de surface requis pour le transport de tensioactifs sur une interface d'épaisseur nulle.

Il existe un besoin spécifique d'une méthode géométrique basée sur le VoF qui prenne en charge le raffinement de maillage adaptatif (AMR), gère les tensioactifs solubles avec un échange dynamique volume-interface, et fonctionne dans un cadre open-source évolutif.

Méthodologie
Les auteurs présentent un cadre numérique hybride implémenté dans le code open-source Basilisk, qui combine une méthode géométrique Volume de Fluide (VoF) avec une approche Champ de Phase diffuse.

1. Représentation de l'interface :

   • VoF ($c$) : Un champ scalaire représentant la fraction volumique d'un fluide, utilisé pour le suivi conservatif de l'interface et la résolution de la quantité de mouvement. Il assure une conservation exacte de la masse et prend en charge la reconstruction géométrique.
   • Champ de Phase ($\phi$) : Une variable auxiliaire évoluant via une formulation d'interface diffuse conservatrice précise basée sur Allen-Cahn (ACDI). Elle sert de support lisse pour le transport de tensioactifs, fournissant des vecteurs normaux, une courbure et des gradients de surface régularisés sans nécessiter de calculs explicites de dérivées de surface sur une interface nette.
   • Couplage : Le Champ de Phase est périodiquement réinitialisé à partir du champ VoF en utilisant une fonction de distance signée pour maintenir la cohérence.

2. Équations gouvernantes :

   • Mouvement du fluide : Équations de Navier-Stokes incompressibles avec densité et viscosité variant spatialement. Les forces de tension de surface sont incorporées comme forces volumiques : une force capillaire isotrope ($\sigma \kappa \mathbf{n} \delta_s$) et une force tangentielle de Marangoni ($\nabla_s \sigma \delta_m$).
   • Transport de tensioactifs : Deux champs scalaires sont résolus : la concentration en volume ($F$) et la concentration interfaciale ($f$, une représentation volumique de $\Gamma$).
     • L'équation de transport interfacial est régularisée pour éviter les dérivées de surface, convertissant la concentration de surface $\Gamma$ en un champ volumique $f = \Gamma \delta_\phi$ localisé au sein de l'interface diffuse.
     • Cinétique : L'adsorption et la désorption sont modélisées par des termes sources/puits régularisés ($J$) localisés à l'interface, reliant les concentrations en volume et interfaciales. Le modèle prend en charge des équations d'état générales (par exemple, isotherme de Langmuir) pour relier la tension de surface à la concentration.

3. Discrétisation numérique :

   • Spatiale : Discrétisation par volumes finis sur des grilles quadtree/octree dynamiquement raffinées (AMR). La résolution est concentrée près de l'interface pour résoudre les gradients raides de concentration et d'écoulement.
   • Temporelle : Une séparation d'opérateurs est utilisée. L'advection est traitée explicitement (schéma géométrique non fractionné Bell–Collela–Glaz), tandis que la diffusion, l'anti-diffusion et les termes sources sont traités implicitement en utilisant un solveur de Poisson multigrille modifié.
   • Stabilité : Une condition Courant–Friedrichs–Lewy (CFL) est imposée sur la base d'une vitesse effective incluant à la fois les contributions convectives et anti-diffusives pour assurer la stabilité en présence de forts gradients d'interface.

Contributions clés

• Cadre hybride : Développement d'une méthode qui conserve les avantages de conservation de masse et d'interface nette du VoF tout en utilisant une variable Champ de Phase uniquement comme support lisse pour le transport de tensioactifs, évitant ainsi les complexités géométriques des dérivées de surface sur des interfaces linéaires par morceaux.
• Capacité pour tensioactifs solubles : Extension des travaux précédents VoF-Champ de Phase (limités aux tensioactifs insolubles) au cas entièrement soluble, intégrant des termes couplés d'advection-diffusion en volume et d'échange cinétique.
• Évolutivité et adaptativité : Le cadre est entièrement adaptatif (AMR) et prend en charge les configurations 2D, axisymétriques et 3D avec une haute évolutivité parallèle, permettant la résolution de gradients multi-échelles dans les champs en volume et interfaciaux.
• Implémentation open-source : La méthode est implémentée dans Basilisk, fournissant un outil reproductible et accessible pour la communauté.

Résultats et validation
L'article valide la méthode par une série de tests de référence et de cas d'application :

1. Validation sans écoulement :

   • Sphère en expansion (3D) : Vérification de la récupération des solutions analytiques pour la dilution de tensioactifs insolubles lors d'une expansion isotrope.
   • Interface en rotation (2D) : Confirmation de la convergence d'ordre deux pour l'advection-diffusion de tensioactifs insolubles sur une interface en mouvement.
   • Adsorption/Désorption sur interface plane : Validation du solveur couplé volume-interface par rapport aux solutions analytiques pour l'adsorption à flux constant, la désorption pure et la cinétique complète de Langmuir. La méthode capture correctement la transition vers l'équilibre et la formation de couches d'appauvrissement/enrichissement en volume.

2. Bulle montante axisymétrique :

   • Limite insoluble : Démonstration de la réduction de la vitesse de montée terminale due aux contraintes de Marangoni, avec des résultats convergeant vers la limite d'interface propre à mesure que la concentration en tensioactif diminue.
   • Dynamique soluble : Variation systématique du nombre de Biot (désorption) et du nombre de Damköhler (adsorption).
     • Biot élevé (désorption rapide) : Le système s'approche de la limite d'interface propre car les tensioactifs sont rapidement retirés de l'interface.
     • Damköhler élevé (adsorption rapide) : Le système s'approche de la limite d'incompressibilité de surface (insoluble) car l'interface se sature rapidement.
   • Structure de l'écoulement : Les visualisations montrent la formation de "caps de stagnation" de tensioactifs à l'arrière de la bulle, qui suppriment la recirculation interne et modifient la structure du sillage.

3. Bulle montante tridimensionnelle :

   • Les simulations d'une bulle montant en 3D avec des tensioactifs solubles ont capturé des phénomènes complexes incluant la déformation de forme, les déviations de trajectoire latérale (mouvement hélicoïdal) et la formation de panaches traînants dus à la désorption.
   • La méthode a récupéré avec succès les limites asymptotiques propres et insolubles en 3D, démontrant sa robustesse face aux régimes d'écoulement instationnaires et aux géométries déformables.

Signification et affirmations
L'article affirme fournir un cadre numérique robuste et polyvalent qui comble le fossé entre les modèles classiques de tensioactifs insolubles et les formulations solubles réalistes. Sa signification réside dans :

• Traitement unifié : Il gère le couplage non linéaire, les interfaces évolutives et les configurations entièrement tridimensionnelles au sein d'une seule plateforme open-source.
• Fidélité physique : Il capture l'interaction complexe entre l'hydrodynamique, le transport volume/interfacial et les contraintes de Marangoni sans recourir à des modèles de fermeture de turbulence ou à des approximations empiriques.
• Applicabilité pratique : La méthode convient aux écoulements à surface libre impliquant une reconnexion d'interface (par exemple, rupture de jet, vagues déferlantes) en présence de tensioactifs solubles pertinents sur les plans industriel et biologique.

Les auteurs notent modestement que le cadre actuel est limité à des nombres de Péclet modérés en raison de la diffusion numérique inhérente à la discrétisation. Ils identifient la résolution précise des gradients interfaciaux raides dans le régime à haut nombre de Péclet (dominé par la convection) comme un défi pour les travaux futurs, visant à développer des formulations minimisant la diffusion numérique.