🔬 physics

A generalized, quasi-universal subgrid model for macroscopic moving contact line flows

Cet article présente un modèle de sous-maille généralisé et quasi universel qui couple la relation asymptotique appariée de Cox à des modèles d'angle de contact empiriques afin de permettre des simulations prédictives et indépendantes du maillage des écoulements de lignes de contact mobiles macroscopiques sans nécessiter de calibration spécifique au cas ou de paramètres phénoménologiques.

Auteurs originaux : Vyankatesh Manoj Mundhada, Manoj Kumar Tripathi

Publié 2026-09-29✓ Author reviewed ⓘ

📖 7 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Vyankatesh Manoj Mundhada, Manoj Kumar Tripathi

Article original sous licence CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Imaginez une goutte d'eau tombant sur une surface plane. À mesure qu'elle s'étale, le bord où le liquide, l'air et le solide se rejoignent ne glisse pas simplement comme le pneu d'une voiture sur le pavé. Au lieu de cela, cette limite, connue sous le nom de ligne de contact, se comporte d'une manière qui défie les règles standards de la mécanique des fluides utilisées pour décrire la plupart des liquides en mouvement. Pendant des décennies, les scientifiques ont lutté pour simuler ce mouvement sur ordinateur car la physique à l'extrême bord de la goutte implique des échelles si infimes — descendant jusqu'à la taille de molécules individuelles — que tenter de calculer chaque détail nécessiterait plus de puissance de calcul que ce qui existe actuellement. Pour contourner cela, les chercheurs ont utilisé des raccourcis, ou « modèles de sous-maille », pour deviner ce qui se passe à ce bord microscopique. Cependant, ces raccourcis présentent un défaut majeur : les résultats qu'ils produisent varient souvent en fonction de la finesse avec laquelle l'ordinateur divise l'espace en une grille, rendant les simulations peu fiables pour des prédictions du monde réel.

Une équipe de chercheurs de l'Indian Institute of Science Education and Research à Bhopal a développé une nouvelle façon de gérer ce problème. Ils ont créé un modèle qui permet aux ordinateurs de simuler les lignes de contact mobiles avec précision sans avoir besoin de connaître les détails impossibles à mesurer du monde microscopique. Leur approche comble le fossé entre la physique minuscule et invisible à l'extrémité d'une goutte et le comportement plus large et visible du fluide. Ce faisant, ils ont produit des simulations qui restent cohérentes et précises, quel que soit le format de la grille de l'ordinateur, une étape cruciale vers la prédiction du comportement des liquides dans tout, de l'impression à jet d'encre à la récupération du pétrole, sans avoir besoin d'ajuster le modèle pour chaque nouvelle situation.

Le cœur du problème réside dans un conflit entre la façon dont les fluides se déplacent et la façon dont ils touchent les surfaces solides. Dans la mécanique des fluides standard, une règle appelée « condition de non-glissement » stipule qu'un fluide touchant une surface solide doit cesser tout mouvement à cette frontière. Bien que cela fonctionne bien pour la plupart des écoulements, cela s'effondre au niveau de la ligne de contact d'une goutte en mouvement, créant une singularité mathématique où les forces deviennent infinies. Pour corriger cela, les scientifiques ont historiquement supposé que le fluide glissait très légèrement à la surface, mais la distance de ce glissement est si petite qu'il est impossible de la mesurer directement. Par conséquent, les modèles informatiques précédents devaient deviner ces valeurs microscopiques ou compter sur des paramètres d'ajustement pour correspondre aux données expérimentales. Cela signifiait que si un chercheur voulait simuler un type de liquide différent ou une vitesse différente, il devait souvent repartir de zéro et recalibrer le modèle, le privant ainsi de sa capacité à prédire de nouveaux scénarios.

Les chercheurs de cette étude ont réalisé que l'incohérence des simulations précédentes provenait d'une mauvaise compréhension de la façon dont l'angle de la goutte change à mesure que l'on s'éloigne de la ligne de contact. Ils se sont appuyés sur une théorie bien établie des années 1980 qui décrit comment l'interface du liquide se courbe en raison de la compétition entre la tension superficielle et la friction interne du fluide. Cette théorie montre que l'angle que vous observez dépend de la distance à laquelle vous vous trouvez du bord. Les chercheurs ont combiné cette compréhension théorique avec des formules empiriques existantes — des règles dérivées de l'observation de gouttes réelles — pour créer un nouveau modèle unifié. Au lieu de forcer l'ordinateur à deviner la longueur de glissement microscopique ou l'angle exact au mur, leur modèle calcule l'angle correct à utiliser en fonction de la taille des cellules de la grille de l'ordinateur.

Pour tester leur idée, l'équipe a effectué une série de simulations informatiques d'une goutte s'étalant sur une plaque plane. Ils ont comparé leur nouvelle méthode aux anciennes techniques en utilisant différentes tailles de grille, allant de grossière à très fine. Avec les anciennes méthodes, le changement de taille de grille modifiait considérablement les résultats, faisant en sorte que la goutte simulée s'étalait à des rythmes différents ou prenait des formes différentes. Avec le nouveau modèle, cependant, les résultats restaient pratiquement identiques pour toutes les tailles de grille. La goutte simulée se comportait de manière cohérente, que l'ordinateur utilise une grille grossière ou une grille fine. Cela a démontré que le modèle avait réussi à éliminer la dépendance artificielle à la grille, un phénomène appelé dépendance à la grille.

L'étude a également exploré comment ce nouveau modèle se comporte sous diverses conditions, telles que différentes vitesses d'écoulement et différents types de liquides. Ils ont testé des scénarios où le liquide s'étale (avance) et où il se retire (recule). Dans chaque cas, le nouveau modèle a produit des résultats stables et indépendants de la grille. Une découverte clé est que le modèle ne nécessite aucun « paramètre d'ajustement » dérivé d'expériences spécifiques pour fonctionner. Une fois que les chercheurs ont établi l'échelle à laquelle les formules empiriques originales ont été dérivées, le modèle pouvait prédire le comportement de l'écoulement sans ajustement supplémentaire. Il s'agit d'un départ significatif par rapport aux méthodes précédentes, qui nécessitaient souvent un calibrage unique pour chaque nouvel ensemble de conditions.

Les chercheurs ont également examiné pourquoi certaines anciennes simulations semblaient bien fonctionner par le passé, même sans cette nouvelle correction. Ils ont découvert que dans certaines situations, particulièrement avec les systèmes gaz-liquide, l'erreur introduite par la taille de la grille était naturellement assez petite pour passer inaperçue. De plus, si la grille de l'ordinateur correspondait par hasard à l'échelle spécifique à laquelle les données expérimentales originales ont été collectées, les résultats s'alignaient de manière fortuite avec la réalité. Cependant, il s'agissait d'une question de chance plutôt que d'une solution robuste. Le nouveau modèle élimine cet élément de hasard, fournissant un cadre fiable qui fonctionne sur une large gamme de conditions physiques, incluant différentes densités et viscosités.

L'un des aspects les plus prometteurs de ce travail est son potentiel de puissance prédictive. Parce que le modèle ne repose pas sur la connaissance de la longueur de glissement microscopique ou de l'angle de contact exact au mur — des quantités notoirement difficiles à mesurer — il peut être appliqué à de nouveaux problèmes où ces valeurs sont inconnues. Les chercheurs ont montré qu'en connaissant simplement l'échelle à laquelle les observations expérimentales originales ont été faites, le modèle pouvait être appliqué universellement. Cela signifie que les ingénieurs et les scientifiques pourraient utiliser ces simulations pour concevoir de meilleurs systèmes de revêtement de surfaces, de gestion du flux de fluides dans les roches poreuses ou pour comprendre les processus biologiques dans les poumons, avec la certitude que les résultats ne sont pas des artefacts de la grille de l'ordinateur.

L'étude conclut que si l'approximation de premier ordre de leur modèle fonctionne bien pour la plupart des angles, des scénarios plus complexes impliquant des angles de contact très importants pourraient nécessiter une correction mathématique légèrement plus détaillée. Cependant, l'approche fondamentale consistant à lier la théorie microscopique aux observations macroscopiques a prouvé son efficacité. En traitant la ligne de contact non pas comme un point fixe aux propriétés mystérieuses, mais comme une caractéristique dynamique qui répond à l'échelle d'observation, les chercheurs ont fourni un outil qui est à la fois précis et largement applicable. Ce travail représente un passage du simple ajustement de données à la compréhension de la physique sous-jacente, offrant une voie vers des simulations capables de réellement prédire comment les liquides se comporteront dans le monde réel.

Résumé technique : Un modèle de sous-maille généralisé et quasi universel pour les écoulements macroscopiques à ligne de contact mobile

Énoncé du problème
La simulation numérique des écoulements à ligne de contact mobile (MCL - Moving Contact Line) est fondamentalement mise au défi par la singularité de contrainte découlant de la condition de non-glissement sur les surfaces solides et par la nature intrinsèquement multi-échelle du phénomène. La résolution des échelles de longueur microscopiques (longueurs de glissement de l'ordre du nanomètre) nécessaires pour éliminer cette singularité rend les simulations macroscopiques numériquement intraitable. Par conséquent, des modèles de sous-maille sont nécessaires pour combler le fossé entre la physique microscopique et la résolution macroscopique. Cependant, les approches existantes souffrent souvent de deux limitations critiques :

Dépendance à la grille : Les résultats numériques, en particulier la position de la ligne de contact et la forme de l'interface, varient considérablement avec la résolution du maillage lorsque des modèles d'angle de contact empiriques sont appliqués directement sans tenir compte de l'échelle à laquelle ils ont été dérivés.
Manque de prédictivité : De nombreux modèles invariants d'échelle nécessitent la prescription de paramètres phénoménologiques, tels que la longueur de glissement microscopique ( $\lambda$ ou $s_m$ ) et l'angle de contact mural microscopique ( $\theta_w$ ), qui sont difficiles à déterminer expérimentalement ou a priori. Cette dépendance à un étalonnage spécifique au cas limite la capacité de prédiction des simulations.

Méthodologie
Les auteurs proposent un modèle de sous-maille quasi universel qui couple la théorie de l'appariement asymptotique de Cox avec des modèles empiriques d'angle de contact dynamique. La méthodologie centrale comprend les étapes suivantes :

Cadre théorique : L'étude s'appuie sur la théorie hydrodynamique de Cox, qui relie l'angle de contact apparent $\theta(s)$ à une distance $s$ de la ligne de contact à l'angle de paroi microscopique $\theta_w$ et au nombre capillaire de la ligne de contact ( $Ca_{CL}$ ). La théorie introduit une dépendance logarithmique vis-à-vis du rapport d'échelle de longueur ( $s/s_m$ ).
Stratégie de couplage : Au lieu de prescrire des paramètres microscopiques inconnus ( $\theta_w$ et $s_m$ ), les auteurs dérivent une relation qui élimine ces variables. En reconnaissant que les modèles empiriques (ex. Kistler, Jiang, Tanner) fournissent l'angle de contact $\theta_0$ à une échelle d'observation spécifique $s_0$ (typiquement la résolution des données expérimentales utilisées pour formuler le modèle), ils établissent un pont entre le modèle empirique et l'échelle de la grille.
Dérivation du modèle modifié : Les auteurs introduisent un facteur de correction dépendant de l'échelle $\chi(Ca_{CL}, s)$ appliqué au nombre capillaire dans la relation empirique. En utilisant un développement de Taylor au premier ordre de la fonction $G(\theta)$ de Cox, ils dérivent un nombre capillaire modifié :
$\chi(s) = Ca_{CL} \left[ 1 + \zeta \ln\left(\frac{s}{s_0}\right) \right]$
où $\zeta$ est un paramètre dérivé analytiquement du modèle empirique spécifique et du rapport de viscosité. Cela permet d'imposer un angle de contact dépendant de la taille de la grille $\theta(s)$ qui satisfait la théorie de Cox sans nécessiter la connaissance de $s_m$ ou $\theta_w$ .
Implémentation numérique : Le modèle est implémenté dans un cadre de Volume de Fluide (VoF) utilisant le code Basilisk. La vitesse de la ligne de contact est calculée simultanément avec l'angle de contact et la position à chaque pas de temps pour assurer la cohérence, évitant ainsi les erreurs associées aux grandeurs temporellement décalées. L'approche est testée sur des scénarios d'étalement de goutte et de retrait de plaque axisymétriques à travers une large gamme de nombres de Reynolds ( $Re$ ) et de Weber ( $We$ ).

Résultats clés

Indépendance de la grille : Les modèles d'angle de contact « modifiés » proposés éliminent avec succès la dépendance à la grille. Les simulations menées sur des grilles allant de $\Delta_5$ à $\Delta_9$ (où $\Delta_i = D_0/2^i$ ) montrent que la position de la ligne de contact et la forme de l'interface de la goutte convergent vers une solution unique. En revanche, l'utilisation des modèles empiriques originaux (non modifiés) entraîne des erreurs significatives (jusqu'à ~1,5 % pour la forme de l'interface et plus pour la position de la ligne de contact) qui ne convergent pas avec l'affinement du maillage.
Réduction de l'erreur : Pour les systèmes liquide-liquide (où la dépendance à la grille est la plus sévère), les erreurs relatives moyennes de la position de la ligne de contact ( $E_r$ ) et de la forme de l'interface ( $E_s$ ) ont été réduites d'environ un ordre de grandeur par rapport aux simulations utilisant les modèles originaux. Par exemple, à $Re=100$ , $We=1$ , et $\theta_e=60^\circ$ , l'erreur de position de la ligne de contact est passée de ~0,69 % (Kistler original) à ~0,05 % (Kistler modifié).
Robustesse à travers les régimes : Le modèle a démontré son efficacité à travers diverses valeurs de $Re$ (1–200), $We$ (0,5–100) et angles de contact d'équilibre ( $30^\circ$ – $135^\circ$ ).
Corrections d'ordre supérieur : L'étude a identifié que pour les angles de contact obtus élevés ( $\ge 120^\circ$ ), l'approximation linéaire du premier ordre du facteur de correction $\zeta$ introduit des erreurs plus importantes. Une correction de second ordre a été dérivée et a montré qu'elle réduisait davantage les erreurs dans ces régimes spécifiques, bien que l'approximation de premier ordre reste suffisante pour les angles aigus et modérés.
Systèmes gaz-liquide : Les auteurs notent que pour les systèmes gaz-liquide, la méthode numérique originale présente souvent une dépendance à la grille négligeable sans modification, probablement parce que la variation de l'angle de contact avec la distance est faible et que la taille de la grille correspond incidemment à l'échelle d'observation expérimentale.

Signification et revendications
L'article affirme présenter le premier cadre capable de réaliser des simulations prédictives de flux à ligne de contact mobile sans nécessiter de calibration spécifique au cas ou la connaissance a priori des paramètres microscopiques ( $\theta_w$ et $s_m$ ).

Quasi-universalité : Le modèle est qualifié de « quasi universel » car il repose uniquement sur la validité de la théorie de Cox et du modèle d'angle de contact empirique choisi, indépendamment des paramètres de flux spécifiques comme $Re$ , $We$ , ou les rapports de densité/viscosité.
Dépendance à l'échelle d'observation : La méthode nécessite seulement la connaissance ponctuelle de l'échelle de longueur d'observation ( $s_0$ ) associée au modèle empirique. Une fois $s_0$ fixé pour un modèle donné, le modèle de sous-maille peut être appliqué à n'importe quel problème de flux macroscopique.
Élimination des paramètres d'ajustement : Contrairement aux approches antérieures qui nécessitaient des paramètres d'ajustement pour correspondre aux données expérimentales, cette méthode dérive les corrections nécessaires de manière analytique, permettant une véritable capacité de prédiction.
Généralisabilité : Bien que démontrée en 2D axisymétrique, la formulation est dérivée de descriptions locales de la ligne de contact et n'est pas spécifique à la géométrie axisymétrique, suggérant une extension potentielle aux lignes de contact en trois dimensions.

En résumé, ce travail traite du verrou fondamental de la dépendance à la grille dans les simulations de MCL en couplant rigoureusement la théorie hydrodynamique aux données empiriques, offrant une voie vers des prédictions numériques fiables et sans paramètres pour les écoulements multiphasiques complexes.

Articles similaires