Astrophysical hydrodynamics with a high-order discontinuous Galerkin scheme and adaptive mesh refinement

Abstract

Solving the Euler equations of ideal hydrodynamics as accurately and efficiently as possible is a key requirement in many astrophysical simulations. It is therefore important to continuously advance the numerical methods implemented in current astrophysical codes, especially also in light of evolving computer technology, which favours certain computational approaches over others. Here we introduce the new adaptive mesh refinement (AMR) code TENET, which employs a high order discontinuous Galerkin (DG) scheme for hydrodynamics. The Euler equations in this method are solved in a weak formulation with a polynomial basis by means of explicit Runge-Kutta time integration and Gauss-Legendre quadrature. This approach offers significant advantages over commonly employed second order finite volume (FV) solvers. In particular, the higher order capability renders it computationally more efficient, in the sense that the same precision can be obtained at significantly less computational cost. Also, the DG scheme inherently conserves angular momentum in regions where no limiting takes place, and it typically produces much smaller numerical diffusion and advection errors than a FV approach. A further advantage lies in a more natural handling of AMR refinement boundaries, where a fall-back to first order can be avoided. Finally, DG requires no wide stencils at high order, and offers an improved data locality and a focus on local computations, which is favourable for current and upcoming highly parallel supercomputers. We describe the formulation and implementation details of our new code, and demonstrate its performance and accuracy with a set of two- and three-dimensional test problems. The results confirm that DG schemes have a high potential for astrophysical applications.

23 pages, 12 figures, published in MNRAS November 2015, a movie may be accessed online: https://youtu.be/cTRQP6DSaqA

References in corpus (5)

• PLUTO: a Numerical Code for Computational Astrophysics
• Athena: A New Code for Astrophysical MHD
• Improving the convergence properties of the moving-mesh code AREPO
• Inviscid SPH
• Solving the relativistic magnetohydrodynamics equations with ADER discontinuous Galerkin methods, a posteriori subcell limiting and adaptive mesh refinement

Cited by in corpus (29)

• The Arepo public code release
• Hydrodynamics of core-collapse supernovae and their progenitors
• SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics
• A Validated Nonlinear Kelvin-Helmholtz Benchmark for Numerical Hydrodynamics
• Well balanced Arbitrary-Lagrangian-Eulerian finite volume schemes on moving nonconforming meshes for the Euler equations of gasdynamics with gravity
• Grand challenges in protoplanetary disc modelling
• Formulation of discontinuous Galerkin methods for relativistic astrophysics
• High-order Magnetohydrodynamics for Astrophysics with an Adaptive Mesh Refinement Discontinuous Galerkin Scheme
• Entropy Stable Discontinuous Galerkin Methods for Balance Laws in Non-Conservative Form: Applications to the Euler Equations with Gravity
• The new discontinuous Galerkin methods based numerical relativity program Nmesh
• CosmosDG: An hp-Adaptive Discontinuous Galerkin Code for Hyper-Resolved Relativistic MHD
• Improved Performances in Subsonic Flows of an SPH Scheme with Gradients Estimated using an Integral Approach
• Simulating magnetized neutron stars with discontinuous Galerkin methods
• A spectral deferred correction method for incompressible flow with variable viscosity
• The role of physical and numerical viscosity in hydrodynamical instabilities
• thornado-hydro: a discontinuous Galerkin method for supernova hydrodynamics with nuclear equations of state
• A hp-adaptive discontinuous Galerkin solver for elliptic equations in numerical relativity
• Planet-disc interactions with Discontinuous Galerkin Methods using GPUs
• High-order Discontinuous Galerkin hydrodynamics with sub-cell shock capturing on GPUs
• Python Radiative Transfer Emission code (PyRaTE): non-LTE spectral lines simulations
• Space-time adaptive ADER-DG finite element method with LST-DG predictor and a posteriori sub-cell ADER-WENO finite-volume limiting for multidimensional detonation waves simulation
• Multidimensional upwind hydrodynamics on unstructured meshes using Graphics Processing Units I. Two-dimensional uniform meshes
• 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
• Future Simulations of Tidal Disruption Events
• Star Formation for Predictive Primordial Galaxy Formation
• Spectral Difference method with a posteriori limiting: II- Application to low Mach number flows
• Neutron star evolution by combining discontinuous Galerkin and finite volume methods