[mesa-users] Deuterium burning in PMS stars

Michal Simon michal.simon at stonybrook.edu
Fri Dec 11 15:09:45 EST 2015


Hi Aaron et al,

Your suggested fixes worked well.  I learned a lot over the last few
days of often frustrating MESAing from you, Josiah, and Robert Farmer.
Thanks to all!

What MESA needs is not a thick text but a "MESA for Dummies".

Cheers, Mike




On Thu, Dec 10, 2015 at 10:35 PM, Aaron Dotter <aaron.dotter at gmail.com>
wrote:

> Hi Mike,
>
> It's safest to be explicit with things like this, i.e., specify the
> initial abundances in your inlist rather than rely on whatever MESA does by
> default.
>
> Default behavior may change over time (though I'm not certain that it has
> in this case).
>
> Aaron
>
> On Fri, Dec 11, 2015 at 2:25 PM, Michal Simon <michal.simon at stonybrook.edu
> > wrote:
>
>> Hi Aaron,  Jayashree and I will try this fix tmw, Fri.  But I am very
>> puzzled.  Mesa 6596 pulled values of center h1, h2, li7 out of its table of
>> abundances, namely 0.7, 2e-5,  etc correctly.   Does Mesa 7624 work
>> differently?  Mike
>>
>> Sent from my iPhone
>>
>> On Dec 10, 2015, at 5:29 PM, Aaron Dotter <aaron.dotter at gmail.com> wrote:
>>
>> Hi Mike,
>>
>> What do you expect for center_h2?
>>
>> In your inlist you have not specified any initial composition information
>> so you get some default values, probably not what you want.  May I suggest
>> you try setting these values in the star_job section of your inlist:
>>
>>
>>
>>       set_uniform_initial_composition = .true.
>>
>>
>>          !### initial_h1
>>          !### initial_h2
>>          !### initial_he3
>>          !### initial_he4
>>          ! if `set_uniform_initial_composition` is true, then
>>          ! set hydrogen and helium mass fractions according to the
>> following:
>>          ! If no h2 in current net, then this will be added to h1.
>>          ! If no he3 in current net, then this will be added to he4.
>>
>> (your preferred values here)
>>
>>       initial_h1 = -1
>>       initial_h2 = -1
>>       initial_he3 = -1
>>       initial_he4 = -1
>>
>>
>>          !### initial_zfracs
>>          ! if `set_uniform_initial_composition` is true, then set metal
>> fractions
>>          ! z fractions -- select one of the options defined in
>> `chem/public/chem_def` :
>>
>>          ! + `AG89_zfracs = 1`
>>          ! + `GN93_zfracs = 2`
>>          ! + `GS98_zfracs = 3`
>>          ! + `L03_zfracs = 4`
>>          ! + `AGS05_zfracs = 5`
>>          ! + `AGSS09_zfracs = 6`
>>          ! + `L09_zfracs = 7`
>>          ! + `A09_Prz_zfracs = 8`
>>
>>
>> (from your inlist it seems you want GS98 for the solar composition so set
>> )
>>
>> initial_zfracs = 3
>>
>>
>>
>>
>> Aaron
>>
>> On Fri, Dec 11, 2015 at 7:00 AM, Michal Simon <
>> michal.simon at stonybrook.edu> wrote:
>>
>>> Dear Colleagues,
>>>
>>> I am running MESA 7624 and trying to determine the time when the  very
>>> brief phase of D burning operates in PMS stars.  I am having the history
>>> file list the value of center_h2 as the star contracts.  The values for the
>>> mass
>>> fraction of deuterium are nonsense.
>>>
>>> What am I doing wrong?  I attach inlist_project and
>>> history.data files.
>>>
>>> Cheers and thanks,
>>> Mike Simon
>>>
>>>
>>> ------------------------------------------------------------------------------
>>>
>>> _______________________________________________
>>> mesa-users mailing list
>>> mesa-users at lists.sourceforge.net
>>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20151211/1dda0d9a/attachment.html>


More information about the Mesa-users mailing list