[mesa-users] helium flash mlt bug

[Jakub Ostrowski]

Hi!

Kevin, thank you for the very interesting paper and clarifying the post-flash behavior.

Josiah, I checked that I indeed had some duplicated model numbers in history file but they seem irrelevant and clearly not related to the problem. Anyway, thank you for pointing me out that such thing could make some troubles. I had no idea about it.

Cheers,
Jakub


On 9 sty 2013, at 19:39, Kevin Moore <kmoore at physics.ucsb.edu> wrote:

> Hi Jakub,
> 
> I'm not sure what's going on with the mixing_length_alpha = 1.8 case, but I can tell you that the post-flash behavior looks like a typical helium flash. As the heat released by he burning moves into the core, it ignites successive subflashes in the degenerate material. There's more detail in sec. 2 of our paper here: http://adsabs.harvard.edu/abs/2012ApJ...744L...6B
> 
> -Kevin
> 
> On Jan 9, 2013, at 9:45 AM, Jakub Ostrowski wrote:
> 
>> Hi!
>> 
>> I'm writing to report a bug I found when I was trying to calculate a 1 M_S model evolution.
>> 
>> When mixing_length_alpha = 1.8 there is a problem with helium flash. As you can see in the attached HR diagram the behavior in the vicinity of the tip of giant branch is totally wrong. After this weird behavior is done, everything seems to look pretty normal.
>> 
>> The problem occurs independently of other parameters. I attached the inlist I used, but when I changed other parameters, like enabling overshooting or choosing different reaction rates, the problem was also present.
>> 
>> The behavior was correct when I set mixing_length_alpha = 2.0, as you can see in the second attached figure. It was also correct with 1.8 but in 2 M_S model.
>> 
>> I would also like to ask about the behavior after the helium flash. There is chaotic behavior in the HR diagram (you can see it in both figures) before the proper onset of helium burning. Do you know if there is some physical reason behind it or is it numerical problem?
>> 
>> I use MESA 4740 with latest sdk on OS X 10.8.2. The problem is also present on Ubuntu 12.10 (v4740, latest sdk). If you need more informations or files let me know.
>> 
>> Cheers,
>> Jakub
>> 
>> 
>> <inlist_1m_mlt18><1m-hr-mlt18.ps><1m-hr-mlt20.ps>------------------------------------------------------------------------------
>> Master Java SE, Java EE, Eclipse, Spring, Hibernate, JavaScript, jQuery
>> and much more. Keep your Java skills current with LearnJavaNow -
>> 200+ hours of step-by-step video tutorials by Java experts.
>> SALE $49.99 this month only -- learn more at:
>> http://p.sf.net/sfu/learnmore_122612 _______________________________________________
>> mesa-users mailing list
>> mesa-users at lists.sourceforge.net
>> https://lists.sourceforge.net/lists/listinfo/mesa-users
> 
> 
> ------------------------------------------------------------------------------
> Master Java SE, Java EE, Eclipse, Spring, Hibernate, JavaScript, jQuery
> and much more. Keep your Java skills current with LearnJavaNow -
> 200+ hours of step-by-step video tutorials by Java experts.
> SALE $49.99 this month only -- learn more at:
> http://p.sf.net/sfu/learnmore_122612 
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>