In astrophysics, the two main methods traditionally in use for solving the
Euler equations of ideal fluid dynamics are smoothed particle hydrodynamics and
finite volume discretization on a stationary mesh. However, the goal to
efficiently make use of future exascale machines with their ever higher degree
of parallel concurrency motivates the search for more efficient and more
accurate techniques for computing hydrodynamics. Discontinuous Galerkin (DG)
methods represent a promising class of methods in this regard, as they can be
straightforwardly extended to arbitrarily high order while requiring only small
stencils. Especially for applications involving comparatively smooth problems,
higher-order approaches promise significant gains in computational speed for
reaching a desired target accuracy. Here, we introduce our new astrophysical DG
code TENET designed for applications in cosmology, and discuss our first
results for 3D simulations of subsonic turbulence. We show that our new DG
implementation provides accurate results for subsonic turbulence, at
considerably reduced computational cost compared with traditional finite volume
methods. In particular, we find that DG needs about 1.8 times fewer degrees of
freedom to achieve the same accuracy and at the same time is more than 1.5
times faster, confirming its substantial promise for astrophysical
applications.Comment: 21 pages, 7 figures, to appear in Proceedings of the SPPEXA
symposium, Lecture Notes in Computational Science and Engineering (LNCSE),
Springe