2024/06/09 by Le, Kim-Ngan, Soenjaya, Agus L., Tran, Thanh · 2 citations
#35K59 #35Q60 #65M12 #65M60 #FOS: Mathematics #Numerical Analysis (math.NA)
paper · doi:10.48550/arxiv.2406.05808
The Landau--Lifshitz--Bloch equation (LLBE) describes the evolution of the magnetic spin field in ferromagnets at high temperatures. In this paper, we study the numerical approximation of the LLBE on bounded polytopal domains in ℝd, where d≤ 3. We first establish the existence and uniqueness of strong solutions to the LLBE and propose a linear, fully discrete, conforming finite element scheme for its approximation. While this scheme is shown to converge, the obtained rate is suboptimal. To address this shortcoming, we introduce a viscous (pseudo-parabolic) regularisation of the LLBE, which we call the ε-LLBE. For this regularised problem, we prove the unique existence of strong solutions and establish a rate of convergence of the solution \boldsymboluε of the ε-LLBE to the solution \boldsymbolu of the LLBE as ε→ 0+. Furthermore, we propose a linear, fully discrete, conforming finite element scheme to approximate the solution of the ε-LLBE. Given sufficiently smooth initial data, error analysis is performed to show stability and uniform-in-time convergence of the scheme. Finally, several numerical simulations are presented to corroborate our theoretical results.