Cahn-Hilliard Phase Field modelling captures nanoscale contact line dynamics on high-friction surfaces

Cette étude démontre qu'un modèle de champ de phase de Cahn-Hilliard, lorsqu'il est systématiquement calibré à l'aide de données de dynamique moléculaire pour capturer la friction de la ligne de contact et la dynamique de l'angle de contact, peut reproduire quantitativement les comportements de mouillage à l'échelle nanométrique sur des surfaces à forte friction, comblant ainsi le fossé entre les processus moléculaires et l'hydrodynamique continue.

L'essentiel

Imaginez que vous observez une goutte d'eau glisser sur une vitre. À l'œil nu, elle paraît lisse. Mais si vous pouviez rétrécir jusqu'à la taille d'une molécule, vous verriez une danse chaotique et saccadée où les molécules d'eau heurtent le verre et s'entrechoquent.

Depuis longtemps, les scientifiques tentent d'écrire des programmes informatiques pour prédire exactement comment ces gouttes se déplacent. Ils disposent de deux outils principaux :

1. La Dynamique Moléculaire (DM) : C'est comme une caméra ultra-microscopique à très haute vitesse. Elle suit chaque molécule individuellement. Elle est incroyablement précise mais nécessite un supercalculateur et prend une éternité à s'exécuter.
2. Les Modèles de Champ de Phase (CHNS) : C'est comme une vidéo lisse et continue. Elle traite le liquide comme une masse fluide plutôt que comme des particules individuelles. Elle est rapide et facile à exécuter, mais elle rate souvent les détails minuscules et désordonnés se produisant juste là où le liquide touche la surface solide (la « ligne de contact »).

Le Problème : Le Bord « Collant »
Lorsqu'une goutte se déplace, le bord où elle touche la surface est la partie la plus importante. Dans le monde réel (et dans la caméra microscopique), ce bord se « colle » ou subit une friction. Les modèles de vidéo lisse ont généralement du mal ici car ils supposent que le liquide glisse parfaitement ou glisse d'une manière qui ne correspond pas à la réalité. Ils obtiennent souvent la forme de la goutte incorrecte car ils ne peuvent pas tenir compte de cette « adhérence » microscopique.

La Solution : Une Approche Hybride
Les auteurs de cet article voulaient corriger le modèle de vidéo lisse pour qu'il se comporte exactement comme la caméra microscopique, sans avoir besoin de suivre chaque molécule individuellement. Ils y sont parvenus en créant un protocole d'étalonnage.

Pensez-y comme à l'accordage d'un instrument de musique. Le modèle lisse est l'instrument, et la simulation microscopique est la justesse parfaite.

1. Le Montage : Ils ont simulé de l'eau et de l'hexane (un type d'huile) glissant l'un sur l'autre entre deux parois mobiles, comme un sandwich qu'on presse et qu'on fait glisser.
2. L'Étalonnage : Ils ont d'abord exécuté la simulation microscopique lente et détaillée. Ils ont mesuré exactement à quel point le « bord » de l'eau résistait au mouvement (la friction de la ligne de contact) et comment la surface se courbait.
3. La Correction : Ils ont injecté ces « chiffres de friction » spécifiques dans le modèle de vidéo lisse. Ils n'ont pas simplement deviné ; ils ont ajusté le « cadran de friction » du modèle jusqu'à ce que le bord du modèle lisse se comporte exactement comme celui du modèle microscopique.

Les Résultats : Une Correspondance Parfaite
Une fois ce « cadran de friction » spécifique accordé, le modèle lisse est devenu incroyablement précis. Il pouvait désormais prédire :

• Comment la goutte se courbe : La courbe de la surface de l'eau près de la paroi.
• Jusqu'où la goutte se déplace : La position stable de la ligne de contact.
• Comment l'eau s'écoule : Les motifs tourbillonnaires à l'intérieur du liquide.

L'article affirme qu'en faisant correspondre simplement la friction de la ligne de contact (la résistance du bord au mouvement) aux données microscopiques, le modèle lisse peut reproduire la physique complexe et désordonnée du monde réel.

La Mise en Garde (Le Secret du « Glissement »)
Il y a un tout petit détail que le modèle lisse ne parvient toujours pas à saisir. Dans le monde microscopique, le tout bord de la ligne de contact « glisse » en réalité un tout petit peu plus que le reste du liquide. Le modèle lisse, même parfaitement accordé, n'inclut pas naturellement ce glissement supplémentaire. Les auteurs suggèrent que, bien que leur méthode constitue une énorme amélioration, les modèles futurs devront peut-être ajouter une règle spécifique pour tenir compte de ce bord « glissant » supplémentaire afin d'être 100 % parfaits.

En Résumé
Cet article traite d'apprendre à un modèle informatique simplifié et rapide à se comporter comme un modèle complexe et lent. Ils ont découvert que si vous dites simplement au modèle rapide exactement à quel point le bord de la goutte est « collant » (sur la base de données moléculaires réelles), il peut prédire avec précision comment la goutte se déplace, se courbe et s'écoule, comblant ainsi le fossé entre le monde microscopique des atomes et le monde macroscopique des fluides.

Résumé technique

1. Énoncé du problème

Une modélisation mésoscopique précise des processus de mouillage est cruciale pour des applications allant de la fabrication additive à la production d'énergie. Cependant, décrire le mouvement des lignes de contact mobiles (où les interfaces fluide-fluide rencontrent une surface solide) reste un défi majeur en dynamique des fluides computationnelle (CFD).

• Le problème de la singularité : Les modèles hydrodynamiques classiques (Navier-Stokes) avec des conditions aux limites d'adhérence (no-slip) prédisent une singularité de contrainte non intégrable à la ligne de contact.
• Le dilemme du glissement : Bien que l'introduction d'un glissement de Navier résolve la singularité, les preuves expérimentales et moléculaires suggèrent que les longueurs de glissement sur des surfaces planes et hydrophiles sont souvent nanométriques (ou effectivement nulles). L'utilisation de longueurs de glissement physiquement irréalistes pour la faisabilité numérique compromet la précision physique.
• Le fossé : Il manque un consensus sur la manière de concilier les modèles continus macroscopiques avec la physique à l'échelle moléculaire, en particulier concernant le frottement de la ligne de contact et les mécanismes de dissipation d'énergie sur les surfaces à fort frottement.

2. Méthodologie

Les auteurs emploient une approche multi-échelle, couplant des simulations de Dynamique Moléculaire (DM) avec un modèle de Champ de Phase (CP) basé sur les équations de Cahn-Hilliard Navier-Stokes (CHNS).

A. Simulations de Dynamique Moléculaire (DM)

• Système : Un système biphasique eau/hexane dans une configuration d'écoulement de Couette à deux phases, confiné entre des parois de type silice.
• Configuration : Le système utilise les champs de force SPC/E pour l'eau et OPLS-AA pour l'hexane. La paroi solide est modélisée comme une surface de type silice avec une mouillabilité ajustable (cas hydrophile et hydrophobe).
• Objectif : Générer des données de « vérité terrain » pour la courbure de l'interface, les angles de contact dynamiques et les profils d'écoulement sous diverses vitesses de paroi ($u_w \in [1.12, 4.84]$ m/s).
• Observables : Les simulations DM permettent la mesure directe des coefficients de frottement de la ligne de contact et l'extraction des propriétés de transport (viscosité, tension superficielle) sans ajustement empirique.

B. Modélisation par Champ de Phase (CHNS)

• Équations gouvernantes : Le modèle résout l'équation couplée de Cahn-Hilliard (pour l'évolution de phase) et les équations de Navier-Stokes (pour la quantité de mouvement du fluide).
• Caractéristiques clés :
  • Interface diffuse : L'interface possède une épaisseur finie ($\epsilon$), permettant un mécanisme de glissement diffusif local qui régularise la singularité de la ligne de contact sans nécessiter de glissement de Navier macroscopique.
  • Conditions aux limites : Inclut une condition de relaxation d'énergie de paroi pour régir la dynamique de l'angle de contact et une longueur de glissement de Navier sub-nanométrique ($\ell_N$) strictement pour la stabilité numérique.
• Protocole d'étalonnage : Les auteurs établissent un protocole systématique pour paramétrer le modèle CHNS directement à partir des données DM :
  1. Épaisseur de l'interface ($\epsilon$) : Fixée à 0,3 nm (correspondant à la largeur intrinsèque de la DM).
  2. Mobilité ($M$) : Contrainte par la condition de Limite d'Interface Aiguë (SIL) ($\epsilon \leq 4\sqrt{M\eta^*}$) pour minimiser les croisements de lignes de courant non physiques.
  3. Frottement de la ligne de contact ($\mu_f$) : Le seul paramètre étalonné empiriquement. Il est extrait en faisant correspondre la relation dérivée de la DM entre la vitesse de la ligne de contact et l'angle de contact dynamique.
  4. Viscosité : Mesurée indépendamment à partir de la DM.

3. Contributions clés

1. Protocole d'étalonnage systématique : L'article démontre qu'un modèle de Champ de Phase peut reproduire quantitativement les résultats de la DM avec un ajustement minimal des paramètres. Crucialement, il identifie le frottement de la ligne de contact comme le seul paramètre physique nécessitant un étalonnage empirique contre les données DM. Tous les autres paramètres (mobilité, épaisseur de l'interface) sont déterminés par des contraintes numériques (limite d'interface aiguë) ou des mesures DM indépendantes.
2. Reproduction quantitative de la dissipation : L'étude sépare et capture avec succès les deux canaux principaux de dissipation d'énergie dans le mouillage :
   • Dissipation visqueuse : Capturée par les équations de Navier-Stokes (courbure de l'interface).
   • Frottement de la ligne de contact : Capturé par le coefficient de frottement étalonné (inclinaison de l'angle de contact).
3. Résolution du croisement de lignes de courant : En adhérant strictement à la condition de limite d'interface aiguë pour le paramètre de mobilité $M$, les auteurs minimisent le croisement non physique des lignes de courant à travers l'interface liquide-liquide, un artefact courant dans les simulations de Champ de Phase.
4. Identification des limites du modèle : Le travail met en évidence une divergence entre le frottement CHNS étalonné et le frottement mesuré par DM, l'attribuant à un glissement localisé à la ligne de contact qui n'est pas entièrement capturé par la diffusion standard du Champ de Phase seule.

4. Résultats clés

• Courbure de l'interface : Le modèle CHNS reproduit quantitativement les profils de courbure de l'interface ($\theta_I(z)$) observés en DM pour les lignes de contact avançantes et récedantes sur une large gamme de vitesses.
• Déplacement de la ligne de contact :
  • Pour les surfaces hydrophobes ($\theta_0 \approx 97,3^\circ$), le modèle montre un excellent accord avec le déplacement à l'état stationnaire de la DM.
  • Pour les surfaces hydrophiles ($\theta_0 \approx 80,9^\circ$), le modèle sous-estime le déplacement à haute vitesse. Cela est attribué à l'omission du terme de couplage de flux diffusif dans l'équation de la quantité de mouvement, qui devient significatif à haute vitesse.
• Structure de l'écoulement : Les champs de vitesse et les motifs de lignes de courant (incluant les zones de recirculation dans la phase de faible viscosité) montrent un accord qualitatif et quantitatif avec la DM.
• Sensibilité des paramètres :
  • Mobilité ($M$) : Affecte la courbure de l'interface loin de la paroi et le déplacement stationnaire, mais n'affecte pas l'angle de contact dynamique.
  • Frottement ($\mu_f$) : Le seul déterminant du comportement de l'angle de contact dynamique.
  • Glissement de Navier ($\ell_N$) : Le glissement sub-nanométrique introduit pour la stabilité a un impact négligeable sur les observables physiques dans la plage testée.

5. Importance et conclusion

Ce travail fournit un cadre robuste pour combler le fossé entre les descriptions moléculaires et macroscopiques du mouillage.

• Validation des modèles mésoscopiques : Il prouve que les modèles de Champ de Phase, lorsqu'ils sont soigneusement étalonnés à l'aide de données DM pour le frottement de la ligne de contact, peuvent capturer avec précision l'hydrodynamique complexe à l'échelle nanométrique sur les surfaces à fort frottement.
• Insight physique : L'étude confirme que le frottement de la ligne de contact est l'entrée moléculaire dominante requise pour une modélisation mésoscopique précise, distincte de la viscosité de volume.
• Perspectives futures : Les auteurs notent que, bien que le modèle actuel soit robuste, il manque un « composant de glissement » spécifique localisé à la ligne de contact. Les travaux futurs devraient explorer les Conditions aux Limites de Navier Généralisées (GNBC) ou des modèles de glissement où la longueur de glissement est une fonction de la variable de phase pour concilier pleinement les modèles continus avec la physique moléculaire.

En résumé, l'article établit que le frottement de la ligne de contact est le paramètre empirique critique pour les modèles de mouillage et démontre un protocole rigoureux et basé sur la physique pour intégrer les données de Dynamique Moléculaire dans les simulations continues, faisant progresser considérablement la fiabilité de la modélisation multi-échelle du mouillage.