Thermodynamics and stability of equilibrium/non-equilibrium steady states in thermodynamically isolated/open systems -- case study for compressible heat conducting fluid

Cet article examine les calculs nécessaires à la construction d'un fonctionnel de type Lyapunov pour analyser la stabilité non linéaire des états stationnaires, à l'équilibre ou non, dans des systèmes isolés ou ouverts composés de fluides compressibles conducteurs de chaleur.

L'essentiel

Le Grand Équilibre : Comment un fluide chaud trouve son calme

Imaginez que vous êtes un chef cuisinier. Vous avez une grande casserole remplie d'eau (un fluide). Parfois, vous la laissez reposer sur le feu éteint (système isolé), et parfois, vous maintenez le fond de la casserole chaud tout en laissant le haut refroidir (système ouvert).

La question que se pose l'auteur, Vít Průša, est simple mais profonde : Si je perturbe cette eau (en la remuant ou en changeant sa température), va-t-elle finir par retrouver son état calme et stable, ou va-t-elle devenir chaotique ?

Pour répondre à cela, il ne suffit pas de dire "oui, c'est logique". Il faut le prouver mathématiquement avec des outils très puissants. Voici comment il y arrive, étape par étape.

1. Le Problème : Le Chaos vs. Le Calme

Dans la vie de tous les jours, si vous versez un peu de lait dans du café, il finit par se mélanger uniformément. Si vous agitez l'eau dans une casserole isolée, elle finit par s'arrêter et atteindre une température uniforme. C'est ce qu'on appelle un état d'équilibre.

Mais les mathématiques sont difficiles. Comment prouver qu'un système complexe (comme l'air dans une pièce ou l'eau dans un tuyau) reviendra toujours au calme, même si on le secoue violemment ?

2. L'Outil Magique : La "Boussole de l'Énergie" (Fonctionnelle de Lyapunov)

Pour prouver la stabilité, les mathématiciens utilisent un outil appelé une fonctionnelle de Lyapunov.

• L'analogie de la colline : Imaginez que l'état de votre fluide est une bille roulant sur une colline.
  • Le bas de la colline est l'état calme (l'équilibre).
  • Le haut de la colline est l'état agité (perturbé).
  • La gravité fait rouler la bille vers le bas.

L'auteur construit une "boussole" mathématique (une fonction spéciale) qui mesure la "hauteur" de la bille.

• Si la bille est en haut, la boussole indique une valeur élevée.
• Si la bille est en bas, la boussole indique zéro.
• Le secret : Cette boussole doit toujours diminuer avec le temps (comme la bille qui descend). Si on peut prouver que cette valeur diminue toujours jusqu'à zéro, alors on sait que le fluide va forcément revenir au calme.

3. Le Cas 1 : La Casserole Fermée (Système Isolé)

Imaginez une boîte parfaitement isolée. Rien n'entre, rien ne sort.

• Le défi : L'auteur doit créer cette "boussole" en combinant deux choses : l'énergie totale (la vitesse du fluide + sa chaleur) et l'entropie (le désordre).
• La découverte : Il découvre que pour que la boussole fonctionne, le fluide doit respecter certaines règles physiques de base (comme le fait que la chaleur se propage toujours du chaud vers le froid, et que la pression augmente quand on compresse un gaz).
• Le résultat : Il prouve que, tant que ces règles physiques sont respectées, la "boussole" descendra toujours vers zéro. Le fluide finira inévitablement par s'arrêter et avoir une température uniforme. C'est la stabilité.

4. Le Cas 2 : La Casserole Ouverte (Système Non-Équilibré)

Maintenant, imaginez que le fond de la casserole est chauffé et le haut refroidi. L'eau ne sera jamais au même endroit ni à la même température. Il y a un courant constant (convection). C'est un état stationnaire non-équilibre.

• Le défi : La bille ne peut pas aller au "bas de la colline" car il n'y a pas de bas unique. La colline est déformée par la chaleur qui rentre et sort.
• L'astuce géniale (Le "Truc de Correction") : L'auteur utilise une astuce mathématique appelée "correction affine".
  • Imaginez que vous prenez la "boussole" du cas fermé (où tout est calme) et que vous la "tordiez" légèrement pour qu'elle s'adapte à la nouvelle forme de la colline (celle où il y a du courant).
  • Il utilise des variables intelligentes (comme la quantité de mouvement au lieu de la vitesse) pour que cette tordure fonctionne parfaitement.
• Le résultat : Même dans ce cas complexe, il montre qu'il existe une "boussole" qui diminue toujours. Cela prouve que si vous perturbez ce courant d'eau, il finira par revenir à son état de flux stable, et non pas devenir chaotique.

5. Pourquoi c'est important ?

Ce papier est comme un manuel de sécurité pour les ingénieurs.

• Si vous concevez un moteur, un réacteur nucléaire ou un système de climatisation, vous voulez être sûr que si vous faites une erreur de calcul ou si un capteur dérape, le système ne va pas exploser ou devenir instable.
• L'auteur nous dit : "Ne vous inquiétez pas. Si les lois de la physique de base sont respectées (ce qui est le cas pour la plupart des fluides réels), la nature a une façon intrinsèque de revenir à l'ordre."

En résumé

L'auteur a construit un outil mathématique de mesure de la "distance" vers le calme.

1. Il a prouvé que cet outil fonctionne pour les systèmes fermés (comme une boîte scellée).
2. Il a utilisé une astuce intelligente pour adapter cet outil aux systèmes ouverts (comme un tuyau avec de l'eau qui coule).
3. Il a démontré que cet outil diminue toujours avec le temps, prouvant ainsi que la stabilité est garantie par les lois de la thermodynamique.

C'est une victoire des mathématiques pour confirmer ce que notre intuition nous dit depuis toujours : le chaos finit toujours par se calmer, à condition que les règles du jeu (la physique) soient respectées.

Résumé technique

1. Problématique

L'article aborde la question fondamentale de la stabilité asymptotique des états stationnaires (d'équilibre ou hors équilibre) pour des fluides compressibles conducteurs de chaleur, régis par les équations de Navier-Stokes-Fourier (NSF).

Le travail se concentre sur deux scénarios distincts :

• Question 1 (Systèmes isolés) : La stabilité d'un état de repos homogène spatiallement dans un récipient thermodynamiquement isolé (pas d'échange de masse, de quantité de mouvement ou de chaleur avec l'extérieur). L'hypothèse est que toute perturbation initiale, pourvu qu'elle conserve la masse et l'énergie totales, doit converger vers l'état d'équilibre homogène.
• Question 2 (Systèmes ouverts) : La stabilité d'un état stationnaire inhomogène spatiallement dans un récipient mécaniquement isolé mais thermiquement ouvert (échange de chaleur via des conditions aux limites de Dirichlet non nulles). Ici, l'objectif est de montrer la convergence vers un état stationnaire non uniforme (gradient de température).

Le défi principal réside dans la construction d'une fonctionnelle de type Lyapunov pour des systèmes non linéaires et continus, capable de contrôler la taille des perturbations (densité, vitesse, température) et de garantir leur décroissance temporelle.

2. Méthodologie

L'auteur utilise une approche rigoureuse combinant la thermodynamique des milieux continus et l'analyse fonctionnelle :

• Analyse linéarisée : Une première analyse est menée sur les équations linéarisées autour de l'état de repos. Cela permet d'identifier les conditions de stabilité thermodynamique nécessaires (positivité de la chaleur spécifique à volume constant et de la dérivée de la pression par rapport à la densité).
• Construction de fonctionnelles de Lyapunov (Non-linéaire) :
  • L'auteur rejette l'entropie nette seule comme fonctionnelle de Lyapunov, car elle ne contrôle pas la taille des perturbations (elle peut être nulle même si le système n'est pas à l'équilibre).
  • Il propose une fonctionnelle basée sur le principe de maximisation de l'entropie sous contraintes (énergie totale et masse constantes). Cette fonctionnelle est construite en introduisant des multiplicateurs de Lagrange ($\lambda_1, \lambda_2$) identifiés via des dérivées de Gâteaux.
  • Technique de correction affine : Pour les systèmes ouverts (état inhomogène), l'auteur applique une « astuce de correction affine » (affine correction trick). Il construit la fonctionnelle pour l'état hors équilibre à partir de la fonctionnelle connue pour l'état d'équilibre homogène, en soustrayant la valeur de l'état stationnaire et son terme linéaire (dérivée première).
• Lien avec la distance de Bregman : Une contribution majeure de la méthodologie est l'identification de la partie thermodynamique de la fonctionnelle comme une distance de Bregman générée par l'énergie interne modifiée (convexe). Cela permet de prouver rigoureusement la positivité de la fonctionnelle.
• Variables naturelles : L'analyse privilégie l'utilisation de la densité ($\rho$) et de l'entropie spécifique par unité de volume ($\eta = \rho s$) ou de l'impulsion ($p = \rho v$) plutôt que la vitesse, car ces variables rendent les termes d'énergie cinétique et interne convexes.

3. Contributions Clés

1. Identification explicite des multiplicateurs de Lagrange : L'article fournit une dérivation complète pour identifier les multiplicateurs $\lambda_1 = 1/\hat{\theta}$ et $\lambda_2 = -(\hat{p}_{th}/\hat{\rho} + \hat{\psi})/\hat{\theta}$, assurant que la fonctionnelle s'annule à l'équilibre et a une dérivée nulle.
2. Preuve de non-négativité via la convexité : L'auteur démontre que la fonctionnelle de Lyapunov est non négative et s'annule si et seulement si le système est à l'état stationnaire. Cette preuve repose sur la convexité stricte de l'énergie interne modifiée $e(\eta, \rho)$ par rapport à ses variables naturelles, ce qui est équivalent aux conditions de stabilité thermodynamique classiques ($c_V > 0$ et $\partial p_{th}/\partial \rho > 0$).
3. Unification des cas isolés et ouverts : La méthode de la « correction affine » permet de traiter les états stationnaires inhomogènes (systèmes ouverts) en utilisant la structure de la fonctionnelle des systèmes isolés, généralisant ainsi l'approche de stabilité.
4. Lien avec le travail de Feireisl : L'article montre que la fonctionnelle construite est équivalente (à une constante près liée à la conservation de la masse) à la fonctionnelle d'entropie relative (ou énergie balistique) utilisée par Feireisl et ses collaborateurs dans l'analyse d'unicité faible-forte des solutions des équations NSF.

4. Résultats Principaux

• Stabilité des systèmes isolés : Il est prouvé que pour un fluide compressible conducteur de chaleur dans un récipient isolé, tout état initial compatible (même masse et énergie) converge vers l'état de repos homogène. La fonctionnelle de Lyapunov $V_{meq}$ décroît strictement dans le temps :
  $$\frac{d}{dt} V_{meq} = -\hat{\theta} \int_{\Omega} \left( \frac{\tilde{\lambda}(\text{div } v)^2 + 2\nu \mathbb{D}^{\delta} : \mathbb{D}^{\delta} + \kappa \frac{|\nabla \theta|^2}{\theta}}{\theta} \right) dv \leq 0$$
  L'égalité n'est atteinte qu'à l'équilibre.

• Stabilité des systèmes ouverts : Pour un système avec des conditions aux limites de température fixes (inhomogènes), la fonctionnelle modifiée $V_{non-eq}$ (obtenue par correction affine) est également non négative et sa dérivée temporelle est non positive, garantissant la stabilité de l'état stationnaire inhomogène.

• Cas du gaz parfait caloriquement parfait : L'article fournit des expressions explicites de la fonctionnelle pour ce cas particulier, confirmant que les termes sont des divergences de Bregman classiques (formes de type $x \ln x - x + 1$).

5. Signification et Impact

Ce travail est significatif car il fournit une justification thermodynamique rigoureuse de la stabilité non linéaire des fluides compressibles, reliant directement les propriétés de stabilité des équations aux lois de la thermodynamique (convexité des potentiels thermodynamiques).

• Rigueur mathématique : Il comble le fossé entre l'intuition physique (le système tend vers l'équilibre) et la preuve mathématique formelle pour des systèmes continus non linéaires.
• Généralité : La méthode de construction de la fonctionnelle via la maximisation de l'entropie sous contraintes et la correction affine est applicable à une large classe de matériaux et de conditions aux limites, au-delà du cas spécifique traité ici.
• Outils pour l'analyse future : L'identification de la fonctionnelle comme une distance de Bregman ouvre la voie à l'utilisation d'outils puissants de l'analyse convexe pour l'étude de la convergence des solutions, de l'unicité et de la régularité des équations de Navier-Stokes-Fourier.

En résumé, l'article établit que les conditions de stabilité thermodynamique classiques ($c_V > 0$, compressibilité positive) sont non seulement nécessaires pour la stabilité linéaire, mais suffisantes pour garantir la stabilité non linéaire globale via la construction explicite d'une fonctionnelle de Lyapunov adaptée.