Suppressing Low-Frequency Cancellation Errors in Far-Field Computation with the MFIE, Müller, and Multitrace Müller Integral Equations

Technische Samenvatting: Onderdrukking van laagfrequente annuleringfouten in verre-veld berekeningen met MFIE, Müller en Multitrace Müller integraalvergelijkingen

Probleemstelling
Randintegraalvergelijkingen (BIE's), zoals de Magnetische Veld Integrale Vergelijking (MFIE) voor Perfect Elektrisch Geleiders (PEC's) en de Müller-vergelijkingen voor diëlektrische objecten, worden breed gebruikt voor elektromagnetische verstrooiingsmodellering. Hoewel vergelijkingen van de tweede soort, zoals de MFIE en de Müller-vergelijkingen, inherent goed geconditioneerd zijn in de dense-mesh en lage-frequentie regimes, lijden ze onder een fundamenteel numeriek probleem dat bekend staat als laagfrequente annuleringfouten (of de lage-frequentie catastrofe) wanneer deze worden geëvalueerd in rekenkunde met eindige precisie.

Dit fenomeen ontstaat door de uiteenlopende asymptotische schaling van de Helmholtz-componenten van oppervlaktestromen (irrotationele/niet-solenoidale, harmonische/globale lus, en solenoidale/lokale lus) wanneer de hoekfrequentie $\omega \to 0$. In standaard drijvende-komma representaties worden de kleinere magnitude componenten (typisch de niet-solenoidale component in PEC-verstrooiing) overstemd door dominante componenten, wat leidt tot een verlies van significante cijfers. Hoewel deze fouten verwaarloosbaar kunnen blijven bij nabijveld- of oppervlaktestroomreconstructies, degraderen ze de nauwkeurigheid van verre-veld evaluaties kritisch. Specifiek hangt het verre-veld patroon af van een combinatie van deze componenten waarbij de niet-solenoidale component, ondanks dat deze asymptotisch kleiner is, bijdraagt op dezelfde orde als de solenoidale delen. Standaard discretisaties falen in het behouden van de noodzakelijke schaleringsrelaties om dit verlies aan precisie te voorkomen, wat resulteert in onbetrouwbare stralingspatronen bij lage frequenties.

Methodologie
De auteurs stellen een robuust numeriek kader voor dat quasi-Helmholtz projectoren combineert met een op maat gemaakte herschalingsprocedure om deze fouten te onderdrukken, terwijl de gunstige conditionering van de onderliggende gemengde discretisaties behouden blijft.

1. Quasi-Helmholtz Decompositie: De oppervlaktestromen worden gescheiden in niet-solenoidale ($j_\Sigma$) en solenoidale ($j_{\Lambda H}$) componenten met behulp van quasi-Helmholtz projectoren ($P_\Sigma$ en $P_{\Lambda H}$). Deze decompositie is stabiel met betrekking tot mesh-verfijning.
2. Rechter Herschaling: Een hulp-onbekende wordt geïntroduceerd via een herschalingsoperator $P = (i\omega c_0^{-1} d) P_\Sigma + P_{\Lambda H}$, waarbij $d$ een karakteristieke lengte is. Dit balanceert de magnituden van de componenten, waardoor ze van vergelijkbare orde ($O(1)$) zijn in de hulp-variabele, wat het verlies van significante cijfers tijdens opslag voorkomt.
3. Stabilisatie in Coëfficiëntruimte: Om de kunstmatige ill-conditioning tegen te gaan die wordt geïntroduceerd door de rechter herscaling en om de excitatievector te balanceren, passen de auteurs een strategie toe die eerder werd gebruikt voor PMCHWT-formuleringen. Het systeem wordt eerst vermenigvuldigd met de inverse Gram-matrix ($G^{-1}$), die de operator terug naar de coëfficiëntruimte brengt. Vervolgens wordt een linker herscalingsoperator ($P^{-1}$) toegepast.
   • Voor de MFIE levert dit het gestabiliseerde systeem op: $(\frac{1}{2}I - P^{-1}G^{-1}K(\kappa_0)P)i = P^{-1}G^{-1}h_{inc}$.
   • Voor de Müller-vergelijking (homogene diëlektrica) en de Multitrace (MT) Müller-vergelijking (composietstructuren), wordt dit uitgebreid met blokdiagonaal uitbreidingen van de projectoren en herschalingsoperatoren.
4. Statische Extractie: Om te waarborgen dat de theoretische $O(\omega^2)$ schaling van specifieke operatorblokken numeriek behouden blijft, extraheren de auteurs analytisch de statische bijdragen van de Green's functie kernen vóór de kwadratuur. Dit voorkomt dat kwadratuurfouten uitvergroten tot spookachtige laagfrequente termen.
5. Verre-veld Evaluatie: Het verre veld wordt direct berekend vanuit de herschale hulp-onbekenden met behulp van een formulering die de annulering-vrije representatie respecteert, waardoor de reconstructie van de fysieke stroom wordt vermeden die de annuleringfouten opnieuw zou introduceren.

Belangrijkste Bijdragen

• Robuuste Formuleringen voor MFIE en Müller: Het artikel presenteert de eerste analoge stabilisatiemethoden voor de MFIE en Müller-vergelijkingen die niet vertrouwen op symmetrisatie of hulp-ladingvergelijkingen. Dit vermijdt het verhogen van het aantal onbekenden of het compliceren van de systeemstructuur.
• Extensie naar Composietstructuren: De methode wordt succesvol uitgebreid naar de globale Multitrace Müller-formulering, wat nauwkeurige verre-veld berekeningen mogelijk maakt voor complexe composiet diëlektrische objecten met juncties, niet-conforme meshes en hoge materiaalkontrasten.
• Behoud van Conditionering: In tegenstelling tot veel stabilisatietechnieken die de conditioneringsgetal verslechteren, behoudt de voorgestelde methode de goed geconditioneerde aard van de gemengde MFIE en Müller discretisaties in de lage-frequentie en dense-mesh regimes.
• Afhandeling van Topologie: De aanpak adresseert meervoudig verbonden geometrieën. Hoewel de statische MFIE-operator een door topologie geïnduceerde nulruimte bezit op dergelijke oppervlakken, zorgt de voorgestelde methode ervoor dat dit de verre-veld nauwkeurigheid niet in gevaar brengt, aangezien de vlakgolf-excitatie niet koppelt aan de nulruimte-modi in de verre-veld evaluatie.

Numerieke Resultaten
Uitgebreide numerieke experimenten valideren de voorgestelde methoden over een breed frequentiebereik (van full-wave tot extreem lage frequenties, bijv. $10^{-4}$ Hz) en diverse geometrieën:

• PEC-verstrooiing: Voor zowel een eenheidssfeer (eenvoudig verbonden) als een vierkante torus (meervoudig verbonden), behoudt de gestabiliseerde MFIE (qHP-MFIE) de nauwkeurigheid waar klassieke en gemengde MFIE formuleringen falen bij lage frequenties.
• Homogene Diëlektrica: De gestabiliseerde Müller-formulering (qHP-Müller) presteert beter dan klassieke en gemengde Müller discretisaties voor sferen en tori met hoge permittiviteit contrasten ($\epsilon_r = 100$).
• Composiet Diëlektrica: De gestabiliseerde MT-Müller formulering (qHP MT-Müller) demonstreert robuustheid voor multiscale, hoog-contrast composiet objecten met juncties, en komt overeen met referentie-oplossingen van commerciële solvers (Altair Feko) tot aan extreem lage frequenties.
• Conditionering: Iteratie-aantallen voor GMRES-solvers blijven begrensd en nemen vaak af met de voorgestelde regularisatie vergeleken met de ongestabiliseerde gemengde formuleringen, zelfs voor meervoudig verbonden geometrieën.

Betekenis en Claims
De auteurs beweren dat hun werk een robuust en relatief eenvoudig alternatief biedt voor elektromagnetische verstrooiingssimulaties die variëren van matig hoge tot extreem lage frequenties. Door zich te concentreren op de MFIE en Müller-vergelijkingen, bieden zij een oplossing die de computationele overhead en complexiteit vermijdt die geassocieerd wordt met het stabiliseren van eerste-orde formuleringen (zoals EFIE of PMCHWT) of het introduceren van hulp-vergelijkingen.

Het artikel benadrukt dat de voorgestelde regularisatie bijzonder voordelig is voor:

• Composiet configuraties bestaande uit meerdere objecten en juncties.
• Tijddomein simulaties, waar bestaande benaderingen tot onwerkbaar complexe systemen kunnen leiden.
• Nauwkeurige verre-veld voorspelling zonder onnodige complexiteit te introduceren, waardoor de essentiële componenten voor verre-veld berekening beschermd worden tegen annulering in eindige precisie.

De auteurs merken op dat hoewel de methode de intrinsieke resonantie-instabiliteit van de MFIE of de door topologie geïnduceerde breakdown bij nul-frequentie op meervoudig verbonden oppervlakken niet elimineert, het deze problemen effectief ontkoppelt van de verre-veld evaluatie, waardoor snelle en nauwkeurige oplossingen met Krylov iteratieve solvers mogelijk zijn. Toekomstig werk wordt gesuggereerd voor het uitbreiden van dit kader naar tijddomein vergelijkingen en gemengde PEC-diëlektrische verstrooiers.