In this talk, we propose a novel exactly divergence-free and well-balanced hybrid finite volume / finite element scheme for the numerical solution of the incompressible viscous and resistive magnetohydrodynamics (MHD) equations on staggered unstructured mixed-element meshes in two and three space dimensions. The algorithm is based on the splitting the equations into several subsystems so that each of them can be discretized with a particular scheme to preserve some fundamental structural features of the MHD system at the discrete level.

The use of face-based staggered grids allows to account for the divergence-free conditions of the velocity field and of the magnetic field in a rather natural manner. In particular, the non-linear convective and the viscous terms in the momentum equation are solved at the aid of an explicit finite volume scheme, while the magnetic field is evolved in an exactly divergence-free manner via an explicit finite volume method based on a discrete form of the Stokes law in the edges/faces of each primary element. The latter method is stabilized by the proper choice of the numerical resistivity in the computation of the electric field in the vertices/edges of the 2D/3D elements. To achieve higher order of accuracy, a piecewise linear polynomial is reconstructed for the magnetic field, which is guaranteed to be exactly divergence-free via a constrained L^2 projection. Finally, a classical continuous finite element approach is employed to compute the pressure in the vertices of the primary mesh. Besides, we take into account the known equilibrium solution at each step of the new algorithm so that the method becomes exactly well-balanced.

The proposed methodology is validated against known exact and numerical reference solutions including an MHD lid-driven cavity benchmark. Furthermore, long-time simulations of Soloviev equilibrium solutions in simplified 3D tokamak configurations are studied to assess the capability of the method to maintain stationary equilibria exactly over very long integration times in general grids.