A posteriori error bounds for finite element approximations of time-dependent mean field games

Technische Samenvatting: A Posteriori Foutgrenzen voor Eindige-Elementenbenaderingen van Tijdafhankelijke Mean Field Games

Probleemformulering
Dit werk behandelt de numerieke benadering van tweede-orde, tijdafhankelijke Mean Field Game (MFG) systemen. Het systeem bestaat uit een gekoppeld paar partiële differentiaalvergelijkingen: een achterwaartse Hamilton–Jacobi–Bellman (HJB) vergelijking voor de waardefunctie $u$ en een voorwaartse Kolmogorov–Fokker–Planck (KFP) vergelijking voor de dichtheid van spelers $m$. Het systeem is gedefinieerd op een begrensd domein $\Omega \subset \mathbb{R}^d$ over een tijdsinterval $(0, T)$, onderworpen aan initiële condities voor $m$ en terminale condities voor $u$. De vergelijkingen bevatten een constante diffusiecoëfficiënt $\nu > 0$, een Hamiltoniaan $H$ die $C^{1,1}$ is met betrekking tot de gradiëntvariabele, en koppelingsoperatoren $F$ en $S$ die Lipschitz-continu en monotoon zijn. Een cruciale vereiste voor de fysieke validiteit van de oplossing is het behoud van de niet-negativiteit van de dichtheid $m$ op het discrete niveau.

Methodologie
De auteurs analyseren een algemene klasse van gestabiliseerde eindige-elementenmethoden (FEM) voor dit systeem. De ruimtelijke discretisatie maakt gebruik van continue stuksgewijze affine eindige elementen, terwijl de temporele discretisatie gebruikmaakt van een impliciet Euler-schema. Om de niet-negativiteit van de dichtheid en de stabiliteit van het schema te waarborgen, worden abstracte stabilisatietermen geïntroduceerd. Deze stabilisaties kunnen massalumping voor tijdderivaten en ruimtelijke stabilisaties voor advectieve termen omvatten.

De kern van de methodologie betreft een rigoureuze a posteriori foutanalyse. De auteurs stellen eerst een continu-niveau equivalentie vast tussen de foutnorm en de duale norm van de residu. Vervolgens breiden zij deze analyse uit naar de discrete setting, wat specifieke uitdagingen met zich meebrengt vanwege het gebrek aan temporele conformiteit in de standaard discontinu Galerkin-achtige ruimtes die voor de dichtheid worden gebruikt. Om deze kloof te overbruggen, definiëren de auteurs uitgebreide normen op somruimtes ($V_i + Y$) en maken zij gebruik van temporele reconstructie-operatoren ($I_i$) om de discrete oplossingen naar de continue oplossingsruimte $Y$ te mappen.

De foutschatter $\eta$ wordt geconstrueerd uit verschillende componenten:

1. Residu-schatters: Gebaseerd op de volumeresidu's en de ruimtelijke sprongresidu's van de HJB- en KFP-vergelijkingen.
2. Temporele sprong-schatters: Die het gebrek aan conformiteit over de tijdstappen heen meten.
3. Stabilisatie-schatters: Die de fout kwantificeren die wordt geïntroduceerd door de stabilisatietermen.
4. Rand-schatters: Die rekening houden met fouten in de initiële en finale tijdcondities.
5. Data-oscillatietermen: Die de benaderingsfout van de probleemdata representeren.

Belangrijkste Bijdragen en Resultaten
Het artikel presenteert drie primaire theoretische resultaten:

1. Equivalentie van Fout- en Residunormen: De auteurs bewijzen dat voor het continue probleem de fout in de oplossing (gemeten in passende Bochner–Sobolev normen) equivalent is aan de duale norm van het residu. Cruciaal is dat deze equivalentie globaal geldt voor elk paar $(u, m)$ met een niet-negatieve dichtheid, en niet slechts in een buurt van de exacte oplossing. Dit breidt eerdere resultaten van stationaire problemen uit naar het tijdafhankelijke geval.

2. Betrouwbaarheid en Efficiëntie van Algemene Schatters: Voor de algemene klasse van gestabiliseerde FEM leiden de auteurs foutgrenzen af die laten zien dat de totale fout van bovenaf begrensd wordt (betrouwbaarheid) en van onderaf (efficiëntie) door de voorgestelde schatter $\eta$, tot aan de data-oscillatietermen.

   • Betrouwbaarheid: De fout wordt begrensd door de schatter plus de temporele data-oscillatie.
   • Efficiëntie: De schatter wordt begrensd door de fout plus de ruimtelijke data-oscillatie.
   • Een belangrijke technische bijdrage is de behandeling van de niet-conforme aard van de discrete dichtheid in de tijd, wat wordt bereikt door de reconstructie-operatoren en uitgebreide normen.

3. Stabilisatie-vrije Grenzen voor Gestructureerde Schema's: Het artikel identificeert een specifieke subklasse van stabilisaties (betreffende massalumping en patchwise affine-preserverende ruimtelijke stabilisaties) waarbij de stabilisatie-schatters begrensd kunnen worden door de residu- en de temporele sprong-schatters. Onder de voorwaarde dat de tijdstapgrootte $\tau$ niet excessief klein is ten opzichte van het kwadraat van de roostergrootte ($h^2 \lesssim \tau$), tonen de auteurs aan dat de stabilisatie-schatters niet expliciet berekend hoeven te worden. Dit levert een lokaal berekenbare en lokaal efficiënte schatter op die uitsluitend steunt op standaard residu- en sprongtermen.

Significantie
Het artikel beweert de eerste a posteriori foutgrenzen te bieden voor een algemene klasse van gestabiliseerde eindige-elementenbenaderingen van tijdafhankelijke MFG-systemen. Hoewel eerdere werken (bijv. [25]) dergelijke grenzen vaststelden voor stationaire MFG's, breidt dit werk de theorie uit naar de tijdafhankelijke setting, waarbij de extra complexiteit van temporele niet-conformiteit en de koppeling van voorwaartse-achterwaartse paraboolvergelijkingen wordt geadresseerd.

De significantie ligt in het mogelijk maken van adaptieve roosterverfijning voor MFG's. Door betrouwbare en efficiënte foutschatters te bieden, stelt dit werk de ontwikkeling van algoritmen mogelijk die het rooster en de tijdstappen lokaal kunnen verfijnen op basis van de berekende fout, wat potentieel aanzienlijke computationele besparingen oplevert vergeleken met uniforme verfijning. De afleiding van stabilisatie-vrije grenzen voor specifiek relevante stabilisatieschema's versterkt het nut van de resultaten verder, omdat het de noodzaak wegneemt om complexe, stabilisatie-specifieke foutindicatoren in de praktijk te berekenen.