On a decoupled linear FEM integrator for Eddy-current-LLG

We propose a numerical integrator for the coupled system of the eddy-current equation with the nonlinear Landau-Lifshitz-Gilbert equation. The considered effective field contains a general field contribution, and we particularly cover exchange, anisotropy, applied field, and magnetic field (stemming from the eddy-current equation). Even though the considered problem is nonlinear, our scheme requires only the solution of two linear systems per time-step. Moreover, our algorithm decouples both equations so that in each time-step, one linear system is solved for the magnetization, and afterwards one linear system is solved for the magnetic field. Unconditional convergence -- at least of a subsequence -- towards a weak solution is proved, and our analysis even provides existence of such weak solutions. Numerical experiments with a micromagnetic benchmark problem underline the performance of the proposed algorithm.Comment: 13 pages, 4 figure

A decoupled and unconditionally convergent linear FEM integrator for the Landau-Lifshitz-Gilbert equation with magnetostriction

To describe and simulate dynamic micromagnetic phenomena, we consider a coupled system of the nonlinear Landau-Lifshitz-Gilbert equation and the conservation of momentum equation. This coupling allows one to include magnetostrictive effects into the simulations. Existence of weak solutions has recently been shown in Carbou et al. (2011) (Global weak solutions for the Landau-Lifschitz equation with magnetostriction. Math. Meth. Appl. Sci., 34, 1274-1288). In our contribution, we give an alternate proof which additionally provides an effective numerical integrator. The latter is based on linear finite elements (FEs) in space and a linear-implicit Euler time-stepping. Despite the nonlinearity, only two linear systems have to be solved per timestep, and the integrator fully decouples both equations. Finally, we prove unconditional convergence—at least of a subsequence—towards, and hence existence of, a weak solution of the coupled system, as timestep size and spatial mesh size tend to zero. We conclude the work with numerical experiments, which study the discrete blow-up of the LLG equation as well as the influence of the magnetostrictive term on the discrete blow-u

A decoupled and unconditionally convergent linear FEM integrator for the Landau-Lifshitz-Gilbert equation with magnetostriction

To describe and simulate dynamic micromagnetic phenomena, we consider a coupled system of the nonlinear Landau-Lifshitz-Gilbert equation and the conservation of momentum equation. This coupling allows one to include magnetostrictive effects into the simulations. Existence of weak solutions has recently been shown in Carbou et al. (2011) (Global weak solutions for the Landau-Lifschitz equation with magnetostriction. Math. Meth. Appl. Sci., 34, 1274-1288). In our contribution, we give an alternate proof which additionally provides an effective numerical integrator. The latter is based on linear finite elements (FEs) in space and a linear-implicit Euler time-stepping. Despite the nonlinearity, only two linear systems have to be solved per timestep, and the integrator fully decouples both equations. Finally, we prove unconditional convergence—at least of a subsequence—towards, and hence existence of, a weak solution of the coupled system, as timestep size and spatial mesh size tend to zero. We conclude the work with numerical experiments, which study the discrete blow-up of the LLG equation as well as the influence of the magnetostrictive term on the discrete blow-u

Spin-polarized transport in ferromagnetic multilayers: An unconditionally convergent FEM integrator

We propose and analyze a decoupled time-marching scheme for the coupling of the Landau-Lifshitz-Gilbert equation with a quasilinear diffusion equation for the spin accumulation. This model describes the interplay of magnetization and electron spin accumulation in magnetic and non-magnetic multilayer structures. Despite the strong nonlinearity of the overall PDE system, the proposed integrator requires only the solution of two linear systems per time-step. Unconditional convergence of the integrator towards weak solutions is proved.Comment: 23 pages, 1 figur

Spin-polarized transport in ferromagnetic multilayers: an unconditionally convergent FEM integrator

23 pages, 1 figureWe propose and analyze a decoupled time-marching scheme for the coupling of the Landau-Lifshitz-Gilbert equation with a quasilinear diffusion equation for the spin accumulation. This model describes the interplay of magnetization and electron spin accumulation in magnetic and non-magnetic multilayer structures. Despite the strong nonlinearity of the overall PDE system, the proposed integrator requires only the solution of two linear systems per time-step. Unconditional convergence of the integrator towards weak solutions is proved