[Mesa-users] The nuclear burning after log T = 9.85 is stop ?
Francis Timmes
fxt44 at mac.com
Wed Sep 21 17:51:13 UTC 2022
please keep all messages on mesa-users.
reaclib rates are fit up to 10e9 k. bad things can (and do) happen when
these polynomial/rational fits are extrapolated beyond their fitting range.
mesa limits the temperature to 10e9 k to avoid such extrapolations and issues a warning.
fxt
> On Sep 18, 2022, at 8:24 PM, 辛文宇 <wyxin at nao.cas.cn> wrote:
>
> Hi, Pro. Timmes,
>
> Thanks very much for your help.
> I tried to uncomment them and re-run the inlist_finish.
> when log T_c = 9.97, the following warning appears:
>
> WARNING: evaluating rates with lgT=10.001623 which is above
> lgT=10.000000, rates have been truncated
>
> I use the reaction rates for jina reaclib. Are the rates in reaclib
> only include the the part for T < 10 GK ?
>
> Thanks,
> Wenyu
>
> > -----Original Messages-----
> > From: "Francis Timmes" <fxt44 at mac.com>
> > Sent Time: 2022-09-18 00:32:48 (Sunday)
> > To: "辛文宇" <wyxin at nao.cas.cn>, xinwenyu16 at mails.ucas.ac.cn
> > Cc: "Frank Timmes" <fxt44 at mac.com>, mesa-users <mesa-users at lists.mesastar.org>
> > Subject: Re: [Mesa-users] The nuclear burning after log T = 9.85 is stop ?
> >
> > > ... nuclear burning in the outer shell seems to stopped.
> >
> > inlist_finish turns off nuclear reactions for this dynamical phase of evolution:
> >
> > eps_nuc_factor = 0d0
> > dxdt_nuc_factor = 0d0
> >
> > https://docs.mesastar.org/en/release-r22.05.1/reference/controls.html#eps-nuc-factor
> > https://docs.mesastar.org/en/release-r22.05.1/reference/controls.html#dxdt-nuc-factor
> >
> > fxt
> >
> >
> >
> >
> > > On Sep 13, 2022, at 8:04 PM, 辛文宇 via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> > >
> > > Hi,
> > >
> > >
> > >
> > > I am using a example '25M_pre_ms_to_core_collapse' in test_suite of MESA 12778
> > >
> > > to evolve a star form ZAMS to Fe core collapse.
> > >
> > > Some of the inlist files are attached. The following three header files are used successively.
> > >
> > > inlist_to_lgT_9.9_header
> > >
> > > inlist_convert_header
> > >
> > > inlist_finish_header
> > >
> > > When log T_c reaches 9.85, the 'inlist_to_lgT_9.9_header' is finished and
> > >
> > > the 'inlist_convert_header' is continue. This file seems evolve the models forwards 100 steps.
> > >
> > > Then the file 'inlist_finish_header' will finish the folloing evolution until the log Rho_c reaches
> > >
> > > 10.4 g cm^-3. However, in this part, the logT and logRho increases, but the nuclear burning
> > >
> > > in the outer shell seems to stopped. I am sure why.
> > >
> > > Here also attached a Kippenhahn diagram for M = 28 Msun.
> > >
> > > When logT_c =9.85, log t_final - t = -5, at the end of the calculation log t_final - t = -8.
> > >
> > > It seems that from the center to M_r=8 Msun, all the nuclear burning is stopped.
> > >
> > > I also checked the energy generation rate of C, O and Si burning in the profileXX.data
> > >
> > > they are 0 in this stage.
> > >
> > >
> > >
> > > Thanks,
> > >
> > > Wenyu Xin
> > >
> > > <inlist.zip><28M_kipp.pdf>_______________________________________________
> > > mesa-users at lists.mesastar.org
> > > https://lists.mesastar.org/mailman/listinfo/mesa-users
> > >
> </inlist.zip></mesa-users at lists.mesastar.org></mesa-users at lists.mesastar.org></fxt44 at mac.com></wyxin at nao.cas.cn></fxt44 at mac.com>
More information about the Mesa-users
mailing list