Stabilization-Free H(curl) and H(div)-Conforming Virtual Element Method

Technische Samenvatting: Stabilisatievrije $H(\text{curl})$- en $H(\text{div})$-conforme Virtuele Element Methode

Probleemstelling
De Virtuele Element Methode (VEM) biedt aanzienlijke voordelen ten opzien van de Eindige Element Methode (FEM) door het hanteren van algemene ster-vormige polytopen, wat het ideaal maakt voor complexe geometrieën en adaptieve roosterverfijning. Echter, standaard VEM-implementaties voor $H(\text{curl})$- en $H(\text{div})$-ruimten vereisen stabilisatietermen om coerciviteit te waarborgen. Deze noodzaak komt voort uit het feit dat de virtuele elementruimte $V(K)$ groter is dan de polynoomruimte $P_k(K)$, waardoor de standaard $L^2$-projectie instabiel is (d.w.z. $\|\Pi_k v\|_K \not\gtrsim \|v\|_K$). Deze stabilisatietermen brengen verschillende nadelen met zich mee: ze hebben een negatieve impact op de numerieke stabiliteit en conditiegetallen, vereisen probleem-specifieke ontwerpen (met name voor niet-lineaire problemen) en vertonen een isotroop gedrag dat moeite heeft met anisotrope coëfficiënten of roosters. Bestaande stabilisatievrije benaderingen vereisen ofwel excessieve interne vrijheidsgraden (DoFs) die de computationele kosten verhogen, of wijzigen de projectieruimte op manieren die de geometrische flexibiliteit van VEM verwateren of essentiële eigenschappen zoals commutativiteit met differentiaaloperatoren niet behouden.

Methodologie
De auteurs stellen een nieuwe, stabilisatievrije VEM voor voor algemene orde $H(\text{curl})$- en $H(\text{div})$-conforme ruimten in $\mathbb{R}^3$. De kern van de methodologie bestaat uit het construeren van nieuwe serendipity-projectoren en bijbehorende serendipity-ruimten met een geminimaliseerd aantal DoFs.

1. Serendipity-projectoren: In tegenstelling tot eerdere werken die vertrouwen op triviale projectoren of interne DoFs, definiëren de auteurs specifieke projectoren ($\Pi^S_k$) voor de vlak-, rand- en knoopruimten. Deze projectoren mappen virtuele functies naar polynoomruimten met behulp van enkel "geconserveerde" rand-DoFs en een specifieke set interne momenten. De projectoren worden gedefinieerd via orthogonaliteitscondities waarbij de normale/tangentiële sporen en interne momenten betrokken zijn (bijv. $(v - \Pi^S_k v, x_K \times p_{kr})_K = 0$).
2. Serendipity-ruimten: De auteurs definiëren serendipity-ruimten ($V^S$) door de standaard virtuele elementruimten te beperken, zodanig dat de "gereduceerde" DoFs worden bepaald door de geconserveerde DoFs via de serendipity-projector. Dit zorgt ervoor dat de $L^2$-projectie op een voldoende hoog-orde polynoomruimte $P_l$ stabiel en berekenbaar is met enkel de rand-DoFs.
3. De Rham-complex: De methode is ontworpen om het volledige De Rham-complex keten in $\mathbb{R}^3$ te behandelen ($H^1 \to H(\text{curl}) \to H(\text{div}) \to L^2$). De constructie waarborgt exacte commutativiteit met differentiaaloperatoren ($\nabla, \nabla\times, \nabla\cdot$), waardoor de topologische structuur van het probleem behouden blijft.
4. Stabiliteitsanalyse: De auteurs bewijzen dat voor een voldoende grote projectieorde $l$, de normequivalentie $\|\Pi_l v\|_K \approx \|v\|_K$ geldt voor functies in de serendipity-ruimte. Dit elimineert de noodzaak voor stabilisatietermen terwijl optimale benaderingseigenschappen behouden blijven.

Belangrijkste Bijdragen

• Volledige De Rham-behandeling: Het werk biedt een verenigde behandeling voor het volledige De Rham-complex in 3D, wat zorgt voor naadloze integratie tussen verschillende differentiaaloperatoren, een kenmerk dat niet aanwezig is in eerdere stabilisatievrije werken die beperkt zijn tot $H^1$.
• Behoud van Wiskundige Eigenschappen: De voorgestelde serendipity-ruimten handhaven rigoureus de randcontinuïteit, operatorconformiteit, exacte commutativiteit met differentiaaloperatoren en B-compatibiliteit.
• Computationele Efficiëntie: Door het aantal DoFs te minimaliseren en gebruik te maken van serendipity-projecties die uitsluitend uit randgegevens worden berekend, reduceert de methode de computationele overhead aanzienlijk. Dit is bijzonder gunstig voor eigenwaardeproblemen waarbij interne DoFs niet geëlimineerd kunnen worden via statische condensatie.
• Rigoureuze Foutanalyse: Het artikel biedt een theoretisch bewijs van optimale benaderingseigenschappen en stabiliteit voor de hoog-orde polynoomprojecties.

Resultaten
De theoretische claims worden gevalideerd door numerieke experimenten op Maxwell eigenwaardeproblemen met Voronoi-polyedrische roosters.

• Convergentiecijfers: De numerieke resultaten bevestigen de theoretische convergentiecijfers. Voor oneindig regelmatige oplossingen bereikt de methode een $O(h^{m+1})$ convergentie in de $L^2$-norm en $O(h^{2m+2})$ voor eigenwaarden. Voor singuliere oplossingen (waar de regulariteit beperkt wordt door $s$), zijn de snelheden respectievelijk $O(h^s)$ en $O(h^{2s})$.
• Computationele Kosten: De assemblage-tijd van de voorgestelde Stabilisatievrije VEM (SFVEM) is slechts marginaal groter dan die van conventionele FEM, wat de efficiëntie van de methode bevestigt. De vereiste projectieorde $l$ blijkt in de praktijk dicht bij de polynoomorde $k$ te liggen.
• Vergelijking: De methode slaagt erin de complicaties gerelateerd aan stabilisatie te elimineren terwijl de nauwkeurigheid van conventionele FEM en gestabiliseerde VEM wordt geëvenaard.

Betekenis en Claims
De auteurs beweren dat dit werk een fundamenteel andere benadering biedt voor stabilisatievrije VEM die de beperkingen van bestaande strategieën (excessieve DoFs of verlies van geometrische flexibiliteit) overwint. Door stabilisatietermen te elimineren, verbetert de methode de numerieke stabiliteit en de conditiegetallen, wat het bijzonder waardevol maakt voor:

• Mixed formuleringen.
• Niet-lineaire problemen.
• Toepassingen die robuust numeriek gedrag vereisen over diverse roosterconfiguraties.

Het artikel concludeert dat hoewel de techniek theoretisch kan worden gegeneraliseerd naar hogere dimensies, de praktische implementatiecomplexiteit en de snelle groei van rand-DoFs de discussie momenteel beperken tot de 3D-setting. De voorgestelde methode vertegenwoordigt een theoretisch rigoureuze en computationeel efficiënte alternatief voor het oplossen van partiële differentiaalvergelijkingen, waarmee de toepasbaarheid van VEM in computationele wiskunde en engineering wordt uitgebreid.