Error estimates for a helicity-preserving finite element discretisation of an incompressible magnetohydrodynamics system
L. Beirao da Veiga, K. Hu, L. Mascotto
TL;DR
The paper addresses the problem of deriving rigorous error estimates for a seven-field, helicity- and energy-preserving finite element discretisation of incompressible MHD. It develops a structure-preserving FE method based on a de Rham sequence with commuting projections, a discrete curl operator, and divergence-free constraints to ensure discrete conservation of energy and magnetic and cross helicities. The main result is an optimal-order a priori error estimate for the semi-discrete scheme, showing convergence of the velocity and magnetic field (and associated variables) with rate $h^{k+1}$ under suitable regularity, while preserving key physical invariants. This work provides theoretical guarantees for the reliability and physical fidelity of helicity-preserving MHD simulations and lays a foundation for fully discrete analyses with time integrators that maintain quadratic invariants.
Abstract
We derive error estimates of a finite element method for the approximation of solutions to a seven-fields formulation of a magnetohydrodynamics model, which preserves the energy of the system, and the magnetic and cross helicities on the discrete level.
