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

[Jelena Petrović]

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
>>