High-Order ADER-DG Hydrodynamics with ExaHyPE: Implementation, Validation, and Astrophysical Benchmarking
Abstract
We describe a high-order ADER-DG solver for the compressible Euler equations within the ExaHyPE framework. The implementation combines a high-order ADER-DG polynomial representation, a local space-time DG predictor, adaptive mesh refinement, and an a posteriori subcell finite-volume limiter. We test the code on a deliberately mixed set of one- and two-dimensional problems: a strong-shock Sod-type problem, the Shu-Osher shock-entropy interaction, the Woodward-Colella blast wave, a contact-driven vortex sheet, and a shock-interface interaction. The one-dimensional cases recover the expected Euler wave patterns and show clear order-dependent gains in smooth and oscillatory regions. The two-dimensional cases probe a different part of the method, namely contact preservation, shear-driven roll-up, baroclinic vorticity deposition, and Richtmyer-Meshkov-type growth. In these tests the high-order update gives the expected resolution away from discontinuities, whereas the subcell limiter keeps the calculation stable near shocks and steep interfaces. The resulting code provides a reproducible ExaHyPE implementation for idealised inviscid, non-relativistic flows in which shocks, contacts, and multidimensional interfaces are the dominant features.
Cite
@article{arxiv.2605.17132,
title = {High-Order ADER-DG Hydrodynamics with ExaHyPE: Implementation, Validation, and Astrophysical Benchmarking},
author = {Andrés Mauricio Suárez Mantilla and Leonardo Castañeda Colorado},
journal= {arXiv preprint arXiv:2605.17132},
year = {2026}
}
Comments
38 pages, 8 figures, benchmark datasets and implementation resources publicly available