A Discontinuous Galerkin Solver in the FLASH Multi-Physics Framework

TL;DR

This paper introduces a discontinuous Galerkin solver integrated into the FLASH multi-physics simulation framework, enabling advanced astrophysical simulations with validated performance and interactions with physics modules.

Contribution

The paper presents a novel DG solver module for FLASH, tailored for astrophysical multi-physics simulations, with detailed implementation and validation.

Findings

The DG solver is successfully integrated into FLASH.

Extensive validation confirms accuracy and stability.

The module interacts effectively with gravity and radiative transfer.

Abstract

In this paper, we present a discontinuous Galerkin solver based on previous work by Markert et al. (2021) 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.