Towards Exascale Simulations of the ICM Dynamo with WENO-Wombat
arXiv:1808.10633
Abstract
In galaxy clusters, modern radio interferometers observe non-thermal radio sources with unprecedented spatial and spectral resolution. For the first time, the new data allows to infer the structure of the intra-cluster magnetic fields on small scales via Faraday tomography. This leap forward demands new numerical models for the amplification of magnetic fields in cosmic structure formation - the cosmological magnetic dynamo. Here we present a novel numerical approach to astrophyiscal MHD simulations aimed to resolve this small-scale dynamo in future cosmological simulations. As a first step, we implement a fifth order WENO scheme in the new code WOMBAT. We show that this scheme doubles the effective resolution of the simulation and is thus less expensive than common second order schemes. WOMBAT uses a novel approach to parallelization and load balancing developed in collaboration with performance engineers at Cray Inc. This will allow us scale simulation to the exaflop regime and achieve kpc resolution in future cosmological simulations of galaxy clusters. Here we demonstrate the excellent scaling properties of the code and argue that resolved simulations of the cosmological small scale dynamo within the whole virial radius are possible in the next years.
References in corpus (7)
- Athena: A New Code for Astrophysical MHD
- Compressible Turbulence in Galaxy Clusters: Physics and Stochastic Particle Re-acceleration
- An Unsplit Godunov Method for Ideal MHD via Constrained Transport in Three Dimensions
- A Divergence-Free Upwind Code for Multidimensional Magnetohydrodynamic Flows
- Plasma instabilities and magnetic-field growth in clusters of galaxies
- A Validated Nonlinear Kelvin-Helmholtz Benchmark for Numerical Hydrodynamics
- Vorticity, Shocks and Magnetic Fields in Subsonic, ICM-like Turbulence