Papers
Topics
Authors
Recent
Gemini 2.5 Flash
Gemini 2.5 Flash
169 tokens/sec
GPT-4o
7 tokens/sec
Gemini 2.5 Pro Pro
45 tokens/sec
o3 Pro
4 tokens/sec
GPT-4.1 Pro
38 tokens/sec
DeepSeek R1 via Azure Pro
28 tokens/sec
2000 character limit reached

IMEX-RK methods for Landau-Lifshitz equation with arbitrary damping (2312.15654v1)

Published 25 Dec 2023 in math.NA, cs.NA, math-ph, and math.MP

Abstract: Magnetization dynamics in ferromagnetic materials is modeled by the Landau-Lifshitz (LL) equation, a nonlinear system of partial differential equations. Among the numerical approaches, semi-implicit schemes are widely used in the micromagnetics simulation, due to a nice compromise between accuracy and efficiency. At each time step, only a linear system needs to be solved and a projection is then applied to preserve the length of magnetization. However, this linear system contains variable coefficients and a non-symmetric structure, and thus an efficient linear solver is highly desired. If the damping parameter becomes large, it has been realized that efficient solvers are only available to a linear system with constant, symmetric, and positive definite (SPD) structure. In this work, based on the implicit-explicit Runge-Kutta (IMEX-RK) time discretization, we introduce an artificial damping term, which is treated implicitly. The remaining terms are treated explicitly. This strategy leads to a semi-implicit scheme with the following properties: (1) only a few linear system with constant and SPD structure needs to be solved at each time step; (2) it works for the LL equation with arbitrary damping parameter; (3) high-order accuracy can be obtained with high-order IMEX-RK time discretization. Numerically, second-order and third-order IMEX-RK methods are designed in both the 1-D and 3-D domains. A comparison with the backward differentiation formula scheme is undertaken, in terms of accuracy and efficiency. The robustness of both numerical methods is tested on the first benchmark problem from National Institute of Standards and Technology. The linearized stability estimate and optimal rate convergence analysis are provided for an alternate IMEX-RK2 numerical scheme as well.

Definition Search Book Streamline Icon: https://streamlinehq.com
References (49)
  1. L. Landau and E. Lifshitz, “On the theory of the dispersion of magnetic permeability in ferromagnetic bodies,” Physikalische Zeitschrift Der Sowjetunion, vol. 63, no. 9, pp. 153–169, 1935.
  2. T. Gilbert, “A Lagrangian formulation of the gyromagnetic equation of the magnetization field,” Physical Review, vol. 100, p. 1243, 1955.
  3. R. An, “Optimal error estimates of linearized Crank–Nicolson Galerkin method for Landau–Lifshitz equation,” Journal of Scientific Computing, vol. 69, no. 1, pp. 1–27, 2016.
  4. R. An and W. Sun, “Analysis of backward Euler projection FEM for the Landau-Lifshitz equation,” IMA Journal of Numerical Analysis, vol. 42, no. 3, pp. 2336–2360, 2022.
  5. S. Bartels and P. Andreas, “Convergence of an implicit finite element method for the Landau–Lifshitz–Gilbert equation,” SIAM journal on numerical analysis, vol. 44, no. 4, pp. 1405–1419, 2006.
  6. S. Bartels, J. Ko, and A. Prohl, “Numerical analysis of an explicit approximation scheme for the Landau–Lifshitz–Gilbert equation,” Mathematics of Computation, vol. 77, no. 262, pp. 773–788, 2008.
  7. W. Chen, Y. Liu, C. Wang, and S. Wise, “An optimal-rate convergence analysis of a fully discrete finite difference scheme for Cahn-Hilliard-Hele-Shaw equation,” Math. Comp., vol. 85, pp. 2231–2257, 2016.
  8. W. Chen, C. Wang, S. Wang, X. Wang, and S. Wise, “Energy stable numerical schemes for ternary Cahn-Hilliard system,” J. Sci. Comput., vol. 84, p. 27, 2020.
  9. K. Cheng and C. Wang, “Long time stability of high order multi-step numerical schemes for two-dimensional incompressible Navier-Stokes equations,” SIAM Journal on Numerical Analysis, vol. 54, pp. 3123–3144, 2016.
  10. S. Conde, S. Gottlieb, Z. Grant, and J. Shadid, “Implicit and implicit-explicit strong stability preserving Runge-Kutta methods with high linear order,” Journal of Scientific Computing, vol. 83, no. 2, pp. 667–690, 2017.
  11. L. Dong, C. Wang, S. Wise, and Z. Zhang, “Optimal rate convergence analysis of a numerical scheme for the ternary Cahn-Hilliard system with a Flory-Huggins-deGennes energy potential,” Journal of Computational and Applied Mathematics, vol. 406, p. 114474, 2022.
  12. C. Duan, C. Liu, C. Wang, and X. Yue, “Convergence analysis of a numerical scheme for the porous medium equation by an energetic variational approach,” Numerical Mathematics: Theory, Methods and Applications, vol. 13, pp. 1–18, 2020.
  13. C. Duan, W. Chen, C. Liu, C. Wang, and X. Yue, “A second order accurate, energy stable numerical scheme for one-dimensional porous medium equation by an energetic variational approach,” Communications in Mathematical Sciences, vol. 20, no. 4, pp. 987–1024, 2022.
  14. W. E and X. Wang, “Numerical methods for the Landau–Lifshitz equation,” SIAM journal on numerical analysis, pp. 1647–1665, 2001.
  15. A. Fuwa, T. Ishiwata, and M. Tsutsumi, “Finite difference scheme for the Landau–Lifshitz equation,” Japan journal of industrial and applied mathematics, vol. 29, no. 1, pp. 83–110, 2012.
  16. S. Gottlieb, C. Shu, and E. Tadmor, “Strong stability-preserving high-order time discretization methods,” SIAM Review, vol. 43, no. 1, pp. 89–112, 2001.
  17. S. Gottlieb and C. Wang, “Stability and convergence analysis of fully discrete Fourier collocation spectral method for 3-D viscous Burgers’ equation,” Journal of Scientific Computing, vol. 53, pp. 102–128, 2012.
  18. Z. Guan, C. Wang, and S. Wise, “A convergent convex splitting scheme for the periodic nonlocal Cahn-Hilliard equation,” Numer. Math., vol. 128, pp. 377–406, 2014.
  19. Z. Guan, J. Lowengrub, and C. Wang, “Convergence analysis for second order accurate schemes for the periodic nonlocal Allen-Cahn and Cahn-Hilliard equations,” Math. Methods Appl. Sci., vol. 40, no. 18, pp. 6836–6863, 2017.
  20. C. Liu, C. Wang, Y. Wang, and S. Wise, “Convergence analysis of the variational operator splitting scheme for a reaction-diffusion system with detailed balance,” SIAM Journal on Numerical Analysis, vol. 60, no. 2, pp. 781–803, 2022.
  21. C. Liu, C. Wang, S. Wise, X. Yue, and S. Zhou, “A positivity-preserving, energy stable and convergent numerical scheme for the Poisson-Nernst-Planck system,” Mathematics of Computation, vol. 90, pp. 2071–2106, 2021.
  22. C. Wang and J.-G. Liu, “Convergence of gauge method for incompressible flow,” Mathematics of Computation, vol. 69, pp. 1385–1407, 2000.
  23. R. Samelson, R. Temam, C. Wang, and S. Wang, “Surface pressure Poisson equation formulation of the primitive equations: Numerical schemes,” SIAM Journal on Numerical Analysis, vol. 41, pp. 1163–1194, 2003.
  24. C. Wang, J.-G. Liu, and H. Johnston, “Analysis of a fourth order finite difference method for incompressible Boussinesq equation,” Numerische Mathematik, vol. 97, pp. 555–594, 2004.
  25. X.-P. Wang, C. Garcıa-Cervera, and W. E, “A Gauss–Seidel projection method for micromagnetics simulations,” Journal of Computational Physics, vol. 171, no. 1, pp. 357–372, 2001.
  26. C. Xie, C. Garcıa-Cervera, C. Wang, Z. Zhou, and J. Chen, “Second-order semi-implicit projection methods for micromagnetics simulations,” Journal of Computational Physics, vol. 404, p. 109104, 2020.
  27. C. Nan and H. Song, “Error estimates of local discontinuous Galerkin method with implicit-explicit Runge Kutta for two-phase miscible flow in porous media,” Applied Numerical Mathematics, vol. 169, pp. 334–350, 2021.
  28. H. Wang, S. Wang, Q. Zhang, and C. Shu, “Local discontinuous Galerkin methods with explicit-implicit time marching for multi-dimensional convection-diffusion problems,” ESAIM: Mathematical Modeling and Numerical Analysis, vol. 50, pp. 1083–1105, 2016.
  29. H. Wang, Q. Zhang, S. Wang, and C. Shu, “Local discontinuous Galerkin methods with explicit–implicit–null time discretizations for solving nonlinear diffusion problems,” Science China Mathematics, vol. 63, no. 1, pp. 183–204, 2020.
  30. H. Wang, J. Zheng, F. You, H. Guo, and Q. Zhang, “Local discontinuous Galerkin method with explicit-implicit time marching for incompressible miscible displacement problem in porous media,” Journal of Scientific Computing, vol. 78, pp. 1–28, 2019.
  31. U. Ascher, S. Ruuth, and R. Spiteri, “Implicit-explicit Runge–Kutta methods for time–dependent partial differential equations,” Applied Numerical Mathematics, vol. 25, no. 2-3, pp. 151–167, 1997.
  32. Y. Zheng and J. Zhu, “Switching field variation in patterned submicron magnetic film elements,” Journal of Applied Physics, vol. 81, no. 8, pp. 5471–5473, 1997.
  33. T. Schrefl, D. Suess, W. Scholz, H. Forster, V. Tsiantos, and J. Fidler, “Finite element micromagnetics,” Computational Electromagnetics, pp. 165–181, 2003.
  34. F. Alouges and P. Jaisson, “Convergence of a finite element discretization for the Landau–Lifshitz equations in micromagnetism,” Mathematical Models and Methods in Applied Sciences, vol. 16, no. 02, pp. 299–316, 2006.
  35. C. García-Cervera and W. E, “Improved Gauss–Seidel projection method for micromagnetics simulations,” IEEE transactions on magnetics, vol. 39, no. 3, pp. 1766–1770, 2003.
  36. I. Cimrák, “Error estimates for a semi–implicit numerical scheme solving the Landau–Lifshitz equation with an exchange field,” IMA journal of numerical analysis, vol. 25, no. 3, pp. 611–634, 2005.
  37. H. Gao, “Optimal error estimates of a linearized backward euler fem for the Landau–Lifshitz equation,” SIAM Journal on Numerical Analysis, vol. 52, no. 5, pp. 2574–2593, 2014.
  38. S. Boscarino, F. Filbet, and G. Russo, “High order semi–implicit schemes for time dependent partial differential equations,” Journal of Scientific Computing, vol. 68, no. 3, pp. 975–1001, 2016.
  39. J. Chen, C. Wang, and C. Xie, “Convergence analysis of a second–order semi–implicit projection method for Landau–Lifshitz equation,” Applied Numerical Mathematics, vol. 168, pp. 55–74, 2021.
  40. Y. Cai, J. Chen, C. Wang, and C. Xie, “A second-order numerical method for Landau–Lifshitz–Gilbert equation with large damping parameters,” Journal of Computational Physics, vol. 451, p. 110831, 2022.
  41. Micromagnetic Modeling Activity Group and Technology, “2000,” https://www.ctcms.nist.gov/~rdm/mumag.org.html.
  42. R. An, H. Gao, and W. Sun, “Optimal error analysis of Euler and Crank–Nicolson projection finite difference schemes for Landau–Lifshitz equation,” SIAM J. Numer. Anal., vol. 59, no. 3, pp. 1639–1662, 2021.
  43. M. Kruzík and A. Prohl, “Recent developments in the modeling, analysis, and numerics of ferromagnetism,” SIAM Rev., vol. 48, no. 3, pp. 439–483, 2006.
  44. I. Cimrák, “A survey on the numerics and computations for the Landau-Lifshitz equation of micromagnetism,” Arch. Comput. Methods Eng., vol. 15, no. 3, pp. 277–309, 2008.
  45. C. J. García-Cervera, “Numerical micromagnetics: a review,” Bol. Soc. Esp. Mat. Apl., vol. 39, pp. 103–135, 2007.
  46. P. Li, C. Xie, R. Du, J. Chen, and X. Wang, “Two improved Gauss-Seidel projection methods for Landau-Lifshitz-Gilbert equation,” J. Comput. Phys., vol. 401, p. 109046, 2020.
  47. G. Akrivis, M. Feischl, B. Kovács, and C. Lubich, “Higher-order linearly implicit full discretization of the Landau-Lifshitz-Gilbert equation,” Math. Comp., vol. 90, pp. 995–1038, 2021.
  48. Y. Cai, J. Chen, C. Wang, and C. Xie, “Error analysis of a linear numerical scheme for the Landau-Lifshitz equation with large damping parameters,” Math. Meth. Appl. Sci., vol. 46, pp. 18 952–18 974, 2023.
  49. C. Gear, “Hybrid methods for initial value problems in ordinary differential equations,” Journal of the Society for Industrial and Applied Mathematics, Series B: Numerical Analysis, vol. 2, no. 1, pp. 69–86, 1965.
Citations (1)
View on Semantic Scholar

Summary

We haven't generated a summary for this paper yet.

Ai Sparkles Streamline Icon: https://streamlinehq.com Summarize Now
X Twitter Logo Streamline Icon: https://streamlinehq.com

Tweets