Convergence of an implicit-explicit midpoint scheme for computational micromagnetics

Abstract

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.

Cited by in corpus (8)

• Second-order semi-implicit projection methods for micromagnetics simulations
• Linear second-order IMEX-type integrator for the (eddy current) Landau-Lifshitz-Gilbert equation
• Convergent tangent plane integrators for the simulation of chiral magnetic skyrmion dynamics
• Computational micromagnetics with Commics
• Iterative solution and preconditioning for the tangent plane scheme in computational micromagnetics
• A Gauss-Seidel projection method with the minimal number of updates for stray field in micromagnetic simulations
• Deterministic strain-control of stability and current-induced motion of skyrmions in chiral magnets
• A second-order numerical method for Landau-Lifshitz-Gilbert equation with large damping parameters