[mesa-users] BV frequency, gradients and mlt problems - bugs?

[Jakub Ostrowski]

Hi Bill,

Thank you for a very quick reply.

The models were created with version 4723 (but I've checked today 4740 calculating a 15 M_S models with Henyey MLT - the problem with gradients is there too, I haven't looked at N^2). I use OS X 10.8.2 and the latest version of SDK.

Presented models have numbers 1080, 1940 and 2080. From the profile files I plotted brunt_N2_dimensionless vs. logT and grada, gradT and gradr vs. logT. The inlist should be enough to recreate the problem because I haven't modified default controls or star_job files.

If you need more information or files let me know. 

Cheers,
Jakub


On 2 sty 2013, at 19:11, Bill Paxton <paxton at kitp.ucsb.edu> wrote:

> Hi Jakub,
> 
> I'll take a look.  Thanks for sending the inlist.  In addition it is helpful if problem reports also include the mesa version number, the operating system you are running, and the compiler.   Also please give me a detailed guide to creating the exact models that you show in the plots so I can reproduce the case.
> 
> Thanks,
> Bill
> 
> 
> 
> 
> 
> 
> 
> On Jan 2, 2013, at 9:58 AM, Jakub Ostrowski wrote:
> 
>> Hi!
>> 
>> I'm doing some tests with 15 M_S blue supergiant models and I found a few things that seriously bother me. I don't know if these are bugs or just problems of inexperienced user.
>> 
>> First of all, I started with plots of dimensionless Brunta-Vaisala frequency in function of logT for three models on different stages of evolution. As you can see in the attached figures, first model doesn't undergo central helium burning while the two following models do so and hence they have convective cores. There have been chosen to have similar T_eff. My first problem are the weird spikes that you can see. Some of them look like they have a bad sign, that means they should go to zero instead of infinity, but there are also peeks which don't seem to be related to presence of any convective zone. Has anyone else experienced similar problem? For me it looks like numerical problem, but of course I might be wrong.
>> 
>> Later on, I was trying to find a cause of this weird behavior of N and then I've found another problem. I plotted gradT, grada and gradr and there seems to be a problem with them. I guess that in core or in any shell that undergoes nuclear burning gradT is arbitrary set to equal grada (so that we have gradT = grada) and in the outer convective zones there is MLT in use. In these zones there are recipes how to set gradT using some combination of both grada and gradr. And here comes the problem: in the outer convective zone related to the Z-bump it seems to work properly (as you can see in the figure) but in the outermost zone related to He opacity bump it looks like the gradT = grada just like in the core. I've checked the default settings of mixing and I can't see any reason of these behavior.
>> 
>> I attached some figures and used inlist file. I would be grateful for any help because I have no idea how to fix these issues.
>> 
>> Cheers,
>> Jakub
>> 
>> 
>> 
>> 
>> <inlist_15m_simple>
>> 
>> 
>> <15m-n2.pdf><15m-grad.pdf>------------------------------------------------------------------------------
>> 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
>