[Mesa-users] Test Suite: 1M_pre_ms_to_wd

[Aaron Dotter]

Hi Erica,

Are you talking about what you showed in the plot a couple of emails back?
If so, then what you are seeing is the evolution from the early pre-main
sequence to the main sequence.

Aaron

On Thu, Nov 18, 2021 at 3:59 PM Erica A Sawczynec <ericasaw at utexas.edu>
wrote:

> Hi Aaron,
> Even with commenting out the max_model_number in the inlist I am still
> having the same issue. The star appears to have reached the end of the MS
> under 1 Gyr. Any other ideas?
> Thanks,
> Erica
>
> On Nov 18, 2021, at 2:34 PM, Aaron Dotter via Mesa-users <
> mesa-users at lists.mesastar.org> wrote:
>
> Hi Erica,
>
> With a max_years_for_timestep = 1d5 and max_model_number = 6500 your
> simulation will end with an age of ~650 Myr.
>
> That does *not* mean it has reached any specific evolutionary point, just
> that it has reached the max_model_number and that is one criterion for MESA
> to terminate the evolution.
>
> A 1 Msun model (at solar Z) should have a main sequence lifetime of ~10
> Gyr so you would need about 10^5 timesteps to reach the end of the main
> sequence with a max_years_for_timestep = 1d5.
>
> Aaron
>
> On Thu, Nov 18, 2021 at 3:24 PM Erica A Sawczynec <ericasaw at utexas.edu>
> wrote:
>
>> Hi Aaron,
>> Will this cause the star to evolve to the MS turn-off in 6500 steps
>> regardless of if it should be there or not? I tried changing this value to
>> be bigger and ran into the same issue. Here is a plot I made using TULIPS
>> of the radius as a function of time. You can see the star is at the MS
>> turn-off at 638 Myr. Trying a run with turning that off now!
>> Thanks,
>> Erica
>>
>> <600Myr_turnoff.png>
>>
>> On Nov 18, 2021, at 2:17 PM, Aaron Dotter <aaron.dotter at gmail.com> wrote:
>>
>> Hi Erica,
>>
>> You have max_model_number = 6500 in the inlist.  Try removing this from
>> the inlist (or set to -1).
>>
>> Aaron
>>
>> On Thu, Nov 18, 2021 at 3:08 PM Erica A Sawczynec <ericasaw at utexas.edu>
>> wrote:
>>
>>> Hi All!
>>> I have been tinkering with the 1M_pre_ms_to_wd test suite to try to look
>>> at the depth of the convection zone boundary as a function of opacity. I
>>> have an altered version that no longer uses the to_end_wd inlist and I set
>>> the envelope mass limit so the to_end_agb inlist doesn’t actually go too
>>> far past the MS turn-off. I want to finely sample the  MS so I tried
>>> turning on the max_years_for_timestep keyword in the debugging section. I
>>> set it to something small like 1d5, which feels reasonable for the MS, but
>>> I’ve noticed whenever I do this, the test suite evolves a 1M_sol star to
>>> the end of the MS in ~600 Myr (which isn’t possible).
>>> I thought that the max_years_for_timestep would keep the time to evolve
>>> to the MS turn-off the same, but increase the sampling of the MS, is that
>>> not the case? I even tried increasing the mass model number (up to 65000),
>>> thinking maybe that imposed some sort of restriction on the model, but even
>>> then the age of the star at the end of the run is way too young. I’ve
>>> attached my changed inlist below, note that it will “fail” with the "failed
>>> to create end_agb.mod when running inlist_to_end_agb_header” error, but
>>> there is still a history.data file where you can look at the results of the
>>> run.
>>> I cannot help but feeling I am missing something silly here, but I just
>>> can’t figure out what it is. Any help is greatly appreciated!
>>> Thanks,
>>> Erica
>>>
>>> _______________________________________________
>>> mesa-users at lists.mesastar.org
>>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>>>
>>>
>> _______________________________________________
> 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/20211118/91512c72/attachment.htm>