[mesa-users] Convective core disappears when use_Ledoux=.true.

Fri Jun 13 18:05:32 EDT 2014

Hi Rich and Bill,
>> I’m playing around with evolving a 5 Msun model across the ZAMS. To my surprise, I find that core convection disappears (leading to a positive Brunt N2) when I set the use_Ledoux flag to .true.
> 
> Just like what happens for the 3 Msun case in Fig 15 of the 2nd mesa paper, assuming you are running with no semiconvection.
> 
> 
>> Can someone explain what’s going on here, please!
> 
> I look forward to hearing the explanation myself!

I suspect the problem is always the same: the Ledoux criterion applied from the 'radiative side' leads (incorrectly) to smaller, eventually vanishing, convective cores. 
I've attached a recent email thread where this issue was discussed in some detail.

hope this helps,

cheers,

Andrea   

>> Da: Arlette Grotsch <Arlette.Noels at ulg.ac.be>
>> Oggetto: Re: Consistently extending the core boundary
>> Data: 09 giugno 2014 14:32:46 GMT+1
>> A: Andrea Miglio <miglioa at bison.ph.bham.ac.uk>
>> Cc: Ehsan Moravveji <Ehsan.Moravveji at ster.kuleuven.be>, David Arnett <wdarnett at gmail.com>, Conny Aerts <Conny.Aerts at ster.kuleuven.be>, Bill Paxton <paxton at kitp.ucsb.edu>, Warrick Ball <wball at astro.physik.uni-goettingen.de>, Diego Bossini <dbossini at bison.ph.bham.ac.uk>, Thoul Anne <anne.thoul at ulg.ac.be>, Pieter Degroote <Pieter.Degroote at ster.kuleuven.be>
>> 
>> Dear all,
>> 
>> Just an example of how the core behaves in a low mass star with the Ldx criterion applied on the convective side and on the radiative side, and with the Schwarzschild criterion. You can see the near identical behavior of the core evolution when the criterions are applied on the convective sides. The big advantage is of course a more coherent Brunt-Väisälä frequency. 
>> 
>> This means that using Ledoux or Schwarzschild criterion is the same for a convective cores, which is obvious since the mixing is assumed to be perfect so the criteria are identical. The problem of the Ledoux criterion arises when a radiative layer (not adjacent to a convective zone) becomes convective. The criterion should be Ledoux in that case. If grad rad is intermediate between grad ad and grad Ldx, the layer is not convective since convection is the result of the non-linear development of the linear instability of dynamically unstable gravity modes. Only some high degree g+ modes are likely found unstable, leading to semi-convection.
>> 
>> I recommend the very interesting recent paper by Georgy, C. et al. 2014, MNRAS 439, L6, which gives arguments in favor of Ledoux.
>> 
>> I am fully aware that the « real » situation is closer to 3D simulations than to our poor LMLT, but, half a loaf is better than no loaf😉 for the time being.
>> 
>> Cheers,
>> Arlette
>> 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mcSch.pdf
Type: application/pdf
Size: 16194 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20140613/03ef6a6c/attachment.pdf>
-------------- next part --------------
>> 
>> 
>> 
>> 
>> Da: Arlette Grotsch <Arlette.Noels at ulg.ac.be>
> Oggetto: Re: Consistently extending the core boundary
> Data: 09 giugno 2014 13:56:33 GMT+1
> A: Andrea Miglio <miglioa at bison.ph.bham.ac.uk>
> Cc: Ehsan Moravveji <Ehsan.Moravveji at ster.kuleuven.be>, David Arnett <wdarnett at gmail.com>, Conny Aerts <Conny.Aerts at ster.kuleuven.be>, Bill Paxton <paxton at kitp.ucsb.edu>, Warrick Ball <wball at astro.physik.uni-goettingen.de>, Diego Bossini <dbossini at bison.ph.bham.ac.uk>, Thoul Anne <anne.thoul at ulg.ac.be>, Pieter Degroote <Pieter.Degroote at ster.kuleuven.be>
> 
> Dear all,
> 
> I of course fully agree with Andrea. Our paper is not at all presenting a new way of treating convective zones and convective boundaries. With the simple LMLT implemented in most codes (including our own CLES code), the boundary is searched for through an interpolation process of the function grad rad - grad ad (or grad Ledoux) in order to locate its zero, even when the function presents a discontinuity, which is the case with Schwarzschild when μ is discontinuous and is the case with Ledoux when gradient μ is discontinuous. 
> 
> We show that this leads to smaller convective cores in massive stars and vanishing convective cores in low mass stars when Ledoux is applied. Some problems also arise with Schwarzschild's criterion especially during core helium burning but H burning can be affected as well.
> 
> If the estimated convective core mass is too small at a given iteration, the interpolation through a discontinuity will inevitably block the boundary and an expansion of the core can only be found when the initial estimation is too large. 
> 
> The remedy is an extrapolation from the layers inside the core which leads to grad rad = grad ad or grad Ldx on the convective side of the boundary. We also advocate the use of a double mesh point at each boundaries and chemical discontinuities, which is a great help for stability analyses.
> 
> I shall be very happy to discuss these points in Geneva. I shall briefly present some of the related problems in my talk, especially those coming from an intermediate convection zone, for which, apart for showing the problems, I do not have any satisfying solutions😥.
> 
> All the best,
> Arlette 
> 

>> Le 9 juin 2014 à 13:48, Andrea Miglio <miglioa at bison.ph.bham.ac.uk> a écrit :
>> 
>>> Dear all,
>>> sorry to get into the discussion quite late.
>>> 
>>> The Gabriel et al paper is not  advocating for a new method to determine convective boundaries, but just re-stating a point which was known decades ago but then got somehow lost.  The risk we are facing is that, thanks to these wonderful asteroseismic constraints, we may end up claiming we are testing convection while what we are doing is merely fitting observations with models in which a simple  - yet tricky to implement - theory is not implemented correctly.
>>> 
>>> This is e.g. the case when computing models with the current implementation of the Ledoux criterion in MESA, which leads to very small convective cores and which  also impacts on the efficiency of any additional process that depends on the definition of a convective boundary (e.g. semiconvection). Another example (which we are focussing on in Birmingham) are the properties of convective cores during the core-He burning phase.
>>> Diego Bossini here is working on the tricky (for any code) definition of the convective core in core-Helium burning stars and, although he made some progress,  he is now focussing on getting a paper out: this is why he has  been very quiet lately (apologies in particular to Warrick and Eshan for the little feedback from Birmingham)
>>> 
>>> This message (and the paper) is by no means a criticism to MESA itself: the vast majority of the codes neglects this aspect. I am sure 3D hydro simulations will show us the way, but a precise implementation of the Schwarzschild criterion I would say is a pre-requisite for trying to introduce a more accurate description inspired by 3D hydro simulations.  
>>> 
>>> ok, after all this blathering, and having unfortunately too little time to get my hands dirty with the code in the next few weeks, what I suggest is that  we set aside some time in Geneva and discuss the details. Arlette will also be in there and I am sure she would be happy to get involved in the discussion (I’ve added her in cc).
>>> 
>>> thanks,
>>> 
>>> regards,
>>> 
>>> Andrea 
> 

Il giorno 13/giu/2014, alle ore 22:31, Bill Paxton <paxton at kitp.ucsb.edu> ha scritto:

> 
> On Jun 13, 2014, at 2:14 PM, Richard Townsend wrote:
> 
>> Hi folks —
>> 
>> I’m playing around with evolving a 5 Msun model across the ZAMS. To my surprise, I find that core convection disappears (leading to a positive Brunt N2) when I set the use_Ledoux flag to .true.
> 
> Just like what happens for the 3 Msun case in Fig 15 of the 2nd mesa paper, assuming you are running with no semiconvection.
> 
> 
>> Can someone explain what’s going on here, please!
> 
> I look forward to hearing the explanation myself!
> 
> 
>> 
>> Also, things look more or less OK when I set use_Ledoux = .false. — but the Brunt is a bit noisy in the mu-gradient zone. What’s considered the best approach to smoothing this noise?
> 
> 
> 
> smoothing is the best approach to smoothing.  try increasing the value for
> 
>      num_cells_for_smooth_brunt_B = 2
>         ! number of cells on either side to use in weighted smoothing of brunt_B
> 
> 
> -b
> ------------------------------------------------------------------------------
> HPCC Systems Open Source Big Data Platform from LexisNexis Risk Solutions
> Find What Matters Most in Your Big Data with HPCC Systems
> Open Source. Fast. Scalable. Simple. Ideal for Dirty Data.
> Leverages Graph Analysis for Fast Processing & Easy Data Exploration
> http://p.sf.net/sfu/hpccsystems
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users