A computer program for the 2-D magnetostatic problem based on integral equations for the field of the conductors and boundary elements

The iterative design of the two-dimensional cross section of a beam transport magnet having infinitely permeable iron boundaries requires a fast means of computing the field of the conductors. Solutions in the form of series expansions are used for rectangular iron boundaries, and programs based on the method of images are used to simulate circular iron boundaries. It is pointed out that a single procedure or program for dealing with an arbitrary iron boundary would be useful. The present program has been tested with rectangular and circular iron boundaries and provision has been made for the use of other curves. It uses complex contour integral equations for the field of the constant-current density conductors and complex line integrals for the field of the piecewise-linear boundary elements. The results obtained show that the quality of the solution at the boundary is dependent on the segmentation, for a roughly constant number of segments