The new discontinuous Galerkin methods based numerical relativity program Nmesh
The new discontinuous Galerkin methods based numerical relativity program Nmesh
复制标题
DOI:
10.1088/1361-6382/acaae7
复制
发表时间:
2023-01-19
影响因子:
3.5
通讯作者:
Pirog,Michal
中科院分区:
文献类型:
--
作者:
Tichy,Wolfgang;Ji,Liwei;Pirog,Michal
Interpreting gravitational wave observations and understanding the physics of astrophysical compact objects such as black holes or neutron stars requires accurate theoretical models. Here, we present a new numerical relativity computer program, called Nmesh, that has the design goal to become a next generation program for the simulation of challenging relativistic astrophysics problems such as binary black hole or neutron star mergers. In order to efficiently run on large supercomputers, Nmesh uses a discontinuous Galerkin method together with a domain decomposition and mesh refinement that parallelizes and scales well. In this work, we discuss the various numerical methods we use. We also present results of test problems such as the evolution of scalar waves, single black holes and neutron stars, as well as shock tubes. In addition, we introduce a new positivity limiter that allows us to stably evolve single neutron stars without an additional artificial atmosphere, or other more traditional limiters.