[mesa-users] off-center Ne and O burning (again…)

Stan Woosley woosley at ucolick.org
Sun Apr 21 15:43:57 EDT 2013


Just a thought. When semi convection is small or turned off in a code
the response is often to create many small alternating convection zones,
often only a few zones thick. These come and go out of phase and 
act as a staircase that can mimic the effects of semiconvective mixing. I think
Icko Iben first saw this many years ago.

Anyway, it may not be  a good scheme to ignore convection anywhere
that the gradients say it should occur, even if its just a few zones.

   Stan



On Apr 21, 2013, at 12:52 AM, Jakub Ostrowski <ostrowski at astro.uni.wroc.pl> wrote:

> Hi!
> 
> I'm working almost exclusively on 15 M_s models and I have found quite a big difference between the v4906 and previous ones too - but somewhere else than you.
> 
> I studied properties of models after depletion of central H, logTeff ~ 4.30, logL ~ 4.63. The behavior of Brunt-Vaisala frequency is much smoother in v4906 than before and finally now there is convection where it should be according to behavior of gradients.
> 
> And this is due to one change in controls.default. Now there is
> 
> prune_bad_cz_min_Hp_height = 0.1 ! lower limit on radial extent of cz
> 
> whereas previously the value was 1.0. I don't know if this is useful for you but for me this previously overlooked parameter made a huge change in case of my model. But it seems there is no dramatic influence on the shape of  track on HRD.
> 
> Take care,
> Jakub
> 
> On 21 kwi 2013, at 06:59, Tuguldur Sukhbold <sukhbold at ucolick.org> wrote:
> 
>> Dear MESA users,
>> 
>> I am getting 2 different results with the exact same inlist for two different versions of the code: 4798 and the current 4906. Since MESA is constantly being developed, I understand that we should expect changes between the versions. However I am little surprised to see such a 'big' change in two releases came out almost one after the other.
>> 
>> The attached plot shows the central T-Rho plane evolution for 15Msun model with the following simple inlist:
>> 
>> &star_job
>> 
>>      mesa_dir = '' 
>>      create_pre_main_sequence_model = .true.
>> 
>>      set_rate_c12ag = 'Buchmann'
>>      num_special_rate_factors = 1
>>      reaction_for_special_factor(:) = 'r_c12_ag_o16'
>>      special_rate_factor(1) = 1.2
>> 
>>      eos_file_prefix = 'mesa'
>>      kappa_file_prefix = 'gs98'
>> 
>> / !end of star_job
>> 
>> &controls
>> 
>>      initial_mass = 15
>>      initial_Z = 0.02
>>      initial_Y = 0.25
>> 
>>      mixing_length_alpha = 2
>> 
>>      RGB_wind_scheme = 'Nieuwenhuijzen'
>>      AGB_wind_scheme = '' 
>>      RGB_to_AGB_wind_switch = 1d-4
>>      Nieuwenhuijzen_wind_eta = 1.0d0
>> 
>>      photostep = 50
>>      profile_interval = 25
>>      history_interval = 1
>>      terminal_interval = 10
>>      write_header_frequency = 10
>> 
>> !hign
>> !      xa_central_lower_limit_species(1) = 'h1'
>> !      xa_central_lower_limit(1) = 0.72
>> !hdep
>> !      xa_central_lower_limit_species(1) = 'h1'
>> !      xa_central_lower_limit(1) = 0.01
>> !heign
>> !      xa_central_upper_limit_species(1) = 'c12'
>> !      xa_central_upper_limit(1) = 0.01
>> !hedep
>> !      xa_central_lower_limit_species(1) = 'he4'
>> !      xa_central_lower_limit(1) = 0.01
>> !cign
>> !      log_center_temp_limit = 8.69897 !0.5d8K
>> !cdep
>>      log_center_temp_limit = 9.07918 !1.2d9K
>> !odep
>> !      xa_central_lower_limit_species(1) = 'o16'
>> !      xa_central_lower_limit(1) = 0.05
>> 
>> /
>> 
>> If you guys remember, I have previously brought up the issue of off-center O burning, in the case of alpha_mlt = 1, in version 4798. This behavior was gone if I set this parameter to alpha_mlt = 2, as it creates a larger CO core. However, when I try to recreate the older results with the current version of MESA - 4906 - I am facing the off-center burning issue again with alpha_mlt = 2. 
>> 
>> In the response to my previous email, Bill has reported that he got the same result (no off-center burning) with alpha_mlt=2 in his version of MESA - 4917.
>> 
>> MESA version	alpha_mlt=1	alpha_mlt=2
>> 4798			y			n
>> 4906			y			y
>> 4917			y			n
>> 
>> (yes/no to off-center Ne and O burning)
>> 
>> I haven't seen anything in the release email from Bill for 4906, that might be relevant to what I see here. 
>> 
>> In general, 15Msun is a funny mass in MESA. It is sitting on the border of two different paths of evolutions of CO core: (1) degeneracy is important and off-center burnings of Ne and O as we've seen (2) degeneracy effects are no longer important and everything ignites in the center first. In other words, 15Msun model essentially with all the default settings, is producing a CO core mass that borderlines these two paths. This is why changing parameters here and there are producing these rather big changes for this 15Msun case. I don't think one will see this behavior in bigger mass models in any version of the code with any alpha_mlt, simply because the CO core will be big enough.
>> 
>> From my experience, this funny mass is usually much lower in other codes, more like in the 8-12Msun range.
>> 
>> Tuguldur
>> 
>> 
>> 
>> <inc.png>
>> ------------------------------------------------------------------------------
>> Precog is a next-generation analytics platform capable of advanced
>> analytics on semi-structured data. The platform includes APIs for building
>> apps and a phenomenal toolset for data science. Developers can use
>> our toolset for easy data analysis & visualization. Get a free account!
>> http://www2.precog.com/precogplatform/slashdotnewsletter_______________________________________________
>> mesa-users mailing list
>> mesa-users at lists.sourceforge.net
>> https://lists.sourceforge.net/lists/listinfo/mesa-users
> 
> 
> ------------------------------------------------------------------------------
> Precog is a next-generation analytics platform capable of advanced
> analytics on semi-structured data. The platform includes APIs for building
> apps and a phenomenal toolset for data science. Developers can use
> our toolset for easy data analysis & visualization. Get a free account!
> http://www2.precog.com/precogplatform/slashdotnewsletter
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
> 
> !DSPAM:212,5173a27a7832272969976!
> 





More information about the Mesa-users mailing list