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

Technische Samenvatting: Over de stabiliteit van oplossingen van niet-Newtoniaanse Navier–Stokes–Fourier-achtige systemen in het superkritische geval

Probleemstelling
Het artikel onderzoekt de wiskundige analyse van een driedimensionaal, homogeen, incompressibel, niet-Newtoniaans, warmtegeleidend fluïdum. Het systeem wordt beheerst door een gemodificeerd Navier–Stokes–Fourier-systeem bestaande uit de impulsvergelijking, de incompressibiliteitsvoorwaarde en de warmtevergelijking met een bronterm die de dissipatie van kinetische energie vertegenwoordigt. De constitutieve relatie voor de Cauchy-spanningstensor is algemeen, waarbij machtswet- en Ladyzhenskaya-modellen worden omvat met een machtswet-exponent $p$ in het bereik $p \in [8/5, 11/5)$. Dit bereik wordt geïdentificeerd als "superkritisch", een regime waarin standaard resultaten voor regulariteit, uniciteit en energiegelijkheid niet bekend zijn. Het probleem is opgesteld op een begrensde Lipschitz-domein $\Omega \subset \mathbb{R}^3$ met no-slip randvoorwaarden voor snelheid en een voorgeschreven niet-uniforme temperatuur (Dirichlet) op de rand. Het primaire doel is het vaststellen van de existentie van globale oplossingen in de tijd voor willekeurige initiële data met eindige energie en het bewijzen van de nietlineaire stabiliteit van de unieke stationaire toestand $(0, \hat{\vartheta})$ naarmate de tijd naar oneindig gaat.

Methodologie en Nieuwe Concepten
De auteurs pakken het gebrek aan regulariteit in het superkritische regime aan door een nieuw concept van oplossing te introduceren, de $b$-zwakke oplossing. Deze formulering wijkt af van klassieke zwakke oplossingen door de interne energiegelijkheid te vervangen door een entropie-ongelijkheid en, cruciaal, door een $b$-gecorrigeerde totale energie-ongelijkheid te introduceren.

1. Entropie-ongelijkheid: In plaats van de warmtevergelijking voldoet de oplossing aan een entropie-ongelijkheid afgeleid door de warmtevergelijking te delen door temperatuur $\vartheta$. Hiermee wordt gebruikgemaakt van de entropievariabele $\eta = \ln \vartheta$, die betere integreerbaarheidseigenschappen bezit dan $\vartheta$ zelf.
2. $\gamma$-Correctiefunctie: Om de niet-homogene Dirichlet-randvoorwaarden voor temperatuur te behandelen, die direct gebruik van standaard totale energie-ongelijkheden verhinderen (door oncontroleerbare randflux-termen), introduceren de auteurs een $\gamma$-correctiefunctie $b(\vartheta, \hat{\vartheta})$. Deze functie maakt de constructie van een gemodificeerde energiefunctionaal mogelijk die rekening houdt met de warmteflux door de rand.
3. $b$-Gecorrigeerde Totale Energie-ongelijkheid: Door de warmtevergelijking te vermenigvuldigen met de correctiefunctie en deze te combineren met de impulsvergelijking, leiden de auteurs een gemodificeerde energiebalans af. Deze ongelijkheid bevat een correctieterm met betrekking tot de primitieve van $b$, aangeduid als $B$, en zorgt ervoor dat de totale energie (kinetisch plus thermisch, aangepast door $B$) een differentiële ongelijkheid vervult die geschikt is voor stabiliteitsanalyse.
4. Existentiebewijs: De existentie van $b$-zwakke oplossingen wordt vastgesteld via Faedo–Galerkin-benaderingen. De auteurs leiden a priori schattingen af die onafhankelijk zijn van de benaderingsparameter $n$, gebruikmakend van de kinetische energie-ongelijkheid, de entropie-ongelijkheid en de $b$-gecorrigeerde totale energie-ongelijkheid. De limietovergang berust op de monotoniciteit van de spanningstensor, de solenoidale Lipschitz-truncatiemethode om de spanninglimiet te identificeren, en compactheidsargumenten (Aubin–Lions lemma).
5. Stabiliteitsanalyse: De stabiliteit van de stationaire toestand wordt geanalyseerd met behulp van een Lyapunov-functionaal $L_{\alpha, \beta}$, die de "afstand" tussen de oplossing en de stationaire toestand meet. Dit functionaal incorporeert de hulpfuncties $H_\alpha$ en $G$ om niet-constante warmtegeleidingscoëfficiënten te behandelen. Het bewijs omvat het renormaliseren van de entropie-ongelijkheid en het combineren ervan met de $b$-gecorrigeerde energie-ongelijkheid om een differentiële ongelijkheid voor het Lyapunov-functionaal af te leiden, wat leidt tot exponentiële verval-schattingen.

Kernresultaten

• Theorem 1.3 (Globale Existentie): Voor elke initiële data met eindige energie ($v_0 \in L^2_{0,div}$, $\vartheta_0 \in L^1$) en randtemperatuur $\hat{\vartheta}$ die aan standaard regulariteit voldoet, bestaat er een globale oplossing in de tijd van het type $b$-zwak. De oplossing voldoet aan een uniforme globale schatting in $L^\infty(0, +\infty; L^1(\Omega))$ voor de temperatuur. Opmerkelijk is dat dit resultaat geldt voor machtswet-exponenten $p > 6/5$, waarmee eerdere resultaten die beperkt waren tot $p > 9/5$ of Neumann-randvoorwaarden werden uitgebreid.
• Theorem 1.4 (Nietlineaire Stabiliteit): De unieke stationaire toestand $(0, \hat{\vartheta})$ is nietlineair stabiel. Specifiek:
  • Het snelheidveld vervalt exponentieel in de tijd: $\|v(\tau)\|_2^2 \leq e^{-\mu \tau} \|v(0)\|_2^2$ (met specifieke modificaties voor grote initiële data wanneer $p < 2$).
  • De volledige oplossing (snelheid en temperatuur) convergeert naar de stationaire toestand in de zin van het Lyapunov-functionaal $L_{\alpha, \beta}$. Voor elke begrensde verzameling van initiële data bestaan er constanten $\mu, \beta > 0$ waarvoor het Lyapunov-functionaal exponentieel vervalt: $\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$.

Betekenis en Claims
Het artikel claimt het eerste resultaat te leveren dat simultaan de existentie van globale oplossingen in de tijd en hun lange-termijn nietlineaire stabiliteit vaststelt voor deze specifieke klasse van niet-Newtoniaanse, warmtegeleidende fluïda in het superkritische regime met niet-homogene Dirichlet-randvoorwaarden.

De auteurs benadrukken dat eerdere werken over soortgelijke systemen ofwel de existentie van gladde oplossingen aannamen, ofwel de analyse beperkten tot subkritische regimes waar klassieke energie-gelijkheden gelden, ofwel vertrouwden op uniciteitsresultaten die niet beschikbaar zijn in het superkritische geval voor schijn-verdunnende (shear-thinning) vloeistoffen. De introductie van het concept van de $b$-zwakke oplossing wordt gepresenteerd als een noodzakelijke technische vooruitgang om het gebrek aan energie-gelijkheid en de uitdagingen van de randvoorwaarden te overwinnen. De stabiliteitsresultaat wordt beschreven als een belangrijke stap voorwaarts, die aantoont dat ondanks de potentiële niet-uniciteit van zwakke oplossingen in dit regime, de fysiek relevante klasse van $b$-zwakke oplossingen een stabiel lange-termijn gedrag vertoont dat convergeert naar het unieke evenwicht. De analyse vermijdt expliciet de Oberbeck–Boussinesq-benadering, waardoor de volledige koppeling wordt geadresseerd waarbij kinetische energie-dissipatie fungeert als een bron in de warmtevergelijking, een kenmerk dat aanzienlijke wiskundige moeilijkheden introduceert vergeleken met vereenvoudigde modellen.