Search

We propose here a novel stabilization strategy for the PMCHWT equation that cures its frequency and conductivity related instabilities and is obtained by leveraging quasi-Helmholtz projectors. The resulting formulation is well-conditioned in the entire low-frequency regime, including the eddy current one, and can be applied to arbitrarily penetrable materials, ranging from dielectric to conductive ones. In addition, by choosing the rescaling coefficients of the quasi-Helmholtz components appropriately, we prevent the typical loss of accuracy occurring at low frequency in the presence of inductive and capacitive type magnetic frill excitations, commonly used in circuit modeling to impose a potential difference. Finally, leveraging on quasi-Helmholtz projectors instead than on the standard Loop-Star decomposition, our formulation is also compatible with most fast solvers and is amenable to multiply connected geometries, without any computational overhead for the search for the global loops of the structure. The efficacy of the proposed preconditioning scheme when applied to both simply and multiply connected geometries is corroborated by numerical examples.

In this contribution we analyze the spectral properties of some commonly used boundary integral operators in computational electromagnetics and of their discrete counterparts, highlighting peculiar features of their spectra. In particular, a comparison with the eigenvalues of the continuous operators will be presented that highlights deviations in the high frequency regime and impacts, in a peculiar way, the accuracy of the numerical solutions of each formulation. A study and a proactive analysis of numerical results from standard boundary element solvers and the predictions from the theoretical analysis will corroborate the analytical framework employed and the validity of our observations.

Surface integral equations (SIEs)-based boundary element methods are widely used for analyzing electromagnetic scattering scenarii. However, after discretization of SIEs, the spectrum and eigenvectors of the boundary element matrices are not usually representative of the spectrum and eigenfunctions of the underlying surface integral operators, which can be problematic for methods that rely heavily on spectral properties. To address this issue, we delineate some efficient algorithms that allow for the computation of matrix square roots and inverse square roots of the Gram matrices corresponding to the discretization scheme, which can be used for revealing the spectrum of standard electromagnetic integral operators. The algorithms, which are based on properly chosen expansions of the square root and inverse square root functions, are quite effective when applied to several of the most relevant Gram matrices used for boundary element discretizations in electromagnetics. Tables containing different sets of expansion coefficients are provided along with comparative numerical experiments that evidence advantages and disadvantages of the different approaches. In addition, to demonstrate the spectrum-revealing properties of the proposed techniques, they are applied to the discretization of the problem of scattering by a sphere for which the analytic spectrum is known., Comment: 8 pages, 9 figures

This paper introduces a time-domain combined field integral equation for electromagnetic scattering by a perfect electric conductor. The new equation is obtained by leveraging the quasi-Helmholtz projectors, which separate both the unknown and the source fields into solenoidal and irrotational components. These two components are then appropriately rescaled to cure the solution from a loss of accuracy occurring when the time step is large. Yukawa-type integral operators of a purely imaginary wave number are also used as a Calderon preconditioner to eliminate the ill-conditioning of matrix systems. The stabilized time-domain electric and magnetic field integral equations are linearly combined in a Calderon-like fashion, then temporally discretized using an appropriate pair of trial functions, resulting in a marching-on-in-time linear system. The novel formulation is immune to spurious resonances, dense discretization breakdown, large-time step breakdown and dc instabilities stemming from non-trivial kernels. Numerical results for both simply-connected and multiply-connected scatterers corroborate the theoretical analysis., Comment: 13 pages

This work presents a fast direct solver strategy allowing full-wave modeling and dosimetry at terahertz (THz) frequencies. The novel scheme leverages a preconditioned combined field integral equation together with a regularizer for its elliptic spectrum to enable its compression into a non-hierarchical skeleton, invertible in quasi-linear complexity. Numerical results will show the effectiveness of the new scheme in a realistic skin modeling scenario.

A fully kinetic 2D axisymmetric Particle-in-Cell (PIC) model is used to examine the effects of background facility pressure on the plasma transport and propulsive efficiency of magnetic nozzles. Simulations are performed for a low-power (150 W class) cathode-less radio-frequency (RF) plasma thruster, operating with xenon, between background pressures up to 10$^{-2}$ Pa and average electron discharge temperatures of 4 - 16 eV. When the electron temperature within the near-plume region reaches 8 eV, a decisive reduction in performance occurs: at 10$^{-2}$ Pa, in-plume power losses surpass 25% of the discharge energy flux. Given that the ionization energy for Xe is 12 eV, the 8 eV threshold indicates that a consistent percentage of electrons has energy enough to trigger ionization. On the other hand, when the temperature is below such threshold, the primary collisions are charge-exchange and inelastic ion scattering, and the power loss remains less than 10%. It is established that losses in the considered HPT are significant if the facility pressure is greater than 10$^{-3}$ Pa, at absorbed powers larger than 130 W. At the nominal 150 W, this results in a 15% thrust reduction. When facility pressure is taken into consideration over ideal vacuum simulations, numerical error is reduced to <30% when compared to experimental thrust measurements at 10$^{-3}$ Pa., Comment: 45 pages, 10 figures, 5 tables. Submitted to Acta Astronautica

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 non-solenoidal contributions, then rescaled in order 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 paper 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.

Recent contributions showed the benefits of operator filtering for both preconditioning and fast solution strategies. While previous contributions leveraged laplacian-based filters, in this work we introduce and study a different approach leveraging the truncation of appropriately chosen spectral representations of operators' kernels. In this contribution, the technique is applied to the operators of the 2D TE- and TM-electric field integral equations (EFIE). We explore two different spectral representations for the 2D Green's function that lead to two distinct types of filtering of the EFIE operators. Numerical results corroborate the effectiveness of the newly proposed approaches, also in the Calder\'on preconditioned EFIE, Comment: 3 pages, 3 figures, to be published in ICEAA 2023

Quasi-Helmholtz decompositions are fundamental tools in integral equation modeling of electromagnetic problems because of their ability of rescaling solenoidal and non-solenoidal components of solutions, operator matrices, and radiated fields. These tools are however incapable, per se, of modifying the refinement-dependent spectral behavior of the different operators and often need to be combined with other preconditioning strategies. This paper introduces the new concept of filtered quasi-Helmholtz decompositions proposing them in two incarnations: the filtered Loop-Star functions and the quasi-Helmholtz Laplacian filters. Because they are capable of manipulating large parts of the operators' spectra, new families of preconditioners and fast solvers can be derived from these new tools. A first application to the case of the frequency and h-refinement preconditioning of the electric field integral equation is presented together with numerical results showing the practical effectiveness of the newly proposed decompositions.

We present a Calder\'on preconditioning scheme for the symmetric formulation of the forward electroencephalographic (EEG) problem that cures both the dense discretization and the high-contrast breakdown. Unlike existing Calder\'on schemes presented for the EEG problem, it is refinement-free, that is, the electrostatic integral operators are not discretized with basis functions defined on the barycentrically-refined dual mesh. In fact, in the preconditioner, we reuse the original system matrix thus reducing computational burden. Moreover, the proposed formulation gives rise to a symmetric, positive-definite system of linear equations, which allows the application of the conjugate gradient method, an iterative method that exhibits a smaller computational cost compared to other Krylov subspace methods applicable to non-symmetric problems. Numerical results corroborate the theoretical analysis and attest of the efficacy of the proposed preconditioning technique on both canonical and realistic scenarios.

This paper focuses on fast direct solvers for integral equations in the low-to-moderate-frequency regime obtained by leveraging preconditioned first kind or second kind operators regularized with Laplacian filters. The spectral errors arising from boundary element discretizations are properly handled by filtering that, in addition, allows for the use of low-rank representations for the compact perturbations of all operators involved. Numerical results show the effectiveness of the approaches and their effectiveness in the direct solution of integral equations.

The inverse source problem in electromagnetics has proved quite relevant for a large class of applications. In antenna diagnostics in particular, Love solutions are often sought at the cost of an increase of the dimension of the linear system to be solved. In this work, instead, we present a reduced-in-size single current formulation of the inverse source problem that obtains one of the Love currents via a stable discretization of the Steklov-Poincar\'e boundary operator leveraging dual functions. The new approach is enriched by theoretical treatments and by a further low-frequency stabilization of the Steklov-Poincar\'e operator based on the quasi-Helmholtz projectors that is the first of its kind in this field. The effectiveness and practical relevance of the new schemes are demonstrated via both theoretical and numerical results.

Inverse source approaches have shown their relevance for several applications in the past years. They rely on the solution of an ill-posed problem where near-field/current data is reconstructed starting from far-field (or less informative field) information. Standard strategies, including the physically constrained ones using Love conditions, result in linear systems to be pseudoinverted which are still ill-conditioned due to the lack of information from the evanescent fields. In this work we present a generalized pseudoinverse for those problems that allows the inclusion of extra constraints from the evanescent field space, when available. This is obtained by dropping some of the standard Moore-Penrose (MP) requirements and using the resulting degrees of freedom to obtain a generalized pseudoinverse that shows favorable performance in several cases of practical interest.

Brain-Computer Interfaces (BCIs) based on Steady State Visually Evoked Potentials (SSVEPs) have proven effective and provide significant accuracy and information-transfer rates. This family of strategies, however, requires external devices that provide the frequency stimuli required by the technique. This limits the scenarios in which they can be applied, especially when compared to other BCI approaches. In this work, we have investigated the possibility of obtaining frequency responses in the EEG output based on the pure visual imagination of SSVEP-eliciting stimuli. Our results show that not only that EEG signals present frequency-specific peaks related to the frequency the user is focusing on, but also that promising classification accuracy can be achieved, paving the way for a robust and reliable visual imagery BCI modality.

This work focuses on a quasi-linear-in-complexity strategy for a hybrid surface-wire integral equation solver for the electroencephalography forward problem. The scheme exploits a block diagonally dominant structure of the wire self block -- that models the neuronal fibers self interactions -- and of the surface self block -- modeling interface potentials. This structure leads to two Neumann iteration schemes further accelerated with adaptive integral methods. The resulting algorithm is linear up to logarithmic factors. Numerical results confirm the performance of the method in biomedically relevant scenarios.

Several strategies are available for solving the inverse source problem in electromagnetics. Among them, many have been focusing in retrieving Love currents by solving, after regularization, for Love's electric and magnetic currents. In this work we present a dual-element discretization, analysis, and stabilization of an inverse source formulation providing Love data by solving for only one current. This results in substantial savings and allows for an effective quasi-Helmholtz projector stabilization of the resulting operator. Theoretical considerations are complemented by numerical tests showing effectiveness and efficiency of the newly proposed method.

Widely employed for the accurate solution of the electroencephalography forward problem, the symmetric formulation gives rise to a first kind, ill-conditioned operator ill-suited for complex modelling scenarios. This work presents a novel preconditioning strategy based on an accurate spectral analysis of the operators involved which, differently from other Calder\'on-based approaches, does not necessitate the barycentric refinement of the primal mesh (i.e., no dual matrix is required). The discretization of the new formulation gives rise to a well-conditioned, symmetric, positive-definite system matrix, which can be efficiently solved via fast iterative techniques. Numerical results for both canonical and realistic head models validate the effectiveness of the proposed formulation.

This work presents a fast direct solver strategy for electromagnetic integral equations in the high-frequency regime. The new scheme relies on a suitably preconditioned combined field formulation and results in a single skeleton form plus identity equation. This is obtained after a regularization of the elliptic spectrum through the extraction of a suitably chosen equivalent circulant problem. The inverse of the system matrix is then obtained by leveraging the Woodbury matrix identity, the low-rank representation of the extracted part of the operator, and fast circulant algebra yielding a scheme with a favorable complexity and suitable for the solution of multiple right-hand sides. Theoretical considerations are accompanied by numerical results both of which are confirming and showing the practical relevance of the newly developed scheme.

This paper extends the concept of Laplacian filtered quasi-Helmholtz decompositions we have recently introduced, to the basis-free projector-based setting. This extension allows the discrete analyses of electromagnetic integral operators spectra without passing via an explicit Loop-Star decomposition as previously done. We also present a fast scheme for the evaluation of the filters in quasi linear complexity in the total number of unknowns. Together with the fact that only a logarithmic number of these filters are required for solving the h-refinement breakdown of electric field integral equation, this results in an effective preconditioner that rivals Calder\'on strategies in performance without relying on barycentric refinements. Numerical results confirm the theoretically predicted behavior and the effectiveness of the approach.

In this paper 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, transcranial magnetic and deep brain stimulation (TMS and DBS, respectively). 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 in 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 re-balancing 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.

The Electric Field Integral Equation (EFIE) is notorious for its ill-conditioning both in frequency and h-refinement. Several techniques exist for fixing the equation conditioning problems based on hierarchical strategies, Calderon techniques, and related technologies. This work leverages on a new approach, based on the construction of tailored spectral filters for the EFIE components which allow the block renormalization of the EFIE spectrum resulting in a provably constant condition number for the equation. This is achieved without the need for a barycentric refinement and with low computational overhead compared with other schemes. In particular, only sparse matrices are required in addition to the EFIE original matrix. Numerical results will show the robustness of our scheme and its application to the solution of realistic problems.

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., Comment: 10 pages, 5 figures

The CFIE used for solving scattering and radiation problems, although a resonance-free formulation, suffers from an ill-conditioning that strongly depends on the frequency and discretization density, both in the low- and high-frequency regime, resulting in slow convergence rates for iterative solvers. This work presents a new preconditioning scheme for the CFIE that cures the low- and the high-frequency as well as the dense discretization breakdown. The new preconditioner for the CFIE is based on a spherical harmonics analysis and the proper regularization with Helmholtz-type operators. Numerical results have been obtained to prove the effectiveness of this new formulation in real scenarios.

Despite its several qualities, the Poggio-Miller-Chang-Harrington-Wu-Tsai (PMCHWT) formulation for simulating scattering by dielectric media suffers from numerical instabilities and severe ill-conditioning at low frequencies. While this drawback has been the object of numerous solution attempts in the standard low-frequency breakdown regime for scattering problems, the eddy-current regime requires a specific treatment. In this contribution, we present an extension of the recently introduced quasi-Helmholtz projectors based preconditioning of the PMCHWT to obtain an equation stable at low frequencies and specifically in three regimes relevant for eddy currents analyses: (i) when the frequency decreases with a constant conductivity, (ii) when the conductivity increases at fixed frequency and (iii) when the frequency decreases while keeping the product of the frequency and the conductivity constant. Being based on quasi-Helmholtz projectors our new strategy does not further degrade the conditioning of the original equation and is compatible with existing fast solvers. The resulting full-wave formulation is capable of smoothly transitioning from simulations at high frequencies to the different low-frequency regimes which encompass, in particular, eddy currents applications. Numerical results demonstrate the validity of our approach in all the regimes with a special emphasis given to eddy currents.

Fast and accurate resolution of electromagnetic problems via the \ac{BEM} is oftentimes challenged by conditioning issues occurring in three distinct regimes: (i) when the frequency decreases and the discretization density remains constant, (ii) when the frequency is kept constant while the discretization is refined and (iii) when the frequency increases along with the discretization density. While satisfactory remedies to the problems arising in regimes (i) and (ii), respectively based on Helmholtz decompositions and Calder\'on-like techniques have been presented, the last regime is still challenging. In fact, this last regime is plagued by both spurious resonances and ill-conditioning, the former can be tackled via combined field strategies and is not the topic of this work. In this contribution new symmetric scalar and vectorial electric type formulations that remain well-conditioned in all of the aforementioned regimes and that do not require barycentric discretization of the dense electromagnetic potential operators are presented along with a spherical harmonics analysis illustrating their key properties.

The surface Boundary Element Method (BEM) is one of the most commonly employed formulations to solve the forward problem in electroencephalography, but the applicability of its classical incarnations is lamentably limited to piece-wise homogeneous media. Several head tissues, however, are strongly anisotropic due to their complex underlying micro-structure. This implies that standard boundary integral formulations oversimplify the electrical properties of the head and produce unrealistic solutions, something that drastically limits the suitability and impact of BEM technologies to brain imaging. This contribution addresses this issue by observing that the brain anisotropy in the white matter is due to the presence of neuronal wire-like structures. We then extend the well known wire integral equations used for high frequency problems to the imperfectly conducting quasi-static case and we propose a new hybrid wire/surface/volume integral equation. When applied on multimodal magnetic resonance images combined with tractography, this new approach can flexibly and realistically handle the conductivity anisotropy in any head compartment providing high level of accuracy and efficiency. The beneficial properties of the new formulation together with its impact on brain imaging is demonstrated via numerical results on both canonical and realistic case scenarios.

Solving the electroencephalography (EEG) forward problem is a fundamental step in a wide range of applications including biomedical imaging techniques based on inverse source localization. State-of-the-art electromagnetic solvers resort to a computationally expensive volumetric discretization of the full head to account for its complex and heterogeneous electric profile. The more efficient, popular in biomedical imaging circles, but unfortunately oversimplifying Boundary Element Method (BEM) relies instead on a piecewise-uniform approximation that severely curbs its application in high resolution EEGs. This contribution lifts the standard BEM contraints by treating the local anisotropies with adequate wire and thin volume integral equations that are tailored to specific structures of the fibrous white matter and the inhomogeneous skull. The proposed hybrid integral equation formulation thereby avoids the full volumetric discretization of the head medium and allows for a realistic and efficient BEM-like solution of the anisotropic EEG forward problem. The accuracy and flexibility of the proposed formulation is demonstrated through numerical experiments involving both canonical and realistic MRI-based head models.

The electroencephalography (EEG) forward problem, the computation of the electric potential generated by a known electric current source configuration in the brain, is a key step of EEG source analysis. In this problem, it is often desired to model the anisotropic conductivity profiles of the skull and of the white matter. These profiles, however, cannot be handled by standard surface integral formulations and the use of volume finite elements is required. Leveraging on the representation theorem using an anisotropic fundamental solution, this paper proposes a modified symmetric formulation for solving the EEG forward problem by a surface integral equation which can take into account anisotropic conductivity profiles. A set of numerical results is presented to corroborate theoretical treatments and to show the impact of the proposed approach on both canonical and real case scenarios., Comment: Physics in medicine and biology (2018)

The symmetric formulation of the electroencephalography (EEG) forward problem is a well-known and widespread equation thanks to the high level of accuracy that it delivers. However, this equation is first kind in nature and gives rise to ill-conditioned problems when the discretization density or the brain conductivity contrast increases, resulting in numerical instabilities and increasingly slow solutions. This work addresses and solves this problem by proposing a new regularized symmetric formulation. The new scheme is obtained by leveraging on Calderon identities which allow to introduce a dual symmetric equation that, combined with the standard one, results in a second kind operator which is both stable and well-conditioned under all the above mentioned conditions. The new formulation presented here can be easily integrated into existing EEG imaging packages since it can be obtained with the same computational technology required by the standard symmetric formulation. The performance of the new scheme is substantiated by both theoretical developments and numerical results which corroborate the theory and show the practical impact of the new technique.

Source localization based on electroencephalography (EEG) has become a widely used neuroimagining technique. However its precision has been shown to be very dependent on how accurately the brain, head and scalp can be electrically modeled within the so-called forward problem. The construction of this model is traditionally performed by leveraging Finite Element or Boundary Element Methods (FEM or BEM). Even though the latter is more computationally efficient thanks to the smaller interaction matrices it yields and near-linear solvers, it has traditionally been used on simpler models than the former. Indeed, while FEM models taking into account the different media anisotropies are widely available, BEM models have been limited to isotropic, piecewise homogeneous models. In this work we introduce a new BEM scheme taking into account the anisotropies of the white matter. The boundary nature of the formulation allows for an efficient discretization and modelling of the fibrous nature of the white matter as one-dimensional basis functions, limiting the computational impact of their modelling. We compare our scheme against widely used formulations and establish its correctness in both canonical and realistic cases., Comment: Update the title of the paper to reflect the title of the published version

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.

Several recent contributions have envisioned the possibility of increasing currently exploitable maximum channel capacity of a free space link, both at optical and radio frequencies, by using vortex waves, i.e. carrying Orbital Angular Momentum (OAM). Our objective is to disprove these claims by showing that they are in contradiction with very fundamental properties of Maxwellian fields. We demonstrate that the Degrees of Freedom (DoF) of the field cannot be increased by the helical phase structure of electromagnetic vortex waves beyond what can be done without invoking this property. We also show that the often-advocated over-quadratic power decay of OAM beams with distance does not play any fundamental role in the determination of the channel DoF., Comment: 8 pages, 7 figures

We present a Calder\'on preconditioner for the electric field integral equation (EFIE), which does not require a barycentric refinement of the mesh and which yields a Hermitian, positive definite (HPD) system matrix allowing for the usage of the conjugate gradient (CG) solver. The resulting discrete equation system is immune to the low-frequency and the dense-discretization breakdown and, in contrast to existing Calder\'on preconditioners, no second discretization of the EFIE operator with Buffa-Christiansen (BC) functions is necessary. This preconditioner is obtained by leveraging on spectral equivalences between (scalar) integral operators, namely the single layer and the hypersingular operator known from electrostatics, on the one hand, and the Laplace-Beltrami operator on the other hand. Since our approach incorporates Helmholtz projectors, there is no search for global loops necessary and thus our method remains stable on multiply connected geometries. The numerical results demonstrate the effectiveness of this approach for both canonical and realistic (multi-scale) problems.