[Mesa-users] [EXT] Re: PISN models for high metallicity
辛文宇
xinwenyu16 at mails.ucas.ac.cn
Sat Jul 6 16:55:28 UTC 2024
Hi EbF,
Thanks for your so many suggestions very much.
I testing all of your comments.
1. turn off the mlt++ in inlist_hydro_on, failure.
2. set Pextra_factor = 1 in inlist_hydro_on, failure.
3. try switching from zoning in logR to zoning to logtau, failure.
4. okay_to_remesh = .false. in inlist_hydro_on.
The Riemann hydro is turn on in step 2205. After some steps of retries, log dt decrease to -6.
Then after ~200 steps, the crash takes place again.
I attached the logRho and entropy near the surface. There are small changes.
The inlist_hydro_on is also attached. Only this file is changed.
Do you have any other comments?
Thanks,
Wenyu
-----原始邮件-----
发件人:"Ebraheem Farag" <ekfarag at asu.edu>
发送时间:2024-07-06 01:19:07 (星期六)
收件人: "Mathieu Renzo" <mrenzo at arizona.edu>
抄送: mesa-users at lists.mesastar.org, xinwenyu16 at mails.ucas.ac.cn
主题: Re: [Mesa-users] [EXT] Re: PISN models for high metallicity
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 and more recently Farmer et al. 2019), although one could invoke B-fields to quench the winds at high Z and bypass this limit (e.g., Georgy et al. 2017).
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>
发送时间: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$
_______________________________________________ mesa-users at lists.mesastar.orghttps://lists.mesastar.org/mailman/listinfo/mesa-users
--
Mathieu Renzo
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/20240707/ebb1fbeb/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: logrho_entropy_surface.pdf
Type: application/pdf
Size: 24075 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240707/ebb1fbeb/attachment.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_hydro_on
Type: application/octet-stream
Size: 675 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240707/ebb1fbeb/attachment.obj>
More information about the Mesa-users
mailing list