Computational Relativistic Astrophysics With Adaptive Mesh Refinement: Testbeds

TL;DR

This paper demonstrates the use of adaptive mesh refinement in numerical simulations of neutron star binary systems, achieving high accuracy comparable to large unigrid simulations on standard workstations.

Contribution

It introduces AMR techniques for relativistic astrophysics simulations, enabling high-accuracy modeling of neutron star mergers on less powerful hardware.

Findings

AMR simulations match $1025^3$ unigrid accuracy

First high-resolution NS binary inspiral simulation on a workstation

Potential for broader applications in relativistic astrophysics

Abstract

We have carried out numerical simulations of strongly gravitating systems based on the Einstein equations coupled to the relativistic hydrodynamic equations using adaptive mesh refinement (AMR) techniques. We show AMR simulations of NS binary inspiral and coalescence carried out on a workstation having an accuracy equivalent to that of a $1025^3$ regular unigrid simulation, which is, to the best of our knowledge, larger than all previous simulations of similar NS systems on supercomputers. We believe the capability opens new possibilities in general relativistic simulations.