BDF2-type integrator for Landau-Lifshitz-Gilbert equation in micromagnetics: unconditional weak convergence to weak solutions

Abstract

We consider the Landau-Lifshitz-Gilbert equation (LLG) that models time-dependent micromagnetic phenomena. We propose a full discretization that employs first-order finite elements in space and a BDF2-type two-step method in time. In each time step, only one linear system of equations has to be solved. We employ linear interpolation in time to reconstruct the discrete space-time magnetization. We prove that the integrator is unconditionally stable and thus guarantees that a subsequence of the reconstructed magnetization converges weakly in $H^1$ towards a weak solution of LLG in the space-time domain. Numerical experiments verify that the proposed integrator is indeed first-order in space and second-order in time.