Numerical Simulation of a Three-Dimensional Strongly Magnetic Material Tensor in Wave Vector Domain

Demagnetization effects are typically ignored in numerical simulations of magnetic material tensors. However, as the magnetization susceptibility increases, demagnetization effects become more pronounced, thereby gradually decreasing the accuracy of magnetic material tensor simulations. An iterative calculation method based on the wave vector domain is proposed to address this issue in the context of strongly magnetic materials. This method uses the integral equations satisfied by the magnetic potential, employs a two-dimensional Fourier transform to convert the three-dimensional integral problem into a one-dimensional problem, and utilizes the shape function method to solve this one-dimensional problem. An iterative tight-binding operator is introduced to ensure the convergence of the algorithm. High-precision numerical simulations of strongly magnetic material tensors are achieved on the basis of this foundation. Comparisons with analytical solutions for strongly magnetic spherical and spherical shell models demonstrate that the proposed method is applicable to a wide range of magnetization susceptibilities. Under normal conditions of strong magnetization (χ = 1SI), the method converges rapidly and closely matches the analytical solutions. At higher magnetization levels (1SI < χ ≤ 100SI), it stabilizes after a finite number of iterations while maintaining high computational accuracy. At extremely strong magnetization levels (100SI < χ ≤ 500SI), the convergence speed and computational accuracy of the method decrease relatively, but it still satisfies the requirements for computational accuracy.