This simulation shows how to evolve a stable, single neutron star, and is a low resolution example of the case presented on our release paper; see Fig. 15 of https://arxiv.org/abs/1111.3344. It uses a TOV solution as initial data and follows the initial (numerical) perturbations, writing the maximum density to disk. For the test described here, we set up a self-gravitating fluid sphere described by a polytropic equation of state \(P = K \rho^\Gamma\) with K = 100 and \(\Gamma = 2\), and an initial central density of \(\rho = 1.28 \times 10^{-3}\). This model can be taken to represent a non-rotating NS with a mass of \(M = 1.4 M_\odot\). The initial spike is due to the perturbation of the solution resulting from the interpolation onto the evolution grid. The remaining oscillations are mainly due to the interaction of the star and the artificial atmosphere and are present during the whole evolution.
Parameter file | tov_gallery_polytrop.par |
---|---|
Thornlist | ET Riemann release (ET_2022_05) |
Setup | 10 cores; 2 threads |
approx. runtime | 45 minutes |
Results (50MB) | tov-20220512.tar.gz |
The plot of the maximum of the density can be plotted using kuibit using this script.