[Mesa-users] modeling a planet

[Francis Timmes]

hi nuno,

per https://docs.mesastar.org/en/release-r22.11.1/reference/star_job.html#relax-core 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.

fxt




> On Feb 14, 2023, at 5:01 AM, Nuno Moedas via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> 
> Hi,
> 
> 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 error.
> 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.
> 
> Sincerely,
> Nuno Moedas<inlist_creat_planet_mod.f><inlist_core_planet_mod.f><inlist_evolve_planet_mod.f>_______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
>