Computational micromagnetics in three dimensions is of increasing interest with the development of magnetostrictive sensors and actuators. In solving the Landau-Lifshitz-Gilbert (LLG) equation, the governing equation of magnetic dynamics for ferromagnetic materials, we need to evaluate the effective field. The effective field consists of several terms, among which the demagnetizing field is of long-range nature. Evaluating the demagnetizing field directly requires work of O(N2) for a grid of N cells and thus it is the bottleneck in computational micromagnetics. A fast hierarchical algorithm using multipole approximation is developed to evaluate the demagnetizing field. We first construct a mesh hierarchy and divide the grid into boxes of different levels. The lowest level box is the whole grid while the highest level boxes are just cells. The approximate field contribution from the cells contained in a box is characterized by the box attributes, which are obtained via multipole approximation. The algorithm computes field contributions from remote cells using attributes of appropriate boxes containing those cells, and it computes contributions from adjacent cells directly. Numerical results have shown that the algorithm requires work of O(NlogN) and at the same time it achieves high accuracy. It makes micromagnetic simulation in three dimensions feasible.
[1]
G. Heber.
W. F. Brown, Magnetoelastic Interactions. (Springer Tracts in Natural Philosophy, Vol. 9) VIII + 155 S. m. 13 Fig. Berlin/Heidelberg/New York 1966. Springer-Verlag. Preis geb. DM 38,–
,
1967
.
[2]
Jian-Gang Zhu,et al.
Micromagnetic studies of thin metallic films (invited)
,
1988
.
[3]
Masud Mansuripur,et al.
Demagnetizing field computation for dynamic simulation of the magnetization reversal process
,
1988
.
[4]
J. L. Blue,et al.
Using multipoles decreases computation time for magnetostatic self-energy
,
1991
.
[5]
S. W. Yuan,et al.
Fast adaptive algorithms for micromagnetics
,
1992
.
[6]
A. Hubert,et al.
Solving Micromagnetic Problems. Towards an Optimal Numerical Method
,
1993
.
[7]
Xiaobo Tan,et al.
Computational micromagnetics for magnetostrictive actuators
,
2000,
Smart Structures.
[8]
L. Greengard.
The Rapid Evaluation of Potential Fields in Particle Systems
,
1988
.