An accurate scalar potential finite element method for linear,two-dimensional magnetostatics problems

We have assessed the accuracy of a commercially available computer software package for finite element method calculations of magnetostatic fields. The computer program, MSC/NASTRAN,

1 Available from the MacNeal-Schwendler Corporation, Los Angeles, CA 90041, U.S.A.

is well known for its wide applicability in structural analysis and heat transfer problems. We exploit the fact that the differential equations of magnetostatics are identical to those for heat transfer if the magnetic field problem is formulated with the reduced scalar potential.¹ Consequently, the powerful, optimized numerical routines of NASTRAN can immediately be applied to two- and three-dimensional linear magneto-statics problems. Application of the NASTRAN reduced scalar potential approach to a ‘worst case’ two-dimensional problem for which an analytic solution is available has yielded much better accuracy than was recently reported² for a reduced scalar potential calculation using a different finite element program. Furthermore, our method exhibits completely satisfactory performance with regard to computational expense and accuracy for a linear electromagnet with an air gap. Our analysis opens the way for large three-dimensional magnetostatics calculations at far greater economy than is possible with the more commonly used vector potential and boundary integral methods.