A divergence-free semi-implicit finite volume scheme for ideal, viscous and resistive magnetohydrodynamics

Abstract

In this paper we present a novel pressure-based semi-implicit finite volume solver for the equations of compressible ideal, viscous and resistive magnetohydrodynamics (MHD). The new method is conservative for mass, momentum and total energy and in multiple space dimensions it is constructed in such a way as to respect the divergence-free condition of the magnetic field exactly, also in the presence of resistive effects. This is possible via the use of multi-dimensional Riemann solvers on an appropriately staggered grid for the time evolution of the magnetic field and a double curl formulation of the resistive terms. The new semi-implicit method for the MHD equations proposed here discretizes all terms related to the pressure in the momentum equation and the total energy equation implicitly, making again use of a properly staggered grid for pressure and velocity. The time step of the scheme is restricted by a CFL condition based only on the fluid velocity and the Alfvén wave speed and is not based on the speed of the magnetosonic waves. Our new method is particularly well-suited for low Mach number flows and for the incompressible limit of the MHD equations, for which it is well-known that explicit density-based Godunov-type finite volume solvers become increasingly inefficient and inaccurate due to the increasingly stringent CFL condition and the wrong scaling of the numerical viscosity in the incompressible limit. We show a relevant MHD test problem in the low Mach number regime where the new semi-implicit algorithm is a factor of 50 faster than a traditional explicit finite volume method, which is a very significant gain in terms of computational efficiency. However, our numerical results confirm that our new method performs well also for classical MHD test cases with strong shocks. In this sense our new scheme is a true all Mach number flow solver.

26 pages, 12 figures,1 table

References in corpus (5)

• Multi-dimensional Numerical Scheme for Resistive Relativistic MHD
• A pressure-based semi-implicit space-time discontinuous Galerkin method on staggered unstructured meshes for the solution of the compressible Navier-Stokes equations at all Mach numbers
• High order ADER schemes for a unified first order hyperbolic formulation of Newtonian continuum mechanics coupled with electro-dynamics
• Semi-implicit discontinuous Galerkin methods for the incompressible Navier-Stokes equations on adaptive staggered Cartesian grids
• An efficient semi-implicit method for three-dimensional non-hydrostatic flows in compliant arterial vessels

Cited by in corpus (13)

• A structure-preserving staggered semi-implicit finite volume scheme for continuum mechanics
• A semi-implicit hybrid finite volume / finite element scheme for all Mach number flows on staggered unstructured meshes
• Simulation of non-Newtonian viscoplastic flows with a unified first order hyperbolic model and a structure-preserving semi-implicit scheme
• A staggered semi-implicit hybrid finite volume / finite element scheme for the shallow water equations at all Froude numbers
• A novel structure preserving semi-implicit finite volume method for viscous and resistive magnetohydrodynamics
• A semi-implicit compressible model for atmospheric flows with seamless access to soundproof and hydrostatic dynamics
• A new thermodynamically compatible finite volume scheme for magnetohydrodynamics
• A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors
• High order asymptotic preserving finite difference WENO schemes with constrained transport for MHD equations in all sonic Mach numbers
• A well-balanced and exactly divergence-free staggered semi-implicit hybrid finite volume/finite element scheme for the incompressible MHD equations
• Numerical viscosity and resistivity in MHD turbulence simulations
• Locally Structure-Preserving div-curl operators for high order Discontinuous Galerkin schemes
• Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows