[mesa-users] Discontinuity

[Chris Mankovich]

Hi Kenny,

What happens if you add an else clause to make sure that the turnover time
is zero in all radiative zones? With the loop written as

         turnover_time = 0.0
         do k = 1, nz
            if (s% mixing_type(k) == convective_mixing) then
               if (k == nz) then
                  dr = s% r(k) - s% R_center
               else
                  dr = s% r(k) - s% r(k+1)
               end if
               ! s% conv_vel in cm/s
               ! dr in cm
               turnover_time = turnover_time + (dr/ s% conv_vel(k))
            end if
         end do

you risk "contaminating" later radiative zones with the turnover time
calculated for the last convective zone. Just a thought.


chris

On Sat, Jul 9, 2016 at 1:29 PM, Francis Timmes <fxt44 at mac.com> wrote:

> 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
>
>
>
> ------------------------------------------------------------------------------
> 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160711/71e1289a/attachment.html>