[Mesa-users] inlist for ~2Ms star terminating on the TP-AGB, mesa v10398
Warrick Ball
W.H.Ball at bham.ac.uk
Wed Sep 5 16:17:35 EDT 2018
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
>
>
More information about the Mesa-users
mailing list