Abstract

Contributed Talk - Splinter Computation

Tuesday, 19 September 2017, 14:30   (HS3)

Towards Fast High-Order Magnetohydrodynamics in the AREPO Code

T. Guillet^1,2, R. Pakmor², V. Springel^1,2
1Zentrum für Astronomie der Universität Heidelberg
²Heidelberg Institute for Theoretical Studies

I present our latest developments on MHD schemes for the cosmological hydrodynamical simulation code AREPO. The code implements finite volume solvers for hydrodynamics and MHD, and can operate on two major types of grids: either a fully dynamic quasi-Lagrangian unstructured Voronoi mesh, or an octree-based adaptive mesh refinement (AMR) grid. The moving mesh solver is based on a second-order unsplit Godunov scheme. For AMR grids, we implemented a higher-order discontinuous Garlerkin (DG) hydrodynamics solver, which we have extended to MHD based on locally divergence-free basis functions. I discuss some challenges to achieve a robust and efficient implementation of high-order methods for MHD that can be applied to astrophysical problems. I also highlight some key aspects of DG schemes for computational efficiency on modern CPU architectures.