Gallery: Single stable neutron star

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
This example was last tested on 12-May-2022.


The plot of the maximum of the density can be plotted using kuibit using this script.

Maximum density over time