A hybrid high-order scheme for the stationary, incompressible magnetohydrodynamics equations

We propose and analyse a hybrid high-order (HHO) scheme for stationary incompressible magnetohydrodynamics equations. The scheme has an arbitrary order of accuracy and is applicable on generic polyhedral meshes. For sources that are small enough, we prove error estimates in energy norm for the velocity and magnetic field, and $L^2$-norm for the pressure; these estimates are fully robust with respect to small faces, and of optimal order with respect to the mesh size. Using compactness techniques, we also prove that the scheme converges to a solution of the continuous problem, irrespective of the source being small or large. Finally, we illustrate our theoretical results through 3D numerical tests on tetrahedral and Voronoi mesh families.