Computational Relativistic Astrophysics With Adaptive Mesh Refinement: Testbeds

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. AMR coalescences of neutron stars can now be simulated with sufficient resolution covering the neutron stars while having the computational domain extend to the local wave zone. We show an AMR simulation carried out with a workstation having an accuracy equivalent to that of a 1025^3 regular uni-grid 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