[Mesa-users] inlist for ~2Ms star terminating on the TP-AGB, mesa v10398

Francis Timmes fxt44 at mac.com
Wed Sep 5 16:14:39 EDT 2018


grab my work directory or just the inlist from the 2018 mesa summer school.
it takes a 2.1 msun from zams to through all the thermal pulses to a wd.
i've also used it on many other masses.

http://cococubed.asu.edu/mesa_market/education.html

it is not minimal as i aim to publish something with it (paper I submitted).

fxt





> On Sep 5, 2018, at 1:05 PM, amber lauer via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> 
> I'm using an inlist created exactly to mimic your minimal inlist (see below), with some file saving settings. Which SDK are you using? I and some of my colleagues are experiencing the same behavior with other models that are known to go through extensive pulses, where it hangs on the second pulse. I'm just trying to check all the options. Thanks.
> 
> &star_job
> 
>   ! begin with a pre-main sequence model
>     create_pre_main_sequence_model = .true.
> 
>   ! save a model at the end of the run
>     save_model_when_terminate = .true.
>     save_model_filename = 'TP_AGB.mod'
> 
>   ! display on-screen plots
>     pgstar_flag = .true.
> 
> / !end of star_job namelist
> 
> 
> &controls
> 
>   ! starting specifications
>     initial_mass = 1.5 ! in Msun units
> 
>     max_model_number=50000    
> 
> 
> / ! end of controls namelist
> 
> 
> On Wed, Sep 5, 2018 at 3:56 PM Warrick Ball <W.H.Ball at bham.ac.uk> wrote:
> Hi Amber,
> 
> In my run, the thermal pulses get started around model number 21000.  I 
> killed the run at model number 107063, which is the point by which I 
> counted 34ish thermal pulses.
> 
> This isn't really my area: I seldom evolve models much past the RGB bump, 
> let alone beyond helium ignition.  I don't know which of the extra 
> controls might be holding you up!  That said, I'm not used to the rotation 
> controls but it looks to me like you're creating a 0.98 Msun star 
> and setting the rotation rate to 204.8 km/s.  Is that not a bit quick?  I 
> might be misinterpreting things.
> 
> Cheers,
> Warrick
> 
> 
> ------------
> Warrick Ball
> Postdoc, School of Physics and Astronomy
> University of Birmingham, Edgbaston, Birmingham B15 2TT
> W.H.Ball at bham.ac.uk
> +44 (0)121 414 4552
> 
> 
> On Wed, 5 Sep 2018, amberlauer at gmail.com wrote:
> 
> > Hey Warrick, just a quick question about your minimal inlist_project.
> > 
> > Can you give me a rough idea how far you ran it before you reached 34 thermal pulses?
> > 
> > I ran the attached inlist, modeled after yours, and at 2.98e9 yrs, 50K models, and a log L/L=~4, log Teff=!~3.6 it is still hung on the second
> > pulse. FYI it stays roughly at that position for about 30
> > k models, see attached images. I'm running v 10398, sdk from 1/18.  Thanks
> > 
> > 
> > 
> > grid1_020000.png grid1_030000.png
> > grid1_020000.png
> > 
> > 
> > grid1_050000.png
> > 
> > 
> > 
> > On Tue, Apr 24, 2018 at 3:17 AM Warrick Ball via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> >       Hi,
> >
> >       I realise this is basically resolved but the question piqued my curiosity.
> >       In my relatively-limited experience of modelling post-RGB stars, I thought
> >       low-mass MESA models usually started going through the TP-AGB without any
> >       special input.  So I spun up a default run in r10398 with the following
> >       almost trivial `inlist_project`:
> > 
> >
> >       &star_job
> >            create_pre_main_sequence_model = .true.
> >       / !end of star_job namelist
> > 
> >
> >       &controls
> >            initial_mass = 1.5 ! in Msun units
> >       / ! end of controls namelist
> > 
> >
> >       Sure, it's slow, but I left it running overnight and after 18 hours on 2
> >       cores I count 34 thermal pulses.
> >
> >       Cheers,
> >       Warrick
> > 
> >
> >       ------------
> >       Warrick Ball
> >       Postdoc, School of Physics and Astronomy
> >       University of Birmingham, Edgbaston, Birmingham B15 2TT
> >       wball at bison.ph.bham.ac.uk
> >       +44 (0)121 414 4552
> >
> >       On Mon, 23 Apr 2018, mesa-users at lists.mesastar.org wrote:
> >
> >       > Hi all,
> >       > Would anyone be able to point me to an inlist which generates a TP-AGB star (terminating during the thermal pulsations) that runs
> >       with MESA v 10398? I had success with the attached inlist ('inlist_tashibu_2017a') used by Tashibu et al (2017) in mesa version 8118,
> >       but it produces a number of parameter name and convergence issues in v10398. I have guessed at the appropriate changes when MESA
> >       points out a naming conflict, but have thus far failed to produce a converged model. One of my attempted adjusted inlists is attached
> >       as 'inlist_tashibu_editing_attempt'
> >       >
> >       > At the moment, I am not so interested in diagnosing why Tashibu's exact physical prescription no longer works, I would just like to
> >       be able to run a generic TP-AGB model *without* using one of MESA's pre-built prems models to start.
> >       >
> >       > Any insight is greatly appreciated,
> >       >
> >       >
> >       > Meridith Joyce
> >       > Ph.D. Candidate, Wilder Lab
> >       > Department of Physics and Astronomy
> >       > Dartmouth College
> >       >
> >       _______________________________________________
> >       mesa-users at lists.mesastar.org
> >       https://lists.mesastar.org/mailman/listinfo/mesa-users
> > 
> > 
> > 
> > --
> > Amber Lauer
> > Postdoctoral Researcher
> > Triangle Universities Nuclear Lab
> > Duke University
> > amber.lauer at tunl.duke.edu
> > 
> >
> 
> 
> -- 
> Amber Lauer
> Postdoctoral Researcher
> Triangle Universities Nuclear Lab
> Duke University
> amber.lauer at tunl.duke.edu
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
> 




More information about the Mesa-users mailing list