A Crank–Nicolson scheme for the Landau–Lifshitz equation without damping

• Published in: Journal of computational and applied mathematics Vol. 234; no. 2; pp. 613 - 623
• Main Authors: Jeong, Darae, Kim, Junseok
• Format: Journal Article
• Language: English
• Published: Kidlington Elsevier B.V 15.05.2010; Elsevier
• Subjects: Crank–Nicolson; Exact sciences and technology; Finite difference method; Landau–Lifshitz equation; Mathematical analysis; Mathematics; Nonlinear algebraic and transcendental equations; Nonlinear multigrid method; Numerical analysis; Numerical analysis. Scientific computation; Partial differential equations; Sciences and techniques of general use; Crank–Nicolson; 65Z05; Nonlinear multigrid method; Landau–Lifshitz equation; Finite difference method; Damping; Transcendental equation; Landau-Lifshitz equation; Crank Nicolson method; Numerical method; Discrete system; Algorithm; Discrete scheme; Convergence; Equation system; Non linear equation; Crank-Nicolson; Parabolic equation; Numerical analysis; Applied mathematics; Numerical solution; Algebraic equation; Nonlinearity; Landau Lifshitz equation
• ISSN: 0377-0427; 1879-1778
• DOI: 10.1016/j.cam.2010.01.002

An accurate and efficient numerical approach, based on a finite difference method with Crank–Nicolson time stepping, is proposed for the Landau–Lifshitz equation without damping. The phenomenological Landau–Lifshitz equation describes the dynamics of ferromagnetism. The Crank–Nicolson method is very popular in the numerical schemes for parabolic equations since it is second-order accurate in time. Although widely used, the method does not always produce accurate results when it is applied to the Landau–Lifshitz equation. The objective of this article is to enumerate the problems and then to propose an accurate and robust numerical solution algorithm. A discrete scheme and a numerical solution algorithm for the Landau–Lifshitz equation are described. A nonlinear multigrid method is used for handling the nonlinearities of the resulting discrete system of equations at each time step. We show numerically that the proposed scheme has a second-order convergence in space and time.