[Mesa-users] modeling a planet
Nuno.Moedas at astro.up.pt
Thu Feb 16 16:55:18 UTC 2023
Thank you for your answer. I will try and see what I obtain.
A 2023-02-15 18:00, Francis Timmes escreveu:
> hi nuno,
> you can experiment with significantly smaller values of
> dlg_core_mass_per_step than the default 1e-3.
>> ... where the planet only is a core and does not present any mass in
>> the envelope.
> see section 2.3 of mesa II. the core is inert and not presently
> evolved in any way. this means a pure core and no envelope will not be
> possible as there is no material to evolve.
> it can be useful to make some plots of key quantities as the core
> fraction increases. this often suggests paths forward when
> non-convergence is reached in a sequence of models.
>> On Feb 14, 2023, at 5:01 AM, Nuno Moedas via Mesa-users
>> <mesa-users at lists.mesastar.org> wrote:
>> I am trying to model a planet of about 2 Jupiter masses using
>> mesa-r12778. I used as starting point the test suit that is available
>> in MESA. I was able to create the model and add the core and evolve
>> it. Now I am trying to create a grid where a change the mass of the
>> core from 0 until reaches the total mass of the planet. However when I
>> reach more or less half the mass of the planet in the core the model
>> does not converge in the step of adding the core. I was trying to
>> change the parameters of dlg_core_mass_per_step and
>> relax_core_years_for_dt, but it seems that still gives me the same
>> stopping because of problems num_backups > max_number_backups
>> This is my first time modelling a planet. I want to reach a point
>> where the planet only is a core and does not present any mass in the
>> envelope. Now I don't know what to do to pass half of the mass of the
>> planet on the core, or if I missing any command in it.
>> I send also the inlists that I use to create the planet, add the core
>> of the planet, and evolve the planet.
>> Thank you.
>> mesa-users at lists.mesastar.org
More information about the Mesa-users