Positivity-preserving entropy filtering for the ideal magnetohydrodynamics equations

Abstract

In this work, we present a positivity-preserving adaptive filtering approach for discontinuous spectral element approximations of the ideal magnetohydrodynamics equations. This approach combines the entropy filtering method (Dzanic and Witherden, J. Comput. Phys., 468, 2022) for shock capturing in gas dynamics along with the eight-wave method for enforcing a divergence-free magnetic field. Due to the inclusion of non-conservative source terms, an operator-splitting approach is introduced to ensure that the positivity and entropy constraints remain satisfied by the discrete solution. Furthermore, a computationally efficient algorithm for solving the optimization process for this nonlinear filtering approach is presented. The resulting scheme can robustly resolve strong discontinuities on general unstructured grids without tunable parameters while recovering high-order accuracy for smooth solutions. The efficacy of the scheme is shown in numerical experiments on various problems including extremely magnetized blast waves and three-dimensional magnetohydrodynamic instabilities.

24 pages, 17 figures

References in corpus (7)

• Vertical structure of a supernova-driven turbulent magnetized ISM
• Nonlinear Evolution of the Magnetohydrodynamic Rayleigh-Taylor Instability
• Positivity-Preserving Entropy-Based Adaptive Filtering for Discontinuous Spectral Element Methods
• Subcell limiting strategies for discontinuous Galerkin spectral element methods
• A positivity preserving strategy for entropy stable discontinuous Galerkin discretizations of the compressible Euler and Navier-Stokes equations
• Three-Dimensional MHD Simulation of Caltech Plasma Jet Experiment: First Results
• Monolithic parabolic regularization of the MHD equations and entropy principles