The accuracy of the electric field integral equation (EFIE) can be substantially improved using high-order discretizations. However, this equation suffers from ill-conditioning and deleterious numerical effects in the low-frequency regime, often jeopardizing its solution. This can be fixed using quasi-Helmholtz decompositions, in which the source and testing elements are separated into their solenoidal and nonsolenoidal contributions, and then rescaled to avoid both the low-frequency conditioning breakdown and the loss of numerical accuracy. However, standard quasi-Helmholtz decompositions require handling discretized differential operators that often worsen the mesh refinement ill-conditioning and require the finding of the topological cycles of the geometry, which can be expensive when modeling complex scatterers, especially in high-order. This article solves these drawbacks by presenting the first extension of the quasi-Helmholtz projectors to high-order discretizations and their application to the stabilization of the EFIE when discretized with high-order basis functions. Our strategy will not require the identification of the cycles and will provide constant condition numbers for decreasing frequencies. Theoretical considerations will be accompanied by numerical results showing the effectiveness of our method in complex scenarios.
High-Order Quasi-Helmholtz Projectors: Definition, Analyses, Algorithms / Bourhis, Johann; Merlini, Adrien; Andriulli, Francesco P.. - In: IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION. - ISSN 0018-926X. - STAMPA. - 72:4(2024), pp. 3572-3579. [10.1109/tap.2024.3369752]
High-Order Quasi-Helmholtz Projectors: Definition, Analyses, Algorithms
Bourhis, Johann;Andriulli, Francesco P.
2024
Abstract
The accuracy of the electric field integral equation (EFIE) can be substantially improved using high-order discretizations. However, this equation suffers from ill-conditioning and deleterious numerical effects in the low-frequency regime, often jeopardizing its solution. This can be fixed using quasi-Helmholtz decompositions, in which the source and testing elements are separated into their solenoidal and nonsolenoidal contributions, and then rescaled to avoid both the low-frequency conditioning breakdown and the loss of numerical accuracy. However, standard quasi-Helmholtz decompositions require handling discretized differential operators that often worsen the mesh refinement ill-conditioning and require the finding of the topological cycles of the geometry, which can be expensive when modeling complex scatterers, especially in high-order. This article solves these drawbacks by presenting the first extension of the quasi-Helmholtz projectors to high-order discretizations and their application to the stabilization of the EFIE when discretized with high-order basis functions. Our strategy will not require the identification of the cycles and will provide constant condition numbers for decreasing frequencies. Theoretical considerations will be accompanied by numerical results showing the effectiveness of our method in complex scenarios.File | Dimensione | Formato | |
---|---|---|---|
High-Order_Quasi-Helmholtz_Projectors_Definition_Analyses_Algorithms.pdf
accesso aperto
Tipologia:
2a Post-print versione editoriale / Version of Record
Licenza:
Creative commons
Dimensione
7.5 MB
Formato
Adobe PDF
|
7.5 MB | Adobe PDF | Visualizza/Apri |
Pubblicazioni consigliate
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.
https://hdl.handle.net/11583/2991303