Linear second-order IMEX-type integrator for the (eddy current) Landau–Lifshitz–Gilbert equation. (28th November 2019)

• Record Type: Journal Article
• Authors: Di Fratta, Giovanni; Pfeiler, Carl-Martin; Praetorius, Dirk; Ruggeri, Michele; Stiftner, Bernhard
• Is Part Of: IMA journal of numerical analysis. Volume 40:Number 4(2020)
• Page Start: 2802
• Page End: 2838
• Publication Date: 2019-11-28
• Subjects: micromagnetism -- finite elements -- linear second-order time integration -- implicit-explicit time-marching scheme -- unconditional convergence; Numerical analysis -- Periodicals; 519.405
• DOI: 10.1093/imanum/drz046
• Languages: English
• ISSNs: 0272-4979

Abstract

Abstract: Combining ideas from Alouges et al. (2014, A convergent and precise finite element scheme for Landau–Lifschitz–Gilbert equation. Numer. Math., 128, 407–430) and Praetorius et al. (2018, Convergence of an implicit-explicit midpoint scheme for computational micromagnetics. Comput. Math. Appl., 75, 1719–1738) we propose a numerical algorithm for the integration of the nonlinear and time-dependent Landau–Lifshitz–Gilbert (LLG) equation, which is unconditionally convergent, formally (almost) second-order in time, and requires the solution of only one linear system per time step. Only the exchange contribution is integrated implicitly in time, while the lower-order contributions like the computationally expensive stray field are treated explicitly in time. Then we extend the scheme to the coupled system of the LLG equation with the eddy current approximation of Maxwell equations. Unlike existing schemes for this system, the new integrator is unconditionally convergent, (almost) second-order in time, and requires the solution of only two linear systems per time step.