Error estimates for completely discrete FEM in energy-type and weaker norms

Technische Samenvatting van "Error estimates for completely discrete FEM in energy-type and weaker norms"

Probleemstelling
Het artikel behandelt de numerieke analyse van eindige elementenmethoden (FEM) voor randwaardeproblemen met betrekking tot lineaire tweede-orde elliptische differentiaalvergelijkingen met een diffusie-convectie-reactiestructuur. De problemen zijn gedefinieerd op begrensde polyedrische Lipschitz-domeinen in $\mathbb{R}^d$ ($d \in \{2, 3\}$) en bevatten gemengde randvoorwaarden van de eerste, tweede en derde soort. Een centraal focuspunt ligt op "volledig discrete" formuleringen, waarbij de discrete eindige elementruimten niet noodzakelijkerwijs conform zijn (dat wil zeggen niet ingebed in de continue oplossingsruimte van de zwakke formulering) en de discrete problemen niet noodzakelijkerwijs aan standaard consistentie-eigenschappen voldoen. De auteurs streven ernaar een verenigd abstract kader te bieden voor het afleiden van foutschattingen in zowel energie-type normen als zwakkere normen (zoals $L^2$), zonder vooraf conformiteit of consistentie te veronderstellen.

Methodologie
De auteurs ontwikkelen een algemeen variationeel kader gebaseerd op twee reële Banachruimten, $U$ en $V$, en hun discrete tegenhangers $U_h$ en $V_h$. De methodologie verloopt via de volgende stappen:

1. Abstracte Foutanalyse:

   • Stabiliteit en Consistentie: Het artikel generaliseert de tweede lemma van Strang om niet-conforme en niet-consistente discretisaties te behandelen. Het definieert consistentie op een deelruimte $W \subset U$ door te eisen dat de discrete bilineaire vorm $a_h$ zodanig kan worden uitgebreid dat de continue oplossing $u$ de discrete vergelijking exact oplost (tot de discrete lineaire vorm $\ell_h$).
   • Energienorm Schattingen: Stelling 3 biedt een foutschatting in de discrete energienorm $\|u - u_h\|_{U_h}$. Deze schatting is afhankelijk van de benaderingsfout (beste benadering in de discrete ruimte), de stabiliteitsconstante (inf-sup conditie), en een consistentiefoutterm die de afwijking van de continue oplossing van de discrete formulering meet.
   • Schattingen in Zwakkere Normen (Gegeneraliseerde Aubin-Nitsche): Het artikel breidt het klassieke Aubin-Nitsche dualiteitsargument uit naar de niet-conforme/niet-consistente setting. Stelling 5 leidt een foutschatting af in een zwakkere norm $\|u - u_h\|_{Z'}$ (waarbij $Z$ een reflexieve Banachruimte is). Deze fout wordt ontbonden in vier termen die betrokken zijn bij:
     • Het product van de primaire en duale benaderingsfouten.
     • Consistentiefouten van de primaire en duale problemen.
     • Benaderingsfouten voortvloeiend uit het verschil tussen continue en discrete bilineaire/lineaire vormen.
   • De analyse steunt op de oplosbaarheid en regulariteit van een adjoint probleem en diens discrete tegenhanger.

2. Toepassing op Specifieke Discretisaties:
   De abstracte theorie wordt toegepast op twee concrete gestabiliseerde discretisaties voor de diffusie-convectie-reactievergelijking:

   • Voorbeeld I: Crouzeix-Raviart (CR1) Discretisatie: Een niet-conforme methode die gebruikmaakt van stuksgewijze lineaire polynomen die continu zijn op de barycentra van de elementvlakken, maar discontinu over de vlakken heen. De methode bevat upwind-stabilisatietermen voor het convectiegedeelte. De auteurs bewijzen coerciviteit en leiden optimale foutschattingen af in de $H^1$-achtige norm en de $L^2$-norm.
   • Voorbeeld II: Interior Penalty Discontinuous Galerkin (IPG) Methoden: Een familie van methoden (inclusief SIPG, NIPG en IIPG) die gebruikmaken van discontinue stuksgewijze polynomen van graad $k$. De formulering bevat penalty-termen voor sprongen in de oplossing en fluxen, evenals upwind-stabilisatie. De auteurs vestigen stabiliteit (coerciviteit) onder specifieke voorwaarden op de penalty-parameter en mesh-regulariteit, en bewijzen convergentiesnelheden.

Kernbijdragen en Resultaten

• Verenigd Kader: Het artikel presenteert een verenigde aanpak die niet vereist dat de discrete ruimten deelruimten zijn van de continue ruimten (niet-conform) noch dat de discrete formulering consistent is met de continue een (niet-consistent). Dit maakt de analyse van een bredere klasse van gestabiliseerde methoden mogelijk.
• Gegeneraliseerde Aubin-Nitsche Techniek: Een significante bijdrage is de uitbreiding van het dualiteitsargument naar volledig discrete settings. Stelling 5 biedt een rigoureuze grens voor fouten in zwakkere normen, waarbij expliciet rekening wordt gehouden met consistentiefouten in zowel het primaire als het adjoint probleem.
• Optimale Foutschattingen:
  • Voor de Crouzeix-Raviart methode bewijst het artikel optimale eerste-orde convergentie in de energienorm ($O(h)$) en tweede-orde convergentie in de $L^2$-norm ($O(h^2)$), uitgaande van voldoende regulariteit van de oplossing en coëfficiënten.
  • Voor IPG methoden worden optimale convergentiesnelheden vastgesteld. In de energienorm is de fout $O(h^m)$ voor polynoomgraad $k$ (met $m \le k$). In zwakkere normen toont het artikel aan dat voor de symmetrische IPG (SIPG, $\theta = -1$), het adjoint probleem consistent blijft, wat optimale foutschattingen in negatieve normen ($O(h^{k+m+1})$) mogelijk maakt. Voor niet-symmetrische varianten ($\theta \neq -1$) wordt aangetoond dat het adjoint probleem niet-consistent is, wat extra schatting van de consistentiefoutterm vereist.
• Stabiliteitsanalyse: De auteurs leveren gedetailleerde bewijzen voor de coerciviteit van de discrete bilineaire vormen, waarbij de noodzakelijke voorwaarden op de penalty-parameters en de mesh-geometrie (vorm- en contactregulariteit) worden geïdentificeerd om stabiliteit te garanderen.

Betekenis en Claims
De auteurs stellen dat hoewel bepaalde aspecten van de uiteenzetting (zoals hun versie van de tweede lemma van Strang en varianten van het dualiteitsargument) niet geheel nieuw zijn, de waarde van het werk ligt in het verenigde abstracte kader. Dit kader verheldert en kwantificeert de wisselwerking tussen stabiliteit, benaderingsfouten en consistentiefouten zonder te vertrouwen op de restrictieve aannames van conformiteit of consistentie.

De auteurs beweren dat deze aanpak:

1. De analyse van meer algemene eindige elementenmethoden mogelijk maakt, inclusief die die inherent niet-conform of niet-consistent zijn.
2. Een systematische manier biedt om optimale foutschattingen in zwakkere normen af te leiden door de Aubin-Nitsche techniek te generaliseren naar deze bredere settings.
3. Richtlijnen biedt voor het opzetten en analyseren van nieuwe stabiele en convergente schema's voor diffusie-convectie-reactieproblemen met gemengde randvoorwaarden.

De auteurs merken bescheiden op dat de theoretische resultaten worden toegepast op twee specifieke voorbeelden (CR1 en IPG) om het nut van het kader te demonstreren, maar dat de abstracte resultaten bedoeld zijn om toepasbaar te zijn op een breder scala aan modellen en numerieke technieken. Het werk bouwt voort op en breidt eerdere literatuur uit (bijv. [DPD18], [DPD21], [DPE12]) door het niet opnemen van interpolatie-operatoren in de definitie van consistentie en door gemengde randvoorwaarden en convectietermen algogener te behandelen.