Abstract

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.

Highlights

  • The paper [1] gives an overview of existing methods and programs for the numerical simulation of magnetic systems

  • The difficulty of applying the integral approach is related to the singularity of the kernel of the integral equations

  • In order to build a discretization for the integral equations (2) the region of calculations should be divided into the tetrahedrons, satisfying the requirements of finite element method (FEM)

Read more

Summary

Introduction

The paper [1] gives an overview of existing methods and programs for the numerical simulation of magnetic systems. A. Sapozhnikov, The volume integral equation method. Integral equation method [2,3,4,5] for magnetic system calculations are discussed in this article. We consider the approach, which allows to construct high order accuracy approximation method for initial problem discretization. The main steps of the modelling process are being discussed: discretization of the initial equations for linear approximation of the magnetic field, algorithm of area division into elements, matrix element calculation procedures, a method of solving a corresponding system of nonlinear equations. The numerical results obtained for the modelling of quadrupole and dipole magnet agree with those obtained by the famous code TOSCA [6]

Volume integral equation method
Finite element mesh generating
Matrix element calculation
Iterative methods for nonlinear systems solving
Magnet system modelling
Conclusion
Full Text
Paper version not known

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call