[Mesa-users] inlist for ~2Ms star terminating on the TP-AGB, mesa v10398
amber lauer
amberlauer at gmail.com
Tue Sep 11 12:54:35 EDT 2018
HI all. I feel a bit sheepish, the source was my pgstar settings. In the
Log_L plot, I was using star_age for the xaxis, without really
understanding how to set it to be useful. Switching the x axis to
model_number revealed that the pulses were there all along. Thanks for the
help and I'm glad it turned out to be something silly.
On Thu, Sep 6, 2018 at 4:30 PM Warrick Ball <W.H.Ball at bham.ac.uk> wrote:
> 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
> >
> >
--
Amber Lauer
Postdoctoral Researcher
Triangle Universities Nuclear Lab
Duke University
amber.lauer at tunl.duke.edu <amber.lauer at duke.edu>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20180911/3b51c65a/attachment.html>
More information about the Mesa-users
mailing list