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

[Warrick Ball]

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.\A0 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.\A0 So I spun up a default run in r10398 with the following
>       almost trivial `inlist_project`:
> 
>
>       &star_job
>       \A0 \A0 \A0create_pre_main_sequence_model = .true.
>       / !end of star_job namelist
> 
>
>       &controls
>       \A0 \A0 \A0initial_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
> 
>