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

[Jelena Petrović]

Thanks both Francis and Rob, I will try this,

Best,
Jelena


On 2020-09-13 21:07, Francis Timmes wrote:
> 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