Periodic micromagnetic finite element method

TL;DR

The paper introduces a periodic micromagnetic finite element method (PM-FEM) for efficiently solving 1D, 2D, and 3D periodic problems using the Landau-Lifshitz-Gilbert equation, extending existing FEM micromagnetic solvers to handle periodicities.

Contribution

It develops a novel PM-FEM that incorporates periodic Green's functions and fast Fourier transform techniques to efficiently simulate periodic micromagnetic structures.

Findings

PM-FEM performs comparably or better than non-periodic codes.

It effectively handles complex 1D, 2D, and 3D periodic structures.

The method is efficiently implemented on CPU and GPU architectures.

Abstract

Periodic micromagnetic finite element method (PM-FEM) is introduced to solve periodic unit cell problems using the Landau-Lifshitz-Gilbert equation. PM-FEM is applicable to general problems with 1D, 2D, and 3D periodicities. PM-FEM is based on a non-periodic FEM-based micromagnetic solver and extends it in several aspects to account for periodicities, including the computation of exchange and magnetostatic fields. For the exchange field, PM-FEM modifies the sparse matrix construction for computing the Laplace operator to include additional elements arising due to the periodicities. For the magnetostatic field, the periodic extensions include modifications in the local operators, such as gradient, divergence, and surface magnetic charges as well as the long-range superposition operator for computing the periodic scalar potential. The local operators are extended to account for the periodicities similar to handling the Laplace operator. For the long-range superposition operator, PM-FEM utilizes a periodic Green's function (PGF) and fast spatial convolutions. The PGF is computed rapidly via exponentially rapidly convergent sums. The spatial convolutions are accomplished via a modified fast Fourier transform based adaptive integral method that allows calculating spatial convolutions with non-uniform meshes in $O(N\log N)$ numerical operations. PM-FEM is implemented on CPU and GPU based computer architectures. PM-FEM allows efficiently handling cases of structures contained withing the periodic unit cell touching or not touching its boundaries as well as structures that protrude beyond the unit cell boundaries. PM-FEM is demonstrated to have about the same or even higher performance than its parent non-periodic code. The demonstrated numerical examples show the efficiency of PM-FEM for highly complex structures with 1D, 2D, and 3D periodicities.