[Mesa-users] inlist for ~2Ms star terminating on the TP-AGB, mesa v10398
Warrick Ball
W.H.Ball at bham.ac.uk
Thu Sep 6 16:30:08 EDT 2018
Hi Amber,
I'm using the January 27, 2018 version (20180127). I don't save PGStar
output but if you want anything from the default history, I've attached
the history file from my run on your inlist. For what it might be worth,
I've attached the terminal output, too.
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 Thu, 6 Sep 2018, amberlauer at gmail.com wrote:
> Hey Warrick, which SDK are you using?
>
> Do you have any pgstar output from this point that you can share?
>
> On Thu, Sep 6, 2018 at 9:39 AM Warrick Ball <W.H.Ball at bham.ac.uk> wrote:
> Hi Amber,
>
> Just to be sure, I copied your inlist into a fresh work folder from an
> unmodified copy of r10398 and it starts going through thermal pulses at
> the same point as my original run did. The terminal output so far (model
> 34845) is exactly the same as I got before.
>
> So I'm really not sure what we've done differently...
>
> Cheers,
> Warrick
>
> PS: Sorry about misreading the inlist you attached to your first email
> yesterday. My confusion was 100% self-inflicted, owing to a mixup on
> where attachments were saved...
>
> ------------
> 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, mesa-users at lists.mesastar.org wrote:
>
> > Oh, I misunderstood. Your email had a different inlist attached (when I
> > looked on the forum website, anyway) [1].
> >
> > I'll check in the morning if my run was actually doing something fishy but
> > that looks like it should work...
> >
> > Cheers,
> > Warrick
> >
> > [1]
> > https://lists.mesastar.org/pipermail/mesa-users/attachments/20180905/cbaceb76/attachment-0001.obj
> >
> > ------------
> > 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:
> >
> >> 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
> >
> >
> >
>
>
>
> --
> Amber Lauer
> Postdoctoral Researcher
> Triangle Universities Nuclear Lab
> Duke University
> amber.lauer at tunl.duke.edu
>
>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: log.txt.gz
Type: application/gzip
Size: 3190337 bytes
Desc:
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20180906/e1b4ad0c/attachment.gz>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: history.data.gz
Type: application/gzip
Size: 2390815 bytes
Desc:
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20180906/e1b4ad0c/attachment-0001.gz>
More information about the Mesa-users
mailing list