[Mesa-users] About nse.net

Thu Dec 30 15:37:50 UTC 2021

Hi
The network is too big for mesa to handle. You need to remove some
isotopes, getting the size down to around 300 should work.

Rob

On Fri, 24 Dec 2021 at 05:49, 辛文宇 <wyxin at nao.cas.cn> wrote:

> Dear Farmer,
>
>
> I calculate the models with mesa-r15140 now. When the Tc reaches 5E9 K
> (log Tc~9.7), I try to change the mesa_128.net to nse.net.
>
> But it shows the errors as follows.
>
>
> __________________________________________________________________
> ______________________________
>
>        step    lg_Tmax     Teff     lg_LH      lg_Lnuc     Mass
>  H_rich     H_cntr     N_cntr     Y_surf   eta_cntr   zones  retry
>    lg_dt_yr    lg_Tcntr    lg_R     lg_L3a     lg_Lneu     lg_Mdot
> He_core    He_cntr    O_cntr     Z_surf   gam_cntr   iters
>      age_yr    lg_Dcntr    lg_L     lg_LZ      lg_Lphoto   lg_Dsurf
>  C_core     C_cntr     Ne_cntr    Si_cntr  v_div_cs       dt_limit
>
> _________________________________________________________________________________________________
>
> free(): corrupted unsorted chunks
> free(): invalid pointer
>
> Program received signal SIGABRT: Process abort signal.
>
> Backtrace for this error:
>
> Program received signal SIGABRT: Process abort signal.
>
> Backtrace for this error:
> #0  0x7f6171abe20f in ???
> #0  0x7f6171abe20f in ???
> #1  0x7f6171abe18b in ???
> #1  0x7f6171abe18b in ???
> #2  0x7f6171a9d858 in ???
> #2  0x7f6171a9d858 in ???
> #3  0x7f6171b083ed in ???
> #3  0x7f6171b083ed in ???
> #4  0x7f6171b1047b in ???
> #4  0x7f6171b1047b in ???
> #5  0x7f6171b11cab in ???
> #5  0x7f6171b121c1 in ???
> #6  0x7f61749b1a64 in ???
> #6  0x7f6171b15a0c in ???
> #7  0x7f6171c72241 in free_stacks
> at /build/glibc-eX1tMB/glibc-2.31/nptl/allocatestack.c:277
> #8  0x7f6171c73521 in queue_stack
> at /build/glibc-eX1tMB/glibc-2.31/nptl/allocatestack.c:311
> #9  0x7f6171c73521 in __deallocate_stack
> at /build/glibc-eX1tMB/glibc-2.31/nptl/allocatestack.c:801
> #10  0x7f6171c73521 in __free_tcb
> at /build/glibc-eX1tMB/glibc-2.31/nptl/pthread_create.c:368
> #11  0x7f6171c73881 in start_thread
> at /build/glibc-eX1tMB/glibc-2.31/nptl/pthread_create.c:573
> #12  0x7f6171b9a292 in ???
> #13  0xffffffffffffffff in ???
> ./rn1: line 6: 2881595 Aborted                 (core dumped) ./star
> DATE: 2021-12-24
> TIME: 12:03:43
>
> Thanks,
>
> Wenyu
>
>
> -----Original Messages-----
> *From:*"Rob Farmer" <robert.j.farmer37 at gmail.com>
> *Sent Time:*2021-12-14 03:01:48 (Tuesday)
> *To:* "Francis Timmes" <fxt44 at mac.com>
> *Cc:* "辛文宇" <wyxin at nao.cas.cn>, mesa-users <mesa-users at lists.mesastar.org>
> *Subject:* Re: [Mesa-users] About nse.net
>
> Hi Wenyu,
>
> The likely problem you have is that with nse.net you are running into the
> rate 55cu(g,p)54ni which is an exothermic photo disintegration.
>
> In mesa 12778 this is not handled properly. We would use what reaclib
> stated was the forward rate to compute the reverse rate in detailed
> balance. However this definition is wrong for exothermic photo
> disintegrations leading to bogus rates which causes the solver issues (the
> rate is ~10^100 at low temperatures and decreases with temperature!).
>
> The solution is to upgrade to mesa 15140 where this is fixed (by ignoring
> what reaclib says and defining the reverse rate as the endothermic
> reaction).
>
> Rob
>
> On Mon, 13 Dec 2021 at 18:30, Francis Timmes via Mesa-users <
> mesa-users at lists.mesastar.org> wrote:
>
>> hi wenyu,
>>
>> 5e9 k is probably too hot to be mapping in a change in the reaction
>> network.
>> perhaps try changing the network when the central temperature is a bit
>> cooler, say 2e9 k.
>>
>> fxt
>>
>>
>>
>>
>>
>> > On Dec 12, 2021, at 9:30 PM, 辛文宇 via Mesa-users <
>> mesa-users at lists.mesastar.org> wrote:
>> >
>> > Dear everyone,
>> >
>> >
>> >
>> > I use the example of "25M_pre_ms_to_core_collapse" in 12778 version to
>> evolve stars
>> >
>> > from ZAMS to Fe core collapse. When the central T reaches 5E+9 (log
>> Tc=9.7), I want
>> >
>> > change the nuclear network into "nse.net". The network in the former
>> stages I used is
>> >
>> > mesa_128.net. The former evolution is smoothly. However, when the
>> network changes
>> >
>> > to nse.net. The following problem is shown:
>> >
>> > ======================================================================
>> >
>> >  netint failed: abs(h).lt.stpmin   3.4753213900951310E-029
>>  9.9999999999999997E-029
>> >  netint failed: abs(h).lt.stpmin   2.8806830949925469E-029
>>  9.9999999999999997E-029
>> >  retry    1012
>> >                            1st model retry log10(dt/yr)
>>  -1.9254117971048657D+01
>> >  netint failed: abs(h).lt.stpmin   3.4753213900951310E-029
>>  9.9999999999999997E-029
>> >  netint failed: abs(h).lt.stpmin   2.8806830949925469E-029
>>  9.9999999999999997E-029
>> >  retry    1012
>> >                                                      dt
>> 8.7894878303463986D-13
>> >                                      min_timestep_limit
>> 9.9999999999999998D-13
>> >
>> >  stopping because of problems dt < min_timestep_limit
>> >
>> >  terminated evolution: cannot find acceptable model
>> >
>> >  termination code: min_timestep_limit
>> >
>> > ======================================================================
>> >
>> > The inlist file the detail output information file are attached.
>> >
>> > I am not sure how to modify the file and makes the evolution
>> >
>> > successfully.
>> >
>> >
>> >
>> > Thanks,
>> >
>> > Wenyu
>> >
>> >
>> >
>> >
>> <inlist_to_lgT_9.7><inlist_to_lgT_9.9><slurm-364461.out>_______________________________________________
>> > mesa-users at lists.mesastar.org
>> > https://lists.mesastar.org/mailman/listinfo/mesa-users
>> >
>>
>> _______________________________________________
>> mesa-users at lists.mesastar.org
>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20211230/7ea69d1f/attachment.htm>