On three-dimensional flows of thermo-viscoelastic fluids of Giesekus type

1. Probleemstelling

Het artikel behandelt de wiskundige analyse van de stroming van een oncompressibele, warmtegeleidende, visco-elastische vloeistof van het Giesekus-type in een begrensd driedimensionaal domein $\Omega \subset \mathbb{R}^3$ over een tijdsinterval $(0, T)$.

Het systeem koppelt de behoudswetten voor massa, lineaire impuls, inwendige energie en de evolutie van de elastische spanning (weergegeven door de linkse Cauchy–Green-tensor $B$ of de vervormingsgradient $F$). De leidende vergelijkingen zijn:

• Oncompressibiliteit: $\text{div } v = 0$.
• Impulsbalans: $\rho \partial_t v + \rho \text{div}(v \otimes v) - \text{div } T = \rho f$.
• Spanningsevolutie: Een evolutievergelijking voor de tensor $B$ (of $F$) die de bovenconvectieve Oldroyd-afgeleide en een dissipatieve term specifiek voor het Giesekus-model bevat.
• Balans van inwendige energie: $\rho \partial_t e + \rho \text{div}(ev) + \text{div } q - T : D = 0$.

Belangrijkste kenmerken van het model:

• Temperatuurafhankelijkheid: Alle materiaaleigenschappen (viscositeit $\nu$, relaxatietijd $\tau$, warmtegeleidingsvermogen $\kappa$ en de functie voor de elastische modulus $g$) zijn afhankelijk van de temperatuur $\theta$.
• Giesekus-dissipatie: Het model bevat een niet-lineaire dissipatieve term in de spanningsevolutie, waardoor het zich onderscheidt van het Oldroyd-B-model.
• Randvoorwaarden: Navier-glijrandvoorwaarden voor de snelheid ($v \cdot n = 0, (Tn)_\tau = -\gamma^* v_\tau$) en nul warmtestroom ($\nabla \theta \cdot n = 0$).
• Initiële data: Er wordt uitgegaan van grote data, wat betekent dat er geen aannames over kleinheid worden opgelegd aan de initiële snelheid, temperatuur of vervormingsgradient.

2. Methodologie

De auteurs hanteren een verfijnd zwak-sterk raamwerk gecombineerd met geavanceerde compactheidstools om de sterke niet-lineariteiten en het gebrek aan regulariteit in de oplossingen te behandelen.

A. Thermodynamische consistentie en herformulering

• Entropie en energie: De auteurs leiden het systeem af vanuit een thermodynamisch consistent raamwerk met behulp van de Helmholtz-vrije energie. Ze vervangen de vergelijking voor inwendige energie door een entropie-ongelijkheid en een totaal-energie-ongelijkheid. Deze aanpak is cruciaal omdat de druk $p$ in de zwakke zin mogelijk niet goed gedefinieerd is voor bepaalde randvoorwaarden, terwijl de energie-ongelijkheid wel geldig blijft.
• Variabele transformatie: Om de evolutievergelijking voor $B$ te behandelen, die termen bevat die niet onder afgeleiden staan (wat de identificatie van limieten bemoeilijkt), introduceren de auteurs een variabele $F$ zodanig dat $B = FF^T$. Ze herformuleren de spanningsevolutievergelijking in termen van $F$ (een vergelijking van het transport-diffusie-type), following de strategie van Masmoudi. Dit zorgt voor betere controle over de niet-lineaire termen.

B. Benaderingsschema

Het bestaanbewijs rust op een meerstaps benaderingsprocedure die het volgende omvat:

1. Galerkin-benadering: Projectie van de snelheid, vervormingsgradient en inwendige energie op eindig-dimensionale deelruimten.
2. Regularisatie: Invoering van kunstmatige diffusietermen ($\varepsilon_4 \Delta F$, $\varepsilon_7 \Delta e$) en afsnijfuncties ($\Lambda_{\varepsilon_3}$, $g_{\varepsilon_1}$) om het bestaan van gladde benaderende oplossingen te garanderen.
3. Iteratieve limieten: De auteurs laten de limiet toe in een specifieke volgorde van parameters ($m, n, l \to \infty$ voor Galerkin-indexen, gevolgd door $\varepsilon_7, \varepsilon_4, \varepsilon_2, \varepsilon_1, \varepsilon_5, \varepsilon_6, \varepsilon_3 \to 0$).

C. Compactheid en convergentietools

De kern van de analyse ligt in het aantonen van de sterke convergentie van de temperatuur $\theta$ en de vervormingsgradient $F$, wat noodzakelijk is om de limiet te nemen in de niet-lineaire constitutieve relaties.

• Convergentie van temperatuur: De auteurs gebruiken het Div-Curl-lemma (Murat-Tartar) in combinatie met de entropievergelijking. Door specifieke testfuncties te construeren die de entropie $\eta$ en de temperatuur betrekken, bewijzen ze dat een deelrij van temperaturen bijna overal convergeert.
• Convergentie van vervormingsgradient: Ze passen de renormalisatietechniek (DiPerna-Lions) toe en een specifieke Gronwall-type ongelijkheid voor de grootheid $|F|^2 - |F|^2$ (waarbij $F$ de limiet is). Dit houdt in dat het snelheidsveld wordt opgesplitst in een "Stokes-achtig" deel en een restterm om de convectieve termen te behandelen, wat uiteindelijk de sterke convergentie van $F$ in $L^2$ bewijst.
• Bijtende limieten: Chacon's bijtende lemma wordt gebruikt om zwakke convergentie in $L^1$-ruimtes te behandelen voor grootheden zoals de dissipatietermen.

3. Belangrijkste bijdragen

1. Globaal-in-tijd bestaan voor grote data: Het artikel vestigt het bestaan van een globale zwakke oplossing voor het volledige 3D Giesekus-model met temperatuurafhankelijkheid, zonder dat er behoefte is aan:
   • Kleinheid van initiële data.
   • Kunstmatige spanningsdiffusie (een gebruikelijke regularisatie in eerdere literatuur).
   • Onfysische groeiveronderstellingen voor materiaaleigenschappen (bijvoorbeeld warmtegeleidingsvermogen mag begrensd zijn, in tegenstelling tot eerdere werken die polynoomgroei vereisten).
2. Thermodynamisch consistente zwakke oplossing: De oplossing is gedefinieerd in een "dissipatief zwak" zin, waarbij de totale energie-ongelijkheid en de entropie-ongelijkheid worden voldaan. Dit zorgt ervoor dat de oplossing voldoet aan de tweede wet van de thermodynamica.
3. Behandeling van de Giesekus-niet-lineariteit: De auteurs behandelen succesvol de specifieke niet-lineaire dissipatieterm van het Giesekus-model ($\alpha(B-I)$ en $\tau(B^2-B)$) in aanwezigheid van temperatuurbekoppeling, een probleem dat open bleef voor grote data in 3D.
4. Algemene materiaaleigenschappen: Het resultaat geldt voor een brede klasse van continue materiaalfuncties (viscositeit, relaxatietijd, enz.) die slechts natuurlijke begrensdheids- en monotoniteitsvoorwaarden voldoen, waardoor het model toepasbaar is op een breed scala aan fysische scenario's.

4. Hoofdresultaten

Stelling: Onder de aanname dat het domein $\Omega$ Lipschitz is en de initiële data $(v_0, F_0, \theta_0)$ behoren tot de natuurlijke energie-entropie functieruimten (specifiek $v_0 \in L^2$, $F_0 \in L^2$, $\theta_0 \in L^1$ met $\det F_0 > 0$), bestaat er een globale-in-tijd dissipatieve zwakke oplossing $(v, \theta, F)$ voor het systeem.

De oplossing voldoet aan:

• Zwakke sequentiële stabiliteit: De oplossing is stabiel onder limieten van benaderende rijen.
• Energie-ongelijkheid: De totale energie (kinetisch + inwendig) voldoet aan een differentiaal-ongelijkheid die rekening houdt met randdissipatie.
• Entropie-ongelijkheid: De entropieproductie is niet-negatief, waardoor aan de tweede wet van de thermodynamica wordt voldaan.
• Regulariteit: De oplossing bezit specifieke regulariteits eigenschappen, waaronder $v \in L^\infty(0,T; L^2) \cap L^2(0,T; W^{1,2})$, $\theta \in L^\infty(0,T; L^1) \cap L^{5/3}(Q_T)$, en $F \in L^\infty(0,T; L^2) \cap L^4(Q_T)$.

5. Betekenis

Dit werk vertegenwoordigt een aanzienlijke doorbraak in de wiskundige theorie van complexe vloeistoffen. Voorheen berustten globale existentie-resultaten voor warmtegeleidende visco-elastische vloeistoffen vaak op:

• Vereenvoudigde modellen (bijvoorbeeld Oldroyd-B met spanningsdiffusie).
• Aannames over kleine data.
• Tweedimensionale settingen.
• Onfysische aannames over de groei van het warmtegeleidingsvermogen.

Door deze beperkingen weg te nemen, bieden Bulíček, Los en Woźnicki een rigoureuze basis voor het modelleren van realistische, grootschalige, driedimensionale stromingen van complexe materialen (zoals polymeren en biologische vloeistoffen) waarbij thermische effecten en elastische spanningen sterk gekoppeld zijn. De ontwikkelde methodologie, met name de combinatie van entropie-gebaseerde compactheid en de specifieke behandeling van de Giesekus-niet-lineariteit zonder kunstmatige diffusie, stelt een nieuwe standaard voor de analyse van visco-elastische vloeistoffen van het type rate-type.