Modelling Capillary Rise with a Slip Boundary Condition: Well-posedness and Long-time Dynamics of Solutions to Washburn's Equation

Cet article étend l'équation de Washburn en y intégrant une condition de glissement, démontrant ainsi l'existence et l'unicité globales d'une solution positive bien posée dont la dynamique à long terme converge vers l'équilibre, que ce soit de manière monotone ou oscillatoire, pour tout paramètre de glissement positif et un rapport initial de hauteur spécifique.

L'essentiel

Le Grand Voyage de la Goutte d'Eau : Quand le Glissement Change la Don

Imaginez que vous trempez une paille fine dans un verre d'eau. Vous savez tous ce qui se passe : l'eau monte toute seule dans la paille, comme par magie, jusqu'à s'arrêter à une certaine hauteur. C'est ce qu'on appelle la capillarité.

Depuis plus d'un siècle, les scientifiques utilisent une équation célèbre (l'équation de Washburn) pour prédire exactement comment l'eau monte. Mais il y a un petit problème : cette vieille équation suppose que l'eau "colle" parfaitement aux parois de la paille (comme si elle était collée avec de la super-glue). En réalité, à l'échelle microscopique, l'eau peut parfois glisser un peu le long des parois, comme un patineur sur la glace.

Ce papier de recherche (par Rapajić, Simić et Süli) fait deux choses principales :

1. Il met à jour la vieille équation pour inclure ce "glissement".
2. Il prouve mathématiquement que cette nouvelle équation fonctionne parfaitement et que l'eau finira toujours par s'arrêter à la bonne hauteur, même si elle oscille un peu avant.

---

1. La Nouvelle Règle du Jeu : Le Patineur sur la Glace

Dans l'ancien modèle, on imaginait que l'eau, au contact de la paroi, s'arrêtait net (vitesse zéro). C'est ce qu'on appelle la condition "sans glissement". Mais les physiciens ont remarqué que cela créait une contradiction bizarre au tout début du mouvement : si l'eau colle à la paroi, comment peut-elle monter ? C'est le "paradoxe".

L'analogie du patineur :
Imaginez que l'eau est un patineur sur une patinoire.

• L'ancien modèle (Sans glissement) : Le patineur est attaché par les chevilles au bord de la glace. Pour avancer, il doit faire un effort énorme pour se décoller.
• Le nouveau modèle (Avec glissement) : Le patineur porte des patins. Il glisse légèrement sur le bord. Cela change la façon dont il accélère et freine, mais il arrive toujours au même endroit.

Les auteurs ont ajouté un paramètre, disons le "coefficient de glissement", à l'équation. Si ce coefficient est nul, on retrouve l'ancienne équation. S'il est positif, l'eau glisse un peu, ce qui modifie la vitesse de montée, surtout au début.

2. La Preuve Mathématique : Pas de Chaos, Juste de l'Ordre

Les mathématiciens adorent prouver que leurs équations sont "saines". Ils veulent être sûrs de trois choses (ce qu'on appelle la "bien-posée" d'un problème) :

1. Existence : Une solution existe-t-elle ? (Oui, l'eau va bien monter).
2. Unicité : Y a-t-il une seule façon pour l'eau de monter ? (Oui, pas de magie ni de bifurcations bizarres).
3. Stabilité : Si on change un tout petit peu la hauteur de départ, est-ce que le résultat change radicalement ? (Non, c'est prévisible).

L'analogie du pont :
Imaginez que vous construisez un pont (votre équation) pour traverser une rivière (le temps).

• Les auteurs ont prouvé que leur pont ne s'effondre jamais, même si vous commencez à marcher dessus avec un pied un peu plus en avant que l'autre (condition initiale).
• Ils ont utilisé des outils mathématiques puissants (comme le théorème de Picard-Lindelöf, qu'on peut voir comme un GPS très précis) pour montrer que, peu importe comment l'eau commence à bouger, elle suivra toujours une trajectoire unique et définie.

Ils ont aussi corrigé une erreur dans un ancien papier (de 2025, mais basé sur des travaux antérieurs) qui avait essayé de prouver cela mais avait laissé des trous dans sa logique. Eux, ils ont comblé tous les trous !

3. Le Voyage vers la Destination Finale

Le but ultime de l'eau dans la paille est d'atteindre une hauteur d'équilibre (appelée hauteur de Jurin). C'est le moment où la force qui tire l'eau vers le haut (la tension de surface) est exactement contrebalancée par le poids de l'eau qui tire vers le bas.

Deux façons d'arriver :

• La montée douce (Monotone) : L'eau monte lentement et s'arrête doucement, comme une voiture qui freine pour se garer.
• La montée en oscillation (Oscillatoire) : L'eau dépasse la hauteur idéale, redescend un peu, remonte, et finit par se stabiliser. C'est comme un ressort qu'on comprime et qui oscille avant de se caler.

La découverte clé :
Les auteurs ont proumé que le "glissement" (le patinage) ne change pas le fait que l'eau arrive à destination. Il change juste comment elle arrive (plus vite, plus lentement, ou avec plus d'oscillations).

• Si le glissement est fort, l'eau peut monter plus vite au début.
• Mais peu importe le glissement, l'eau finira toujours par s'arrêter à la bonne hauteur.

Ils ont même dessiné une "carte de sécurité" (le bassin d'attraction) qui garantit que, tant que l'eau commence à une hauteur raisonnable (pas trop haute, pas trop basse), elle ne s'échappera pas et finira par se stabiliser.

En Résumé

Ce papier est comme une mise à jour du manuel d'instructions pour comprendre comment les liquides montent dans les tuyaux fins.

• Avant : On pensait que l'eau collait toujours aux murs.
• Maintenant : On sait qu'elle peut glisser, et on a mis à jour les maths pour en tenir compte.
• Le résultat : Même avec ce glissement, les lois de la physique restent solides. L'eau monte, oscille peut-être un peu, mais finit toujours par trouver sa place au bon endroit. C'est une victoire pour la rigueur mathématique et pour notre compréhension de la nature !

Résumé technique

1. Problématique et Contexte

L'objectif principal de ce travail est d'étendre l'équation classique de Washburn, qui décrit la montée capillaire d'un liquide dans un tube vertical, en y intégrant une condition de glissement (slip condition) à la paroi du tube.

• Contexte physique : L'équation de Washburn (1921) combine les effets d'inertie, de viscosité, de gravité et de tension superficielle. Cependant, l'hypothèse classique de "non-glissement" (vitesse nulle à la paroi) crée un paradoxe (paradoxe de Huh–Scriven) lors du mouvement de la ligne de contact, car la vitesse du liquide à la paroi ne peut être nulle si la colonne de liquide monte.
• Problème mathématique : Des travaux antérieurs (notamment [9]) ont tenté de prouver l'existence et l'unicité globale des solutions, mais leurs preuves présentaient des lacunes critiques, notamment concernant la vérification des hypothèses des théorèmes de point fixe (Schauder et Banach) et l'application incorrecte du théorème de Picard-Lindelöf en raison de la non-lipschitzianité du terme non linéaire à $t=0$.
• Objectif : Dérivation rigoureuse de l'équation à partir des principes fondamentaux, preuve de l'existence et de l'unicité globale pour des données initiales plus générales, et analyse de la stabilité à long terme.

2. Méthodologie

A. Dérivation du modèle (Principes fondamentaux)

Les auteurs dérivent l'équation de Washburn non pas en appliquant directement la seconde loi de Newton à un volume de contrôle, mais en partant des lois de conservation de la masse et de la quantité de mouvement (équations de Navier-Stokes) pour un fluide newtonien incompressible.

• Hypothèses : Écoulement axisymétrique, profil de vitesse de Poiseuille modifié par un paramètre de glissement $L$ (longueur de glissement).
• Résultat : L'introduction du paramètre de glissement $\beta = (1 + 4L/R)^{-1}$ dans le terme visqueux. L'équation dimensionnelle obtenue est :
  $$\frac{d}{dt} \left[ \rho h(t) \frac{dh(t)}{dt} \right] + \frac{8\mu\beta}{R^2} h(t) \frac{dh(t)}{dt} + \rho g h(t) = \frac{2\gamma \cos \theta}{R}$$
  où $h(t)$ est la hauteur de la colonne liquide.

B. Analyse dimensionnelle et réduction de modèle

Une nouvelle mise à l'échelle est introduite pour obtenir une forme adimensionnelle :
$$\omega(HH')' + \beta HH' + H = 1$$
où $H$ est la hauteur adimensionnelle, $T$ le temps adimensionnel, et $\omega$ un paramètre lié aux nombres d'Ohnesorge et de Bond.

• Les auteurs analysent systématiquement quatre régimes d'écoulement (réductions du modèle) en fonction de la taille de $\omega$ (négligeant la gravité, l'inertie, ou la viscosité), illustrant comment le paramètre de glissement $\beta$ influence la dynamique (pente plus raide avec glissement).

C. Preuve d'existence et d'unicité (Analyse mathématique)

C'est la contribution majeure du papier. Pour contourner les problèmes de régularité à $t=0$ (terme $\sqrt{u}$ non lipschitzien), les auteurs utilisent une approche différente de celle de [9] :

1. Transformation : Introduction de la variable $u(s) = \frac{1}{2}H(T)^2$ pour transformer l'équation en une forme où la non-linéarité est mieux contrôlée.
2. Régularisation : Résolution d'un problème régularisé avec un paramètre $\epsilon > 0$ pour assurer la condition de Lipschitz.
3. Compacité : Utilisation du théorème d'Arzelà-Ascoli pour extraire une sous-suite convergente lorsque $\epsilon \to 0$, prouvant ainsi l'existence d'une solution globale.
4. Unicité : Démonstration de l'unicité pour des données initiales $\alpha = h_0/h_e \in [0, 3/2]$ en utilisant des inégalités de Grönwall et le théorème de point fixe de Precup (Théorème 11.3 de [10]) pour le cas singulier $\alpha=0$.
5. Bien-posé : Preuve de la dépendance continue des solutions par rapport aux données initiales (sens de Hadamard).

D. Analyse de stabilité

• Stabilité linéaire : Analyse des valeurs propres autour du point critique (hauteur d'équilibre). Le paramètre de glissement $\beta$ modifie la valeur critique $\omega^* = \beta^2/4$ qui sépare les comportements monotones ($\omega < \omega^*$) des comportements oscillatoires ($\omega > \omega^*$).
• Stabilité globale : Construction d'une fonction de Lyapunov adaptée pour prouver la stabilité asymptotique globale.
• Bassin d'attraction : Utilisation du principe d'invariance de LaSalle pour identifier le bassin d'attraction, garantissant que toutes les solutions partant des conditions initiales physiques ($0 \le h_0 \ll h_e$) convergent vers l'équilibre.

3. Résultats Clés

1. Existence et Unicité Globale : Il est prouvé qu'il existe une unique solution positive et bornée pour l'équation de Washburn avec glissement sur l'intervalle de temps $[0, \infty)$, pour tout $\beta > 0$ et pour un large éventail de données initiales ($\alpha \in [0, 3/2]$).
2. Correction des preuves antérieures : Les lacunes identifiées dans la littérature précédente (notamment dans [9]) concernant l'application des théorèmes de point fixe sont comblées par une approche rigoureuse basée sur la régularisation et la compacité.
3. Influence du glissement : Le paramètre de glissement $\beta$ n'affecte pas la hauteur d'équilibre (qui reste déterminée par la loi de Jurin), mais il modifie la dynamique transitoire :
   • Il accélère la montée initiale.
   • Il déplace la frontière entre les régimes monotones et oscillatoires (via $\omega^*$).
4. Convergence vers l'équilibre : Toute solution issue d'une hauteur initiale positive (même très petite) converge asymptotiquement vers la hauteur d'équilibre, que ce soit de manière monotone ou oscillatoire. Le bassin d'attraction inclut les conditions initiales physiques réalistes.

4. Signification et Impact

• Rigueur Mathématique : Ce travail fournit la première preuve mathématiquement complète et rigoureuse de l'existence et de l'unicité globale pour l'équation de Washburn avec conditions de glissement, corrigeant des erreurs subtiles dans les travaux antérieurs.
• Validité Physique : En intégrant la condition de glissement, le modèle résout le paradoxe de Huh-Scriven et offre une description plus fidèle de la dynamique initiale de la montée capillaire, en accord avec les observations expérimentales.
• Généralité : La méthodologie développée (transformation de variables, régularisation, principe de LaSalle) est robuste et pourrait être appliquée à d'autres problèmes de mécanique des fluides avec des singularités ou des conditions aux limites complexes.
• Perspectives : Les auteurs suggèrent des recherches futures sur l'analyse de valeurs négatives du paramètre de glissement (liées à des changements d'angle de contact) et le développement de méthodes numériques stables pour reproduire ces dynamiques à long terme.

En résumé, ce papier établit une fondation théorique solide pour l'étude de la montée capillaire avec glissement, combinant une dérivation physique rigoureuse à une analyse mathématique avancée des équations différentielles non linéaires.