Virtual element approximation of eigenvalue problems: is the stabilization of the right hand side necessary?

Titel: Virtuele Elementbenadering van Eigenwaardeproblemen: Is Stabilisatie van de Rechterkant Noodzakelijk?

Auteurs: Daniele Boffi, Francesca Gardini, en Lucia Gastaldi.

---

1. Het Probleem

Het artikel behandelt de numerieke benadering van eigenwaardeproblemen die geassocieerd zijn met partiële differentiaalvergelijkingen (PDE's), specifiek voor elliptische zelfgeadjungeerde operatoren, met behulp van de Virtuele Elementmethode (VEM).

In een standaard variatieformulering zoekt men naar eigenwaarden $\lambda$ en eigenfuncties $u$ zodanig dat:
$$a(u, v) = \lambda b(u, v) \quad \forall v \in V$$
Waarbij $a(\cdot, \cdot)$ de bilineaire vorm van de operator (bijv. de Laplace-operator) is en $b(\cdot, \cdot)$ de bilineaire vorm die overeenkomt met het inproduct (meestal de $L^2$-inproduct).

Bij de discretisatie met VEM worden de bilineaire vormen $a_h$ en $b_h$ benaderd via projecties op polynoomruimtes. Een bekend probleem bij VEM is dat deze projecties alleen consistentie garanderen, maar niet altijd stabiliteit. Om dit op te lossen, wordt doorgaans een stabilisatieterm toegevoegd aan zowel de linkerzijde (matrix $A$, geassocieerd met $a$) als de rechterzijde (matrix $B$, geassocieerd met $b$).

De kernvraag: Is de stabilisatie van de matrix $B$ (de rechterkant) echt noodzakelijk voor de convergentie van de eigenwaarden, of kan men deze weglaten zonder de kwaliteit van de oplossing te schaden? Dit is kritiek omdat een onjuiste stabilisatie van $B$ kan leiden tot "spurious eigenvalues" (kunstmatige eigenwaarden) die het spectrum vervuilen, vooral in de buurt van de fundamentele modi.

2. Methodologie

De auteurs analyseren het probleem binnen een abstract raamwerk voor niet-compacte operatoren, gebaseerd op de theorie van Descloux–Nassif–Rappaz.

• Discretisatie: Ze gebruiken standaard VEM-ruimtes $V_h^k$ met polynoomgraad $k$.
• Stabilisatie:
  • De matrix $A$ (linkerkant) wordt wel gestabiliseerd (via een term $S_E$) om coerciviteit te garanderen.
  • De matrix $B$ (rechterkant) wordt niet gestabiliseerd. De discrete vorm $b_h$ wordt uitsluitend gedefinieerd via de $L^2$-projectie $\Pi^0_k$.
• Kern van $B$: Omdat $b_h$ niet gestabiliseerd is, kan de matrix $B$ een niet-triviale kern hebben ($K_{b,h} \neq \{0\}$). Dit betekent dat er functies bestaan waarvoor $b_h(w_h, v_h) = 0$. In het algebraïsche probleem $Ax = \lambda Bx$ zouden dit oneindige eigenwaarden zijn.
• Aanpak: De auteurs beperken het probleem tot de orthogonale complementaire ruimte $K_{b,h}^\perp$ (ten opzichte van $a_h$). Ze definiëren een discrete oplossingoperator $T_h$ op deze ruimte.
• Convergentiecriteria: Om convergentie te bewijzen, moeten twee eigenschappen worden aangetoond:
  1. P1: De operator $T_h$ convergeert naar de continue operator $T$ (in de operator-norm).
  2. P2: De discrete ruimte kan de continue oplossing goed benaderen (approximatie-eigenschap).

3. Belangrijkste Bijdragen en Theoretische Resultaten

Het artikel levert een theoretisch bewijs dat stabilisatie van de rechterkant niet noodzakelijk is voor lage-orde VEM-ruimtes:

• Theorema voor $k=1$ en $k=2$: De auteurs bewijzen dat voor polynoomgraden $k=1$ en $k=2$, de eigenschap P2 geldt, zelfs als de matrix $B$ een niet-triviale kern heeft.
  • Redenering: Voor $k=1$ en $k=2$ geldt een specifieke eigenschap in de VEM-definitie: de gradiëntprojectie $\Pi^\nabla_k$ en de $L^2$-projectie $\Pi^0_k$ gedragen zich zodanig dat de fouten in de benadering van de kern van $B$ voldoende snel verdwijnen naarmate de meshgrootte $h$ afneemt.
  • Dit impliceert dat de discrete eigenwaarden convergeren naar de ware eigenwaarden, zelfs zonder stabilisatieterm in $b_h$.
• Afwezigheid van Spurious Modes: Het bewijs garandeert dat er geen kunstmatige eigenwaarden in het spectrum verschijnen die de resultaten vervuilen, zolang men zich beperkt tot de relevante deelruimte.
• Hoogere Orde: Hoewel het theoretische bewijs specifiek is voor $k=1, 2$, suggereren de numerieke resultaten dat het ook werkt voor hogere graden ($k=3, 4$), hoewel dit theoretisch nog niet volledig is afgedekt.

4. Numerieke Resultaten

De auteurs voeren uitgebreide numerieke tests uit op het Laplace-eigenwaardeprobleem op het eenheidsvierkant, gebruikmakend van vier verschillende mesh-types:

1. T: Driehoeken (Triangular)
2. S: Vierkanten (Square)
3. V: Voronoi-polygoonmeshes
4. H: Hexagonale meshes
5. D: Dyadische meshes (een minder standaard type)

Vergelijking van resultaten:

• Kern-dimensie: Voor $k=1$ en $k=2$ is de dimensie van de kern van $B$ ($K_{b,h}$) overal nul of zeer klein (afhankelijk van de mesh), wat betekent dat $B$ vaak positief definiet is of slechts een zeer kleine singulariteit heeft. Voor $k \ge 3$ is de kern vaak niet-triviaal.
• Convergentie: Ondanks de afwezigheid van stabilisatie in $B$ en het bestaan van een niet-triviale kern voor hogere graden, tonen de tabellen aan dat:
  • De berekende eigenwaarden optimaal convergeren naar de exacte waarden.
  • De convergentieorden komen overeen met de theoretische verwachtingen ($O(h^{k+1})$ of $O(h^k)$ afhankelijk van de norm).
  • Er geen spurious eigenvalues worden waargenomen in het spectrum.
• Mesh-onafhankelijkheid: De resultaten zijn robuust voor alle geteste mesh-types (inclusief onregelmatige Voronoi- en dyadische meshes).

5. Betekenis en Conclusie

De belangrijkste conclusie van dit werk is dat voor standaard VEM-ruimtes met lage orde ($k=1, 2$) de stabilisatie van de rechterkant (matrix $B$) overbodig is voor het oplossen van elliptische eigenwaardeproblemen.

Significantie:

1. Vereenvoudiging: Het elimineert de noodzaak om een extra stabilisatieparameter te kiezen voor de massamatrix ($B$), wat een lastige en gevoelige taak is. Foutieve keuze van deze parameter kan leiden tot vervuiling van het spectrum.
2. Robuustheid: Het bewijst dat VEM-methoden robuust zijn tegenover het negeren van stabilisatie in de massamatrix, zolang men de correcte functionele ruimte gebruikt.
3. Praktische toepassing: Numerieke tests tonen aan dat zelfs voor hogere orden ($k > 2$), waar de theorie nog niet volledig geldt, de methode in de praktijk uitstekend werkt zonder stabilisatie van $B$. Dit opent de deur voor efficiëntere en eenvoudigere implementaties van VEM voor eigenwaardeproblemen.

Kortom, de paper weerlegt de noodzaak van "over-stabilisatie" en biedt een theoretische onderbouwing voor het gebruik van een "stabilisatie-vrije" rechterkant in VEM voor eigenwaardeproblemen.