quasi-Helmholtz Projectors Research Articles

On the Use of Quasi-Helmholtz Projectors for Modeling Structures With Junctions at Arbitrarily Low Frequencies

On the Use of Quasi-Helmholtz Projectors for Modeling Structures With Junctions at Arbitrarily Low Frequencies

A Stabilized Time-Domain Combined Field Integral Equation Using the Quasi-Helmholtz Projectors

A Stabilized Time-Domain Combined Field Integral Equation Using the Quasi-Helmholtz Projectors

On a Low-Frequency and Contrast-Stabilized Full-Wave Volume Integral Equation Solver for Lossy Media

In this article, we present a new regularized electric flux volume integral equation (D-VIE) for modeling high-contrast conductive dielectric objects in a broad frequency range. This new formulation is particularly suitable for modeling biological tissues at low frequencies, as it is required by brain epileptogenic area imaging, but also at higher ones, as it is required by several applications, including, but not limited to, deep brain stimulation (DBS). When modeling inhomogeneous objects with high complex permittivities at low frequencies, the traditional D-VIE is ill-conditioned and suffers from numerical instabilities that result in slower convergence and less accurate solutions. In this work, we address these shortcomings by leveraging a new set of volume quasi-Helmholtz projectors. Their scaling by the material permittivity matrix allows for the rebalancing of the equation when applied to inhomogeneous scatterers and, thereby, makes the proposed method accurate and stable even for high complex permittivity objects until arbitrarily low frequencies. Numerical results, canonical and realistic, corroborate the theory and confirm the stability and the accuracy of this new method both in the quasi-static regime and at higher frequencies.

Improved A-EFIE System for Electromagnetic Simulation in Low Frequency Regime

In this letter, the quasi-Helmholtz projectors are applied to address numerical cancellation of augmented electric field integral equation (A-EFIE). The projectors rescale the term <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"><tex-math notation="LaTeX">$\overline{\bf V} \cdot i{k}_0{\bf j}$</tex-math></inline-formula> and right-hand side to guarantee that the solution is stable at very low frequency. The improved system is full rank without removing the singularity caused by charge neutrality. Compared to the EFIE with quasi-Helmholtz projectors, this letter presents a better convergence performance and takes less computational cost. The numerical results show that the solution of both scattering problems and circuit problems can ensure the accuracy in the low-frequency regime.

Magnetic and Combined Field Integral Equations Based on the Quasi-Helmholtz Projectors

Boundary integral equation methods for analyzing electromagnetic scattering phenomena typically suffer from several of the following problems: (i) ill-conditioning when the frequency is low; (ii) ill-conditioning when the discretization density is high; (iii) ill-conditioning when the structure contains global loops (which are computationally expensive to detect); (iv) incorrect solution at low frequencies due to current cancellations; (v) presence of spurious resonances. In this paper, quasi-Helmholtz projectors are leveraged to obtain a magnetic field integral equation (MFIE) formulation that is immune to drawbacks (i)-(iv). Moreover, when this new MFIE is combined with a regularized electric field integral equation, a new quasi-Helmholtz projector combined field integral equation is obtained that also is immune to (v). Numerical results corroborate the theory and show the practical impact of the newly proposed formulations.

Large Time Step and DC Stable TD-EFIE Discretized With Implicit Runge–Kutta Methods

The Time Domain-Electric Field Integral Equation (TD-EFIE) and its differentiated version are widely used to simulate the transient scattering of a time dependent electromagnetic field by a Perfect Electrical Conductor (PEC). The time discretization of the TD-EFIE can be achieved by a space-time Galerkin approach or, as it is considered in this contribution, by a convolution quadrature using Implicit Runge-Kutta methods. The solution is then computed using the Marching-On-in-Time (MOT) algorithm. The differentiated TD-EFIE has two problems: (i) the system matrix suffers from ill-conditioning when the time step increases (low frequency breakdown) and (ii) it suffers from the DC instability, i.e. the formulation allows for the existence of spurious solenoidal currents that grow slowly in the solution. In this work, we show that (i) and (ii) can be alleviated by leveraging quasi-Helmholtz projectors to separate the Helmholtz components of the induced current and rescale them independently. The efficacy of the approach is demonstrated by numerical examples including benchmarks and real life applications.

Eddy Current Modeling in Multiply Connected Regions via a Full-Wave Solver Based on the Quasi-Helmholtz Projectors

Eddy currents are central to several industrial applications and there is a strong need for their efficient modeling. Existing eddy current solution strategies are based on a quasi-static approximation of Maxwell’s equations for lossy conducting objects and thus their applicability is restricted to low frequencies. On the other hand, available full-wave solvers such as the Poggio-Miller-Chang-Harrington-Wu-Tsai (PMCHWT) equation become highly ill-conditioned and inaccurate in eddy current settings. This work presents a new well-conditioned and stable full-wave formulation which encompasses the simulation of eddy currents. Our method is built upon the PMCHWT equation and thus remains valid over the entire frequency range. Moreover, our scheme is also compatible with structures containing holes and handles (multiply connected geometries). The effectiveness of quasi-Helmholtz projectors is leveraged to obtain a versatile solver, which is computationally efficient and allows for a seamless transition between low and high frequencies. The stability and accuracy of the new method are demonstrated both theoretically and through numerical experiments on canonical and realistic structures.

On a Low-Frequency and Refinement Stable PMCHWT Integral Equation Leveraging the Quasi-Helmholtz Projectors

Classical Poggio–Miller–Chan–Harrington–Wu–Tsai (PMCHWT) formulations for modeling radiation and scattering from penetrable objects suffer from ill-conditioning when the frequency is low or when the mesh density is high. The most effective techniques to solve these problems, unfortunately, either require the explicit detection of the so-called global loops of the structure, or suffer from numerical cancellation at extremely low frequency. In this contribution, a novel regularization method for the PMCHWT equation is proposed, which is based on the quasi-Helmholtz projectors. This method not only solves both the low frequency and the dense mesh ill-conditioning problems of the PMCHWT, but it is immune from low-frequency numerical cancellations and it does not require the detection of global loops. This is obtained by projecting the range space of the PMCHWT operator onto a dual basis, by rescaling the resulting quasi-Helmholtz components, by replicating the strategy in the dual space, and finally, by combining the primal and the dual equations in a Calderon-like fashion. Implementation-related treatments and details alternate the theoretical developments in order to maximize impact and practical applicability of the approach. Finally, numerical results corroborate the theory and show the effectiveness of the new schemes in real case scenarios.

A DC Stable and Large-Time Step Well-Balanced TD-EFIE Based on Quasi-Helmholtz Projectors

The marching-on-in-time (MOT) solution of the time-domain electric field integral equation (TD-EFIE) has traditionally suffered from a number of issues, including the emergence of spurious static currents (dc instability) and ill-conditioning at large-time steps (low frequencies). In this contribution, a space-time Galerkin discretization of the TD-EFIE is proposed, which separates the loop and star components of both the equation and the unknown. Judiciously integrating or differentiating these components with respect to time leads to an equation which is free from dc instability. By choosing the correct temporal basis and testing functions for each of the components, a stable MOT system is obtained. Furthermore, the scaling of these basis and testing functions ensure that the system remains well conditioned for large-time steps. The loop–star decomposition is performed using quasi-Helmholtz projectors to avoid the explicit transformation to the unstable bases of loops and stars (or trees), and to avoid the search for global loops, which is a computationally expensive operation.