A Third-order Implicit-Explicit Runge-Kutta Method for Landau-Lifshitz Equation with Arbitrary Damping Parameters

Yan Gui,Rui Du,Cheng Wang
2024-07-07
Abstract:A third-order accurate implicit-explicit Runge-Kutta time marching numerical scheme is proposed and implemented for the Landau-Lifshitz-Gilbert equation, which models magnetization dynamics in ferromagnetic materials, with arbitrary damping parameters. This method has three remarkable advantages:~(1) only a linear system with constant coefficients needs to be solved at each Runge-Kutta stage, which greatly reduces the time cost and improves the efficiency; (2) the optimal rate convergence analysis does not impose any restriction on the magnitude of damping parameter, which is consistent with the third-order accuracy in time for 1-D and 3-D numerical examples; (3) its unconditional stability with respect to the damping parameter has been verified by a detailed numerical study. In comparison with many existing methods, the proposed method indicates a better performance on accuracy and efficiency, and thus provides a better option for micromagnetics simulations.
Numerical Analysis
What problem does this paper attempt to address?
The problem this paper attempts to address is the development of an efficient and accurate numerical method to solve the Landau-Lifshitz equation with arbitrary damping parameters. The Landau-Lifshitz equation is used to describe the behavior of magnetization dynamics in ferromagnetic materials and is a vector nonlocal nonlinear system with non-convex constraints and potential degeneracy. Traditional numerical methods face issues such as low efficiency, insufficient accuracy, or limitations on the values of damping parameters when dealing with such equations. To this end, the authors propose a third-order implicit-explicit Runge-Kutta time discretization method (IMEX-RK3), which has the following advantages: 1. **Efficient solution of linear systems**: Each Runge-Kutta stage only requires solving a linear system with constant coefficients, significantly reducing computation time and improving efficiency. 2. **Optimal convergence rate without damping parameter restrictions**: This method maintains third-order time accuracy in both 1D and 3D numerical examples, regardless of the size of the damping parameter. 3. **Unconditional stability**: Detailed numerical studies have verified the unconditional stability of this method with respect to the damping parameter. Compared to many existing methods, this method demonstrates better performance in terms of accuracy and efficiency, providing a better choice for micromagnetic simulations.