By Tikhonov’s regularization method and Richardson’s method in mathematical modeling magnetostatics problems
ЛИТ, ОИЯИ, г. Дубна, Московской области
A nonlinear magnetostatics inverse problem is investigated for a case when it is needed to create a required magnetic field using conductors which coordinates vary on condition that current value is equal in all the conductors. It is known that the same problems fall under the category of non correct problems. Mathematical statement of this type of nonlinear magnetostatics inverse problems is given. The proposed numerical algorithm using the regularization method by A.N.Tikhonov permits to overcome comparatively easily the difficulties connected with badly conditioned equation systems to which usually the magnetostatics inverse problems reduce. This algorithm allows one to calculate the existing winding geometry of superconducting dipole ironless magnet which ensures a magnetic field homogeneity up to within a rectangular aperture. Magnetic field distribution which is created by wires of an infinitely thin cross section being positioned near mag¬netic field surfaces and providing magnetic field index correc¬tion in the working region of accelerator is calculated by the mirror reflection method. Some expressions by which it is possible to calculate magnetic field gradient of arbitrary set of correcting wires are given. Current magnitude in correcting wires for 10 GeV JINR synchrophasotron magnetic field index correction under field level about 0.023 T are calculated and optimized by the regularization method. This work gives results of numerical simulation of a superconducting magnetic focusing system. When modeling this system, an additional control was done over the condition approximation accuracy by using Richardson method. Design, optimization and construction of magnetic systems are a complex engineering and physical issue. Mathematical modeling of the magnetostatics problems allows us to make easy and accelerate the process of research on physical properties of a specific magnetic system. In particular, the mathematical simulation provides a way for studying those segments of the magnet where it is quite difficult and even impossible for making measurements (for example, in the core or very thin constructions). Thus, it is the mathematical modeling that provides a way for finding an optimal configuration of the magnet in each particular case. It should be noted that in the course of modeling a specific magnetic system some problems related to the peculiarities of its design arise. This is connected particularly with the choice of an optimal calculated grid for the boundary magnetostatics problem. It is the problem that arose for the authors of this work, because the magnetic system considered is of small size and has a very thin gap. As a result, various variants of the magnetic system with a round aperture have been investigated numerically, various configurations of the system have been proposed for realization.
The work is dedicated to the memory of the doctor of physical and mathematical sciences I. A. Shelaeva.