[Mesa-users] ppisn: neutrino luminosity as a threshold for turning on Riemann solver

Jing-Ze Ma mjz18 at mails.tsinghua.edu.cn
Sat May 8 02:47:12 UTC 2021

Hi Francis and Pablo,

Thank you for your comments! So far I've tried a 60 solar mass and a 90 solar mass model with x_ctrl(3) = 12.5d0, and their final fates show no significant difference than the runs with x_ctrl(3) = 9d0. But I'll try more runs to see if it works for all masses.


-----Original Messages-----
From:"Pablo Marchant" <pamarca at gmail.com>
Sent Time:2021-05-07 01:13:15 (Friday)
To: "Francis Timmes" <fxt44 at mac.com>
Cc: "Jing-Ze Ma" <mjz18 at mails.tsinghua.edu.cn>, mesa-users <mesa-users at lists.mesastar.org>
Subject: Re: [Mesa-users] ppisn: neutrino luminosity as a threshold for turning on Riemann solver

And just as a quick comment, many of those options that regulate the turning on/off of hydro in the ppisn test case have evolved to attempt to make the test as reliable as possible. In my experience at least, shifting the threshold for neutrino luminosity too high can cause problems in pulses after the first one, where the star can be very bloated owing to the energy injection in the first pulse. Without hydro, the envelope of the model can struggle and drive timesteps down.

Still, would be pretty interesting for someone to study more in detail how much these numerical options affect the outcome. I did a few tests and did not find significant differences, other than a varying reliability of the models. But still, I never studied it systematically.


On Thu, May 6, 2021 at 6:58 PM Francis Timmes via Mesa-users <mesa-users at lists.mesastar.org> wrote:

dear jing-ze,

this is an excellent experiment for you to run!
then you'll be able to quantify the differences.


> On May 5, 2021, at 10:38 PM, Jing-Ze Ma via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> Hi everyone,
> I've been experimenting with the ppisn test_suite. I noticed that the Riemann solver will be turned on when neutrino luminosity reaches a certain value before dynamical instability. In Merchant 2019 paper it's set to be 10d0 (10^10 solar luminosity), whereas in the test_suite it's 9d0. If I want to delay the turning on of Riemann solver and set this value to be 12d0, will the results be much different (e.g., give very different BH masses or make a PPISN event become a direct iron core collapse)?
> Thank you in advance.
> Cheers,
> Jing-Ze
> --
> Jing-Ze Ma
> Undergraduate Student,
> Tsinghua University
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users

mesa-users at lists.mesastar.org


Pablo Marchant Campos
M.Sc on Astrophysics, Universidad Católica de Chile
PhD on Astrophysics, Argelander-Institut für Astronomie, Universität Bonn
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210508/2ccae57b/attachment.htm>

More information about the Mesa-users mailing list