A Reynolds- and Hartmann-semirobust hybrid method for magnetohydrodynamics

TL;DR

This paper introduces a new hybrid numerical method for magnetohydrodynamics that is convection-semirobust, improves convergence in diffusion regimes, and reduces computational complexity, validated by theoretical analysis and numerical experiments.

Contribution

It presents a novel hybrid approach for MHD equations that is convection-semirobust, with enhanced convergence and computational efficiency features.

Findings

Method is convection-semirobust for smooth solutions.

Achieves improved convergence in diffusion-dominated regimes.

Reduces algebraic problem size via static condensation.

Abstract

We propose and analyze a new method for the unsteady incompressible magnetohydrodynamics equations on convex domains with hybrid approximations of both vector-valued and scalar-valued fields. The proposed method is convection-semirobust, meaning that, for sufficiently smooth solutions, one can derive a priori estimates for the velocity and the magnetic field that do not depend on the inverse of the diffusion coefficients. This is achieved while at the same time providing relevant additional features, namely an improved order of convergence for the (asymptotic) diffusion-dominated regime, a small stencil (owing to the absence of inter-element penalty terms), and the possibility to significantly reduce the size of the algebraic problems through static condensation. The theoretical results are confirmed by a complete panel of numerical experiments.