[Mesa-users] ppisn: evolve >150 Ms massive stars in 15140

[Jing-Ze Ma]

Hi mesa users,


I am having trouble evolving massive stars with mass above 150 Ms past He main sequence in ppisn test_suite. During relax_to_radiative_core or helium burning, I keep getting retries like


 retry: max residual jumped -- give up in solver
 retry: logT < hydro_mtx_min_allowed_logT
 retry: logRho > hydro_mtx_max_allowed_logRho


which slow things down and make it hard to get past He main sequence stage or even enter the main sequence stage. I am keeping everything as defaults in ppisn test_suite except with rotation_flag turned off. I am not sure whether this is caused by different mesa versions or not, since I am using 15140.




A possible solution I am experimenting is to add


 okay_to_reduce_gradT_excess = .true.
 gradT_excess_lambda1 = -1


and then turn off okay_to_reduce_gradT_excess when hydro is on, which would evolve a 150 Ms to the end.


By further adding


 change_mass_years_for_dt = 1d-2


I manage to evolve a 200 Ms to the end.
But for a 250 Ms model, I have not found a satisfying solution yet. Tweaking pre_ms_T_c or mesh_delta_coeff or turning off conv_vel_flag is not helping much.


I suppose that either this is a new problem from new version, or this has been solved in Marchant 2019 paper but not written in ppisn test_suite. Any suggestion would be deeply appreciated.


Thanks in advance!


Cheers,
Jing-Ze
 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210521/c4dc5430/attachment.htm>