Numerical solution of the forward magnetic field problem for models with irregular polyhedron discretization taking into account the "demagnetization effect"

A performance-effective numerical method for magnetic field calculation is proposed. The method works with irregular polyhedron discretization which enables us to construct models with magnetic objects of arbitrary shape. As a case study, a model of a well in plane-parallel layer is considered. The model is approximated with dense irregular grid, elements of which are polyhedrons. With the help of conjugate gradient method, we solve "demagnetization effect" equation and calculate total magnetic field on the plane above the well. For a well of 0.25m radius and 8m height "demagnetization effect" is of order of 2% relative to the field induced by the object placed in equivalent of Earth magnetic field.