Accurate numerical methods for micromagnetics simulations with general geometries

Carlos J. García-Cervera, Zydrunas Gimbutas, E. Weinan

Research output: Contribution to journalArticlepeer-review

55 Scopus citations


In current FFT-based algorithms for micromagnetics simulations, the boundary is typically replaced by a staircase approximation along the grid lines, either eliminating the incomplete cells or replacing them by complete cells. Sometimes the magnetizations at the boundary cells are weighted by the volume of the sample in the corresponding cell. We show that this leads to large errors in the computed exchange and stray fields. One consequence of this is that the predicted switching mechanism depends sensitively on the orientation of the numerical grid. We present a boundary-corrected algorithm to efficiently and accurately handle the incomplete cells at the boundary. We show that this boundary-corrected algorithm greatly improves the accuracy in micromagnetics simulations. We demonstrate by using A. Arrott's example of a hexagonal element that the switching mechanism is predicted independently of the grid orientation.

Original languageEnglish (US)
Pages (from-to)37-52
Number of pages16
JournalJournal of Computational Physics
Issue number1
StatePublished - Jan 1 2003

All Science Journal Classification (ASJC) codes

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • General Physics and Astronomy
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics


  • Cartesian grid
  • Landau-Lifshitz equation
  • Micromagnetics
  • Neumann problems
  • Stray field


Dive into the research topics of 'Accurate numerical methods for micromagnetics simulations with general geometries'. Together they form a unique fingerprint.

Cite this