A new thermodynamically compatible finite volume scheme for magnetohydrodynamics

Abstract

In this paper we propose a novel thermodynamically compatible finite volume scheme for the numerical solution of the equations of magnetohydrodynamics (MHD) in one and two space dimensions. As shown by Godunov in 1972, the MHD system can be written as overdetermined symmetric hyperbolic and thermodynamically compatible (SHTC) system. More precisely, the MHD equations are symmetric hyperbolic in the sense of Friedrichs and satisfy the first and second principles of thermodynamics. In a more recent work on SHTC systems, \cite{Rom1998}, the entropy density is a primary evolution variable, and total energy conservation can be shown to be a \textit{consequence} that is obtained after a judicious linear combination of all other evolution equations. The objective of this paper is to mimic the SHTC framework also on the discrete level by directly discretizing the \textit{entropy inequality}, instead of the total energy conservation law, while total energy conservation is obtained via an appropriate linear combination as a \textit{consequence} of the thermodynamically compatible discretization of all other evolution equations. As such, the proposed finite volume scheme satisfies a discrete cell entropy inequality \textit{by construction} and can be proven to be nonlinearly stable in the energy norm due to the discrete energy conservation. In multiple space dimensions the divergence-free condition of the magnetic field is taken into account via a new thermodynamically compatible generalized Lagrangian multiplier (GLM) divergence cleaning approach. The fundamental properties of the scheme proposed in this paper are mathematically rigorously proven. The new method is applied to some standard MHD benchmark problems in one and two space dimensions, obtaining good results in all cases.

References in corpus (11)

• Second Order Accurate Schemes for Magnetohydrodynamics With Divergence-Free Reconstruction
• Multidimensional HLLE Riemann solver; Application to Euler and Magnetohydrodynamic Flows
• A provably entropy stable subcell shock capturing approach for high order split form DG for the compressible Euler Equations
• A general framework to construct schemes satisfying additional conservation relations. Application to entropy conservative and entropy dissipative schemes
• Ideal GLM-MHD: About the entropy consistent nine-wave magnetic field divergence diminishing ideal magnetohydrodynamics equations
• A divergence-free semi-implicit finite volume scheme for ideal, viscous and resistive magnetohydrodynamics
• Rarefaction Shocks, Shock Errors and Low Order of Accuracy in ZEUS
• A high-order nonconservative approach for hyperbolic equations in fluid dynamics
• A semi-implicit hybrid finite volume / finite element scheme for all Mach number flows on staggered unstructured meshes
• On thermodynamically compatible finite volume schemes for continuum mechanics
• A new continuum model for general relativistic viscous heat-conducting media

Cited by in corpus (3)

• A well-balanced and exactly divergence-free staggered semi-implicit hybrid finite volume/finite element scheme for the incompressible MHD equations
• An Eulerian hyperbolic model for heat transfer derived via Hamilton's principle: analytical and numerical study
• An exactly curl-free finite-volume scheme for a hyperbolic compressible barotropic two-phase model