$hp$-a posteriori error estimates for hybrid high-order methods applied to biharmonic problems

Probleemstelling
Dit artikel behandelt de a posteriori foutschatting voor Hybrid High-Order (HHO) methoden toegepast op het biharmonische probleem ($\Delta^2 u = f$) met ingeklemde randvoorwaarden ($u = \partial_n u = 0$ op $\partial\Omega$). Het probleem is geformuleerd op twee- en driedimensionale polytopale Lipschitz-domeinen met een algemene topologie (mogelijk met gaten). Het doel is om residu-gebaseerde $hp$-foutschatters af te leiden die robuust zijn ten opzichte van de maashoogte ($h$), de polynoomgraad ($p$ of $k$) en de domeintopologie.

Methodologie
De auteurs ontwikkelen een raamwerk om de totale fout te begrenzen door deze te ontbinden in conforme en niet-conforme componenten.

1. Foutontbinding: De fout wordt gesplitst in een conform deel (gemeten in de energiemorm van een gereconstrueerde $H^2_0$-functie) en een niet-conform deel (de afstand tussen de discrete oplossing en de conforme ruimte).
2. Niet-conforme foutcontrole:
   • Om de niet-conforme fout te begrenzen, construeren de auteurs een conforme $H^2_0$-reconstructie van de discrete oplossing met behulp van een $C^1$-partitie van eenheid.
   • Deze partitie wordt opgebouwd met behulp van $C^1$-conforme eindelelementruimten (specifiek de constructie uit [40] gebaseerd op Alfeld-splitsingen van het rooster).
   • Cruciaal is dat de auteurs lokale Helmholtz-ontbindingen gebruiken op vertex-patches (sterren). In tegen tegenstelling tot eerdere benaderingen die vertrouwden op globale Helmholtz-ontbindingen (die constanten introduceren die afhankelijk zijn van het aantal gaten in het domein), zorgt deze lokale benadering ervoor dat de stabiliteitsconstante alleen afhankelijk is van de vorm-regulariteit van het rooster.
3. Conforme foutcontrole:
   • Twee verschillende bovengrenzen voor de conforme fout worden afgeleid, elk gebruikmakend van een specifieke interpolatieoperator:
     • Setting 1: Gebruikt de canonieke hybride eindelemen-interpolatieoperator ($C^{k+2}_T$). Deze benadering steunt op de eigenschap dat de compositie van deze operator met de HHO-reconstructie een $H^2$-elliptische projectie oplevert. De resulterende bovengrens bevat alleen de stabilisatieterm en data-oscillatie, maar is afhankelijk van een numerieke aanname over de interpolatiefout (Aanname 3.4).
     • Setting 2: Gebruikt de gemodificeerde Babuška–Suri $hp$-interpolatieoperator ($I^{k+2}_{mBS}$). Deze benadering vereist geen aanvullende aannames, maar introduceert extra foutindicatortermen, specificaal bulk-residuen en normale/tangentiële flux-sprongen.
4. Residu-schatters: De tekst definieert lokale foutindicatoren bestaande uit stabilisatietermen, bulk-residuen, normale flux-sprongen en tangentiële sprongen, geschaald door passende $hp$-factoren.

Belangrijkste Bijdragen

• Eerste $hp$-schatter voor HHO op biharmonische problemen: Dit werk levert de eerste residu-gebaseerde $hp$-$a$ posteriori foutschatter voor HHO-methoden toegepast op biharmonische problemen op algemene 2D en 3D domeinen.
• Lokale Helmholtz-ontbinding met topologische robuustheid: Het artikel stelt een nieuwe stabiliteitsschatting vast voor lokale Helmholtz-ontbindingen op vertex-sterren. Het bewijst dat de stabiliteitsconstante uitsluitend afhangt van de vorm-regulariteit van het rooster, waardoor de degradatie van constanten geassocieerd met het aantal gaten in het domein wordt vermeden, een bekend probleem bij globale ontbindingsbenaderingen.
• Duale Bovengrenzen: De afleiding van twee verschillende bovengrenzen voor de conforme fout biedt flexibiliteit. De ene bovengrens is "schoon" (afhankelijk van enkel stabilisatie en oscillatie) maar voorwaardelijk op een numeriek geverifieerde aanname; de andere is onvoorwaardelijk maar bevat complexere residu-termen.
• Extensie naar Algemene Topologie: De methodologie is expliciet ontworpen om domeinen met een willekeurige topologie (meerdere gaten) te kunnen verwerken, een kenmerk dat vaak moeilijk aan te pakken is in $hp$-schattingen voor vierde-orde problemen.

Resultaten

• Theoretische Grenzen: Het hoofdbewijs (Theorem 5.6) stelt vast dat de totale fout (gebroken Hessiaan-norm plus stabilisatie) begrensd wordt door de som van de voorgestelde foutindicatoren. De grenzen zijn getoond als $h$-optimaal en $p$-suboptimaal (met maximaal één orde, of potentieel een halve orde onder de geverifieerde aanname).
• Numerieke Verificatie:
  • Verificatie van Aanname: Numerieke experimenten in 2D en 3D bevestigen de geldigheid van Aanname 3.4 met betrekking tot de interpolatiefout van de canonieke hybride operator.
  • Convergentiecijfers: Tests op gladde oplossingen in 3D laten zien dat zowel de energie-fout als de schatter convergeren met de optimale snelheid $O(\text{DoFs}^{-(k+1)/d})$.
  • Effectiviteitsindices: De schatters leveren matige effectiviteitsindices (variërend tussen ongeveer 1.2 en 2.0) over diverse polynoomgraden en roosterverfijningen.
  • Topologische Robuustheid: Experimenten op 2D-domeinen met een toenemend aantal gaten (0, 1, 4, 8) tonen aan dat de effectiviteitsindices stabiel blijven, wat de theoretische claim bevestigt dat de schatter robuust is ten opzichte van de domeintopologie.
  • Adaptieve Prestaties: $h$-adaptieve algoritmen, gestuurd door de schatter, slagen erin om optimale convergentiecijfers te herstellen voor singuliere oplossingen in zowel 2D (L-vormig domein) als 3D (punt-singulariteit), waarbij ze superieur zijn aan uniforme verfijning.

Betekenis
De auteurs positioneren dit werk als het eerste dat een $C^1$-partitie van eenheid combineert met lokale Helmholtz-ontbindingen op vertex-sterren om niet-conforme fouten te begrenzen in de $hp$-setting voor biharmonische problemen. De primaire betekenis ligt in het vermogen om de niet-conforme fout te controleren zonder dat de stabiliteitsconstanten verslechteren door de topologische complexiteit van het domein (aantal gaten). Bovendien merken de auteurs op dat hoewel de benadering is ontwikkeld voor HHO-methoden, de techniek voor het begrenzen van niet-conforme fouten via lokale Helmholtz-ontbindingen potentieel toepasbaar is op andere niet-conforme methoden, zoals Interior Penalty Discontinuous Galerkin (IPDG) methoden. Het artikel claimt geen lagere foutgrenzen te bieden, waarbij wordt opgemerkt dat dergelijke grenzen zouden volgen uit aanpassingen van bestaande argumenten, maar waarschijnlijk een pessimistische $p$-suboptimaliteit zouden vertonen vanwege de beperkingen van huidige $hp$-technieken voor het construeren van specifieke bubble-functies en extensie-operatoren.