Abstract
Reasons were found for the low accuracy of the demagnetizing field in the neighborhood of the corner regions calculated by the previously developed numerical scheme based on the Fredkin-Koehler method, and improved methods were developed. The improvement includes redivision of the computing cells along the boundary of the calculation region, numerical integration of the magnetic potential on the boundary using additional calculation points obtained by Lagrange interpolation, and use of a functional form which approximates the analytic solution of the magnetic potential in the neighborhood of the corners of the calculation region. The improved method was found to decrease the numerical error of the demagnetizing field at the calculation points close to the corners from about 10% to 2.5% in the case of a square prism of an infinite length magnetized along an edge of its cross-section. The effect of improving the demagnetizing field was confirmed by simulating a two-dimensional reversal process in a square prism based on the Landau-Lifshitz-Gilbert equation.
