Maintaining Pressure Positivity in Magnetohydrodynamic Simulations

Dinshaw S Balsara,Daniel Spicer
DOI: https://doi.org/10.1006/jcph.1998.6108
IF: 4.645
1999-01-01
Journal of Computational Physics
Abstract:Higher order Godunov schemes for solving the equations of magnetohydrodynamics (MHD) have recently become available. Because such schemes update the total energy, the pressure is a derived variable. In several problems in laboratory physics, magnetospheric physics, and astrophysics the pressure can be several orders of magnitude smaller than either the kinetic energy or the magnetic energy. Thus small discretization errors in the total energy can produce situations where the gas pressure can become negative. In this paper we design a linearized Riemann solver that works directly on the entropy density equation. We also design switches that allow us to use such a Riemann solver safely in conjunction with a normal Riemann solver for MHD. This allows us to reduce the discretization errors in the evaluation of the pressure variable. As a result we formulate strategies that maintain the positivity of pressure in all circumstances. We also show via test problems that the strategies designed here work.
physics, mathematical,computer science, interdisciplinary applications
What problem does this paper attempt to address?