Re: [mesa-users] Eddington Limit

[Pablo M. <pa...@gm...>]

Sorry for the late reply.

The default mdot_edd is meant for black holes, so you should adjust it for
anything else! This probably merits a few comments in the defaults file so
people don't get confused. If you want to change mdot_edd you don't need to
modify directly the code in binary/private, I created a hook for this so
you can just define it in the src/run_binary_extras.f of your work
directory. For instructions, read the comments in

$MESA_DIR/binary/other/mod_other_mdot_edd.f90

On Mon, Dec 7, 2015 at 11:38 PM, Francis Timmes <fx...@ma...> wrote:

> consider applying for the 2016 summer school when it opens;
> we’ll make the announcement on this forum.
>
> fxt
>
>
>
> > On Dec 7, 2015, at 2:34 PM, Kenny Van <kv...@ua...> wrote:
> >
> > Hey,
> >
> > Wow thanks, those summer school pdfs are very helpful. They go over most
> of what I'm trying to look at.
> >
> > Thanks
> >
> >
> >
> > On Mon, Dec 7, 2015 at 1:16 PM, Josiah Schwab <jws...@be...>
> wrote:
> > Hi Kenny,
> >
> > > I was going through the code to try and set a different maximum
> mdot_edd
> > > and I noticed inside of binary_mdot.f90 that mdot_edd seems to be
> defined
> > > oddly
> > >
> > > if (.not. b% use_es_opacity_for_mdot_edd) then
> > >             mdot_edd = 4*pi*b% s1% cgrav(1)*b% m(b% a_i)/(clight*b% s1%
> > > opacity(1))
> > >          else
> > >             mdot_edd = 4*pi*b% s1% cgrav(1)*b% m(b%
> a_i)/(clight*0.2*(1+b%
> > > s1% surface_h1))
> > >
> > > Which seem to correspond to:
> > >
> > > Mdot = 4 pi G M / (c K)
> > >
> > > and
> > >
> > > Mdot = 4 pi G M / (c (1 + X))
> > >
> > > Is this equivalent to the normal representation of
> > >
> > > Mdot = 4 pi c R / K?
> > >
> > > While the units appear to match, I would assume that the two are not
> > > necessarily equivalent. Unless I'm misunderstanding the code here. Does
> > > anyone have some insight to this?
> >
> > Well, it is equivalent if you set R = GM/c^2. :)
> >
> > There's not a universally useful definition of Mdot_edd, given that MESA
> > can evolve both point mass and non-point mass companions.
> >
> > Users should make use of the other_mdot_edd hook and define the Mdot_edd
> > that is appropriate for your problem.
> >
> > You may be interested to look at Sterl Phinney's labs from the most
> > recent MESA summer school.
> >
> >
> http://mesastar.org/teaching-materials/2015-mesa-summer-school/phinney-lecture-2-advanced-mass-transfer-and-binary-extras
> >
> > Josiah
> >
> >
> >
> ------------------------------------------------------------------------------
> > Go from Idea to Many App Stores Faster with Intel(R) XDK
> > Give your users amazing mobile app experiences with Intel(R) XDK.
> > Use one codebase in this all-in-one HTML5 development environment.
> > Design, debug & build mobile apps & 2D/3D high-impact games for multiple
> OSs.
> >
> http://pubads.g.doubleclick.net/gampad/clk?id=254741911&iu=/4140_______________________________________________
> > mesa-users mailing list
> > mes...@li...
> > https://lists.sourceforge.net/lists/listinfo/mesa-users
>
>
>
> ------------------------------------------------------------------------------
> Go from Idea to Many App Stores Faster with Intel(R) XDK
> Give your users amazing mobile app experiences with Intel(R) XDK.
> Use one codebase in this all-in-one HTML5 development environment.
> Design, debug & build mobile apps & 2D/3D high-impact games for multiple
> OSs.
> http://pubads.g.doubleclick.net/gampad/clk?id=254741911&iu=/4140
> _______________________________________________
> mesa-users mailing list
> mes...@li...
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>



-- 
Pablo Marchant Campos
M.Sc on Astrophysics, Universidad Católica de Chile
PhD student, Argelander-Institut für Astronomie