[Mesa-users] inlist for ~2Ms star terminating on the TP-AGB, mesa v10398
amber lauer
amberlauer at gmail.com
Thu Sep 6 13:54:05 EDT 2018
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 at duke.edu>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20180906/707a0428/attachment.html>
More information about the Mesa-users
mailing list