Second-order semi-implicit projection methods for micromagnetics simulations

Micromagnetics simulations require accurate approximation of the magnetiza- tion dynamics described by the Landau-Lifshitz-Gilbert equation, which is non- linear, nonlocal, and has a non-convex constraint, posing interesting challenges in developing numerical methods. In this paper, we propose two s...

Full description

Bibliographic Details
Authors: García-Cervera, C.J., Wangd, C., Zhoue, Z., Chena, J., Xie, C.
Format: article
Status:Published version
Publication Date:2019
Country:España
Institution:Basque Center for Applied Mathematics (BCAM)
Repository:BIRD. BCAM's Institutional Repository Data
OAI Identifier:oai:bird.bcamath.org:20.500.11824/1164
Online Access:http://hdl.handle.net/20.500.11824/1164
https://doi.org/10.1016/j.jcp.2019.109104
Access Level:Open access
Keyword:Micromagnetics simulation
Landau-Lifshitz-Gilbert equation
Backward differentiation formula
Second-order accuracy
Hysteresis loop
Description
Summary:Micromagnetics simulations require accurate approximation of the magnetiza- tion dynamics described by the Landau-Lifshitz-Gilbert equation, which is non- linear, nonlocal, and has a non-convex constraint, posing interesting challenges in developing numerical methods. In this paper, we propose two second-order semi-implicit projection methods based on the second-order backward differen- tiation formula and the second-order interpolation formula using the informa- tion at previous two temporal steps. Unconditional unique solvability of both methods is proved, with their second-order accuracy verified through numerical examples in both 1D and 3D. The efficiency of both methods is compared to that of another two popular methods. In addition, we test the robustness of both methods for the first benchmark problem with a ferromagnetic thin film material from National Institute of Standards and Technology.