Electric Potential and Field Calculation of Charged Bem Triangles and Rectangles by Gaussian Cubature

Abstract

Abstract-It is a widely held view that analytical integration is more accurate than the numerical one. In some special cases, however, numerical integration can be more advantageous than analytical integration. In our paper we show this benefit for the case of electric potential and field computation of charged triangles and rectangles applied in the boundary element method (BEM). Analytical potential and field formulas are rather complicated (even in the simplest case of constant charge densities), they have usually large computation times, and at field points far from the elements they suffer from large rounding errors. On the other hand, Gaussian cubature, which is an efficient numerical integration method, yields simple and fast potential and field formulas that are very accurate far from the elements. The simplicity of the method is demonstrated by the physical picture: the triangles and rectangles with their continuous charge distributions are replaced by discrete point charges, whose simple potential and field formulas explain the higher accuracy and speed of this method. We implemented the Gaussian cubature method for the purpose of BEM computations both with CPU and GPU, and we compare its performance with two different analytical integration methods. The ten different Gaussian cubature formulas presented in our paper can be used for arbitrary high-precision and fast integrations over triangles and rectangles.

Towards model-driven reconstruction in atom probe tomography

“…where F ij and G ij are elements of the the k × k problem matrices, y defines a point on the boundary element Γ j , and nj is this boundary element's normal. The matrix elements defined in equation (10) and equation (11) are calculated by numerically integrating [22] over surface panels (see figure 2(a)), except for diagonal matrix elements for which slower analytical solutions are required to treat singularities. Unlike in 2D [17], the exact form of these integrals is complicated and can be found in [23].…”

“…This starts at the leaf nodes of the tree and moves upward until reaching the root node. All series expansion integrals, given in appendix A, are solved numerically [22]. Such an operation has a combined computational complexity of O(k • log(k)), where k is the total number of surface boundary elements.…”

Improved determination of the $β$-$\overlineν_e$ angular correlation coefficient $a$ in free neutron decay with the $a$SPECT spectrometer

“…At that point the applied currents, voltages (see Table I) as well as the measured work functions of the particular electrode segments have to be set as input parameters. The different methods used for charge density and field calculation aredescribed in [60][61][62]. The calculated fields together with the geometrical arrangement are then used for the particle tracking, performed with the KASSIOPEIA package [57].…”

Electric Field Mitigation Near Over Head Power Lines Using Voltage Unbalance

“…When designing a transmission line this limit is not crossed. Also, proper care has been taken in order to keep minimum clearance between transmission lines [5] .…”

“…Thus this method is very accurate at positions far from the elements. However, for positions very close to the elements, the Gaussian cubature method is not precise enough, therefore analytical integration has to be used there [5,7,8].…”