A Kernel-Free Boundary Integral Method for 2-D Magnetostatics Analysis

Performing magnetostatic analysis accurately and efficiently is very important for multi-objective optimization of electromagnetic device designs. In this research, a kernel-free boundary integral method (KFBIM) has been introduced for solving magnetic computations in a toroidal core geometry in 2-D. This study is very relevant in the design and optimization of toroidal inductors or transformers used in electrical systems, where lighter weight, higher inductance, higher efficiency, and lower leakage flux are required. The governing partial differential equations (PDEs) have been formulated as a system of the Fredholm integral equations of the second kind. Unlike traditional boundary integral methods or boundary element methods, KFBIM does not require an analytical form of Green’s function for evaluating integrals via numerical quadrature. Instead, KFBIM computes integrals by solving an equivalent interface problem on a Cartesian mesh. Compared with traditional finite difference methods for solving the governing PDEs directly, KFBIM produces a well-conditioned linear system. Therefore, the KFBIM requires only a fixed number of iterations when an iterative method [e.g., generalized minimal residual method (GMRES)] is applied for solving the linear system, and the numerical solution is not sensitive to computer round-off errors. The obtained results are then compared with a commercial finite element solver (ANSYS), which shows excellent agreement. It should be noted that, compared with FEM, the KFBIM does not require a body-fit mesh and can achieve high accuracy with a coarse mesh. In particular, the calculations of the magnetic potential and the tangential field intensity on the boundaries are more stable and exhibit almost no oscillations.