Indirect Galerkin–Bubnov boundary element method for solving integral equations in electromagnetics (original) (raw)
Related papers
Application of the boundary elements method to 3-D radiation problems
1993
In order to analyze the electromagnetic fields generated by perfectly conducting bodies composed by interconnected wires and solid objects, whichever its geometry, the Boundary Element Method is applied. The basic formulation and the associated items are commented. Results of problems of radiation are given and compared with available analytic solutions.
International Journal for Research in Applied Science and Engineering Technology -IJRASET, 2020
This work is dedicated to the numerical results and the implementation of the method coupling a discontinuous Galerkin with an integral representation (CDGIR). The originality of this work lies in the choice of discretization by discontinuous Galerkin element and a mixed form for Maxwell's equations. The numerical tests justify the effectiveness of the proposed approach. Keywords: Finite element method, Maxwell equations, Discontinuous Galerkin method, fictitious domain, integral representation, time-harmonic. I. INTRODUCTION Mathematically, the phenomenon of the electromagnetic waves propagation is generally modeled by the system of equations known as the Maxwell equations. There are two modes of the Maxwell equations to be treated, a first mode that is known by the time domain Maxwell equations in which the evolution of electromagnetic fields is studied as a function of time and the second mode that is known by the frequency domain Maxwell equations where one studies the behavior of electromagnetic fields when the source term follows a harmonic dependence in time. Numerical modeling has become the most important and widely used tool in various fields such as scientific research. The finite-difference methods (FDM), the finite element methods (FEM) and the finite volume methods (FVM) are the three classes of methods known for the numerical resolution of the problems of electromagnetic waves propagation. In 1966, Yee cited the first efficient method in [42] which is the finite-difference methods in the time domain (FDMTD). When diffraction problems are posed in unbounded domain, the use of these methods induces a problem. In order to solve it, two techniques are used. The first consists in reducing to a bounded domain by truncating the computational domain, then it is necessary to impose an artificial condition on the boundary on the truncation boundary. The second technique consists in writing an equivalent problem posed on the boundary of the obstacle, it is therefore what is called the theory of integral equations. The numerical resolution can then be done by discretizing the problem by collocation (method of moments, method of singularity) or by a finite element discretization of the boundary. In 1980, Nedelec introduces the edge finite element method developed in [31] which is also available in [29, 30]. With the conservation of energy, this method also possesses several advantages; it allows to treat unstructured meshes (complex geometries) as it can be used with high orders (see [41, 24, 29]). In recent years, research has revealed a new technique known as Discontinuous Galerkin Methods (GDM); this strategy is based on combining the advantages of FEM and FVM methods since it approaches the field in each cell by a local basis of functions by treating the discontinuity between neighboring cells by approximation FVM on the flows. Initially, these methods have been proposed to treat the scalar equation of neutron transport (see [35]). In the field of wave propagation, precisely for the resolution of the Maxwell equations in the time domain, many schemes are based on two forms of formulations: a concentrated flux formulation (see [16, 34]) and an upwind flux formulation (see [22, 12]). Discontinuous Galerkin methods have shown their effectiveness in studying the problem with discrete eigenvalues (see [23]). In frequency domain, for the resolution of Maxwell equations, the majority consider the second order formulation (see [25, 32, 33]), as others study the formulation of the first order as in [6, 20]. This strategy of the CDGIR method allows us to write a problem in an unbounded domain into an equivalent problem in a domain bounded by a fictitious boundary where a transparent condition is imposed. This transparent condition is based on the use of the integral form of the electric and magnetic fields using the Stratton-chu formulas (see [7]). This process has been studied, in the
International Journal of Numerical Modelling: Electronic Networks, Devices and Fields, 2020
In the framework solving electromagnetic scattering problems with a Leontovich impedance boundary condition (LIBC), the numerical resolution requires the evaluation of singular integrals appearing in the discretization of the variational formulation. Our main interest is to pricisely evaluate these integrals. Thus, we propose an analytic method to approximate it. The performance of this method will be evaluated by calculating the radar cross section (RCS). Then, we compare RCS to the Mie series solution for the unit sphere in various configurations.
2010
Boundary element methods (BEMs) are an increasingly popular approach to the modeling of electromagnetic scattering both by perfect conductors and dielectric objects. Several mathematical, numerical, and computational techniques pullulated from the research into BEMs, enhancing its efficiency. The Fast Multipole Method (FMM) and its descendants accelerate the matrix-vector product that constitutes the BEM's computational bottleneck. In particular, dedicated FMMs have been conceived for the computation of the electromagnetic scattering at complex metallic and/or dielectric objects in free space and in layered background media. Calderón preconditioning of the BEM's system matrix lowers the number of matrix-vector products required to reach an accurate solution, and thus the time to reach it. Parallelization distributes the remaining workload over a battery of affordable computational nodes, diminishing the wall-clock computation time. In honor of our former colleague and mentor, Prof. F. Olyslager, an overview of some dedicated BEMs for large and complex EM problems developed within the Electromagnetics Group at Ghent University is presented. Recent results that ramified from Prof. Olyslager's scientific endeavors are included in the survey.
IEEE Transactions on Antennas and Propagation, 2013
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.
2022 Photonics & Electromagnetics Research Symposium (PIERS)
In this paper, we implement the high order boundary condition (HOIBC) in the integral equations to solve the scattering problem. We study the HOIBC operators that appear in the boundary integral equations and we propose a discretization of these operators in spaces H(div) and H(curl). Then, we give validations of the HOIBC by calculating the radar cross section (RCS), that show the improvement in accuracy over the standard impedance boundary condition (SIBC). The computation code implements the H-matrix approach with Adaptive Cross Approximation (ACA) and MPI parallelization.
Surface integral equation methods for multi-scale electromagnetic problems
2014 XXXIth URSI General Assembly and Scientific Symposium (URSI GASS), 2014
This work investigates the efficient and robust integral equation based solution of large multi-scale electromagnetic problems. The major technical ingredients in the proposed work include: (i) a scalable domain decomposition method for surface integral equations via a novel multi-trace formulation, (ii) a discontinuous Galerkin boundary element method, which employs discontinuous trial and testing functions without continuity requirements across element boundaries, and (iii) an optimized multiplicative Schwarz algorithm using complete second order transmission condition. The results obtained through this research greatly simplify the model preparation and mesh generation for complex electromagnetic simulation. Moreover, It provide an effective preconditioning scheme for the integral equation based solution of multi-scale problems. The strength and flexibility of the proposed method will be illustrated by means of several challenge real-world applications.
Boundary Elements and Other Mesh Reduction Methods XXXIII, 2011
The paper deals with the space-time Galerkin-Bubnov scheme of the Indirect Boundary Element Method for the solution of time domain integral equations for thin wires in the presence of lossy media of the Pocklington and Hallen type, respectively. The presence of a dissipative half-space is taken into account via the corresponding space-time reflection coefficients. Some illustrative computational examples related to the overhead wires and grounding electrodes are presented.
A novel integral formulation for the solution of maxwell equations
IEEE Transactions on Magnetics, 2003
This paper presents a novel integral formulation for the solution of the full Maxwell equations in the frequency domain. The solenoidality of the current density inside a homogeneous conductor is exploited to introduce a two-component electric vector potential, which is numerically expanded in terms of edge elements.