Abstract
This paper proposes a method to generate adaptively 2D and 3D finite element meshes taking into account the continuity requirements of the magnetic field at the interface between two neighboring elements. First, this paper proposes a new error estimator that includes the Zienkiewicz and Zhu error norm estimator and the boundary rules in the electromagnetic field. Using a 2D simple model, this paper decides two parameters of the proposed estimator. Next, this paper presents a 3D mesh generation method based on the Voronoi-Delaunay theory, which ensures that the bounding surface of the domain is contained in the triangulation. The method has the capability to decrease the amount of information on the connectivity of boundary nodes by generating nodes not only in the interior of the domain but also on its surface. Two simple magnetostatic field problems are provided to illustrate the usefulness of the proposed method.