An Unsplit Godunov Method for Ideal MHD via Constrained Transport in Three Dimensions

Abstract

We present a single step, second-order accurate Godunov scheme for ideal MHD which is an extension of the method described by Gardiner & Stone (2005) to three dimensions. This algorithm combines the corner transport upwind (CTU) method of Colella for multidimensional integration, and the constrained transport (CT) algorithm for preserving the divergence-free constraint on the magnetic field. We describe the calculation of the PPM interface states for 3D ideal MHD which must include multidimensional ``MHD source terms'' and naturally respect the balance implicit in these terms by the $\nabla \cdot \mathbf{B}=0$ condition. We compare two different forms for the CTU integration algorithm which require either 6- or 12-solutions of the Riemann problem per cell per time-step, and present a detailed description of the 6-solve algorithm. Finally, we present solutions for test problems to demonstrate the accuracy and robustness of the algorithm.

Extended version of the paper accepted for publication in JCP

Cited by in corpus (7)

• Athena: A New Code for Astrophysical MHD
• Simulating Magnetohydrodynamical Flow with Constrained Transport and Adaptive Mesh Refinement; Algorithms & Tests of the AstroBEAR Code
• Anisotropic Thermal Conduction and the Cooling Flow Problem in Galaxy Clusters
• The Magnetohydrodynamics of Shock-Cloud Interaction in Three Dimensions
• Density Probability Distribution Functions in Supersonic Hydrodynamic and MHD Turbulence
• Simulations of Magnetorotational Turbulence with a Higher-Order Godunov Scheme
• The Magnetothermal Instability in the Intracluster Medium