High-Order-Accurate Continuity Enforcing Nyström Discretization of 3D Maxwell Combined Field Integral Equations
Abstract
In Nyström-collocation discretizations of the electric field integral equation (EFIE), the surface divergence acts on surface densities that may be discontinuous across patch boundaries, which degrades accuracy and convergence. We show that this not only affects the EFIE but every formulation in which the operator occurs, either in the equation itself or in the scattered field computation, and propose a high-order-accurate continuity-enforcing scheme for smooth surfaces as a remedy for the direct and indirect EFIEs, magnetic field integral equations (MFIEs), and regularized combined field integral equations (CFIEs) alike. The scheme comprises two ingredients: i) We show how to discretize the equations via a Chebyshev-based Nyström scheme, which admits closed quadrature rules. ii) Since unknowns and test vectors are in terms of patch-local curvilinear bases, continuity is enforced by a change of basis: we construct sparse mapping matrices assembled solely from the curvilinear geometry description. In doing so, we restore the accuracy of the EFIE such that it can be combined with the MFIEs with equal weights to form CFIEs. Numerical studies for the scattering from canonical and realistic geometries show that all considered formulations individually and combined benefit from the continuity enforcement in terms of better conditioning, reduced iterations of an iterative solver, and several more digits of accuracy in the scattered fields, despite reducing the total number of unknowns.
Related
- A dual-field structure-preserving mixed finite element discretization for incompressible Hall MHD equations
- An efficient coupled block preconditioner for the Marker-and-Cell discretization of the Stokes-Darcy equations
- Regularization and Hölder continuity for complex Hessian equations on Hermitian manifolds
- Numerical Methods for Dynamical Low-Rank Approximations of Stochastic Differential Equations -- Part II: Stochastic discretization
- Parameter-robust preconditioners for a stabilized control-volume discretization of Biot's equations
BibTeX
@misc{hofmann2026highorderaccurate,
title = {{High-Order-Accurate Continuity Enforcing Nyström Discretization of 3D Maxwell Combined Field Integral Equations}},
author = {Bernd Hofmann and Reza Molavi and Constantine Sideris},
howpublished = {arXiv},
year = {2026},
month = aug,
doi = {10.48550/arxiv.2608.28876},
eprint = {2608.28876},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2608.28876},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - High-Order-Accurate Continuity Enforcing Nyström Discretization of 3D Maxwell Combined Field Integral Equations
AU - Bernd Hofmann
AU - Reza Molavi
AU - Constantine Sideris
PY - 2026
DA - 2026/08/28
PB - arXiv
DO - 10.48550/arxiv.2608.28876
UR - https://arxiv.org/abs/2608.28876
AB - In Nyström-collocation discretizations of the electric field integral equation (EFIE), the surface divergence acts on surface densities that may be discontinuous across patch boundaries, which degrades accuracy and convergence. We show that this not only affects the EFIE but every formulation in which the operator occurs, either in the equation itself or in the scattered field computation, and propose a high-order-accurate continuity-enforcing scheme for smooth surfaces as a remedy for the direct and indirect EFIEs, magnetic field integral equations (MFIEs), and regularized combined field integral equations (CFIEs) alike. The scheme comprises two ingredients: i) We show how to discretize the equations via a Chebyshev-based Nyström scheme, which admits closed quadrature rules. ii) Since unknowns and test vectors are in terms of patch-local curvilinear bases, continuity is enforced by a change of basis: we construct sparse mapping matrices assembled solely from the curvilinear geometry description. In doing so, we restore the accuracy of the EFIE such that it can be combined with the MFIEs with equal weights to form CFIEs. Numerical studies for the scattering from canonical and realistic geometries show that all considered formulations individually and combined benefit from the continuity enforcement in terms of better conditioning, reduced iterations of an iterative solver, and several more digits of accuracy in the scattered fields, despite reducing the total number of unknowns.
ER -