Iterative solution and preconditioning for the tangent plane scheme in computational micromagnetics

Abstract

The tangent plane scheme is a time-marching scheme for the numerical solution of the nonlinear parabolic Landau-Lifshitz-Gilbert equation (LLG), which describes the time evolution of ferromagnetic configurations. Exploiting the geometric structure of LLG, the tangent plane scheme requires only the solution of one linear variational form per time-step, which is posed in the discrete tangent space determined by the nodal values of the current magnetization. We develop an effective solution strategy for the arising constrained linear systems, which is based on appropriate Householder reflections. We derive possible preconditioners, which are (essentially) independent of the time-step, and prove that the preconditioned GMRES algorithm leads to linear convergence. Numerical experiments underpin the theoretical findings.2010 Mathematics Subject Classification. 35K35, 65M60, 65F08, 65F10, 65M22, 65Z05.

Gamma-Convergent Projection-Free Finite Element Methods for Nematic Liquid Crystals: The Ericksen Model

“…Altogether, both (4.1) and (4.3) are linear symmetric positive definite systems in the unknowns t i,ℓ h and s i+1 h . The orthogonality constraint in (4.1) can be imposed at the linear algebraic level by introducing a Lagrange multiplier associated with it (see, e.g., the discussion in [8, section 7.2.5]) or via a null-space method as done, e.g., in [25,23,20]. In this work, we implement it using a Lagrange multiplier.…”

Finite element simulation of nonlinear bending models for thin elastic rods and plates

“…A further discussion of related preconditioners in the context of micromagnetics can be found in [41]. We illustrate the construction of the bases of the null space and the performance of different solution strategies in the context of harmonic maps into spheres.…”

“…We review the numerical treatment of rods following [9,18] and plates as proposed in [8,10]. For the efficient iterative solution we adopt ideas from [45,41]. We discuss the treatment of bilayer plates following [14,13], illustrate a method that enforces injectivity of deformations in the case of rods following [19,17], and propose methods for the numerical solution of bending deformations with shearing effects following ideas from [12].…”

Computational micromagnetics with Commics

“…In a variant from [13,38], the formal convergence order in time has been increased from one to two. Effective solution strategies and preconditioning for the resulting constrained linear system have recently been proposed in [55].…”