[Mesa-users] PISN models for high metallicity
辛文宇
xinwenyu16 at mails.ucas.ac.cn
Fri Jul 5 08:15:02 UTC 2024
Hi, fxt,
Thank you very much for giving me this paper. In fact I have noticed this paper before this email.
And I also use the similar resolution settings with this paper.
As I have mentioned, for the models with Z<1d-3, the models works well.
The crashes happen only for Z>.1d-3 case.
> maybe in future postings remove the star executable before assembling a zip archive of a mesa work directory (wink).
OK, thanks for your reminder.
Hi, EbF,
Also thanks for you suggestions. I attached the logRho and entropy as a function of Mr near the surface before crashes,
as well as the mass fraction of several isotopes. I also show the kippenhan diagram.
There is a sharply increase (decrease) in logRho (entropy). Do the crashes occur because this？
I guess this sharply increase in logRho is because I set Pextra_factor = 2, and I need change is as defacult at the end
of He burning?
Thanks,
Wenyu
-----原始邮件-----
发件人:"Ebraheem Farag" <ekfarag at asu.edu>
发送时间:2024-07-05 06:26:30 (星期五)
收件人: "Frank Timmes" <fxtimmes at gmail.com>
抄送: mesa-users <mesa-users at lists.mesastar.org>, "辛文宇" <xinwenyu16 at mails.ucas.ac.cn>
主题: Re: [Mesa-users] PISN models for high metallicity
I can’t run 11701 or 12778 on my current machine, but i would suggest looking at the surface regions of the trho profile to figure out if it it’s just the surface going haywire, especially at higher metallicity. It’s good to see what’s going on before the model crashes.
-EbF
On Thu, Jul 4, 2024 at 4:00 PM Frank Timmes via Mesa-users <mesa-users at lists.mesastar.org> wrote:
hi wenyu,
my first guess is mass or time resolution. there might be some useful resolution settings to consider in https://urldefense.com/v3/__https://zenodo.org/records/6930577__;!!IKRxdwAv5BmarQ!bQd8zAP4sa19edBT_5ot5y-ETo-NtUApGFIVwnCdA6mSRTndGfv0e4ehJKYCLY2Tpmgdy2zyO907bq5vgdSWJggObPg5ig$ , with an accompanying article https://urldefense.com/v3/__https://ui.adsabs.harvard.edu/abs/2022ApJ...937..112F/abstract__;!!IKRxdwAv5BmarQ!bQd8zAP4sa19edBT_5ot5y-ETo-NtUApGFIVwnCdA6mSRTndGfv0e4ehJKYCLY2Tpmgdy2zyO907bq5vgdSWJgiZuHIOyg$ . these models were run with r11701, somewhat close to your r12778, but at a smaller metallicity.
maybe in future postings remove the star executable before assembling a zip archive of a mesa work directory (wink).
fxt
> On Jul 2, 2024, at 3:10 AM, 辛文宇 <xinwenyu16 at mails.ucas.ac.cn> wrote:
>
> Hi everyone,
>
> I am calculating the PISN models with ppisn example in test_suite and in mesa-r12778.
> For Z<1E-3, the models work well. However, when I use Z>1d-3, the models stop at the end of He burning.
> I set okay_to_reduce_gradT_excess = .true. to make the models go through pre-MS, which is turn off at
> power_he_burn > 0.9 L_surf. I also use Pextra_factor = 2 to ensure the model go through He burning.
> However, when T_c reaches 8.988, there are some retries like:
> hydro_mtx: logRho too large 2208 1258 1 1.2040579673786977D+01 4.2901366099600260D+00
> hydro_mtx: logRho too large 2208 1335 1 1.2588366322584088D+01 4.3060459950999634D+00
> hydro_mtx: logRho too large 2208 1258 1 1.2040579673786977D+01 4.2901366099600260D+00
> hydro_mtx: logRho too large 2208 1335 1 1.2588366322584088D+01 4.3060459950999634D+00
> set_vars_for_solver failed: model, nz 2208 1533
> eval_equations: set_newton_vars returned ierr -1
> Finally, the error are as follows:
> Error termination. Backtrace:
>
> Could not print backtrace: unrecognized DWARF version in .debug_info at 6
> #0 0xa90697 in __utils_lib_MOD_mesa_error
> at ../public/utils_lib.f:1021
> #1 0x4f47e2 in get_t_tau
> at ../private/atm_support.f90:310
> #2 0x4f70cc in __atm_support_MOD_get_atm_pt
> at ../private/atm_support.f90:127
> #3 0x50ab78 in __hydro_vars_MOD_get_surf_pt
> at ../private/hydro_vars.f90:1459
> #4 0x50c068 in __hydro_vars_MOD_set_teff_info_for_eqns
> at ../private/hydro_vars.f90:936
> #5 0x8599e5 in pt_eqns_surf
> at ../private/hydro_eqns.f90:2372
> #6 0x8599e5 in eval_equ_for_solver
> at ../private/hydro_eqns.f90:376
> #7 0x85b5bf in __hydro_eqns_MOD_eval_equ
> at ../private/hydro_eqns.f90:64
> #8 0x81f640 in __hydro_newton_procs_MOD_eval_equations
> at ../private/hydro_newton_procs.f90:143
> #9 0x7f82d7 in adjust_correction
> at ../private/star_newton.f90:955
> #10 0x7f82d7 in do_newton
> at ../private/star_newton.f90:631
> #11 0x7f82d7 in __star_newton_MOD_newton
> at ../private/star_newton.f90:183
> #12 0x799795 in newt
> at ../private/solve_hydro.f90:830
> #13 0x799f0d in hydro_newton_step
> at ../private/solve_hydro.f90:766
> #14 0x79bd7e in do_hydro_newton
> at ../private/solve_hydro.f90:163
> #15 0x79bd7e in __solve_hydro_MOD_do_hydro_converge
> at ../private/solve_hydro.f90:163
> #16 0x79eb54 in __struct_burn_mix_MOD_do_struct_burn_mix
> at ../private/struct_burn_mix.f90:207
> #17 0x5d793e in __evolve_MOD_do_evolve_step_part2
> at ../private/evolve.f90:599
> #18 0x4220fd in __star_lib_MOD_star_evolve_step
> at ../public/star_lib.f90:550
> #19 0x43d091 in __run_star_support_MOD_run1_star
> at ../job/run_star_support.f90:130
> #20 0x41b67f in __run_star_MOD_do_run_star
> at /home/wyxin/work/mesa/12778/mesa-r12778/star/job/run_star.f:26
> #21 0x41b746 in run
> at ../src/run.f:16
> #22 0x41b788 in main
> at ../src/run.f:2
> 2.5691E+05 21.505411 7.756651 7.400769 7.643802 -5.045201 0.000000 0.050090 0.087833 0.000177 9.254E-01 -0.221E-01 hold
> log_rel_run_E_err -4.7322703837780633
> rel_E_err -8.2532431125648629D-07
>
> check gamma integral 2.6093959149115335E-002
> newton solver rejected trial model
> s% model_number 2210
> s% hydro_call_number 2953
> nz 2396
> s% num_retries 4
> s% num_backups 0
> s% number_of_backups_in_a_row 0
> log dt/secyer -3.3598229247300733D+00
>
> retry 2210
> retry log10(dt/yr), log10(dt), retry_factor 1 2210 -3.961883 3.537229 0.500000
> get_T_tau: L <= 0 2210 -5.6720755514084073D+42
> File: ../private/atm_support.f90, Line: 310
> DATE: 2024-07-02
> TIME: 08:57:57
>
> I am not sure how to make the models go through the further evolution.
> Is there anyone can give me some suggestions?
> The output file and inlist files are attached.
>
> Thanks,
> Wenyu
> <ppisn2.zip>_______________________________________________
> mesa-users at lists.mesastar.org
> https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!IKRxdwAv5BmarQ!bQd8zAP4sa19edBT_5ot5y-ETo-NtUApGFIVwnCdA6mSRTndGfv0e4ehJKYCLY2Tpmgdy2zyO907bq5vgdSWJgjNV7oIZg$
>
_______________________________________________
mesa-users at lists.mesastar.org
https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!IKRxdwAv5BmarQ!bQd8zAP4sa19edBT_5ot5y-ETo-NtUApGFIVwnCdA6mSRTndGfv0e4ehJKYCLY2Tpmgdy2zyO907bq5vgdSWJgjNV7oIZg$
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240705/21aee521/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: kipp_M120.pdf
Type: application/pdf
Size: 782036 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240705/21aee521/attachment-0004.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: logrho_entropy_surface.pdf
Type: application/pdf
Size: 21151 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240705/21aee521/attachment-0005.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mass_frac_surface.pdf
Type: application/pdf
Size: 22791 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240705/21aee521/attachment-0006.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: logrho_entropy.pdf
Type: application/pdf
Size: 21374 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240705/21aee521/attachment-0007.pdf>
More information about the Mesa-users
mailing list