Discrete and Continuous Models and Applied Computational ScienceDiscrete and Continuous Models and Applied Computational Science2658-46702658-7149Peoples' Friendship University of Russia named after Patrice Lumumba (RUDN University)2219710.22363/2658-4670-2019-27-1-60-69Research ArticleThe volume integral equation method in magnetostatic problemAkishinPavel G<p>Doctor of Physical and Mathematical Sciences, Deputy Head of Scientific Department of Computational Physics, Laboratory of Information Technologies</p>akishin@jinr.ruSapozhnikovAndrey A<p>Junior Researcher of Scientific Department of Computational Physics, Laboratory of Information Technologies</p>asap@jinr.ruJoint Institute for Nuclear Research15122019271606920112019Copyright © 2019, Akishin P.G., Sapozhnikov A.A.2019<p>This article addresses the issues of volume integral equation method application to magnetic system calculations. The main advantage of this approach is that in this case finding the solution of equations is reduced to the area filled with ferromagnetic. The difficulty of applying the method is connected with kernel singularity of integral equations. For this reason in collocation method only piecewise constant approximation of unknown variables is used within the limits of fragmentation elements inside the famous package GFUN3D. As an alternative approach the points of observation can be replaced by integration over fragmentation element, which allows to use approximation of unknown variables of a higher order.In the presented work the main aspects of applying this approach to magnetic systems modelling are discussed on the example of linear approximation of unknown variables: discretisation of initial equations, decomposition of the calculation area to elements, calculation of discretised system matrix elements, solving the resulting nonlinear equation system. In the framework of finite element method the calculation area is divided into a set of tetrahedrons. At the beginning the initial area is approximated by a combination of macro-blocks with a previously constructed two-dimensional mesh at their borders. After that for each macro-block separately the procedure of tetrahedron mesh construction is performed. While calculating matrix elements sixfold integrals over two tetrahedra are reduced to a combination of fourfold integrals over triangles, which are calculated using cubature formulas. Reduction of singular integrals to the combination of the regular integrals is proposed with the methods based on the concept of homogeneous functions. Simple iteration methods are used to solve non-linear discretized systems, allowing to avoid reversing large-scale matrixes. The results of the modelling are compared with the calculations obtained using other methods.</p>finite element methodmagnetostaticsvolume integral equationssystems of nonlinear equationscubature formulaeiterative processметод конечных элементов, магнитостатика, объемные интегральные уравнения, системы нелинейных уравнений, кубатурные формулы, итерационный процесс[C. W. Trowbridge, J. K. Sykulski, Some key developments in computational electromagnetics and their attribution, IEEE transactions on magnetics 42 (4) (2006) 503-508. doi:10.1109/TMAG.2006.872491.][A. G. Armstrong, GFUN3D User Guide. RL-76-029/A, CERN (1976).][A. A. Halacsy, Three-dimensional analysis of magnetic fields, in: Proceedings 3rd International Conference on Magnet Technology, Hamburg, 1970, pp. 113-128.][M. J. Friedman, Mathematical study of the nonlinear singular integral magnetic field equation I, SIAM Journal on Applied Mathematics 39 (1) (1980) 14-20. doi:10.1137/0139003.][M. J. Friedman, Mathematical study of the nonlinear singular integral magnetic field equation II, SIAM Journal on Numerical Analysis 18 (4) (1981) 644-653. doi:10.1137/0718042.][J. Simkin, C. W. Trowbridge, Three dimensional non-linear electro-magnetic field computations using scalar potentials, IEE Proceedings B - Electric Power Applications 127 (6) (1990) 368-374. doi:10.1049/ip-b.1980.0052.][J. A. Stratton, Electromagnetic theory, MCgraw-hill, 1941.][J. D. Jackson, Classical electrodynamics, 2nd Edition, John Wiley & Sons, 1975.][O. C. Zienkiewicz, The finite element method in engineering science, MCgraw-hill, 1971.][J. T. Oden, Finite elements of nonlinear continua, McGraw-Hill, New York, 1971.][W. G. Strang, G. J. Fix, An analysis of the finite element method, Prentice-Hall, 1973.][J. P. Aubin, Approximation of elliptic boundary-value problems, Wiley-Interscience, 1972.][P. G. Akishin, A. A. Sapozhnikov, Linear approximation of volume integral equations for the problem of magnetostatics, in: EPJ Web Conferences. Mathematical Modeling and Computational Physics (MMCP 2017), Vol. 173, 2018, Article Number 03001. doi:10.1051/epjconf/201817303001.][Z. Cendes, Magnetic field computation using Delaunay triangulation and complementary finite element methods, IEEE Transactions on Magnetics 19 (1983) 2551-2554. doi:10.1109/TMAG.1983.1062841.][E. K. Buratynski, A three-dimensional unstructured mesh generator for arbitrary internal boundaries, in: Numerical Grid Generation in Computational Fluid Mechanics: Proceedings, Pineridge Press, Swansea, 1988, pp. 621-631.][M. Berzins, Mesh quality: a function of geometry, error estimates or both?, Engineering with Computers 15 (1999) 236-247. doi:10.1007/s003660050019.][L. Durbeck, Evaporation: a technique for visualizing mesh quality, in: 8th International Meshing Roundtable, Sandia National Laboratories, South Lake Tahoe, 1999, pp. 259-265.][P. G. Akishin, A. A. Sapozhnikov, Automatic generation of three-dimensional grids, JINR, Dubna, 2015. URL http://www1.jinr.ru/Preprints/2015/058(P11-2015-58).pdf][P. G. Akishin, A. A. Sapozhnikov, 3DFEMMesh - program for automatic generation of three-dimensional Mesh. URL http://wwwinfo.jinr.ru/programs/jinrlib/3dfemmesh/ indexe.html][M. Abramowitz, I. Stegun, Handbook of mathematical functions with functions, graphs, and mathematical tables, 1964.][I. P. Mysovskikh, Interpolation cubature formulas [Interpolyatsionnyye kubaturnyye formuly], Nauka, Moscow, 1981, in Russian.][P. G. Akishin, The integral equation method in magnetostatic problems: abstract of a PhD thesis, Ph.D. thesis, JINR, Dubna, in Russian (1983).][J. A. Meijerink, H. A. van der Vorst, An iterative solution method for linear systems of which the coefficient matrix is a symmetric M-matrix, Mathematics of Computation 31 (137) (1977) 148-162. doi:10.2307/2005786.]