Spectral coarse spaces based on indefinite operators: the $H_k$-GenEO method

Technische Samenvatting: Spectrale Coarse Ruimtes Gebaseerd op Indefiniete Operatoren: De Hk-GenEO Methode

Probleemstelling
Het artikel behandelt de numerieke oplossing van sterk indefiniete randwaardeproblemen voortvloeiend uit de heterogene Helmholtz-vergelijking:
$$-\nabla \cdot (A\nabla u) - k^2 n u = f \quad \text{in } \Omega, \quad u = 0 \quad \text{op } \partial\Omega,$$
waarbij $k > 0$ een grote golfgetal (wavenumber) is, $A$ een positief definiete matrix is met $L^\infty$-elementen, en $n$ een positieve $L^\infty$-functie is. De focus ligt op het regime waar $k$ groot is, wat het probleem sterk indefiniet en uitdagend maakt voor standaard iteratieve oplossers. Hoewel klassieke domeindecompositie-methoden (DDM) met spectrale coarse ruimtes (bijv. GenEO) succesvol zijn geweest voor positief semi-definiete problemen, falen ze vaak voor sterk indefiniete problemen wanneer ze worden geconstrueerd op basis van positief semi-definiete lokale eigenproblemen.

Methodologie
De auteurs stellen een nieuwe twee-niveau additive Schwarz-preconditioner voor, genaamd Hk-GenEO, die specifiek is ontworpen voor indefiniete operatoren. De methodologie bestaat uit de volgende kerncomponenten:

1. Lokale Indefiniete Eigenproblemen: In tegenstelling tot eerdere spectrale methoden die vertrouwen op positief semi-definiete lokale operatoren (zoals de Laplaciaan), construeert Hk-GenEO de coarse ruimte door lokale gegeneraliseerde eigenwaardeproblemen (GEVPs) op te lossen gebaseerd op de volledige indefiniete operator van het globale probleem. Op elk coarse subdomein $\Omega^c_i$ is de GEVP gedefinieerd als:
   $$-\nabla \cdot (A\nabla p) - k^2 n p = \lambda \left[ -\nabla \cdot (A\nabla (\Xi_i p)) + k^2 n (\Xi_i p) \right] \quad \text{op } \Omega^c_i,$$
   waarbij $\Xi_i$ een partitie van eenheid (partition of unity) operator is. De rechterzijde is aangepast om positief definiet te zijn (dicht bij de norm gebruikt in de analyse), terwijl de linkerzijde de indefiniete aard van de oorspronkelijke Helmholtz-operator behoudt.

2. Constructie van de Coarse Ruimte: De coarse ruimte $V_0$ wordt gespannen door de eigenfuncties die overeenkomen met de kleinste eigenwaarden (inclusief negatieve en kleine positieve waarden) van deze lokale GEVPs. De methode staat expliciet toe dat de lokale subdomeinen groot genoeg zijn zodat de lokale GEVPs indefiniet blijven, waardoor essentiële oscillerende modi worden gevangen voor hoogfrequente problemen.

3. Theoretisch Kader: De auteurs ontwikkelen een abstracte spectrale theorie voor indefiniete GEVPs. Zij stellen voorwaarden vast waaronder de bilineaire vorm geassocieerd met de GEVP positief definiet is op het complement van de subruimte gespannen door de behouden eigenfuncties. Deze theorie wordt toegepast om voldoende voorwaarden af te leiden voor de robuustheid van de gepreconditioneerde GMRES-oplosser.

4. Robuustheidsvoorwaarden: De analyse biedt expliciete voorwaarden voor robuuste convergentie in termen van de lokale subdomein-diameter $H_\ell$, de coarse subdomein-diameter $H_c$, en de eigenwaarde-tolerantie $\tau$ (de drempelwaarde voor het behouden van eigenmodi). Robuustheid wordt gegarandeerd indien:

   • $H_\ell \lesssim k^{-1}$ (om lokale oplosbaarheid te garanderen).
   • $\tau \gtrsim (1+C_{stab})^2 k^2$ (om voldoende modi te behouden).
   • Cruciaal is dat de theorie geen expliciete bovengrens oplegt aan de coarse diameter $H_c$, waardoor coarse subdomeinen groot kunnen zijn en de lokale GEVPs indefiniet kunnen zijn.

Belangrijkste Bijdragen

• Hk-GenEO Formulering: Introductie van een spectrale coarse ruimte gebouwd direct van indefiniete lokale operatoren, in contrast met eerdere $\Delta$-GenEO of $H$-GenEO benaderingen die positief semi-definiete operatoren gebruikten.
• Abstracte Indefiniete Spectrale Theorie: Ontwikkeling van een algemene theorie voor indefiniete GEVPs, inclus\text{ bij de} bestaan van een basis van eigenvectoren en de constructie van een projectie-operator die stabiliteit garandeert op het complement van de coarse ruimte.
• Expliciete Robuustheidsschattingen: Afleiding van $k$-expliciete voldoende voorwaarden voor de convergentie van GMRES. Het artikel bewijst dat de methode robuust is ten opzichte van het golfgetal $k$, het aantal subdomeinen en de coëfficiënt-heterogeniteit, mits de eigenwaarde-tolerantie $\tau$ gepast wordt gekozen.
• Ontkoppelde Decompositie: Theoretische en numerieke demonstratie dat de coarse decompositie (gebruikt voor eigenproblemen) kan worden ontkoppeld van de lokale decompositie (gebruikt voor één-niveau oplossingen), wat verschillende subdomeingroottes ($H_c$ versus $H_\ell$) toestaat.

Numerieke Resultaten
De auteurs presenteren numerieke experimenten in 2D voor zowel homogene als heterogene (gelaagde media) gevallen met behulp van het FreeFEM-framework.

• Robuustheid t.o.v. $k$: Voor matige tot grote drempelwaarden ($\tau = 0,4, 0,6$) blijven de GMRES-iteratieaantallen stabiel en in essentie onafhankelijk van $k$ en het aantal subdomeinen $N$.
• Gevoeligheid van de Drempelwaarde: Een kleine drempelwaarde ($\tau = 0,2$) leidt tot een significante toename van iteraties naarmate $k$ groeit, wat wijst op een verlies van robuustheid.
• Heterogeniteit: De methode vertoont veerkracht tegen hoge variaties in de coëfficiënt $A$, waarbij de iteratieaantallen slechts een milde afhankelijkheid laten zien van de contrastratio.
• Dimensie van de Coarse Ruimte: De dimensie van de coarse ruimte groeit matig met $k$ (empirisch waargenomen als $O(k^{1,1} \text{ tot } k^{1,4})$), wat aanzienlijk langzamer is dan de groei van de fine ruimte dimensie ($O(k^2)$).
• Ontkoppelde Decompositie: Experimenten tonen aan dat het gebruik van een grovere decompositie voor de coarse ruimte (minder subdomeinen $N$) vergeleken met de lokale oplosser ($Q$) slechts een bescheiden toename in iteratieaantallen oplevert, wat potentieel geheugenbesparing biedt.

Betekenis en Claims
Het artikel claimt dat Hk-GenEO een praktische en schaalbare preconditioner biedt voor sterk indefiniete Helmholtz-type problemen. De primaire betekenis ligt in:

1. Het Overwinnen van de Beperkingen van Positief Semi-Definiete Coarse Ruimtes: Door de indefiniete aard van de operator te gebruiken bij de constructie van de coarse ruimte, vangt de methode de noodzakelijke oscillerende modi op die positief semi-definiete benaderingen bij hoge frequenties missen.
2. Theoretische Rigor: Het biedt de eerste substantiële theoretische analyse voor spectrale coarse ruimtes gebaseerd op indefiniete operatoren, waarmee het verder gaat dan de empirische observaties in eerder werk (bijv. [5, 7]).
3. Flexibiliteit: De methode staat grote coarse subdomeinen en ontkoppelde decomposities toe, wat flexibiliteit biedt in algoritme-ontwerp die voorheen beperkt werd door de vereiste van kleine, coërcieve lokale problemen.

De auteurs merken op dat hoewel hun theoretische bovengrens voor de eigenwaarde-tolerantie $\tau$ (die schaalt met $k^2$) pessimistisch lijkt vergeleken met numerieke observaties (waar een constante $\tau$ vaak volstaat), de methode toch robuust blijft onder deze zwakkere praktische condities. Het werk legt een fundament voor het oplossen van indefiniete PDE's waarbij coëfficiënt-heterogeniteit en hoge golfgetallen voorheen significante uitdagingen vormden voor iteratieve oplossers.