A Constrained Transport Scheme for MHD on Unstructured Static and Moving Meshes

Abstract

Magnetic fields play an important role in many astrophysical systems and a detailed understanding of their impact on the gas dynamics requires robust numerical simulations. Here we present a new method to evolve the ideal magnetohydrodynamic (MHD) equations on unstructured static and moving meshes that preserves the magnetic field divergence-free constraint to machine precision. The method overcomes the major problems of using a cleaning scheme on the magnetic fields instead, which is non-conservative, not fully Galilean invariant, does not eliminate divergence errors completely, and may produce incorrect jumps across shocks. Our new method is a generalization of the constrained transport (CT) algorithm used to enforce the $\nabla \cdot \mathbf{B}=0$ condition on fixed Cartesian grids. Preserving $\nabla \cdot \mathbf{B}=0$ at the discretized level is necessary to maintain the orthogonality between the Lorentz force and $\mathbf{B}$. The possibility of performing CT on a moving mesh provides several advantages over static mesh methods due to the quasi-Lagrangian nature of the former (i.e., the mesh generating points move with the flow), such as making the simulation automatically adaptive and significantly reducing advection errors. Our method preserves magnetic fields and fluid quantities in pure advection exactly.

13 pages, 9 figures, accepted to MNRAS. Animations available at http://www.cfa.harvard.edu/~pmocz/research.html

References in corpus (4)

• Athena: A New Code for Astrophysical MHD
• A discontinuous Galerkin method for solving the fluid and MHD equations in astrophysical simulations
• Linear Structures in the Core of the Coma Cluster of Galaxies
• A switch to reduce resistivity in smoothed particle magnetohydrodynamics

Cited by in corpus (20)

• Improving the convergence properties of the moving-mesh code AREPO
• Accurate, Meshless Methods for Magneto-Hydrodynamics
• Magnetohydrodynamic Simulations of Active Galactic Nucleus Disks and Jets
• A Constrained-Gradient Method to Control Divergence Errors in Numerical MHD
• The large-scale properties of simulated cosmological magnetic fields
• DISCO: A 3D Moving-Mesh Magnetohydrodynamics Code Designed for the Study of Astrophysical Disks
• A moving mesh unstaggered constrained transport scheme for magnetohydrodynamics
• Stellar orbit evolution in close circumstellar disc encounters
• Local simulations of MRI turbulence with meshless methods
• Magnetohydrodynamics with Physics Informed Neural Operators
• Non-ideal magnetohydrodynamics on a moving mesh
• Simulating the magnetorotational instability on a moving mesh with the shearing box approximation
• A Moving Mesh Hydrodynamic Solver for ChaNGa
• Time Dependent Radiation Hydrodynamics on a Moving Mesh
• General Relativistic Hydrodynamics on a Moving-mesh I: Static Spacetimes
• A new first-order formulation of the Einstein equations exploiting analogies with electrodynamics
• Python Radiative Transfer Emission code (PyRaTE): non-LTE spectral lines simulations
• Multidimensional upwind hydrodynamics on unstructured meshes using Graphics Processing Units I. Two-dimensional uniform meshes
• Meshless Methods for Magnetohydrodynamics with Vector Potential
• Numerical Methods for Simulating Star Formation