  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.
