A constrained transport scheme for MHD on unstructured static and moving meshes

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.

[1]  Y. Levin,et al.  MAGNETICALLY LEVITATING ACCRETION DISKS AROUND SUPERMASSIVE BLACK HOLES , 2012, 1201.4873.

[2]  P. Roe,et al.  A Solution-Adaptive Upwind Scheme for Ideal Magnetohydrodynamics , 1999 .

[3]  D. Balsara,et al.  MAGNETOHYDRODYNAMIC MODELING OF SOLAR SYSTEM PROCESSES ON GEODESIC GRIDS , 2013 .

[4]  E. Parker The dynamical state of the interstellar gas and field. II. , 1966 .

[5]  A. Strong,et al.  Propagation of Cosmic-Ray Nucleons in the Galaxy , 1998, astro-ph/9807150.

[6]  Dario Grasso,et al.  Magnetic Fields in the Early Universe , 2001 .

[7]  S. Orszag,et al.  Small-scale structure of two-dimensional magnetohydrodynamic turbulence , 1979, Journal of Fluid Mechanics.

[8]  Mark Vogelsberger,et al.  A discontinuous Galerkin method for solving the fluid and magnetohydrodynamic equations in astrophysical simulations , 2013, 1305.5536.

[9]  S. Sridhar,et al.  Toward a theory of interstellar turbulence. 2. Strong Alfvenic turbulence , 1994 .

[10]  C. Munz,et al.  Hyperbolic divergence cleaning for the MHD equations , 2002 .

[11]  Evghenii Gaburov,et al.  Astrophysical Weighted Particle Magnetohydrodynamics , 2010, 1006.4159.

[12]  Gopal-Krishna,et al.  Was the Cosmic Web of Protogalactic Material Permeated by Lobes of Radio Galaxies During the Quasar Era? , 2001, astro-ph/0108117.

[13]  J. Hawley,et al.  Instability, turbulence, and enhanced transport in accretion disks , 1997 .

[14]  J. Hawley,et al.  Simulation of magnetohydrodynamic flows: A Constrained transport method , 1988 .

[15]  P. Teuben,et al.  Athena: A New Code for Astrophysical MHD , 2008, 0804.0402.

[16]  V. Springel E pur si muove: Galilean-invariant cosmological hydrodynamical simulations on a moving mesh , 2009, 0901.4107.

[17]  Daniel J. Price,et al.  A switch to reduce resistivity in smoothed particle magnetohydrodynamics , 2013, 1309.5437.

[18]  D. Balsara,et al.  A Staggered Mesh Algorithm Using High Order Godunov Fluxes to Ensure Solenoidal Magnetic Fields in Magnetohydrodynamic Simulations , 1999 .

[19]  G. Tóth The ∇·B=0 Constraint in Shock-Capturing Magnetohydrodynamics Codes , 2000 .

[20]  R. Blandford,et al.  Electromagnetic extraction of energy from Kerr black holes , 1977 .

[21]  Andreas Bauer,et al.  Magnetohydrodynamics on an unstructured moving grid , 2011, 1108.1792.

[22]  R. Blandford,et al.  Hydromagnetic flows from accretion discs and the production of radio jets , 1982 .

[23]  K. Kusano,et al.  A multi-state HLL approximate Riemann solver for ideal magnetohydrodynamics , 2005 .

[24]  V. Springel,et al.  Simulations of magnetic fields in isolated disc galaxies , 2012, 1212.1452.

[25]  A. Fabian,et al.  Linear Structures in the Core of the Coma Cluster of Galaxies , 2013, Science.

[26]  Paul C. Duffell,et al.  TESS: A RELATIVISTIC HYDRODYNAMICS CODE ON A MOVING VORONOI MESH , 2011, 1104.3562.

[27]  E. N. Parker THE DYNAMICAL STATE OF THE INTERSTELLAR GAS AND FIELD. III. TURBULENCE AND ENHANCED DIFFUSION. , 1966 .

[28]  J. Michael Picone,et al.  Evolution of the Orszag-Tang vortex system in a compressible medium , 1991 .

[29]  J. Landstreet,et al.  Magnetic Fields of Nondegenerate Stars , 2009, 0904.1938.