[Mesa-users] Trouble w. convective stability criteria

[Rob Farmer]

Please keep replies on the mesa-users mailing list.

What happens if you use mesa's mixing window to show what's going on?
Add pgstar_flag=.true. to star_job and this to a pgstar inlist

  mixing_win_flag=.true.
  Mixing_xaxis_name = 'r_div_R'

  mixing_file_flag=.true.
  mixing_file_interval=10

    Mixing_xmin = 0.92
    Mixing_xmax = 1.0

The attached mixing window shows that mesa is putting overshoot below the
convective boundary. But I did not run the model for long, so you should
check what happens at whatever point in time you think its not working.

You could also try adding log_D_ovr to your mixing plot to show the
overshoot diffusion coefficient

Rob


On Fri, 9 Apr 2021 at 11:56, Peter Juelsgaard <pjuels at gmail.com> wrote:

> H
> Thanks for the suggestion, but sadly, no effect
>  - Peter
>
> On Fri, Apr 9, 2021 at 10:45 AM Rob Farmer <r.j.farmer at uva.nl> wrote:
>
>> Hi,
>> Can you try setting your overshoot as:
>>
>>     overshoot_scheme(1) = 'exponential'
>>     overshoot_zone_type(1) = 'any'
>>     overshoot_zone_loc(1) = 'any'
>>     overshoot_bdy_loc(1) = 'any'
>>     overshoot_f(1) = 0.030
>>     overshoot_f0(1) = 0.008
>>
>> i.e using an array index of (1) instead of (:)
>>
>> Rob
>>
>>
>> On Fri, 9 Apr 2021 at 10:08, Peter Juelsgaard via Mesa-users <
>> mesa-users at lists.mesastar.org> wrote:
>>
>>> Hi all.
>>> I am having trouble with an area of a stellar model, where the
>>> Schwarzvhild criteria fails (ie. the area should be convective), but the
>>> Diffusion coefficient is effectively 0.
>>> The system I work on is:
>>> Debian 3.16.43-2 (as far as I can figure out from what the server login
>>> tells me)
>>> MESA ver. mesa-r12778
>>> mesasdk-20200325
>>>
>>> The trouble arises in 0.927 < r/R < 0.930
>>> [image: MixingTrouble2.png]
>>> It should be noted that the f_OV in the plot title is the effective
>>> overshoot, which I calculate from f_effective = overshoot_f(:) -
>>> overshoot_f0(:). Also, the area where log_D_mix = -3 is an artificial
>>> cutoff I set, log_D_mix < -3 gets set to -3. There is an attached figure,
>>> where I remove this cutoff, and you can see that log_D_mix=-100 in this
>>> area. I have also included the inlist and the loaded model, for
>>> reproducibility.
>>>
>>>
>>> Kind regards
>>>  - Peter Juelsgaard
>>>
>>>
>>>
>>>
>>> I have attached the inlist files, a png. detailing
>>> _______________________________________________
>>> mesa-users at lists.mesastar.org
>>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>>>
>>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210409/19dc0aae/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MixingTrouble2.png
Type: image/png
Size: 36534 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210409/19dc0aae/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mixing_000890.png
Type: image/png
Size: 3948 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210409/19dc0aae/attachment-0001.png>