A discontinuous Galerkin solver in theflashmultiphysics framework

Abstract

In this paper, we present a discontinuous Galerkin solver based on previous work by the authors for magneto-hydrodynamics in form of a new fluid solver module integrated into the established and well-known multi-physics simulation code FLASH. Our goal is to enable future research on the capabilities and potential advantages of discontinuous Galerkin methods for complex multi-physics simulations in astrophysical settings. We give specific details and adjustments of our implementation within the FLASH framework and present extensive validations and test cases, specifically its interaction with several other physics modules such as (self-)gravity and radiative transfer. We conclude that the new DG solver module in FLASH is ready for use in astrophysics simulations and thus ready for assessments and investigations.

High-order discontinuous Galerkin hydrodynamics with sub-cell shock capturing on GPUs

“…Despite these advantages, DG methods have only recently begun to be considered in astrophysics. First implementations and applications include Mocz et al ( 2014 ); Schaal et al ( 2015 ); Kidder et al ( 2017 ); Velasco Romero et al ( 2018 ); Guillet et al ( 2019 ), as well as more recently Lombart & Laibe ( 2021 ); Deppe et al ( 2022 ); Markert, Walch & Gassner ( 2022 ). We here focus on exploring a new implementation of DG that we developed from the ground up for use with graphical processing units (GPUs).…”