Application of the domain-integrated field relations method to the solution of large scale static and stationary magnetic field problems

A new efficient implementation of the domain-integrated field relations approach to computing static and stationary magnetic fields is presented. The philosophy behind this implementation is to combine the robustness of this method with a computationally efficient discretization technique that uses first-order Cartesian expansion functions on a simplicial mesh. The accuracy of the implementation is illustrated by solving a two-dimensional field problem. The ability of the method to accurately model the behavior of the field at interfaces between highly contrasting media and at geometrical singularities is emphasized. A detailed efficiency analysis of the present implementation is performed