Galerkin Surface Integral Equation Method Research Articles

An Efficient and Effective Preconditioner for the Discontinuous Galerkin Domain Decomposition Method of Surface Integral Equation

This paper proposes a new preconditioning technique based on a restricted additive Schwarz (RAS) approach to improve the efficiency of the discontinuous Galerkin surface integral equation method. The RAS preconditioning is implemented efficiently using an oct-tree structure derived from the multilevel fast multipole algorithm (MLFMA) and is constructed using near-field matrices associated with boxes at the finest level. Compared to existing domain decomposition preconditioning methods based on subdomain block matrices, the proposed RAS preconditioning significantly reduces computation time and memory requirements, while providing scalable convergence for iterative solutions. Numerical experiments are presented to demonstrate the performance of the proposed cost-effective preconditioning technique.

An Efficient and Error-Controllable Angular Sweep Algorithm for Electromagnetic Wave Scattering and its Analysis

The angular sweep of electromagnetic wave scattering is formulated as a matrix equation with multiple right-hand sides (RHSs). Although the low-rank approximation of an RHS matrix is a popular choice for reducing the computational costs of multiple RHSs, only a small amount of research has been conducted to explore how this approximation impacts the solution quality. Furthermore, there has not been sufficient research on the quality of the solution as a function of the accuracy of the iterative solver. We present an error analysis of the approximated solution considering both the reduced number of RHSs and the tolerance of the iterative solver. Based on the error analysis, a new angular sweep algorithm is proposed with fine-tuned tolerances of the iterative solver for individual singular vectors. The different tolerances for each singular vector increase the efficiency of the proposed algorithm. Another benefit of the proposed algorithm is that the error can be bounded by a user-defined global tolerance. In addition, a variant of the generalized conjugate residual method for multiple RHSs is introduced to accelerate iterative solvers. Finally, numerical validation is conducted with three examples in which the discontinuous Galerkin surface integral equation method is applied. The experiments support two conclusions: tight upper and lower bounds of the solution error exist, and fine-tuning the tolerances reduces the computational costs.

Massively Parallel Discontinuous Galerkin Surface Integral Equation Method for Solving Large-Scale Electromagnetic Scattering Problems

In this communication, we present flexible and efficient solutions of large-scale electromagnetic scattering problems using the parallel discontinuous Galerkin (DG) surface integral equation (SIE) method. A simplified formulation of the conventional DG is used by removing the interior penalty term, thereby avoiding computational burden of the contour integral in the conventional DG. It is demonstrated numerically that the simplification causes convergence deficiency and changes the accuracy convergence rate of the solution for objects with sharp corners. The former is overcome by constructing preconditioners using the near-field matrix of the whole region, and the latter is proved to have a neglectable influence on the accuracy of the final solutions as long as normal mesh density requirement is satisfied. The well-scaling ternary parallelization approach of the multilevel fast multipole algorithm (MLFMA) is incorporated into the DG domain decomposition framework to reduce the complexity of the solution and strength its capability for large-scale problems, with the distribution of the near-field-related calculation adjusted to adapt the nonuniform meshes for a better workload balance. Numerical results are included to validate the accuracy and demonstrate the scalability and versatility of the proposed method. In addition, we demonstrated the effectiveness of our algorithm by solving a high-definition complicated aircraft model with inlets and exhausts involving over one billion unknowns.

A discontinuous Galerkin surface integral equation method for scattering from IBC targets

AbstractA discontinuous Galerkin (DG) surface integral equation method is proposed for electromagnetic scattering from targets with the impedance boundary condition (IBC). We present electric field integral equation (EFIE), magnetic field integral equation (MFIE), and self‐dual integral equation formulations for the problem and study their numerical performance. On the basis of the results of these experiments, DG is developed for EFIE and MFIE (DG‐EFIE and DG‐MFIE). The convergence of the iterative solutions and the solution accuracy of DG‐EFIE and DG‐MFIE are further investigated for a wide range of surface impedances. The numerical performance of these formulations is found to be nearly complementary with respect to the surface impedance. The capability of the proposed DG solution strategy for calculating scattering from large multiscale IBC targets is demonstrated.

A Discontinuous Galerkin Surface Integral Equation Method for Scattering From Multiscale Homogeneous Objects

A discontinuous Galerkin (DG) surface integral equation approach is proposed for scattering from homogeneous dielectric objects. The formulation of DG for homogeneous bodies is derived from the combined tangential field integral equation. The differences of DG for penetrable and nonpenetrable objects are presented by numerical experiments to demonstrate the numerical mechanism of DG. Numerical experiments demonstrate the great advantages of our presented formulation of DG for homogeneous objects in efficiency, flexibility, and scalability. A series numerical results are presented to show the capability of the presented DG solution for homogeneous bodies, especially for multiscale homogeneous bodies.

A Discontinuous Galerkin Surface Integral Equation Method for Electromagnetic Wave Scattering From Nonpenetrable Targets

We present a discontinuous Galerkin surface integral equation method, herein referred to as IEDG, for time harmonic electromagnetic wave scattering from nonpenetrable targets. The proposed IEDG algorithm allows the implementation of the combined field integral equation (CFIE) using square-integrable, , trial and test functions without any considerations of continuity requirements across element boundaries. Due to the local characteristics of basis functions, it is possible to employ nonconformal surface discretizations of the targets. Furthermore, it enables the possibility to mix different types of elements and employ different order of basis functions within the same discretization. Therefore, the proposed IEDG method is highly flexible to apply adaptation techniques. Numerical results are included to validate the accuracy and demonstrate the versatility of the proposed IEDG method. In addition, a complex large-scale simulation is conducted to illustrate the potential benefits offered by the proposed method for modeling multiscale electrically large targets.