Two improved Gauss-Seidel projection methods for Landau-Lifshitz-Gilbert equation

Panchi Li, Changjian Xie, Rui Du*, Jingrun Chen*, Xiao Ping Wang*

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

9 Citations (Scopus)

Abstract

Micromagnetic simulation is an important tool to study various dynamic behaviors of magnetic order in ferromagnetic materials. The underlying model is the Landau-Lifshitz-Gilbert equation, where the magnetization dynamics is driven by the gyromagnetic torque term and the Gilbert damping term. Numerically, considerable progress has been made in the past decades. One of the most popular methods is the Gauss-Seidel projection method developed by Xiao-Ping Wang, Carlos García-Cervera, and Weinan E in 2001. It first solves a set of heat equations with constant coefficients and updates the gyromagnetic term in the Gauss-Seidel manner, and then solves another set of heat equations with constant coefficients for the damping term. Afterwards, a projection step is applied to preserve the length constraint in the pointwise sense. This method has been verified to be unconditionally stable numerically and successfully applied to study magnetization dynamics under various controls. In this paper, we present two improved Gauss-Seidel projection methods with unconditional stability. The first method updates the gyromagnetic term and the damping term simultaneously and follows by a projection step. The second method introduces two sets of approximate solutions, where we update the gyromagnetic term and the damping term simultaneously for one set of approximate solutions and apply the projection step to the other set of approximate solutions in an alternating manner. Compared to the original Gauss-Seidel projection method which has to solve heat equations 7 times at each time step, the improved methods solve heat equations 5 times and 3 times, respectively. First-order accuracy in time and second-order accuracy in space are verified by examples in both 1D and 3D. In addition, unconditional stability with respect to both the grid size and the damping parameter is confirmed numerically. Application of both methods to a realistic material is also presented with hysteresis loops and magnetization profiles. Compared with the original method, the recorded running times suggest that savings of both methods are about 2/7 and 4/7 for the same accuracy requirement, respectively.

Original languageEnglish
Article number109046
JournalJournal of Computational Physics
Volume401
DOIs
Publication statusPublished - 15 Jan 2020
Externally publishedYes

Keywords

  • Gauss-Seidel projection method
  • Landau-Lifshitz-Gilbert equation
  • Micromagnetic simulation
  • Unconditional stability

Fingerprint

Dive into the research topics of 'Two improved Gauss-Seidel projection methods for Landau-Lifshitz-Gilbert equation'. Together they form a unique fingerprint.

Cite this