On the measurements of numerical viscosity and resistivity in Eulerian MHD codes

Abstract

We propose a simple ansatz for estimating the value of the numerical resistivity and the numerical viscosity of any Eulerian MHD code. We test this ansatz with the help of simulations of the propagation of (magneto)sonic waves, Alfven waves, and the tearing mode instability using the MHD code Aenus. By comparing the simu- lation results with analytical solutions of the resistive-viscous MHD equations and an empirical ansatz for the growth rate of tearing modes we measure the numerical viscosity and resistivity of Aenus. The comparison shows that the fast-magnetosonic speed and wavelength are the characteristic velocity and length, respectively, of the aforementioned (relatively simple) systems. We also determine the dependance of the numerical viscosity and resistivity on the time integration method, the spatial reconstruction scheme and (to a lesser extent) the Riemann solver employed in the simulations. From the measured results we infer the numerical resolution (as a function of the spatial reconstruction method) required to properly resolve the growth and saturation level of the magnetic field amplified by the magnetorotational instability in the post-collapsed core of massive stars. Our results show that it is to the best advantage to resort to ultra-high order methods (e.g., 9th-order Monotonicity Preserving method) to tackle this problem properly, in particular in three dimensional simulations.

arXiv comment36 pages, 35 figures, accepted for publication in ApJS

References in corpus (24)

• On the Nature of Core Collapse Supernova Explosions
• The Magnetorotational Instability in Core Collapse Supernova Explosions
• A large scale dynamo and magnetoturbulence in rapidly rotating core-collapse supernovae
• The Role of Turbulence in Neutrino-Driven Core-Collapse Supernova Explosions
• Neutrino-driven convection versus advection in core collapse supernovae
• Semi-global simulations of the magneto-rotational instability in core collapse supernovae
• The Dominance of Neutrino-Driven Convection in Core-Collapse Supernovae
• Deceleration of arbitrarily magnetized GRB ejecta: the complete evolution
• Magnetic field amplification and magnetically supported explosions of collapsing, non-rotating stellar cores
• Axisymmetric simulations of magnetorotational core collapse: Approximate inclusion of general relativistic effects
• Fast reconnection in relativistic plasmas: the magnetohydrodynamics tearing instability revisited
• MRI channel flows and their parasites
• Neutrino-driven Turbulent Convection and Standing Accretion Shock Instability in Three-Dimensional Core-Collapse Supernovae
• On the Saturation of the Magnetorotational Instability via Parasitic Modes
• General relativistic simulations of pasive-magneto-rotational core collapse with microphysics
• Angular Momentum Transport in Protoplanetary and Black-Hole Accretion Disks: The Role of Parasitic Modes in the Saturation of MHD Turbulence
• Improved Accuracy of High-Order WENO Finite Volume Methods on Cartesian Grids
• On the maximum magnetic field amplification by the magnetorotational instability in core-collapse supernovae
• Numerical simulations of the magnetorotational instability in protoneutron stars: I. Influence of buoyancy
• Global Simulations of Magnetorotational Instability in The Collapsed Core of A Massive Star
• Termination of the magnetorotational instability via parasitic instabilities in core-collapse supernovae
• Implicit large eddy simulations of anisotropic weakly compressible turbulence with application to core-collapse supernovae
• The Evolution and Impacts of Magnetorotational Instability in Magnetized Core-Collapse Supernovae
• An Efficient Implementation of Flux Formulae in Multidimensional Relativistic Hydrodynamical Codes

Cited by in corpus (30)

• Hydrodynamics of core-collapse supernovae and their progenitors
• Magnetorotational core collapse of possible GRB progenitors. I. Explosion mechanisms
• Plasmoid formation in global GRMHD simulations and AGN flares
• Stellar transits across a magnetized accretion torus as a mechanism for plasmoid ejection
• Magnetorotational core collapse of possible GRB progenitors. II. Formation of protomagnetars and collapsars
• Magnetic reconnection and plasmoid formation in three-dimensional accretion flows around black holes
• A three-phase amplification of the cosmic magnetic field in galaxies
• The impact of non-dipolar magnetic fields in core-collapse supernovae
• Magnetic fields in the Galactic halo restrict fountain-driven recycling and accretion
• Striped Blandford/Znajek jets from advection of small scale magnetic field
• Heavy sterile neutrinos in stellar core-collapse
• On The Development of Multidimensional Progenitor Models For Core-collapse Supernovae
• 3D Simulations of Oxygen Shell Burning with and without Magnetic Fields
• Magnetar bursts due to Alfvén wave nonlinear breakout
• An HLLC Riemann Solver for Resistive Relativistic Magnetohydrodynamics
• Non-ideal magnetohydrodynamics on a moving mesh
• Numerical viscosity and resistivity in MHD turbulence simulations
• Tayler-Spruit dynamo in stably stratified rotating fluids: Application to proto-magnetars
• Post-dynamical inspiral phase of common envelope evolution. The role of magnetic fields
• Methodology for estimating the magnetic Prandtl number and application to solar surface small-scale dynamo simulations
• Tayler-Spruit dynamo in binary neutron star merger remnants
• Relativistic gas accretion onto supermassive black Hole binaries from inspiral through merger
• The impact of resistivity on the variability of black hole accretion flows
• General-relativistic resistive-magnetohydrodynamics simulations of self-consistent magnetized rotating neutron stars
• Numerical viscosity in simulations of the two-dimensional Kelvin-Helmholtz instability
• Computational General Relativistic Force-Free Electrodynamics: I. Multi-Coordinate Implementation and Testing
• General relativistic force-free electrodynamics with a discontinuous Galerkin-finite difference hybrid method
• Axisymmetric flows on the torus geometry
• Time integration for neutrino radiation transport using minimally implicit Runge-Kutta methods
• On the Impact of the Numerical Method on Magnetic Reconnection and Particle Acceleration -- I. The MHD case