[Mesa-users] [mesa-users] any recommendation for running 200 Msun at Zsun?

Fri Aug 25 12:44:52 UTC 2023

Thanks Ebraheem for the hint. Idk what's about varcontrol_targe, but it
worked like it was some sort of black magic. Super useful!

Just a follow-up question here: I encountered in some models that I cannot
get out of the pre-MS because I get (inlist attached)

retry: logT < hydro_mtx_min_allowed_logT    1184
                                            num_retries         301
                                     max_number_retries         300

 stopping because of problems -- too many retries
 failed in relax_to_radiative_core
 failed in do_relax_to_radiative_core
 star_create_pre_ms_model ierr          -1
 do_load1_star ierr          -1
 before_evolve_loop ierr

Do you guys know what's up with this?

Cheers,

Amedeo

On Thu, 10 Aug 2023 at 17:06, Ebraheem Farag <ekfarag at asu.edu> wrote:

> Pextra_factor =2 allows Natalie's model to work at 1000 Msun.
>
> For you, Amedeo (I would also rec. Pextra_factor =2), but reduce
> varcontrol_target = 1d-4 -> 1d-3, and reduce your 'mesh_delta_coeff' and
> 'time_delta_coeff' back to default. For some reason varcontrol_target =
> 1d-4 can cause very massive models to evolve extremely slowly. I prefer to
> relax this control and use other timestep controls on the evolution of the
> core to ensure convergence:
>
> like time_delta_coeff = <1d0, or
>
>       !          'lgRho'                delta_lgRho_limit
>
>       !          'lgRho_cntr'           delta_lgRho_cntr_limit
>
>       !          'lgT'                  delta_lgT_limit
>
>       !          'lgT_cntr'             delta_lgT_cntr_limit
>
>       !          'lgT_max'              delta_lgT_max_limit
>
>
> p.s. be careful mixing overshoot with convective_premix.
>
>
> -EbF
>
>
> On Thu, Aug 10, 2023 at 7:30 AM Amedeo Romagnolo <
> amedeoromagnolo at gmail.com> wrote:
>
>> Hi Ebraheem and all of others,
>>
>> Yeah, the Pextra factor was one of the first things I did. No workaround
>> there, unfortunately...
>>
>> One of the things I noticed now that I implemented some of the things I
>> was recommended is that I get this error, which I cannot find the cause the
>> email archive, despite I didn't implement any new EOS thing at all. Any
>> recommendations here?
>>
>> failed while reading
>>> /work/chuck/aromagnolo/MESA/mesa/data/eosFreeEOS_data/mesa-FreeEOS_00z00x.data
>>>  Load1_eosDT_Table ierr          -1           1           1
>>>           NaN                       NaN
>>>  Load1_eosDT_Table ierr        5014           1           1
>>>           NaN                       NaN
>>>
>>
>> On Wed, 9 Aug 2023 at 19:17, Ebraheem Farag <ekfarag at asu.edu> wrote:
>>
>>> Hello Amedeo and Natalie,
>>>
>>> Your models are likely super-eddington at the surface.
>>> Try adopting "Pextra_factor = 2" to stabilize the surface.
>>>
>>> Let me know if this helps.
>>>
>>> -EbF
>>>
>>> On Wed, Aug 9, 2023 at 4:11 AM Natalie Rees via Mesa-users <
>>> mesa-users at lists.mesastar.org> wrote:
>>>
>>>> Hi Amedeo, Jared and all MESA users,
>>>>
>>>> I thought I would jump on this thread because I have a very similar
>>>> issue. I too am trying to evolve very massive stars at Zsun (0.02) up to
>>>> 1000 Msun. I had no problem doing so when I was using version 15140 with
>>>> use_superad_reduction however now I have updated to 23.05.1 my models with
>>>> M> 200 Msun can't evolve through the MS within 5000 timesteps. In
>>>> particular the 1000 Msun only reaches Xc = 0.6 after 5000 timesteps even
>>>> when superad_reduction_diff_grads_limit = 1d-4. I am aware that a lot has
>>>> changed between these two versions but any insight in to what I can do to
>>>> get these massive stars running faster again would be greatly appreciated.
>>>> I have attached my inlists also.
>>>>
>>>> Many thanks,
>>>> Natalie
>>>> ------------------------------
>>>> *From:* Mesa-users <mesa-users-bounces at lists.mesastar.org> on behalf
>>>> of Jared Goldberg via Mesa-users <mesa-users at lists.mesastar.org>
>>>> *Sent:* 08 August 2023 13:52
>>>> *To:* Amedeo Romagnolo <amedeoromagnolo at gmail.com>
>>>> *Cc:* mesa-users <mesa-users at lists.mesastar.org>
>>>> *Subject:* Re: [Mesa-users] [mesa-users] any recommendation for
>>>> running 200 Msun at Zsun?
>>>>
>>>> Hi Amedeo,
>>>>
>>>> The default work directory is not a good starting point for massive
>>>> stars (not to mention very likely not numerically resolved even at high
>>>> mesh and time delta resolution). Especially in metal-rich environments, the
>>>> convection zone around the Fe- and He- opacity peaks can become very
>>>> superadiabatic, leading to huge density inversions. I also see that you're
>>>> using MLT++ (okay_to_reduce_gradT_excess = .true.). There are a bunch of
>>>> tricks to deal with the superadiabatic surface (and other even dirtier
>>>> numerical problems). Check out the controls surrounding
>>>> use_superad_reduction, described in Section 7.2 of MESA VI  (
>>>> https://ui.adsabs.harvard.edu/abs/2023ApJS..265...15J/abstract
>>>> <https://urldefense.com/v3/__https://ui.adsabs.harvard.edu/abs/2023ApJS..265...15J/abstract__;!!IKRxdwAv5BmarQ!fDWOwg2C65ZkJYk7AM7SX_U9Eu27Xe5gwne99EcM6J1vcClntvQnJPKXfh6hi1ePscW-HeavCAH8oKIuVpkzO1M_$>),
>>>> which are more stable using the implicit method. The MESA 6 inlists for the
>>>> Convection section also have a bunch of other tricks which help get through
>>>> various stages of very massive star evolution.
>>>>
>>>> Hope this helps!
>>>> ~Jared
>>>>
>>>>
>>>>
>>>> On Tue, Aug 8, 2023 at 5:52 AM Amedeo Romagnolo via Mesa-users <
>>>> mesa-users at lists.mesastar.org> wrote:
>>>>
>>>> Hi all,
>>>>
>>>> I am trying a thought experiment with very massive stars in highly
>>>> metallic environments with MESA, but I am seriously struggling to make any
>>>> working setup for a 200 Msun star that either doesn't take ages or jumps
>>>> into the HR diagram like crazy.
>>>>
>>>> Here attached my inlist. As you can see I have been already messing
>>>> with the delta parameters quite a bunch, but I don't quite think this is
>>>> the right direction.
>>>>
>>>> Does anyone have any recommendations here?
>>>>
>>>> Cheers,
>>>>
>>>> Amedeo
>>>>
>>>> _______________________________________________
>>>> mesa-users at lists.mesastar.org
>>>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>>>> <https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!IKRxdwAv5BmarQ!fDWOwg2C65ZkJYk7AM7SX_U9Eu27Xe5gwne99EcM6J1vcClntvQnJPKXfh6hi1ePscW-HeavCAH8oKIuVj_8lDpA$>
>>>>
>>>> _______________________________________________
>>>> mesa-users at lists.mesastar.org
>>>>
>>>> https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!IKRxdwAv5BmarQ!fDWOwg2C65ZkJYk7AM7SX_U9Eu27Xe5gwne99EcM6J1vcClntvQnJPKXfh6hi1ePscW-HeavCAH8oKIuVj_8lDpA$
>>>>
>>>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20230825/8c5a1f5e/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_project
Type: application/octet-stream
Size: 6903 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20230825/8c5a1f5e/attachment.obj>