[Mesa-users] [EXT] Re: PISN models for high metallicity
Ebraheem Farag
ekfarag at asu.edu
Fri Jul 5 17:19:07 UTC 2024
Hi Wenyu,
Looking at your inlists and your "120M.log" file I have a few comments.
It appears your model is crashing immediately after Riemann hydro turns on.
Riemann hydro turns on at timestep 2207, and the model crashes at timestep
2208. Check that mlt++ is turn off entirely when Riemann hydro turns on
(perhaps by placing controls in 'inlist_hydro_on') as you are still reading
inlist_before at logTcen ~ 9d0.
Try returning Pextra_factor to 1 when Riemann hydro turns on if you think
that will help.
Try turning off remeshing for some few time steps to see if that helps
"okay_to_remesh = .false.", this will smooth the transition from cell face
to cell center velocities.
As a last resort, try switching from zoning in logR to zoning to logtau -->
set :
split_merge_amr_log_zoning = .false.
split_merge_amr_logtau_zoning = .true. ! adjust
'split_merge_amr_nz_baseline' accordingly.
Happy debugging :)
-EbF
On Fri, Jul 5, 2024 at 2:47 AM Mathieu Renzo via Mesa-users <
mesa-users at lists.mesastar.org> wrote:
> Hi Wenyu,
>
> I agree with Eb that this smells like a surface issue -- likely opacity
> related since it appears at large Z but not at lower Z. I suppose (also
> based on the Kippenhahn you sent) the issue is occurring during the first
> PPISN pulse.
>
> As a test, have you tried running a model with a H-rich envelope? This
> would change quite a bit the surface properties (unless your wind setup is
> such that at these mass and metallicity the H-rich envelope is removed by
> radiatively driven winds before the pulses). If you find that models get
> farther in their evolution (maybe though the first pulse) with an H rich
> envelope, it would strengthen the suspicion that it's boundary layers that
> are the problem.
>
> In that case, there are *many* things one could try to play with, from
> outer boundary conditions, to wind mass loss rates pre-pulse, to doing the
> re-mapping to a new star earlier. But you would need to carry out extensive
> tests of any of these, which can impact the evolution significantly,
> introducing numerical artifacts.
>
> Also, you probably are aware that at high Z PPISN is not expected (cf. Langer
> et al. 2007
> <https://urldefense.com/v3/__https://ui.adsabs.harvard.edu/abs/2007A*26A...475L..19L/abstract__;JQ!!IKRxdwAv5BmarQ!eJty_qjigdgXhcYFpBCRby3xzoKN9OJCcn86gVf0wSVpe_e_Kb0JU25_YIHhQqkaapsl5pq1BAmkWVxes_6o8S27tiFJVA$>
> and more recently Farmer et al. 2019
> <https://urldefense.com/v3/__https://doi.org/10.3847/1538-4357/ab518b__;!!IKRxdwAv5BmarQ!eJty_qjigdgXhcYFpBCRby3xzoKN9OJCcn86gVf0wSVpe_e_Kb0JU25_YIHhQqkaapsl5pq1BAmkWVxes_6o8S2MagqY9A$>),
> although one could invoke B-fields to quench the winds at high Z and bypass
> this limit (e.g., Georgy et al. 2017
> <https://urldefense.com/v3/__https://ui.adsabs.harvard.edu/abs/2017A*26A...599L...5G/abstract__;JQ!!IKRxdwAv5BmarQ!eJty_qjigdgXhcYFpBCRby3xzoKN9OJCcn86gVf0wSVpe_e_Kb0JU25_YIHhQqkaapsl5pq1BAmkWVxes_6o8S0yHF-5Og$>).
>
>
> Hope this helps, and good luck!
> Mathieu
>
> On 7/5/24 10:15, 辛文宇 wrote:
>
>
>
> * External Email *
> ------------------------------
> 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> <ekfarag at asu.edu>
> *发送时间:*2024-07-05 06:26:30 (星期五)
> *收件人:* "Frank Timmes" <fxtimmes at gmail.com> <fxtimmes at gmail.com>
> *抄送:* mesa-users <mesa-users at lists.mesastar.org>
> <mesa-users at lists.mesastar.org>, "辛文宇" <xinwenyu16 at mails.ucas.ac.cn>
> <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$
>>
>>
> _______________________________________________mesa-users at lists.mesastar.orghttps://lists.mesastar.org/mailman/listinfo/mesa-users <https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!IKRxdwAv5BmarQ!eJty_qjigdgXhcYFpBCRby3xzoKN9OJCcn86gVf0wSVpe_e_Kb0JU25_YIHhQqkaapsl5pq1BAmkWVxes_6o8S1vjR60vg$>
>
>
> --
> Mathieu Renzo
> <https://urldefense.com/v3/__https://users.flatironinstitute.org/*mrenzo/__;fg!!IKRxdwAv5BmarQ!eJty_qjigdgXhcYFpBCRby3xzoKN9OJCcn86gVf0wSVpe_e_Kb0JU25_YIHhQqkaapsl5pq1BAmkWVxes_6o8S0IR04qPQ$>
> Assistant Professor
> University of Arizona & Steward Observatory
> _______________________________________________
> mesa-users at lists.mesastar.org
>
> https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!IKRxdwAv5BmarQ!eJty_qjigdgXhcYFpBCRby3xzoKN9OJCcn86gVf0wSVpe_e_Kb0JU25_YIHhQqkaapsl5pq1BAmkWVxes_6o8S1vjR60vg$
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240705/83164174/attachment.htm>
More information about the Mesa-users
mailing list