[mesa-users] Discontinuity

[Francis Timmes]

hi kenny,

it might be worthwhile contacting konstantin and/or natalia
since they have already implemented this in mesa. per the 
open-source, open-knowledge manifesto of mesa I, i'm hopeful
they would be willing to share. 

fxt

 



> On Jul 9, 2016, at 12:29 AM, Kenny Van <kvan at ualberta.ca> wrote:
> 
> Hi,
> 
> So I've been working on incorporating turnover time and stellar wind loss into magnetic braking. The equation that I'm basing this work off of is from the following paper:
> 
> http://arxiv.org/pdf/1511.08847v1.pdf
> 
> in equation 16.
> 
> I've been able to add these additional scaling factors using the "run_binary_extras.f" module which I've attached. The issue that I'm running into is there are some discontinuities when calculating the turnover time in these systems. I've attached a plot of turnover time in relation to model number
> 
> <tt_plot.png> 
> 
> Turnover time in this case only depends on convective velocity and the radius of the convective zone. I do a quick check of the convective velocity and radius of the convective zone using the profile.data files and I can produce the following plots:
> 
> <conv_vel.png><mixing_types.png>
> 
> The convective velocity appears to remain largely the same, while the radius shifts around 0.65 Rsun.
> 
> Does anyone have an insight on why there seems to be this sharp increase in the on the turnover time? I believe it is tied to the change in the radius of convective zone but it might be something else.
> 
> I've attached all the relevant files below. 
> 
> Thanks for any advice anyone can give
> <inlist><inlist1><binary_run.f><run_binary_extras.f><run_star_extras.f>------------------------------------------------------------------------------
> Attend Shape: An AT&T Tech Expo July 15-16. Meet us at AT&T Park in San
> Francisco, CA to explore cutting-edge tech and listen to tech luminaries
> present their vision of the future. This family event has something for
> everyone, including kids. Get more information and register today.
> http://sdm.link/attshape_______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users