Strong convergence of finite element schemes for the stochastic Landau--Lifshitz--Bloch equation
Agus L. Soenjaya
TL;DR
The paper addresses the stochastic Landau--Lifshitz--Bloch equation for high-temperature micromagnetics and develops strong convergence results for semi-implicit and implicit finite element schemes in 1D and 2D, including a regularised 2D formulation. It introduces localisation-based error analysis and new exponential moment bounds for the exact solution, yielding mean-square exponential stability and, in 1D with small noise, uniqueness of the invariant measure. The authors derive sharp convergence rates in $L^2(\Omega)$ and, in 1D, optimal probabilistic rates in $H^1$ for various schemes, supported by numerical experiments validating the theory. The work provides rigorous justification and practical guidance for stochastic FEM discretisations in micromagnetics and sheds light on stability properties under thermal noise.
Abstract
The dynamics of magnetisation in a bounded ferromagnet in $\mathbb{R}^d$ ($d=1,2$) at high temperatures can be described by the stochastic Landau--Lifshitz--Bloch (sLLB) equation, which is a vector-valued quasilinear stochastic partial differential equation. In this paper, assuming adequate regularity of the initial data, we establish strong convergence in $L^2(Ω)$ of several semi-implicit and implicit fully discrete finite element schemes for the sLLB equation, together with explicit convergence rates. The analysis relies on localised error estimates and new exponential moment bounds for the exact solution. As a by-product, these moment bounds yield mean-square exponential stability of solutions and uniqueness of the invariant measure in one spatial dimension under a small noise assumption. We also sharpen existing convergence-in-probability results for the numerical schemes. Numerical experiments are presented to illustrate and support the theoretical findings.
