[Mesa-users] About nse.net

[辛文宇]

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/20211224/c74ad460/attachment.htm>