[Mesa-users] TOV approximation, other_cgrav

[Bill Paxton]

> On Nov 26, 2020, at 10:51 AM, Nicholas Herrington via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> 
> Hi,
> 
> From a previous thread related to modelling a star with a neutron star in its core, I saw Bill provided a other_cgrav routine for the TOV approximation. I have been using the other_cgrav routine to implement the TOV approximation and looking at the other_cgrav.inc provided, I wanted to ask what is the significance of using face centred pressure and density instead of cell centred (I used cell centred variables in my version).

i have no idea what is or is not significant for your application!   But it is a fact of mesa that the momentum equation that uses cgrav applies to cells faces where r and m and cgrav are defined.   Pressures are defined at cell “centers” (actually they are cell average values, but the distinction is blurred in mesa).  That’s convenient for the momentum equation since it uses the pressure difference across the face.



> Additionally the version of the TOV corrections to cgrav in the script follows Gnew = G*f1*f2/f3, would you be able to explain where this comes from? I am not knowledgeable in the solutions to arrive at G, I have been following the correction from https://ui.adsabs.harvard.edu/abs/2020MNRAS.494.2236W/abstract <https://ui.adsabs.harvard.edu/abs/2020MNRAS.494.2236W/abstract> (equation 11 in section 2).

i’m no expert on GR or TOV either, and that’s why this is in a user-accessible hook so you can fix it!

the comments in other_cgrav.inc in src directories for the current, pre-release version, include comments that summarize my entire (lack of) understanding of this issue.

      ! use the Tolman–Oppenheimer–Volkoff (TOV) equation. 
      ! See first equation in https://en.wikipedia.org/wiki/Tolman%E2%80%93Oppenheimer%E2%80%93Volkoff_equation.  
      ! want to replace -G*m/r^2 by -G*m/r^2*(1 + P/(rho c^2))(1 + 4 pi r^3 P /(m c^2))/(1 - 2 G m/(r c^2))

Those factors in ()’s that multiply gravity in TOV in the comment simply multiply standard_cgrav in the code.
	    G = standard_cgrav
            f1 = 1d0 + P/(rho*clight**2)
            f2 = 1d0 + 4d0*pi*r**3*P/(m*clight**2)
            f3 = 1d0 - 2d0*G*m/(r*clight**2)
            s% cgrav(k) = G*f1*f2/f3

-bill







> 
> Thanks
> Nick
> 
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20201126/de624a0a/attachment.htm>