PIER B | |

Progress In Electromagnetics Research B | ISSN: 1937-6472 |

Home > Vol. 74 > pp. 1-21
## ELECTRIC POTENTIAL AND FIELD CALCULATION OF CHARGED BEM TRIANGLES AND RECTANGLES BY GAUSSIAN CUBATUREBy F. Gluck and D. Hilk
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.
2. Segerlind, L. J., "Applied Finite Element Analysis," 3. Cook, R. D., et al., 4. Gupta, O. P., 5. Kythe, P. K., 6. Gaul, L., M. Kogl, and M. Wagner, 7. Beer, G., I. Smith, and Ch. Duenser, 8. Brebbia, C. A., J. C. F. Telles, and L. C. Wrobel, 9. Brebbia, C. A. and J. Dominguez, 10. Corona, T. J., 11. Szilagyi, M., 12. Hawkes, P. W. and E. Kasper, 13. Kraus, Ch., et al., "Final results from phase II of the Mainz neutrino mass search in tritium β decay," 14. Lobashev, V. M., "The search for the neutrino mass by direct method in the tritium beta-decay and perspectives of study it in the project KATRIN," 15. Gluck, F., et al., "The neutron decay retardation spectrometer aSPECT: Electromagnetic design and systematic effects," 16. Angrik, J., et al., " (KATRIN Collaboration), KATRIN Design Report 2004,", Wissenschaftliche Berichte FZKA 7090, http://bibliothek.fzk.de/zb/berichte/FZKA7090.pdf. 17. Prall, M., et al., "The KATRIN pre-spectrometer at reduced filter energy," 18. Arenz, M., et al., "Commissioning of the vacuum system of the KATRIN main spectrometer," 19. Corona, T. J., 20. Furse, D., 21. Groh, S., 22. Chari, M. V. K. and S. J. Salon, 23. Ueberhuber, Ch. W., "Numerical computation," 24. Goldberg, D., "What every computer scientist should know about floating-point arithmetic," 25. Formaggio, J., et al., "Solving for micro- and macro-scale electrostatic configurations using the Robin Hood algorithm," 26. Hanninen, I., M. Taskinen, and J. Sarvas, "Singularity subtraction integral formulae for surface integral equations with rwg, rooftop and hybrid basis functions," 27. Hilk, D., 28. Rao, S. M., et al., "A simple numerical solution procedure for statics problems involving arbitraryshaped surfaces," 29. Okon, E. E. and R. F. Harrington, "The potential due to a uniform source distribution over a triangular domain," 30. Davey, K. and S. Hinduja, "Analytical integration of linear three-dimensional triangular elements in BEM," 31. Tatematsu, A., S. Hamada, and T. Takuma, "Analytical expressions of potential and electric field generated by a triangular surface charge with a high-order charge density distribution," 32. Lopez-Pena, S., A. G. Polimeridis, and J. R. Mosig, "On the analytic-numeric treatment of weakly singular integrals on arbitrary polygonal domains," 33. Carley, M. J., "Analytical formulae for potential integrals on triangles," 34. Durand, E., 35. Eupper, M., 36. Birtles, A. B., B. J. Mayo, and A. W. Bennett, "Computer technique for solving 3-dimensional electron-optics and capacitance problems," 37. Evans, G., 38. Kythe, P. K. and M. R. Schaferkotter, 39. Stroud, A. H., 40. Engels, H., 41. Engeln-Mullges, G. and F. Uhlig, "Numerical Algorithms with C," 42. Radon, J., "Zur mechanischen Kubatur," 43. Hammer, P. C. and A. H. Stroud, "Numerical integration over simplexes," 44. Hammer, P. C. and A. H. Stroud, "Numerical evaluation of multiple integrals II," 45. Hammer, P. C., O. J. Marlowe, and A. H. Stroud, "Numerical integration over simplexes and cones," 46. Gatermann, K., "The construction of symmetric cubature formulas for the square and the triangle," 47. Lyness, J. N. and D. Jespersen, "Moderate degree symmetric quadrature rules for the triangle," 48. Papanicolopulos, S., "Computation of moderate-degree fully-symmetric cubature rules on the triangle using symmetric polynomials and algebraic solving," 49. Albrect, J. and L. Collatz, "Zur numerischen auswertung mehrdimensionaler integrale," 50. Tyler, G. W., "Numerical integration of functions of several variables," 51. Moller, H. M., "Kubaturformeln mit minimaler Knotenzahl," 52. Cools, R. and A. Haegemans, "Another step forward in searching for cubature formulae with a minimal number of knots for the square," 53. Rabinowitz, P. and N. Richter, "Perfectly symmetric two-dimensional integration formulas with minimal numbers of points," 54. Omelyan, I. P. and V. B. Solovyan, "Improved cubature formulae of high degrees of exactness for the square," 55. Dunavant, D. A., "High degree efficient symmetrical gaussian quadrature rules for the triangle," 56. Haegemans, A. and R. Piessens, "Construction of cubature formulas of degree eleven for symmetric planar regions, using orthogonal polynomials," 57. Haegemans, A. and R. Piessens, "Construction of cubature formulas of degree seven and nine symmetric planar regions, using orthogonal polynomials," 58. Wandzura, S. and H. Xiao, "Symmetric quadrature rules on a triangle," 59. Zhang, L., T. Cui, and H. Liu, "A set of symmetric quadrature rules on triangles and tetrahedra," 60. Cools, R. and P. Rabinowitz, "Monomial cubature rules since Stroud: A compilation," 61. Lyness, J. N. and R. Cools, "A survey of numerical cubature over triangles," 62. Cools, R., "Monomial cubature rules since Stroud: A compilation — Part 2," 63. Cools, R., "Advances in multidimensional integration," 64. Cools, R., "An encyclopaedia of cubature formulas," 65. Cools, R., 66. Duffy, M. G., "Quadrature over a pyramid or cube of integrands with a singularity at a vertex," 67. Valerius, K., "The wire electrode system for the KATRIN main spectrometer," 68. Stern, S., 69. Combe, R., 70. Barrett, J., 71. Gosda, W., 72. Barrett, J., J. Formaggio, and T. J. Corona, "A method to calculate the spherical multipole expansion of the electrostatic charge distribution on a triangular boundary element," 73. Lazic, P. H., H. Stefancic, and H. Abraham, "The robin hood method — A novel numerical method for electrostatic problems based on a non-local charge transfer," 74. Walker, D. W. and J. J. Dongarra, "MPI: A standard message passing interface," 75. Stone, E., D. Gohara, and G. Shi, "OpenCL: A parallel programming standard for heterogeneous computing systems," 76. Kahan, W., "Further remarks on reducing truncation errors," 77. ROOT Data Analysis Framework, C++ program package, developed in CERN, https://root.cern.ch/, . 78. Hwu, W.-M. W., (Editor), |

© Copyright 2010 EMW Publishing. All Rights Reserved