Evaluation of resonances: adaptivity and AAA rational approximation of randomly scalarized boundary integral resolvents

Technische Samenvatting: Evaluatie van Resonanties via Adaptieve AAA Rationale Approximatie

Probleemstelling
Het artikel behandelt de computationele uitdaging van het evalueren van akoestische en elektromagnetische resonanties in open en gesloten holtes, evenals in andere verstrooiingsstructuren. Wiskundig gezien komen deze resonanties overeen met oplossingen $(u, k)$ van het Helmholtz-eigenwaardeprobleem $\Delta u + k^2 u = 0$, onderworpen aan homogene randvoorwaarden (Dirichlet, Neumann of Zaremba) op een domein $\Omega$.

• Voor interne problemen representeert $k$ reële eigenwaarden.
• Voor externe en open-holteproblemen representeert $k$ complexe "verstrooiingspolen" of "complexe eigenwaarden" die zich in de onderste helft van het complexe vlak bevinden ($\Im k \le 0$).

Het probleem wordt geherformuleerd als een Nietlineair Eigenwaardeprobleem (NEP) voor een matrix-waardige functie $F(k)$, afgeleid van de discretisatie van randintegraaloperatoren. Het doel is om complexe waarden $k$ te vinden waarvoor $F(k)w = 0$ geldt voor een niet-nul vector $w$.

Methodologie
Het voorgestelde algoritme combineert drie primaire componenten: een willekeurig gescalariseerde resolvent, het AAA rationale approximatiealgoritme, en een adaptieve zoekstrategie met een secant-methode verfijningsfase.

1. Gescalariseerde Resolvent: In plaats van het volledige matrix-NEP direct op te lossen, evalueert de methode een scalaire functie $S(k) = u^* F(k)^{-1} v$, waarbij $u$ en $v$ vaste willekeurige vectoren in $\mathbb{C}^d$ zijn. De polen van $S(k)$ komen overeen met de eigenwaarden van het oorspronkelijke probleem.
2. AAA Rationale Approximatie: De functie $S(k)$ wordt benaderd door een rationale functie $r(k)$ met behulp van het AAA-algoritme. Dit algoritme selecteert gulzig interpolatiepunten en gewichten om een barycentrische rationale approximant te construeren. De polen van deze rationale functie $r(k)$ dienen als initiële benaderingen voor de resonanties.
3. Adaptieve Zoekstrategie: Om ervoor te zorgen dat alle resonanties binnen een zoekdomein worden gevangen, maakt het algoritme gebruik van een recursieve dyadische onderverdeling van het domein (rechthoeken in het complexe vlak of intervallen op de reële lijn).
   • Algoritme 1 (Basis): Evalueert $S(k)$ op een set punten, construeert de AAA-approximant en geeft polen binnen het domein terug.
   • Algoritme 2 (Adaptief): Onderverdeelt het domein recursief. Als het aantal gevonden polen in een subregio verandert bij verdere onderverdeling, wordt het proces herhaald. Dit zorgt ervoor dat geen resonanties worden gemist door een onvoldoende steekproefdichtheid.
4. Verfijning en Verificatie:
   • Secant-methode: Zodra initiële polen zijn geïdentificeerd, wordt een secant-methode terminatiefase toegepast op $1/S(k)$ om exponentiële convergentie naar bijna machineprecisie te bereiken. Deze stap biedt ook een betrouwbare foutinschatting en helpt bij het opsporen van spookeigenwaarden die het AAA-algoritme incidenteel kan produceren.
   • Eigenfuncties: De bijbehorende eigenfuncties worden gerecupereerd door de nulruimte van $F(k)$ te berekenen bij de verfijnde eigenwaarden met behulp van Gaussische eliminatie met pivoting.

Belangrijkste Bijdragen

• Directe Polenextractie: In tegenstelling tot eerdere AAA-gebaseerde benaderingen die rationale approximatie gebruiken om het NEP te lineariseren, gebruikt deze methode direct de polen van de rationale approximant van de gescalariseerde resolvent als de eigenwaarde-schattingen.
• Geometrische Flexibiliteit: De afhankelijkheid van rationale approximatie in plaats van contourintegratie maakt het gebruik van willekeurige zoekdomeinen mogelijk, inclusief rechthoekige regio's en reële intervallen. Dit faciliteert het ontwerp van de adaptieve onderverdelingsstrategie, wat moeilijk te implementeren is met traditionele contourintegratiemethoden (bijv. de trapeziumregel op rechthoeken).
• Robuustheid en Efficiëntie: De combinatie van adaptieve zoekstrategie en secant-verfijning garandeert een hoge nauwkeurigheid (nabij machineprecisie) voor zowel laagfrequente als hoogfrequente toestanden. De methode vereist alleen de evaluatie van de gescalariseerde resolvent bij verschillende golfgetallen, wat extreem goed paralleliseerbaar is.
• Afhandeling van Open Bogen: Het artikel breidt het theoretische kader voor het relateren van resolventpolen aan eigenwaarden uit naar open-boogproblemen, waarbij specifieke sprongrelaties voor de normale afgeleide van de single-layer potentiaal worden gebruikt om randsingulariteiten te behandelen.

Numerieke Resultaten
Het artikel presenteert een verscheidenheid aan numerieke experimenten die de prestaties van de methode demonstreren:

• Algemene NEP's: Het algoritme identificeerde succesvol alle eigenwaarden in de "CD player" (60 eigenwaarden) en "Butterfly" (256 eigenwaarden) problemen uit de NLEVP-collectie.
• Laagfrequente Verstrooiing: Nauwkeurige berekening van eigenwaarden voor vliegervormige en cirkelvormige holtes met openingen. De methode handelde zowel reële eigenwaarden (intern) als complexe polen nabij de reële as (open holtes) af.
• Vergelijking met Contourintegratie: Vergelijkingen met de block SS-methode en Beyns "Integral Algorithm 1" lieten zien dat hoewel Beyns methode sneller kan convergeren zonder verfijning, de voorgestelde methode een vergelijkbare nauwkeurigheid bereikt met lokale verfijning. Cruciaal is dat de voorgestelde methode de noodzaak voor meerdere matrix-oplossingen per frequentiepunt vermijdt (vereist door block SS/Beyn wanneer $L > 1$), wat potentiële voordelen biedt voor grootschalige 3D-problemen waarbij LU-factorisatie niet beschikbaar is.
• Gap Size Asymptotics: De methode werd gebruikt om de convergentie van verstrooiingspolen te bestuderen wanneer een open boog een gesloten curve nadert. De resultaten toonden aan dat vervalrates ($\Im k$) convergeren met snelheden van $O(\theta^4)$ of $O(\theta^8)$ afhankelijk van de modus-uitlijning met de opening, terwijl ruimtelijke frequenties ($\Re k$) convergeren met $O(\theta^2)$ of $O(\theta^4)$.
• Hoogfrequente Prestaties: Het algoritme berekende succesvol 1.244 onderscheidende interne eigenwaarden voor een eenheidscirkel in het interval $[1, 100]$ met bijna machineprecisie. Het visualiseerde ook hoogfrequente eigenfuncties voor vlieger- en raketvormige domeinen, wat een sterke correspondentie toont tussen open-holte eigenfuncties en verstrooiingsoplossingen onder vlakgolfverlichting.

Betekenis en Claims
De auteurs beweren dat de voorgestelde aanpak eenvoudig, gemakkelijk te implementeren en snel convergerend is. De primaire betekenis ligt in het vermogen om zeer nauwkeurige benaderingen van verstrooiingsresonanties en NEP-oplossingen te bieden over een breed frequentiebereik zonder de geometrische beperkingen die worden opgelegd door contourintegratiemethoden. De adaptieve strategie garandeert exhaustieve zoekmogelijkheden, terwijl de secant-verfijningsfase hoge precisie en betrouwbaarheid garandeert. Het artikel positioneert deze methode als een flexibel alternatief voor bestaande contourintegratie-gebaseerde eigensolvers, wat bijzonder voordelig is voor problemen met complexe geometrieën of waar computationele middelen voor herhaalde matrix-factorisaties beperkt zijn.