Convergence of an implicit-explicit midpoint scheme for computational micromagnetics
This work provides a provably convergent numerical method for micromagnetics simulations, addressing the challenge of efficiently handling lower-order field terms.
The authors developed a second-order implicit-explicit midpoint scheme for the Landau-Lifschitz-Gilbert equation and proved unconditional convergence to a weak solution, with numerical experiments confirming the theory.
Based on lowest-order finite elements in space, we consider the numerical integration of the Landau-Lifschitz-Gilbert equation (LLG). The dynamics of LLG is driven by the so-called effective field which usually consists of the exchange field, the external field, and lower-order contributions such as the stray field. The latter requires the solution of an additional partial differential equation in full space. Following Bartels and Prohl (2006) (Convergence of an implicit finite element method for the Landau-Lifschitz-Gilbert equation. SIAM J. Numer. Anal. 44), we employ the implicit midpoint rule to treat the exchange field. However, in order to treat the lower-order terms effectively, we combine the midpoint rule with an explicit Adams-Bashforth scheme. The resulting integrator is formally of second-order in time, and we prove unconditional convergence towards a weak solution of LLG. Numerical experiments underpin the theoretical findings.