Large-scale magnetostatic field calculation in finite element micromagnetics with H2 -matrices
Abstract
Magnetostatic field calculations in micromagnetic simulations can be numerically expensive, particularly in the case of large-scale finite element simulations. The established finite element/boundary element method (FEM/BEM) by Fredkin and Koehler (1990) involves a densely populated matrix with unacceptable numerical costs for problems with a large number of degrees of freedom N. By using hierarchical matrices of H2 type, we show that the memory requirements of this FEM/BEM method can be reduced dramatically, effectively converting the quadratic complexity O (N2) of the problem to a linear one O (N) . We obtain matrix size reductions of nearly 99% in test cases with more than 106 degrees of freedom, and we test the computed magnetostatic energy values by means of comparison with analytic values. The efficiency of the H2 -matrix compression opens the way to large-scale magnetostatic field calculations in micromagnetic modeling, all while preserving the accuracy of the established FEM/BEM formalism.
- Publication:
-
Journal of Magnetism and Magnetic Materials
- Pub Date:
- May 2019
- DOI:
- 10.1016/j.jmmm.2018.12.103
- arXiv:
- arXiv:1811.05731
- Bibcode:
- 2019JMMM..477..118H
- Keywords:
-
- Magnetostatics;
- Boundary element method;
- Hierarchical matrix;
- Micromagnetism;
- Mathematics - Numerical Analysis;
- 78A30;
- 65N38;
- 65F30
- E-Print:
- doi:10.1016/j.jmmm.2018.12.103