[Users] Conservative to primative conversion fails in TOV collapse

[Shoup, Anthony]

Hi Again,

Just to let you know, I changed my resolution from 3.2 (cell size) to 1.6 and the problem went away. I will now run the full simulation and see if I can match that paper's Figure 18.

Thanks, Tony...

[https://email.osu.edu/owa/attachment.ashx?id=RgAAAAAb%2fHy0wVvTSoHQx8OJXAaLBwCiA5IZrwRKTqiLVNbt4xWyAAAAAAFUAADltPc25wRDT4tJbW9en2wXAHKArVn%2fAAAJ&attcnt=1&attid0=EAD8Cse5Lj5uQ6ZJWk98Q%2blj]
Anthony Shoup PhD, Senior Lecturer
College of Arts & Sciences, College of Engineering Departments of Physics, Astronomy, EEIC
315 Science Bldg. | 4250 Campus Dr. Lima, OH 45807
419-995-8018 Office | 419-516-2257 Mobile
shoup.31 at osu.edu<https://email.osu.edu/owa/redir.aspx?C=j5WpnJiBk0W5oVlCbtvB-xiCkA_lbdEIi9hlk7ByHiG7ARrxjwDFmAW8S_XespJbMLJRblY5JKc.&URL=mailto%3ashoup.31%40osu.edu> osu.edu<https://email.osu.edu/owa/redir.aspx?C=j5WpnJiBk0W5oVlCbtvB-xiCkA_lbdEIi9hlk7ByHiG7ARrxjwDFmAW8S_XespJbMLJRblY5JKc.&URL=http%3a%2f%2fosu.edu>
________________________________
From: users-bounces at einsteintoolkit.org <users-bounces at einsteintoolkit.org> on behalf of Shoup, Anthony <shoup.31 at osu.edu>
Sent: Friday, March 15, 2019 10:01 PM
To: Einstein Toolkit Users
Subject: [Users] Conservative to primative conversion fails in TOV collapse

Hi All,

I have been working on reproducing the results in the paper:

The Einstein Toolkit: A Community Computational
Infrastructure for Relativistic Astrophysics  (arXiv: 1111.3344v1 {gr-qc] 14 Nov 2011)

and others as a way to teach myself how to use the Einstein Toolkit and to test out my home-built cluster computer.

I am currently working on the non-rotating TOV neutron star collapse into a BH.  It works well up until it fails trying to convert conserved variables into primitive variables in grhydro (see attached tovCollapse.out file).  I have tried adjusting the atmosphere parameters (specifically GRHydro::rho_abs_min  = 1.e-8, up from 1.e-10)., as per the suggestions in the grhydro documentation, which has helped, but I think the conversion is failing now very near where the AH will/is forming.

I have included my parameter file for the simulation (I lowered evolve K from 98 to 80 just to accelerate the collapse process).

Any help anyone can give would be greatly appreciated.

Sincerely, Tony...

[https://email.osu.edu/owa/attachment.ashx?id=RgAAAAAb%2fHy0wVvTSoHQx8OJXAaLBwCiA5IZrwRKTqiLVNbt4xWyAAAAAAFUAADltPc25wRDT4tJbW9en2wXAHKArVn%2fAAAJ&attcnt=1&attid0=EAD8Cse5Lj5uQ6ZJWk98Q%2blj]
Anthony Shoup PhD, Senior Lecturer
College of Arts & Sciences, College of Engineering Departments of Physics, Astronomy, EEIC
315 Science Bldg. | 4250 Campus Dr. Lima, OH 45807
419-995-8018 Office | 419-516-2257 Mobile
shoup.31 at osu.edu<https://email.osu.edu/owa/redir.aspx?C=j5WpnJiBk0W5oVlCbtvB-xiCkA_lbdEIi9hlk7ByHiG7ARrxjwDFmAW8S_XespJbMLJRblY5JKc.&URL=mailto%3ashoup.31%40osu.edu> osu.edu<https://email.osu.edu/owa/redir.aspx?C=j5WpnJiBk0W5oVlCbtvB-xiCkA_lbdEIi9hlk7ByHiG7ARrxjwDFmAW8S_XespJbMLJRblY5JKc.&URL=http%3a%2f%2fosu.edu>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20190316/02d260dc/attachment.html