🔬 physics

The Statistical physics of unsaturated soil water: kinetic theory and non commutative pore water dynamics

Cet article présente un cadre de mécanique statistique multi-échelle qui unifie la dynamique de l'eau dans les sols non saturés en une équation de continuum cinétique unique, démontrant que les modèles classiques tels que l'équation de Richards et des phénomènes tels que le flux préférentiel et l'hystérésis émergent comme des limites spécifiques ou des conséquences géométriques d'un nombre de Damköhler résolu à l'échelle des pores régissant l'occupation non équilibrée des pores.

Auteurs originaux : Riccardo Rigon

Publié 2026-07-13

📖 7 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Riccardo Rigon

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 ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Imaginez le sol non pas comme une masse boueuse, mais comme un gigantesque labyrinthe en 3D composé de milliards de minuscules tubes creux de différentes largeurs. Certains sont aussi fins qu'un cheveu humain ; d'autres sont aussi larges qu'une paille à boire. À l'intérieur de ce labyrinthe, l'eau cherche son chemin.

Pendant plus d'un siècle, les scientifiques ont tenté de décrire comment cette eau se déplace en utilisant une règle unique et célèbre appelée l'équation de Richards. Considérez cette équation comme un livre de règles au "ralenti". Elle suppose que si vous versez de l'eau sur le sol, l'eau trouve instantanément l'endroit le plus confortable et le plus économe en énergie pour chaque goutte. Elle suppose que l'eau est toujours parfaitement satisfaite, toujours en équilibre, et que la seule chose qui importe est la quantité d'eau présente dans le sol.

La Grande Découverte : Le Sol a une Mémoire (et une Limite de Vitesse)

Cet article, par Riccardo Rigon, soutient que l'ancien livre de règles est incomplet. C'est comme essayer de décrire une fosse de mosh chaotique en ne regardant que la densité moyenne de la foule. On passe à côté du fait que les gens courent, se bousculent et restent coincés dans des endroits spécifiques.

L'article introduit une vision plus détaillée. Au lieu de simplement demander « À quel point le sol est-il humide ? », il demande : « Quels tubes spécifiques sont pleins et lesquels sont vides ? »

Les auteurs proposent une nouvelle théorie « cinétique ». Ils traitent l'eau dans le sol comme une foule de personnes dans un bâtiment comprenant de nombreuses pièces.

L'échelle micro (le couloir) : L'eau se déplace entre les tubes en fonction de la physique (taux de Hagen–Poiseuille). C'est comme de l'eau s'écoulant dans un tuyau.
L'échelle méso (la pièce) : Lorsque vous regardez une pièce entière (un minuscule fragment de sol), l'eau échange constamment sa place, essayant de s'installer dans l'arrangement le plus confortable.
L'échelle macro (le bâtiment) : Lorsque vous zoomez sur tout le champ, vous voyez l'eau passer d'une pièce à l'autre.

Le « Nombre de Damköhler » : Le Compteur de Vitesse du Sol

La partie la plus excitante de cet article est un nombre magique unique qu'ils appellent le nombre de Damköhler (Da). Considérez cela comme un « Compteur de vitesse » pour le sol. Il compare deux vitesses :

La vitesse à laquelle l'eau veut se réorganiser pour trouver la place parfaite (relaxation).
La vitesse à laquelle la pluie frappe le sol (force imposée).

Si la pluie est lente (Da est petit) : L'eau a tout le temps de se réorganiser. Elle trouve les endroits parfaits, tout comme le prédisait l'ancienne équation de Richards. Le sol se comporte comme un fluide calme et lent.
Si la pluie est une averse soudaine (Da est grand) : L'eau est déversée plus vite qu'elle ne peut se réorganiser. Elle n'a pas le temps de trouver les minuscules tubes confortables. Au lieu de cela, elle se bloque dans les gros tubes rapides, créant un « écoulement préférentiel » (comme de l'eau jaillissant d'une lance à incendie).

Le Mystère de l'« Écoulement Préférentiel » Résolu ?

Pendant longtemps, les scientifiques ont pensé que l'« écoulement de Richards » (lent, régulier) et l'« écoulement préférentiel » (rapide, chaotique) étaient deux choses totalement différentes nécessitant deux livres de règles différents.

Cet article soutient que NON. Ce sont le même équation, se comportant différemment selon la vitesse de la pluie.

L'écoulement de Richards est ce qui se passe quand le sol a le temps de reprendre son souffle.
L'écoulement préférentiel est ce qui se passe quand le sol est submergé.

Ce n'est pas que le sol change d'avis ; c'est que la vitesse de l'événement change le résultat. L'article suggère que si vous connaissez la vitesse de la pluie et la taille des tubes, vous pouvez prédire exactement quand le sol passera de l'état « calme » à l'état « chaotique ».

La Boucle d'« Hystérésis » : Un Fantôme Géométrique

Vous avez probablement remarqué que le sol retient l'eau différemment lorsqu'il s'assèche par rapport à lorsqu'il s'humidifie. C'est ce qu'on appelle l'hystérésis.

Vieille idée : Les scientifiques pensaient que cela était dû au fait que le sol possédait des tubes « bistables » — comme des interrupteurs qui restent bloqués sur « on » ou « off » selon l'historique.
Nouvelle idée : L'article soutient que cela n'est pas dû à des interrupteurs bloqués. C'est une question de géométrie et de temps.

Imaginez marcher sur un sentier dans une forêt. Si vous marchez lentement, vous pouvez enjamber chaque branche. Si vous courez, vous pourriez trébucher sur une branche et prendre un chemin différent. Si vous marchez en cercle (humidification puis assèchement), vous pourriez ne pas revenir exactement là où vous étiez parti parce que le chemin de la « course » était différent de celui de la « marche ».

Les auteurs appellent cela une phase géométrique (ou holonomie). C'est comme un « fantôme géométrique ». Le sol se souvient du chemin qu'il a parcouru pour atteindre son humidité actuelle, et pas seulement de l'humidité elle-même. L'article suggère que la « taille » de cette mémoire (la boucle d'hystérésis) croît avec le carré de l'intensité de la pluie. Si vous doublez la vitesse de la pluie, l'effet de mémoire est multiplié par quatre. C'est une prédiction spécifique et testable que les anciens modèles d'« interrupteurs bloqués » ne peuvent expliquer.

Ce qui est Prouvé, ce qui est Suggéré et ce qui est Exclu

Exclu : L'article s'oppose explicitement à l'idée que l'« écoulement préférentiel » et l'« écoulement de Richards » soient deux phénomènes distincts et sans lien ; ils constituent une seule et même transition continue. Il rejette également l'idée que l'hystérésis soit causée par des pores individuels restant « bloqués » dans un état binaire ; c'est plutôt le résultat de la géométrie du réseau global et de la vitesse de la force imposée.
Prouvé (dans les mathématiques de l'article) : Les auteurs ont prouvé mathématiquement que leur nouvelle équation respecte les lois de la physique (comme la conservation de la masse) et que le système tend toujours vers un état d'énergie plus basse (le théorème H). Ils ont montré que leur nouvelle équation peut mathématiquement se transformer en l'ancienne équation de Richards si l'on ralentit suffisamment les processus.
Suggéré / Simulé : La prédiction spécifique selon laquelle la surface de la boucle d'hystérésis croît avec le carré de l'intensité de la pluie ( $H \propto I^2$ ) est une conjecture basée sur leur théorie. Ils suggèrent que cela pourrait être testé en laboratoire. Ils suggèrent également que la « capacité au champ » (le point où le sol cesse de drainer) est en réalité un « seuil de percolation » (un point où le réseau d'eau se brise), une idée plausible soutenue par leur modèle mais qui nécessite une vérification plus concrète sur le terrain.
Simulé uniquement : L'idée que la « non-commutativité » (le fait que l'humidification suivie de l'assèchement donne un résultat différent de l'assèchement suivi de l'humidification) est la cause profonde de l'hystérésis est démontrée dans leurs calculs et dans des simulations informatiques d'un petit réseau de tubes.

L'Essentiel

Cet article ne propose pas seulement une nouvelle formule ; il offre une nouvelle façon de voir le monde. Il suggère que le sol n'est pas une éponge statique, mais un labyrinthe dynamique et sensible au temps. Que l'eau s'écoule lentement et uniformément ou qu'elle jaillisse en courants chaotiques dépend entièrement de la vitesse de la pluie par rapport à la vitesse de réorganisation interne du sol.

Les auteurs sont convaincus que leur nouvelle équation est la « clé maîtresse » qui déverrouille à la fois l'écoulement lent et régulier et l'écoulement rapide et chaotique, les unifiant en un récit unique, beau et mathématiquement rigoureux. Ils invitent les scientifiques à tester leur prédiction : si l'on mesure la boucle d'hystérésis à différentes vitesses de pluie, elle devrait croître de manière quadratique. Si c'est le cas, le « fantôme géométrique » est réel. Si ce n'est pas le cas, la théorie doit être réécrite.

Résumé Technique : Physique Statistique de l'Eau dans les Sols Non Saturés

Énoncé du Problème
La dynamique des milieux poreux non saturés est intrinsèquement irréversible, caractérisée par des configurations hors équilibre qui dépendent de l'historique du forçage (mouillage vs assèchement). Les modèles classiques, principalement l'équation de Richards, supposent un équilibre thermodynamique local et une variable d'état scalaire unique (teneur en eau $\theta$ ou potentiel matriciel $\psi$ ). Ce cadre ne parvient pas à expliquer :

L'hystérésis : Des états macroscopiques identiques correspondant à des configurations internes différentes selon le chemin de mouillage ou d'assèchement.
Le Flux Préférentiel : La coexistence d'une redistribution lente, contrôlée par capillarité, et d'un flux de contournement rapide à travers les grands pores lors de forçages intenses, ce qui viole les hypothèses de domaine unique de l'équation de Richards.
La Dépendance au Taux : L'incapacité des modèles indépendants du taux à capturer la dépendance de l'hystérésis et des régimes d'écoulement vis-à-vis de la vitesse du forçage.

Les modèles existants à double domaine ou à double perméabilité traitent ces problèmes de manière phénoménologique en introduisant des variables d'état supplémentaires et des termes d'échange sans les dériver d'une description physique commune.

Méthodologie
L'article développe une théorie de mécanique statistique en passant par trois échelles, en traitant la fraction d'occupation des pores $g(r, x, t)$ (la fraction des pores de rayon $r$ remplis d'eau à la position $x$ et au temps $t$ ) comme la variable d'état fondamentale.

Échelle Micro (Taux Cinétiques) :
- Les transferts d'eau inter-pores sont régis par des taux de Hagen–Poiseuille $\kappa(r)$ et un potentiel moteur $\Phi(r, r')$ , défini comme la différence de potentiels chimiques entre les classes de pores (incorporant les termes capillaires et gravitationnels).
- Une constante de taux dépendante de l'état $\kappa(r, g)$ est dérivée, transitant de la dynamique de Lucas–Washburn (invasion initiale rapide) vers un écoulement de Hagen–Poiseuille (état stationnaire) à mesure que la fraction de remplissage $g$ augmente.
Échelle Méso (Équation Maître) :
- Le passage à un Volume Élémentaire Représentatif (VER) produit une équation maître pour $g$ .
- L'équation inclut des termes de gain ( $G$ ) et de perte ( $L$ ) représentant l'échange d'eau entre les classes de pores. Ces termes utilisent un « Stosszahlansatz » (hypothèse du chaos moléculaire), supposant l'indépendance statistique des occupations du donneur et du receveur, conduisant à une forme bilinéaire $g(r')[1-g(r)]$ .
- Le système est démontré être un flux de gradient d'Onsager, descendant un fonctionnel d'énergie libre de Gibbs $F[g]$ avec une mobilité dépendante de l'état. Cela garantit le théorème H ( $dF/dt \leq 0$ ) pour les systèmes isothermes et non forcés.
Échelle Continuum (Équation de Champ) :
- La contraction du VER en un point promeut $g$ au statut de champ continu.
- La loi de bilan résultante est :
  $\partial_t g + \nabla \cdot \mathbf{F} = C[g] - E - T$
  où $C[g]$ est l'opérateur de redistribution intra-VER, $\nabla \cdot \mathbf{F}$ est la divergence de transport (flux inter-VER), et $E, T$ sont des puits (évaporation, absorption racinaire).
- Le flux $\mathbf{F}$ est non local en rayon de pore, couplant les différentes classes de pores via la matrice de connectivité $C(r, r')$ .

Contributions Clés et Résultats

Unification des Régimes d'Écoulement : La théorie introduit un paramètre de contrôle dimensionnel unique, le nombre de Damköhler résolu par pore $Da(r, x) = \tau_{redis}(r) / \tau_{forcing}(x)$ , qui compare le temps de relaxation interne au temps d'échelle du forçage.
- $Da \ll 1$ (Quasi-statique) : Le système suit l'équilibre local. La réduction de Chapman–Enskog de l'équation cinétique recupère l'équation de Richards, avec le potentiel matriciel $\psi$ et la conductivité hydraulique $K(\theta)$ émergeant comme des limites.
- $Da \gtrsim 1$ (Cinétique) : La redistribution ne peut plus suivre le rythme du forçage. L'eau occupe les voies géométriquement accessibles (grands pores) plutôt que les voies thermodynamiquement optimales, reproduisant naturellement le flux préférentiel sans hypothèses de double domaine ad hoc.
Origine Géométrique de l'Hystérésis : L'hystérésis n'est pas attribuée à une bistabilité par pore (ex: effets de type "goulot d'encre"), mais à l'holonomie d'un « faisceau de forçage ».
- Les opérateurs de mouillage ( $L_W$ ) et de drainage ( $L_D$ ) ne commutent pas ( $[L_W, L_D] \neq 0$ ) en raison de l'asymétrie entre le remplissage non sélectif (précipitation) et le vidage sélectif (drainage/évaporation contraint par la topologie).
- Cette non-commutativité crée une phase géométrique (état résiduel) après un cycle de forçage fermé.
- Prédiction : L'aire de la boucle d'hystérésis $H$ suit une loi quadratique par rapport à l'intensité du forçage ( $H \propto I^2$ ) pour les cycles lents, une signature absente des modèles indépendants du taux.
Lois Constitutives issues de la Microstructure :
- La conductivité hydraulique $K(\theta)$ est dérivée comme un coefficient de transport (limite spectrale du noyau de mobilité) plutôt que comme un paramètre ajusté.
- La théorie retrouve les limites classiques : modèles de faisceaux capillaires (limite diagonale), modèles de chemin critique (limite spectrale près de la percolation), et modèles de pression capillaire dynamique (approximation du premier moment).
- La capacité au champ est identifiée comme le seuil de percolation où le réseau d'eau se fragmente ( $\theta_c$ ), provoquant la disparition de la conductivité macroscopique.
Variable d'État Irréductible : En dehors de la limite quasi-statique, la teneur en eau scalaire $\theta$ est insuffisante pour caractériser le système. La distribution complète $g(r, x, t)$ est la variable d'état irréductible.

Signification et Revendications
L'article affirme fournir un cadre constructif et fondé sur la thermodynamique qui unifie des phénomènes auparavant distincts (écoulement de Richards, écoulement préférentiel, hystérésis) en une seule équation cinétique.

Rigueur Théorique : La théorie est présentée comme un analogue des milieux poreux de l'équation de Boltzmann, où $g(r)$ correspond à la distribution des vitesses $f(v)$ , et la réduction de Chapman–Enskog reflète la dérivation des équations de Navier–Stokes.
Pouvoir Prédictif : Toutes les entrées (distribution de taille de pores $f(r)$ , connectivité $C(r, r')$ , tortuosité $\bar{\tau}$ ) sont des propriétés géométriques mesurables par imagerie micro-CT, ne nécessitant aucun étalonnage par rapport aux données de rétention ou de conductivité macroscopiques.
Falsifiabilité : La théorie propose des prédictions spécifiques et testables, notamment la mise à l'échelle $H \propto I^2$ de l'aire de la boucle d'hystérésis avec l'intensité du forçage, ainsi que l'existence d'un seuil d'intensité dérivé pour l'apparition du flux préférentiel basé sur la structure du sol et les propriétés du fluide.
Portée : Le cadre est généralisable à d'autres systèmes présentant des transitions contraintes dépendant de l'historique (ex: déformation plastique, hystérésis magnétique) et peut être étendu pour inclure des termes d'adsorption, d'osmose ou thermiques en affinant le potentiel moteur $\Phi$ sans modifier l'algèbre des opérateurs de base.

L'article conclut que la dichotomie entre écoulement à l'équilibre et hors équilibre n'est pas une question de choix de modèles différents, mais la reconnaissance d'un passage continu contrôlé par le nombre de Damköhler au sein d'une théorie cinétique unique et unifiée.

Articles similaires