Vol. 16, No. 1, 2021

Download this article
Download this article For screen
For printing
Recent Issues
Volume 16, Issue 1
Volume 15, Issue 2
Volume 15, Issue 1
Volume 14, Issue 2
Volume 14, Issue 1
Volume 13, Issue 2
Volume 13, Issue 1
Volume 12, Issue 1
Volume 11, Issue 2
Volume 11, Issue 1
Volume 10, Issue 2
Volume 10, Issue 1
Volume 9, Issue 2
Volume 9, Issue 1
Volume 8, Issue 1
Volume 7, Issue 2
Volume 7, Issue 1
Volume 6, Issue 1
Volume 5, Issue 2
Volume 5, Issue 1
Volume 4, Issue 1
Volume 3, Issue 1
Volume 2, Issue 1
Volume 1, Issue 1
The Journal
About the Journal
Editorial Board
Submission Guidelines
Submission Form
Policies for Authors
Ethics Statement
ISSN: 2157-5452 (e-only)
ISSN: 1559-3940 (print)
Author Index
To Appear
Other MSP Journals
Globally divergence-free DG scheme for ideal compressible MHD

Dinshaw S. Balsara, Rakesh Kumar and Praveen Chandrashekar

Vol. 16 (2021), No. 1, 59–98
DOI: 10.2140/camcos.2021.16.59

The high-accuracy solution of the MHD equations is of great interest in various fields of physics, mathematics, and engineering. Higher-order DG schemes offer low dissipation and dispersion as well as the ability to model complex geometries, which is very desirable in various applications. Numerical solution of the MHD equations is made challenging by the fact that the PDE system has an involution constraint. Therefore, we construct high-order, globally divergence-free DG schemes for compressible MHD. The modes of the fluid variables are collocated at the zones of the mesh; the magnetic field components and their higher-order modes are collocated at the faces of the mesh. The fluid equations are evolved using classical DG, while the magnetic fields are evolved using a novel DG-like approach, first proposed by Balsara and Käppeli (J. Comput. Phys. 336 (2017), 104–127). This DG-like method ensures the globally divergence-free evolution of the magnetic field.

The method is built around three building blocks. The first building block consists of a divergence-free reconstruction of the magnetic field. The second building block consists of a DG-like formulation of Faraday’s law that provides a weak-form interpretation of Stokes’ law (as opposed to traditional DG, which relies on Gauss’s law). To provide a physically consistent electric field for the update of Faraday’s law, we use the third building block, which consists of a multidimensional Riemann solver that is evaluated at the edges of the mesh. We recognize that the limiting of facial variables makes the design of the MHD limiter very different from the usual DG limiter. As a result, a limiter strategy is presented for DG schemes which retains the traditional DG limiting approach while building into it a positivity-enforcement step and a step that updates the facial modes in a constraint-preserving fashion. This limiter is crucial to the robust and physically consistent operation of our DG scheme for MHD even at high orders.

It is shown that our schemes meet their design accuracies at second, third, and fourth orders on smooth test problems. Several stringent test problems with complex flow features are presented, which are robustly handled by our DG method.

ideal compressible MHD, divergence-free, discontinuous Galerkin method
Mathematical Subject Classification
Primary: 35L65, 35L67, 65M12, 65M60, 76W05
Received: 11 June 2020
Revised: 8 December 2020
Accepted: 16 December 2020
Published: 19 January 2021
Dinshaw S. Balsara
Department of Physics
University of Notre Dame
Notre Dame, IN
United States
Rakesh Kumar
Centre for Applicable Mathematics
Tata Institute of Fundamental Research
Praveen Chandrashekar
Centre for Applicable Mathematics
Tata Institute of Fundamental Research