Convergence analysis for a finite-volume scheme for the Euler- and Navier-Stokes-Korteweg system via energy-variational solutions

1. Probleemstelling en Context

Het artikel richt zich op de numerieke simulatie van de beweging van samendrukbare vloeistoffen met capillaire effecten (diffuse interfaces), beschreven door de Euler-Korteweg (EK) en Navier-Stokes-Korteweg (NSK) vergelijkingen.

• De Uitdaging: Klassieke oplossingen voor deze hyperbolische en parabolisch-hyperbolische systemen bestaan vaak slechts lokaal in tijd. Bij zwakke oplossingen (weak solutions) treedt het probleem van niet-uniekheid op, wat de voorspellende kracht van de vergelijkingen in gevaar brengt.
• Het Doel: Het bewijzen van de convergentie van een specifiek, structuurbehoudend finite-volume schema naar een goed gedefinieerd concept van een zwakke oplossing, zonder de aanwezigheid van een sterke (klassieke) oplossing vooraf te hoeven aannemen.

2. Methodologie

A. Het Numerieke Schema

De auteurs analyseren een semidiscreet finite-volume schema (geïntroduceerd in eerdere werk [25]) voor het tweedimensionale torus-domein ($\mathbb{T}^2$).

• Discretisatie: Het domein wordt opgedeeld in een Cartesiaans rooster. De variabelen dichtheid ($\rho$) en impuls ($m = \rho u$) worden benaderd door stuksgewijs constante functies per cel.
• Structuurbehoud: Het schema is ontworpen om fundamentele fysische eigenschappen te behouden:
  • Behoud van totale massa en impuls.
  • Dissipatie van de totale energie (voor NSK) of behoud van energie (voor EK, afhankelijk van viscositeit).
• Korteweg-spanning: De derde-orde termen (capillariteit) worden geïntroduceerd via een specifieke discretisatie van de Korteweg-spanningstensor, wat essentieel is voor de stabiliteit.
• Numerieke Dissipatie: Een parameter $\lambda$ wordt toegevoegd voor stabiliteit, gekozen als een functie van de lokale snelheid en geluidssnelheid (Rusanov-type dissipatie).

B. Het Oplossingsconcept: Energie-Variationale Oplossingen

In plaats van te werken met maat-waarde oplossingen (measure-valued solutions) of Young-maten, maken de auteurs gebruik van het relatief nieuwe concept van energie-variational solutions.

• Definitie: Een oplossing wordt gedefinieerd als een triplet $(\rho, m, E)$, waarbij $E$ een variabele is die fungeert als een bovengrens voor de fysieke energie.
• Variationale Ongelijkheid: De oplossing moet voldoen aan een specifieke variationale ongelijkheid die de zwakke vorm van de vergelijkingen combineert met een energie-dissipatie-inegaliteit.
• Voordeel: Dit concept maakt het mogelijk om te werken in standaard energie-ruimtes en vermijdt de complexiteit van defect-maten (defect measures) die nodig zijn bij andere generalisaties. De formulering is stabiel onder zwakke convergentie.

C. Convergentieanalyse

Het bewijs van convergentie verloopt in drie hoofdstappen:

1. A-priori Schattingen: Afgeleid uit de discrete energie-dissipatie-inegaliteit. Dit levert uniforme bounds op voor de dichtheid, impuls, snelheid en de discrete energie in geschikte Banach-ruimtes (o.a. $L^\infty(L^\gamma)$, $L^2(H^1)$).
2. Sterke Convergentie: Gebruikmakend van de Aubin-Lions lemma en discrete Rellich-Kondrachov theorema's, wordt bewezen dat er een sterk convergente deelrij bestaat voor de dichtheid $\rho_h$. Dit is cruciaal om de niet-lineaire termen (zoals $p(\rho)$ en convectie) correct te kunnen behandelen.
3. Limietovergang: Door de limiet $h \to 0$ te nemen in de discrete formulering, wordt aangetoond dat de limiet voldoet aan de definitie van een energie-variational oplossing. Hierbij wordt gebruik gemaakt van de lage halfcontinuïteit van de energie en de specifieke keuze van de regulariteitsgewicht $K(\phi)$.

3. Belangrijkste Resultaten

• Convergentiestelling (Theorem 7): De auteurs bewijzen dat de numerieke oplossingen $(\rho_h, u_h)$ van het finite-volume schema, onder mesh-verfijning ($h \to 0$), convergeren (in een geschikte zwakke en sterke topologie) naar een energie-variational oplossing van het EK- of NSK-systeem.
• Uniekheid van het Concept: De limiet is een geldige oplossing volgens het nieuwe raamwerk, zelfs in afwezigheid van een sterke oplossing.
• Sterke Convergentie van Snelheid: Voor het NSK-systeem (met viscositeit $\mu, \eta > 0$) en onder de aanname dat de dichtheid uniform van nul gescheiden blijft, wordt ook de sterke convergentie van de snelheid $u_h$ naar $u$ bewezen.
• Initiële Waarden: De constructie garandeert dat de energie van de numerieke oplossing op $t=0$ overeenkomt met de initiële energie van het continue systeem.

4. Significantie en Bijdrage

• Eerste Bewijs zonder Sterke Oplossing: Dit is, naar de kennis van de auteurs, het eerste artikel dat de convergentie van een numeriek schema voor EK- en NSK-vergelijkingen bewijst zonder de aanwezigheid van een sterke oplossing als voorwaarde te stellen. Dit is een grote stap voorwaarts in de numerieke analyse van niet-lineaire hyperbolische systemen.
• Vermijden van Young-maten: Door het gebruik van energie-variational oplossingen vermijden de auteurs de introductie van complexe Young-maten of defect-maten, wat de analyse robuuster en directer toepasbaar maakt voor numerieke methoden.
• Structuurbehoud: Het werk benadrukt het belang van structuurbehoudende schema's (die massa, impuls en energie respecteren) als de sleutel tot het aantonen van convergentie voor complexe fysische systemen.
• Toepasbaarheid: De resultaten zijn relevant voor de ontwikkeling van betrouwbare simulaties van vloeistoffen met capillaire effecten, zoals in fase-overgangen, vloeistof-damp grensvlakken en andere complexe stromingsproblemen.

Conclusie

Het artikel levert een fundamentele theoretische onderbouwing voor het gebruik van een specifiek finite-volume schema voor de Euler- en Navier-Stokes-Korteweg systemen. Door de brug te slaan tussen numerieke discretisatie en het concept van energie-variational oplossingen, bieden de auteurs een rigoureuze garantie dat de numerieke resultaten fysisch betekenisvolle oplossingen benaderen, zelfs in regimes waar klassieke theorie faalt.