A Stabilized Multilevel B-Spline-Based Fast Integral Method for the Solution of the Electric Field Integral Equation

Technische Samenvatting: Een Gestabiliseerde Multilevel B-spline-gebaseerde Snelle Integrale Methode voor de EFIE

Probleemstelling
De Elektrische Veld Integrale Vergelijking (EFIE) is een standaardinstrument voor het oplossen van elektromagnetische verstrooiings- en stralingsproblemen, waarbij enkel oppervlakte-discretisatie vereist is. Echter, standaard Boundary Element Method (BEM) discretisaties resulteren in $O(N^2)$ complexiteit voor matrixassemblage, opslag en matrix-vectorproducten (MVP), waarbij $N$ het aantal onbekenden is. Snelle integrale methoden, zoals de Multilevel Fast Multipole Method (MLFMM) en FFT-versnelde benaderingen zoals de Adaptive Integral Method (AIM), zijn ontwikkeld om deze complexiteit te reduceren naar $O(N)$ of bijna-$O(N)$.

Er bestaat een specifieke beperking in FFT-versnelde methoden die vertrouwen op globale Lagrange-interpolatie. Deze methoden vereisen equidistante interpolatiepunten om gebruik te maken van FFT-structuren. Echter, equidistante Lagrange-interpolatie lijdt aan Runge-type instabiliteiten bij hoge polynomiale orden, wat leidt tot een afname in nauwkeurigheid en het voorkomen van robuuste hogere-orde compressie. Hoewel B-splines een oplossing bieden voor Runge-instabiliteiten op equidistante roosters, is hun toepassing in multilevel snelle integrale methoden niet triviaal. In tegen tegenstelling tot Lagrange-polynomen komen de B-spline coëfficiënten niet overeen met de functiewaarden op interpolatiepunten, wat de inversie van bemonsteringsmatrices noodzakelijk maakt. Deze matrices kunnen zeer slecht geconditioneerd raken bij hoge orden, wat de numerieke stabiliteit in gevaar brengt. Bovendien is het handhaven van exacte interlevel-overdrachten (vereist voor multilevel efficiëntie) terwijl stabiliteit wordt gewaarborgd, een aanzienlijke uitdaging gebleken.

Methodologie
De auteurs stellen een gestabiliseerde multilevel B-spline-gebaseerde snelle integrale methode voor die FFT-compatibele kernel-interpolatie combineert met robuuste hogere-orde interpolatie. De kernmethodologie bestaat uit drie primaire componenten:

1. B-spline Kernel-interpolatie:
   In plaats van globale Lagrange-polynomen gebruikt de methode B-splines gedefinieerd op equidistante knoopvectoren. Deze keuze mitigeert Runge-type instabiliteiten terwijl de compatibiliteit met de voor FFT vereiste equidistante bemonstering behouden blijft. De Green-functie wordt gefactoriseerd met behulp van tensorproduct-polynomen. Voor B-splines zijn de interpolatiegewichten geen directe samples van de Green-functie, maar worden ze afgeleid via een bemonsteringsmatrix $B$. De tensorproductstructuur van de B-splines maakt het mogelijk om de inversie van de bemonsteringsmatrix efficiënt uit te voeren in $O(s^3)$ operaties (waarbij $s$ het aantal splines is) in plaats van de naïeve $O(s^6)$, waardoor de totale $O(s^3 \log s)$ complexiteit van de MVP behouden blijft.

2. Stabilisatie via Knoopverwijdering:
   Om de slecht geconditioneerde B-spline bemonsteringsmatrices bij hoge orden aan te pakken, introduceren de auteurs een "knoopverwijderings"-stabilisatiestrategie. Empirische analyse toont aan dat het verwijderen van specifieke knopen nabij de grenzen van de equidistante knoopvector de stabiliteitsconstante (de norm van de inverse bemonsteringsmatrix) aanzienlijk vermindert zonder de equidistante structuur die nodig is voor het multilevel schema te vernietigen.

   • Bufferknopen: Om basisfuncties te behandelen die buiten hun bijbehorende boxen uitsteken (een veelvoorkomend probleem in BEM), breidt de methode de knoopvector uit met "bufferknopen". Dit zorgt ervoor dat het ondersteuningsgebied van alle basisfuncties binnen het interpolatiedomein valt.
   • Exacte Interlevel-overdrachten: Een cruciale vereiste voor multilevel methoden is het vermogen om momenten tussen ouder- en kind-boxen exact over te dragen. De voorgestelde stabilisatiestrategie behoudt de geneste structuur van de knoopvectoren. Door "knoopinvoeging" (een standaard spline-verfijningstechniek) te gebruiken, construeert de methode exacte transfermatrices tussen niveaus. Dit zorgt ervoor dat momenten berekend op het gestabiliseerde (met verwijderde grensknopen) rooster exact kunnen worden gemapt naar het oorspronkelijke equidistante rooster dat wordt gebruikt voor kernel-interpolatie, waardoor machine-precisie nauwkeurigheid over de niveaus heen wordt behouden.

3. Factorisatiestrategieën:
   Het artikel introduceert een nieuwe "ster-gebaseerde" factorisatie voor de gediscretiseerde EFIE, naast bestaande directe en dyadische factorisaties.

   • Ster-gebaseerde Factorisatie: Deze benadering maakt gebruik van de representatie van de scalaire potentiaaloperator met behulp van een loop-star decompositieconcept. Het factoriseert de scalaire potentiaaloperator in een ijle ster-matrix en een single-layer operator. Deze strategie behoudt de nulruimte van de scalaire potentiaaloperator tot op machine-precisie, wat cruciaal is voor de effectiviteit van lage-frequentie preconditioning technieken (bijv. Calderón preconditioning). Dit vereist het berekenen van zowel RWG als piecewise-constant momenten.
   • Directe en Dyadische Factorisaties: Deze worden besproken als alternatieven. De directe factorisatie interpoleert de kernel direct, terwijl de dyadische formulering een dyadische kernel gebruikt. De ster-gebaseerde benadering wordt uitgelicht vanwege haar superieure stabiliteitseigenschappen met betrekking tot de nulruimte.

Belangrijkste Bijdragen

• Gestabiliseerde B-spline Interpolatie: De introductie van een knoopverwijderingsstrategie die het Runge-fenomeen en de slecht geconditioneerde aard van hoge-orde equidistante B-spline interpolatie elimineert, wat robuuste hogere-orde compressie mogelijk maakt.
• Exacte Multilevel Overdrachten: De ontwikkeling van een mechanisme om exacte interlevel-overdrachten uit te voeren met behulp van knoopinvoeging, zelfs wanneer stabilisatie (knoopverwijdering) wordt toegepast. Dit lost het conflict op tussen de behoefte aan stabiele interpolatie en de behoefte aan exacte hiërarchische transfers.
• Ster-gebaseerde Factorisatie: De voorstel van een nieuwe factorisatiestrategie voor de EFIE die de nulruimte van de scalaire potentiaaloperator tot op machine-precisie behoudt, wat de bruikbaarheid voor lage-frequentie problemen verhoogt.
• Complexiteitsbehoud: Demonstratie dat de voorgestelde methode $O(N)$ complexiteit behoudt voor MVP, assemblage en opslag, door gebruik te maken van de compacte ondersteuning van B-splines en de tensorproductstructuur van de bemonsteringsmatrices.

Resultaten
Numerieke resultaten werden gepresenteerd voor zowel canonieke (bol) als realistische (vliegtuig) geometrieën over diverse elektrische groottes (van $0.2\lambda$ tot $4\lambda$).

• Stabiliteit en Nauwkeurigheid: De B-spline-gebaseerde methode vertoonde stabiele hogere-orde interpolatie zonder de instabiliteit die bij Lagrange-benaderingen wordt waargenomen. Foutanalyse toonde convergentie die consistent is met theoretische grenzen naarmate de polynomiale graad $p$ en het aantal splines $s$ toenamen. De stabilisatiestrategie slaagde erin fouten te beheersen, zelfs bij hoge orden ($p=5, 6$), terwijl Lagrange-methoden instabiliteit vertoonden.
• Factorisatie Prestaties: De ster-gebaseerde factorisatie vertoonde het meest robuuste foutgedrag onder de drie geteste strategieën.
• Schaalbaarheid: Prestatieanalyse bevestigde $O(N)$ schaling voor opslag, assemblage tijd en MVP-tijd voor zowel het bol- als het vliegtuigmodel. Voor een bol met $N \approx 3$ miljoen onbekenden varieerde het geheugengebruik van 91 GB tot 182 GB afhankelijk van $p$, en de MVP-tijden varieerden van 7,6s tot 54s. Vergelijkbare lineaire schaling werd waargenomen voor het vliegtuigmodel met tot $\approx 4$ miljoen onbekenden.

Betekenis
Het artikel beweert dat de voorgestelde gestabiliseerde multilevel B-spline methode een praktische, hoog-nauwkeurige alternatief biedt voor Lagrange-interpolatie voor snelle integrale vergelingssolvers. Door het Runge-fenomeen en de slecht geconditioneerde aard van bemonsteringsmatrices te overwinnen, maakt de methode robuuste hogere-orde compressie mogelijk die voorheen onbereikbaar was met equidistante FFT-versnelde schema's. Het behoud van de nulruimte van de scalaire potentiaal via de ster-gebaseerde factorisatie verhoogt verder de bruikbaarheid van de methode voor lage-frequentie problemen. De auteurs merken op dat voorlopige resultaten op conferenties zijn gepresenteerd en dat de methode gecombineerd kan worden met hogere-frequentie technieken zoals MLFMM voor breedbandformuleringen, waarbij de primaire bijdrage de stabilisatie van de lage- tot middenfrequentie kernel-interpolatie raamwerk is.