Color-gradient lattice Boltzmann modeling of wetting boundary condition on curved solid boundaries

Cet article introduit une condition de bord de mouillage pour les surfaces solides courbes dans la méthode de Lattice Boltzmann à gradient de couleur en mettant à jour les paramètres d'ordre sur des nœuds fantômes, un schéma validé sur matériel GPU pour gérer efficacement les contrastes de densité et de viscosité importants tout en minimisant les courants parasites et en reproduisant avec précision les comportements de ligne de contact statiques et dynamiques.

L'essentiel

Imaginez que vous essayez de simuler le comportement d'une goutte d'eau lorsqu'elle frappe une surface courbe, comme une goutte de pluie tombant sur une feuille ou une bulle glissant le long d'un verre incurvé. Pour faire cela sur un ordinateur, les scientifiques utilisent une méthode appelée la Méthode de Lattice Boltzmann. Imaginez cette méthode comme une grille géante et invisible de petits carreaux recouvrant l'écran de l'ordinateur. Chaque carreau contient un peu d'information sur le « fluide », et l'ordinateur met à jour ces carreaux étape par étape pour voir comment le fluide se déplace.

La partie délicate est la condition aux limites — spécifiquement, comment le fluide se comporte lorsqu'il touche une paroi solide. Dans le monde réel, l'eau ne s'arrête pas net contre un mur ; elle forme un angle spécifique (appelé angle de contact) selon que la surface est mouillante (comme du verre propre) ou sèche (comme une voiture cirée).

Le Problème : Le « Fantôme » dans la Machine

Dans la simulation informatique, la paroi solide n'est pas une ligne lisse ; elle est dentelée car elle est composée de carreaux de grille carrés. Pour que les mathématiques fonctionnent, l'ordinateur doit savoir ce que fait le fluide à l'intérieur de la paroi solide, même s'il n'y a pas de fluide là. Ces emplacements imaginaires à l'intérieur de la paroi sont appelés « nœuds fantômes » (ghost nodes).

Les méthodes précédentes pour dire à ces nœuds fantômes quoi faire présentaient des défauts :

• Elles créaient parfois de faux « courants fantômes » (vitesses parasites) où le fluide semblait bouger de lui-même sans aucune force.
• Elles avaient du mal avec les surfaces courbes, agissant souvent comme si elles n'avaient été conçues que pour des parois plates.
• Elles nécessitaient parfois des mathématiques spéciales et compliquées pour gérer un angle neutre (où l'eau ne s'étale ni ne perle).

La Solution : Une Nouvelle Règle pour les Fantômes

Les auteurs de cet article ont introduit une nouvelle règle plus simple pour ces nœuds fantômes.

L'Analogie : Imaginez que le fluide a une « humeur » (représentée par une couleur, de 0 pour le gaz à 1 pour le liquide). Dans le monde réel, cette humeur change de manière fluide entre le gaz et le liquide au moment où l'on traverse la surface.

• Ancienne Méthode : C'était comme essayer de deviner l'humeur d'une personne se tenant derrière un mur en criant une supposition au hasard.
• Nouvelle Méthode : Les auteurs ont réalisé que si vous connaissez l'« humeur » de la personne se tenant juste à l'extérieur du mur (dans le fluide), vous pouvez mathématiquement prolonger cette courbe d'humeur fluide à travers le mur jusqu'au nœud fantôme. Ils demandent simplement : « Si le fluide veut former un angle de 45 degrés ici, quelle doit être l'humeur du nœud fantôme pour que cela se produise ? »

Cette nouvelle règle est comme un pont sans couture. Elle prolonge la forme naturelle de la goutte de fluide jusqu'au bord et légèrement à l'intérieur de la paroi solide, garantissant que l'angle que la goutte fait avec la paroi est exactement celui demandé par le scientifique.

Ce Qu'ils Ont Testé

Pour prouver que leur nouvelle règle fonctionne, ils ont lancé plusieurs simulations sur une puce informatique très puissante (un GPU NVIDIA A100) :

1. La Goutte Statique : Ils ont placé une goutte d'eau sur une plaque plate et sur un cylindre courbe. Ils ont vérifié si la goutte se stabilisait à l'angle exact qu'ils avaient demandé.
   • Résultat : Leur nouvelle règle était plus précise que la meilleure méthode précédente, surtout lorsque l'angle était très aigu (comme une goutte qui perle) ou très plat (comme une goutte qui s'étale).
2. La Particule Flottante : Ils ont simulé un cylindre flottant à la limite entre l'huile et l'eau.
   • Résultat : Leur méthode calculait la position de la ligne de flottaison plus précisément qu'auparavant.
3. La Goutte en Chute Libre : Ils ont simulé une goutte tombant et frappant un cylindre, observant comment elle éclabousse et s'étale.
   • Résultat : La goutte s'est comportée de manière réaliste, et la nouvelle règle n'a provoqué aucun mouvement étrange ou artificiel dans le fluide.

Points Clés à Retenir

• Précision : La nouvelle méthode gère beaucoup mieux les surfaces courbes que les anciennes méthodes, maintenant l'angle correct du fluide, que la paroi soit plate ou arrondie.
• Stabilité : Elle crée très peu de « bruit artificiel » (courants parasites) dans la simulation, ce qui signifie que le fluide semble plus naturel.
• Simplicité : Elle évite d'avoir recours à des mathématiques spéciales et compliquées lorsque l'angle de contact est exactement de 90 degrés (neutre), ce qui était un casse-tête pour les méthodes précédentes.
• Vitesse : En utilisant des puces informatiques modernes (GPU) et un style de programmation spécifique, ils ont rendu les simulations très rapides. Ils ont constaté que l'utilisation d'un format de nombre légèrement moins précis (simple précision) permettait à l'ordinateur de fonctionner deux fois plus vite sans altérer les résultats pour la plupart des tests.

En résumé, les auteurs ont construit un meilleur « livre de règles » pour la façon dont les simulations informatiques gèrent la limite où le liquide rencontre une paroi solide, rendant les gouttes numériques plus semblables et plus réactives aux vraies gouttes, même sur des surfaces courbes.

Résumé technique

Résumé technique : Modélisation de la méthode de Lattice Boltzmann à gradient de couleur pour les conditions limites de mouillage sur des frontières solides courbes

Énoncé du problème
La modélisation précise des phénomènes de mouillage est cruciale pour simuler des processus technologiquement pertinents tels que l'écoulement en milieux poreux, la microfluidique et le séchage de films colloïdaux. Dans les cadres de la méthode de Lattice Boltzmann (LBM) utilisant des interfaces diffuses, la mise en œuvre de conditions limites de mouillage sur des frontières solides courbes présente des défis importants. Les approches existantes souffrent souvent de limitations incluant l'introduction de sources ou de puits de masse artificiels, la génération de vitesses parasites (spureuses) près des interfaces, des restrictions aux géométries planes, et des incohérences dans l'application simultanée des conditions de non-glissement et de mouillage. De plus, de nombreuses méthodes nécessitent un traitement spécial pour les angles de contact neutres ($\theta = 90^\circ$) ou reposent sur des extrapolations linéaires qui introduisent des erreurs à proximité des frontières solides.

Méthodologie
Les auteurs introduisent une condition limite de mouillage spécifiquement conçue pour les frontières solides courbes au sein d'un cadre LBM à gradient de couleur en deux dimensions. Le cœur de la méthodologie repose sur la mise à jour du paramètre d'ordre (champ de couleur, $\phi$) sur les « nœuds fantômes » situés à l'intérieur de la phase solide.

• Règle de mise à jour des nœuds fantômes : Au lieu de prescrire une valeur constante ou d'utiliser une extrapolation linéaire, le schéma étend le profil d'équilibre de la couleur dans la phase solide. La valeur du champ de couleur à un nœud fantôme ($x_G$) est calculée à partir du champ de couleur à un point correspondant dans le domaine fluide ($x_F$) situé à une distance $\Delta x$ le long de la normale à la paroi solide.
• Formulation mathématique : En supposant que la distribution de couleur est proche de l'équilibre, la dérivée du paramètre d'ordre dans la direction de la normale à la paroi est dérivée du profil d'équilibre et de l'angle de contact $\theta$ prescrit. L'intégration de cette relation produit une règle de mise à jour non linéaire :
  $$\phi_G = \frac{\phi_F}{\phi_F + (1 - \phi_F) \exp(\varepsilon)}$$
  où $\varepsilon = -4\Delta x \cos \theta / W$, et $W$ est la largeur de l'interface.
• Caractéristiques clés : Cette formulation garantit que le champ de couleur du nœud fantôme reste borné entre 0 et 1, empênant ainsi tout transfert de masse artificiel. Crucialement, la règle ne nécessite pas de traitement spécial lorsque $\theta = 90^\circ$ (où $\varepsilon = 0$), contrairement aux schémas précédents. La méthode est implémentée à l'aide d'un framework basé sur JAX s'exécutant sur des GPU NVIDIA A100, utilisant la fusion de noyaux (kernel fusion) et la différenciation automatique pour des performances élevées.

Contributions clés et résultats
Le schéma proposé a été validé par rapport à des solutions analytiques et au schéma de référence de Fakhari et al. [20] à travers plusieurs cas de test :

1. Goutte statique sur une plaque plane : Le schéma a systématiquement produit des erreurs plus faibles dans les angles de contact mesurés par rapport à Fakhari et al. sur une plage de $30^\circ$ à $150^\circ$. Il a également généré des courants parasites plus faibles ($5,74 \times 10^{-7}$ à $7 \times 10^{-7}$ unités de réseau) par rapport au schéma de référence.
2. Goutte sur un cylindre solide : Pour une goutte posée sur un cylindre courbe, la méthode proposée a démontré une meilleure précision dans la prédiction de la hauteur de l'apex de la goutte ($h_{max}$), particulièrement pour les angles de contact extrêmes ($30^\circ$ et $150^\circ$), où les erreurs ont été réduites de $\sim 0,76\%$ à $\sim 0,35\%$ par rapport au schéma de référence.
3. Particule fixe à l'interface liquide-gaz : Dans les simulations d'un réseau périodique de cylindres piégés à une interface, le schéma proposé a réduit les erreurs de déplacement vertical de l'interface à $\sim 9\%$ pour les angles extrêmes, contre $\sim 14\%$ pour le schéma de référence.
4. Impact de goutte dynamique : Le schéma a simulé avec succès l'impact de gouttes sur un cylindre à des nombres de Reynolds et de Bond élevés, capturant des comportements dynamiques complexes tels que l'étalement de la goutte sur des surfaces hydrophiles et la division sur des surfaces hydrophobes. Les résultats concordent bien avec les lois d'échelle analytiques pour l'épaisseur du film.
5. Performance de calcul : L'implémentation GPU a montré que les calculs en simple précision (FP32) étaient au moins deux fois plus rapides que ceux en double précision (FP64) tout en maintenant la précision pour les cas statiques et dynamiques modérés. Cependant, le FP32 est devenu instable pour le test d'impact de goutte à haut nombre de Reynolds/Weber, nécessitant le recours au FP64 pour ces scénarios à haute énergie.

Signification et revendications
L'article affirme que la condition limite de mouillage proposée offre une alternative robuste et précise pour simuler des écoulements diphasiques immiscibles sur des frontières solides courbes. En étendant le profil d'équilibre dans la phase solide via une règle spécifique de mise à jour des nœuds fantômes, la méthode :

• Préserve la capacité du modèle à gérer de grands contrastes de densité et de viscosité.
• Produit des courants parasites relativement faibles (ordres de grandeur inférieurs aux modèles de pseudo-potentiel populaires).
• Évite les singularités numériques associées aux angles de contact neutres présentes dans d'autres méthodes.
• Démontre une précision supérieure aux schémas existants (spécifiquement Fakhari et al.) pour les lignes de contact statiques et dynamiques sur des géométries courbes.

Les auteurs notent que, bien que l'implémentation actuelle soit restreinte à une exécution sur un seul GPU et en deux dimensions, la physique sous-jacente est directement applicable aux modèles de champ de phase et d'énergie libre et peut être étendue à trois dimensions. Ce travail constitue une base pour de futures simulations impliquant des particules solides à des interfaces liquide-gaz, notamment lorsqu'elles sont couplées à des avancées dans les calculs de forces capillaires de la ligne de contact et les schémas d'hystérésis.