High-order Magnetohydrodynamics for Astrophysics with an Adaptive Mesh Refinement Discontinuous Galerkin Scheme

Abstract

Modern astrophysical simulations aim to accurately model an ever-growing array of physical processes, including the interaction of fluids with magnetic fields, under increasingly stringent performance and scalability requirements driven by present-day trends in computing architectures. Discontinuous Galerkin methods have recently gained some traction in astrophysics, because of their arbitrarily high order and controllable numerical diffusion, combined with attractive characteristics for high performance computing. In this paper, we describe and test our implementation of a discontinuous Galerkin (DG) scheme for ideal magnetohydrodynamics in the AREPO-DG code. Our DG-MHD scheme relies on a modal expansion of the solution on Legendre polynomials inside the cells of an Eulerian octree-based AMR grid. The divergence-free constraint of the magnetic field is enforced using one out of two distinct cell-centred schemes: either a Powell-type scheme based on nonconservative source terms, or a hyperbolic divergence cleaning method. The Powell scheme relies on a basis of locally divergence-free vector polynomials inside each cell to represent the magnetic field. Limiting prescriptions are implemented to ensure non-oscillatory and positive solutions. We show that the resulting scheme is accurate and robust: it can achieve high-order and low numerical diffusion, as well as accurately capture strong MHD shocks. In addition, we show that our scheme exhibits a number of attractive properties for astrophysical simulations, such as lower advection errors and better Galilean invariance at reduced resolution, together with more accurate capturing of barely resolved flow features. We discuss the prospects of our implementation, and DG methods in general, for scalable astrophysical simulations.

38 pages, 29 figures, 2 tables, accepted for publication in MNRAS

References in corpus (19)

• PLUTO: a Numerical Code for Computational Astrophysics
• Athena: A New Code for Astrophysical MHD
• A High Order Godunov Scheme with Constrained Transport and Adaptive Mesh Refinement for Astrophysical MHD
• Moving mesh cosmology: tracing cosmological gas accretion
• Efficient, High Accuracy ADER-WENO Schemes for Hydrodynamics and Divergence-Free Magnetohydrodynamics
• A robust numerical scheme for highly compressible magnetohydrodynamics: Nonlinear stability, implementation and tests
• Magnetohydrodynamic Simulations of Disk Galaxy Formation: the Magnetization of The Cold and Warm Medium
• Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence
• Kinematic Dynamos using Constrained Transport with High Order Godunov Schemes and Adaptive Mesh Refinement
• Constrained Hyperbolic Divergence Cleaning for Smoothed Particle Magnetohydrodynamics
• A New Numerical Scheme for Cosmic Ray Transport
• Ideal GLM-MHD: About the entropy consistent nine-wave magnetic field divergence diminishing ideal magnetohydrodynamics equations
• ADER discontinuous Galerkin schemes for general-relativistic ideal magnetohydrodynamics
• Constrained hyperbolic divergence cleaning in smoothed particle magnetohydrodynamics with variable cleaning speeds
• A small-scale dynamo in feedback-dominated galaxies - III. Cosmological simulations
• Magnetized Moving Mesh Merger of a Carbon-Oxygen White Dwarf Binary
• Non-ideal magnetohydrodynamics on a moving mesh
• CosmosDG: An hp-Adaptive Discontinuous Galerkin Code for Hyper-Resolved Relativistic MHD
• Planet-disc interactions with Discontinuous Galerkin Methods using GPUs

Cited by in corpus (14)

• The Arepo public code release
• Sphenix: Smoothed Particle Hydrodynamics for the next generation of galaxy formation simulations
• Cosmic magnetic fields with MASCLET: an application to galaxy clusters
• High-order Discontinuous Galerkin hydrodynamics with sub-cell shock capturing on GPUs
• Spectral Difference method with a posteriori limiting: Application to the Euler equations in one and two space dimensions
• Monolithic parabolic regularization of the MHD equations and entropy principles
• WENO-Wombat: Scalable Fifth-Order Constrained-Transport Magnetohydrodynamics for Astrophysical Applications
• A Discontinuous Galerkin Solver in the FLASH Multi-Physics Framework
• Supersonic turbulence simulations with GPU-based high-order Discontinuous Galerkin hydrodynamics
• High order well-balanced Arbitrary-Lagrangian-Eulerian ADER discontinuous Galerkin schemes on general polygonal moving meshes
• Spectral Difference method with a posteriori limiting: II- Application to low Mach number flows
• Numerical Methods for Simulating Star Formation
• DISPATCH methods: an approximate, entropy-based Riemann solver for ideal magnetohydrodynamics
• Constraint preserving discontinuous Galerkin method for ideal compressible MHD on 2-D Cartesian grids