Adaptively meshed boundary integral equation method for nonlinear magnetostatics

A general approach to the treatment of nonlinear problems in magnetostatics using boundary integral equation methods is discussed. Adaptive meshing is used to obtain a more optimal distribution of boundary elements using the dependent variable as the criterion. Permanent magnets are shown to be accurately modeled using this method. Example solutions are presented for a ring magnet in the vicinity of a ferromagnetic sphere. Forces are easily calculated using Maxwell stress and are seen to compare very well against measurement