Computation of magnetostatic field using second order edge elements in 3D

Several second order edge elements have been applied to solving magnetostatic problems. The performances of these elements are compared through an example of magnetic circuit. In order to ensure the compatibility of the system equations and hence the convergence, the current density is represented by the curl of a source field. This avoids an explicit gauge condition which is cumbersome in the case of high order elements.