Neutron star evolution by combining discontinuous Galerkin and finite volume methods
Feb 10, 2025Citations per year
Abstract: (arXiv)
We present here a new hybrid scheme that combines a discontinuous Galerkin (DG) method with finite volume (FV) and finite difference (FD) methods. The computational mesh is divided into smaller elements that touch but do not overlap. Like a pure DG method, our new hybrid scheme requires information exchange only at the surface of neighboring elements. This avoids the need for ghostzones that are usually many points deep in traditional FV implementations. Furthermore, unlike traditional FV implementations, that require information exchange between each element and its 26 surrounding neighbors on non-cuboid meshes, our new hybrid method exchanges information only between each element and its six nearest neighbors. Through this reduction in communication, we aim to retain the high scalability of DG when using large supercomputers. The goal is to use DG in elements with smooth matter fields and to fall back onto the more robust FV/FD method in elements that contain non-smooth shocks or star surfaces. For this we devise trouble criteria to decide whether an element should be evolved with DG or FV/FD. We use the Nmesh program to implement and test the new scheme. We successfully evolve various single neutron star cases. These include the challenging cases of a neutron star initially in an unstable equilibrium migrating to a stable configuration and a boosted neutron star. These cases are simulated for the first time here in full 3D with general relativistic hydrodynamics using DG methods. We also describe additional numerical methods, such as the limiters and the atmosphere treatment we need for our simulations.Note:
- 34 pages, 20 figures
References(24)
Figures(17)
- [1]
- [2]
- [3]
- [4]
- [5]
- [6]
- [6]
- [7]
- [8]
- [9]
- [10]
- [11]
- [12]
- [13]
- [14]
- [15]
- [16]
- [17]
- [18]
- [19]
- [20]
- [21]
- [22]
- [23]