[Mesa-users] PISN models for high metallicity
Ebraheem Farag
ekfarag at asu.edu
Thu Jul 4 22:26:30 UTC 2024
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/20240704/1fcf83c9/attachment.htm>
More information about the Mesa-users
mailing list