[Mesa-users] Massive star evolution after CO core formation

[Francis Timmes]

hi jelana,

thanks. this is a learning opportunity, so i am going to echo rob's suggestion 

> I would suggest setting up the pgstar inlist from
> https://zenodo.org/record/3986828 which will then contain alot of
> information about your star ...

one can see this pgstar inlist in action at

https://www.youtube.com/watch?v=JZFa4WURztI

while this alone isn't going to directly alleviate the issues with your model,
it will give you visual information on what your model is doing (or not).
this can suggest (along with the accompanying project_inlist) things to improve upon.

fxt





> On Sep 10, 2020, at 11:16 PM, Jelena Petrović <jpetrovic at aob.rs> wrote:
> 
> Hi Francis,
> 
> Below are the inlist_project and inlist for ~10Ms star,
> 
> Thanks!
> Jelena
> --------------------------
> &star_job
> 
>  ! begin with a pre-main sequence model
>    create_pre_main_sequence_model = .true.
> 
>  ! save a model at the end of the run
>    save_model_when_terminate = .true.
>    save_model_filename = '10M_at_TAMS.mod'
> 
>  ! display on-screen plots
>    pgstar_flag = .false.
> 
> / !end of star_job namelist
> 
> 
> &controls
> 
>  ! starting specifications
>    initial_mass = 10.710456d0 ! in Msun units
>    initial_y = 0.98d0
>    initial_z = 0.02d0
> 
>     hot_wind_scheme = 'Dutch'
>     Dutch_scaling_factor = 0d0
> 
>     profile_interval = 100
>     max_num_profile_models = 10000
> 
>    photo_interval = 500
>    photo_digits = 8
> / ! end of controls namelist
> 
> ------------------------
> 
> &star_job
> 
>    read_extra_star_job_inlist1 = .true.
>    extra_star_job_inlist1_name = 'inlist_project'
>    create_pre_main_sequence_model = .true.
> / ! end of star_job namelist
> 
> 
> &controls
> 
>    read_extra_controls_inlist1 = .true.
>    extra_controls_inlist1_name = 'inlist_project'
>    initial_y = 0.98d0
>    initial_z = 0.02d0
> 
> / ! end of controls namelist
> 
> 
> &pgstar
> 
>    read_extra_pgstar_inlist1 = .false.
> !   extra_pgstar_inlist1_name = 'inlist_pgstar'
> 
> / ! end of pgstar namelist
> -------------------------
> 
> 
> On 2020-09-11 05:48, Francis Timmes wrote:
>> hi jelena,
>> what is your inlist?
>> fxt
>>> On Sep 10, 2020, at 7:35 AM, Jelena Petrović <jpetrovic at aob.rs> wrote:
>>> Hi,
>>> Sorry, I forgot, the string swaps between 'newton iters' and 'max increase'.
>>> Many thanks,
>>> Jelena
>>> On 2020-09-10 14:51, Rob Farmer wrote:
>>>> Hi
>>>> Firstly, the terminal output contains a string in the bottom right
>>>> corner that tells you what is currently setting mesa's timestep. You
>>>> can then search for that string in the controls.default file to find
>>>> an appropriate control that you could tweak (or atleast use to
>>>> understand what physics is setting the current limit).
>>>> It's hard to say what is limiting your timesteps without seeing any
>>>> plots. I would suggest setting up the pgstar inlist from
>>>> https://zenodo.org/record/3986828 which will then contain alot of
>>>> information about your star and will make it much much easier to debug
>>>> why a problem is occurring.
>>>> Rob
>>>> On Thu, 10 Sep 2020 at 11:56, Jelena Petrović <jpetrovic at aob.rs>
>>>> wrote:
>>>>> Dear all,
>>>>> I am evolving massive stars (above 30Ms) and I have noticed that
>>>>> most of
>>>>> my models have problems evolving further than CO core,
>>>>> time steps get really small and models do not progress anymore.
>>>>> Only sometimes I manage to get the evolution all the way to Fe core.
>>>>> It
>>>>> seems that lower mass models do slightly better.
>>>>> I though this is because of extended hydrogen envelope, so I have
>>>>> tried
>>>>> to evolve helium stars. However, the same happens.
>>>>> When the core is about 25% C and 75% O, the models run, but with
>>>>> extremely small time steps.
>>>>> Just to note that stellar wind is also switched off to ensure that
>>>>> this
>>>>> is not the reason for smaller time steps.
>>>>> Does anybody know what could be the problem?
>>>>> Many thanks!
>>>>> Jelena
>>>>> ----------------------------------------
>>>>> Dr. Jelena Petrovic
>>>>> Senior Research Associate
>>>>> Astronomical Observatory Belgrade
>>>>> Serbia
>>>>> ----------------------------------------
>>>>> _______________________________________________
>>>>> mesa-users at lists.mesastar.org
>>>> https://eur04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.mesastar.org%2Fmailman%2Flistinfo%2Fmesa-users&data=02%7C01%7Cr.j.farmer%40uva.nl%7C0f0fabf818fc486e188c08d8556fe482%7Ca0f1cacd618c4403b94576fb3d6874e5%7C1%7C0%7C637353286344198406&sdata=xzV97QRzUQmqWXDcPfcQSQDW1YMBLdMckeTPdovtyOE%3D&reserved=0
>>> _______________________________________________
>>> mesa-users at lists.mesastar.org
>>> https://lists.mesastar.org/mailman/listinfo/mesa-users