Discrete and Continuous Models and Applied Computational ScienceDiscrete 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)8383ArticlesСтатьиResearch ArticleThe Volume Integral Equations Method in Magnetostatics ProblemsМетод объёмных интегральных уравнений в задачах магнитостатикиAkishinP GАкишинПавел ГригорьевичLaboratory of Information TechnologiesЛаборатория информационных технологийakishin@jinr.ruSapozhnikovA AСапожниковАндрей АлександровичLaboratory of Information TechnologiesЛаборатория информационных технологийasap@jinr.ruJoint Institute for Nuclear ResearchОбъединённый институт ядерных исследований150220142NO2 (2014)№2 (2014)31031508092016Copyright ©; 2014, Акишин П.Г., Сапожников А.А.2014Акишин П.Г., Сапожников А.А.http://creativecommons.org/licenses/by/4.0https://journals.rudn.ru/miph/article/view/8383In this article, use of volume integral equations method for calculations of magnetic systems is considered. GFUN program package based on integral approach to magnetostatics applies the method of collocations and piecewise constant approximations of unknown variables within the elements to discretization of equations. The limitation of this approach is related to singularity of the integral equations kernel. An alternative to collocation method, replacing observation point by integration over discretization elements, is considered in this article. This approach enables one to use higher order approximations for unknown variables. In the context of the finite element method, the piecewise constant and piecewise linear approximations of unknown variables are considered. The problem of computing matrix elements for discretized systems of equations can be reduced to evaluation of sixth-order integrals, singular ones in the general case, over two different elements of the computational region. Possible methods are proposed for calculating this kind of integrals. Iterative processes for solving the arising nonlinear systems of discretized equation are discussed. The proposed approach enables one to build discretizations with higher precision of approximation for the initial volume integral equations of magnetostatics. The proposed method was used for 3D modeling of a dipole magnet. Comparison of results obtained for simulation of the dipole magnet using different versions of integral magnetostatics problem discretization are given.В работе рассматриваются вопросы применения метода объёмных интегральных уравнений для расчёта магнитных систем. В пакете программ GFUN, использующем интегральную постановку задачи магнитостатики, для дискретизации уравнений используются метод коллокаций и кусочно-постоянная аппроксимация неизвестных в пределах элемента разбиения области. Ограниченность применимости данного подхода связана с сингулярностью ядра интегральных уравнений. В данной работе для дискретизации рассматривается альтернативный методу коллокаций подход, в котором точка наблюдения заменяется интегрированием по элементу разбиения. Это позволяет использовать приближения для неизвестных более высокого порядка. В рамках метода конечных элементов рассматриваются кусочно-постоянная и кусочно-линейная аппроксимации неизвестных. Проблема вычисления матричных элементов дискретизованных систем уравнений сводится к вычислению шестикратных, в общем случае сингулярных интегралов по двум различным элементам разбиения расчётной области. Предлагаются методы вычисления интегралов подобного типа. Обсуждаются итеративные методы решения возникающих нелинейных дискретизованных систем уравнений. Данный подход позволяет построить дискретизации исходных объёмных интегральных уравнений магнитостатики с более высоким порядком аппроксимации. Предлагаемый метод использовался для трёхмерного моделирования дипольного магнита. В работе приводится сравнение результатов моделирования дипольного магнита с использованием различных вариантов дискретизации интегральной постановки задачи магнитостатики.magnetostaticsvolume integral equations methodthe finite element methoddiscretizationnonlinear problemiterative processмагнитостатикаметод объёмных интегральных уравненийметод конечных элементовдискретизациянелинейная проблемаитерационный процесс1.Armstrong A. G. et al., 1976. - GFUN3D User Guide. - RL-76-029/A.2.Акишин П. Г. Метод интегральных уравнений в задачах магнитостатики: Автореф. дис. канд. физ.-мат. наук. - Дубна: ОИЯИ, 1983. - 11-83-558.3.Pasciak J. E. An Iterative Algorithm for the Volume Integral Method for Magnetostatics Problems // Comp. Maths. with Appls. - 1982. - Vol. 8, No 8. - Pp. 283-290.4.Meijerink J. A., van der Vorst H. A. An Iterative Solution Method for Linear Systems, of Which the Coefficient Matrix is a Symmetric M-matrix // Math. Comput. - 1977. - Vol. 31, No 137. - Pp. 148-162.