Tetrahedral meshing is the critical element of finite element modeling. Recently, adaptive meshing has been commonly used. In such meshing, according to the Delaunay method, mesh density is connected with the curvature of modeled object’s edge. Such a method is especially efficient during the modeling of mechanical systems. However, the efficiency of commonly used meshing algorithms is strongly limited in surface-focused phenomena, such as eddy current induced by magneto-dynamic processes. The paper proposes the improved method of Delaunay meshing considering the specific requirements of magnetodynamic systems. In the proposed method, tetrahedral mesh density may be flexibly modified according to the needs of modeled physical phenomena, such as eddy current density. As a result, physical effects may be efficiently and accurately described in finite element models. The paper presents the example of implementing the proposed solution for cylindrical wire. The complete source code is available as open-source software for further practical use and development.