On the stability of solutions to non-Newtonian Navier--Stokes--Fourier-like systems in the supercritical case

Technische Zusammenfassung: Zur Stabilität von Lösungen zu nicht-newtonschen Navier–Stokes–Fourier–ähnlichen Systemen im superkritischen Fall

Problemstellung
Die vorliegende Arbeit untersucht die mathematische Analyse eines dreidimensionalen, homogenen, inkompressiblen, nicht-newtonschen, wärmeleitenden Fluids. Das System wird durch ein modifiziertes Navier–Stokes–Fourier-System beschrieben, welches die Impulsgleichung, die Inkompressibilitätsbedingung und die Wärmegleichung mit einem Quellterm umfasst, der die Dissipation kinetischer Energie repräsentiert. Die Konstitutivrelation für den Cauchy-Spannungstensor ist allgemein gehalten und umfasst Potenzgesetze sowie Ladyzhenskaya-Modelle mit einem Potenzgesetz-Exponenten $p$ im Bereich $p \in [8/5, 11/5)$. Dieser Bereich wird als „superkritisch“ identifiziert – ein Regime, in dem Standardergebnisse zur Regularität, Eindeutigkeit und Energiegleichung nach aktuellem Kenntnisstand nicht gelten. Das Problem ist auf einer beschränkten Lipschitz-Domäne $\Omega \subset \mathbb{R}^3$ mit No-Slip-Randbedingungen für die Geschwindigkeit und einer vorgegebenen nicht-homogenen Temperatur (Dirichlet) am Rand formuliert. Das primäre Ziel besteht darin, die Existenz globaler Lösungen in der Zeit für beliebige Anfangsdaten mit endlicher Energie zu etablieren und die nichtlineare Stabilität des eindeutigen stationären Zustands $(0, \hat{\vartheta})$ zum Zeitpunkt $t \to \infty$ zu beweisen.

Methodik und neuartige Konzepte
Die Autoren adressieren den Mangel an Regularität im superkritischen Regime durch die Einführung eines neuartigen Lösungsbegriffs, der als $b$-schwache Lösung bezeichnet wird. Diese Formulierung weicht von klassischen schwachen Lösungen ab, indem sie die interne Energiegleichung durch eine Entropieungleichung ersetzt und entscheidend eine $b$-korrigierte totale Energieungleichung einführt.

1. Entropieungleichung: Anstelle der Wärmegleichung erfüllt die Lösung eine Entropieungleichung, die durch Division der Wärmegleichung durch die Temperatur $\vartheta$ gewonnen wird. Hierbei wird die Entropievariable $\eta = \ln \vartheta$ verwendet, welche über bessere Integrabilitäts-Eigenschaften als $\vartheta$ selbst verfügt.
2. $\gamma$-Korrekturfunktion: Um die nicht-homogenen Dirichlet-Randbedingungen für die Temperatur zu handhaben, welche die direkte Verwendung standardmäßiger totaler Energieungleichungen (aufgrund unkontrollierbarer Randfluss-Terme) verhindern, führen die Autoren eine $\gamma$-Korrekturfunktion $b(\vartheta, \hat{\vartheta})$ ein. Diese Funktion ermöglicht die Konstruktion eines modifizierten Energiefunktionals, das den Wärmefluss durch den Rand berücksichtigt.
3. $b$-korrigierte totale Energieungleichung: Durch Multiplikation der Wärmegleichung mit der Korrekturfunktion und Kombination mit der Impulsgleichung leiten die Autoren eine modifizierte Energiebilanz ab. Diese Ungleichung enthält einen Korrekturterm unter Einbeziehung des Primitivs von $b$, bezeichnet als $B$, und stellt sicher, dass die totale Energie (kinetisch plus thermisch, angepasst durch $B$) eine Differentialungleichung erfüllt, die für die Stabilitätsanalyse geeignet ist.
4. Existenzbeweis: Die Existenz von $b$-schwachen Lösungen wird mittels Faedo–Galerkin-Approximationen etabliert. Die Autoren leiten A-priori-Abschätzungen her, die unabhängig vom Approximationsparameter $n$ sind, unter Verwendung der kinetischen Energieungleichung, der Entropieungleichung und der $b$-korrigierten totalen Energieungleichung. Der Grenzübergang stützt sich auf die Monotonie des Spannungstensors, die solenoidale Lipschitz-Abschneidungsmethode zur Identifizierung des Spannungslimits sowie Kompaktheitsargumente (Aubin–Lions-Lemma).
5. Stabilitätsanalyse: Die Stabilität des stationären Zustands wird unter Verwendung eines Lyapunov-Funktionals $L_{\alpha, \beta}$ analysiert, welches den „Abstand“ zwischen der Lösung und dem stationären Zustand misst. Dieses Funktional bezieht die Hilfsfunktionen $H_\alpha$ und $G$ ein, um den nicht-konstanten Wärmeleitwert zu handhaben. Der Beweis erfolgt durch Renormierung der Entropieungleichung und Kombination mit der $b$-korrigierten Energieungleichung, um eine Differentialungleichung für das Lyapunov-Funktional abzuleiten, was zu exponentiellen Zerfallsschätzungen führt.

Zentrale Ergebnisse

• Theorem 1.3 (Globale Existenz): Für beliebige Anfangsdaten mit endlicher Energie ($v_0 \in L^2_{0,div}$, $\vartheta_0 \in L^1$) und eine Randtemperatur $\hat{\vartheta}$ mit Standardregularität existiert eine global-in-der-Zeit $b$-schwache Lösung. Die Lösung erfüllt eine uniforme globale Abschätzung in $L^\infty(0, +\infty; L^1(\Omega))$ für die Temperatur. Bemerkenswerterweise gilt dieses Ergebnis für Potenzgesetz-Exponenten $p > 6/5$, was vorangegangene Ergebnisse erweitert, die auf $p > 9/5$ beschränkt waren oder Neumann-Randbedingungen erforderten.
• Theorem 1.4 (Nichtlineare Stabilität): Der eindeutige stationäre Zustand $(0, \hat{\vartheta})$ ist nichtlinear stabil. Speziell gilt:
  • Das Geschwindigkeitsfeld zerfällt exponentiell in der Zeit: $\|v(\tau)\|_2^2 \leq e^{-\mu \tau} \|v(0)\|_2^2$ (mit spezifischen Modifikationen für große Anfangsdaten bei $p < 2$).
  • Die vollständige Lösung (Geschwindigkeit und Temperatur) konvergiert im Sinne des Lyapunov-Funktionals $L_{\alpha, \beta}$ gegen den stationären Zustand. Für jede beschränkte Menge von Anfangsdaten existieren Konstanten $\mu, \beta > 0$, sodass das Lyapunov-Funktional exponentiell zerfällt: $\int_\Omega L_{\alpha, \beta}(v(\tau), \vartheta(\tau), \hat{\vartheta}) dx \leq e^{-\mu(\tau-\sigma)} \int_\Omega L_{\alpha, \beta}(v(\sigma), \vartheta(\sigma), \hat{\vartheta}) dx$.

Bedeutung und Behauptungen
Das Paper beansprucht, das erste Ergebnis zu liefern, das gleichzeitig die Existenz globaler Lösungen in der Zeit und deren langfristige nichtlineare Stabilität für diese spezifische Klasse von nicht-newtonschen, wärmeleitenden Fluiden im superkritischen Regime mit nicht-homogenen Dirichlet-Randbedingungen etabliert.

Die Autoren betonen, dass vorangegangene Arbeiten zu ähnlichen Systemen entweder die Existenz glatter Lösungen voraussetzten, die Analyse auf subkritische Regime beschränkten, in denen klassische Energiegleichheiten gelten, oder auf Eindeutigkeitsergebnissen beruhten, die im superkritischen Fall für scherverdünnende Fluide nicht verfügbar sind. Die Einführung des Konzepts der $b$-schwachen Lösung wird als notwendige technische Weiterentwicklung präsentiert, um den Mangel an Energiegleichheit und die Herausforderungen durch die Randbedingungen zu bewältigen. Das Stabilitätsresultat wird als signifikanter Fortschritt beschrieben, der zeigt, dass trotz der potenziellen Nicht-Eindeutigkeit schwacher Lösungen die physikalisch relevante Klasse der $b$-schwachen Lösungen ein stabiles Langzeitverhalten aufweist, das gegen das eindeutige Gleichgewicht konvergiert. Die Analyse vermeidet explizit die Oberbeck–Boussinesq-Approximation und adressiert die volle Kopplung, bei der die kinetische Energiedissipation als Quelle in der Wärmegleichung wirkt – ein Merkmal, das im Vergleich zu vereinfachten Modellen erhebliche mathematische Schwierigkeiten einführt.