[Mesa-users] Small timesteps in binary evolution towards LMXB formation
Ndiogou Niang
niang at sabanciuniv.edu
Sat Mar 23 16:27:27 UTC 2024
Dear MESA users,
I am evolving a 13 M_sun with a 1.6 M_sun system in
.../binary/test_suite/evolve_both_stars, hoping to simulate the formation
of a low mass X-ray binary.
However, I keep getting small timesteps (down to log dt < -8 ). The
timesteps starts to decrease following the burning of Oxygen.
I am attaching the inlists as well as a file showing the last line of the
terminal before the run stopped.
I would really appreciate it if I could receive some guidance on this
problem.
Thank you.
Regards
--
*Ndiogou NIANG *
*Physics, PhD Candidate*
*Faculty of Engineering and Natural SciencesSabancı University*
*Istanbul, Turkey*
*+90 536 918 06 92*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240323/e68823ec/attachment.htm>
-------------- next part --------------
retry 875
retry: max resid -- give up in solver 875
1 875 7.283440 7524.544 0.809386 0.809386 1.600000 1.600000 0.692985 0.005108 0.280000 -2.288336 915 59
-1.0992E+01 7.283440 0.174419 -39.574814 -0.637617 -99.000000 0.000000 0.286460 0.009260 0.020000 0.060537 2
1.5446E+07 1.898417 0.809452 -99.000000 -99.000000 -7.736750 0.000000 0.000026 0.002085 0.020555 0.000E+00 retry
2 875 9.225483 3593.530 -6.686357 10.948636 13.000000 9.561765 0.000000 0.000000 0.280000 7.740749 2025 59
-1.0992E+01 9.036058 2.751667 4.929041 8.249909 -99.000000 3.438235 0.000000 0.508312 0.020000 1.812912 10
1.5446E+07 7.349569 4.680112 10.948636 11.304865 -8.756789 1.696966 0.000000 0.452633 0.000537 0.000E+00 retry
bin 875 14.600000 1.687E+03 2.100E+03 0.000E+00 8.125000 0 1 0.000E+00 1.000000 1.351E+54 -1.337E+27 0.000E+00
-10.992126 1.600000 1.494235 0.000000 0.000E+00 36.182142 3.703E+02 -9.960E-01 0.000E+00 1.000E+99 0.000E+00 -1.337E+27 0.000E+00
1.5446E+07 13.000000 5.645E+02 0.000000 -2.339E+46 4.453187 9.471E+02 -4.039E-01 0.000E+00 0.000E+00 0.000E+00 0.000E+00 1
max_residual > tol_max_residual 876 2.8735439470310622D-05 1.0000000000000001D-05
solver rejected trial model
s% model_number 876
s% solver_call_number 935
nz 2022
s% num_retries 59
dt 3.2134990335632081D-04
log dt/secyer -1.0992125793449739D+01
retry 876
retry: max resid -- give up in solver 876
1 876 7.283440 7524.544 0.809386 0.809386 1.600000 1.600000 0.692985 0.005108 0.280000 -2.288336 915 60
-1.1293E+01 7.283440 0.174419 -39.574814 -0.637617 -99.000000 0.000000 0.286460 0.009260 0.020000 0.060537 2
1.5446E+07 1.898417 0.809452 -15.051500 -99.000000 -7.736750 0.000000 0.000026 0.002085 0.020555 0.000E+00 retry
2 876 9.225483 3593.530 -6.686361 10.948636 13.000000 9.561765 0.000000 0.000000 0.280000 7.740749 2022 60
-1.1293E+01 9.036058 2.751667 4.929041 8.249909 -99.000000 3.438235 0.000000 0.508312 0.020000 1.812912 7
1.5446E+07 7.349569 4.680112 10.948636 11.304865 -8.756789 1.696966 0.000000 0.452633 0.000537 0.000E+00 retry
bin 876 14.600000 1.687E+03 2.100E+03 0.000E+00 8.125000 0 1 0.000E+00 1.000000 1.351E+54 -1.337E+27 0.000E+00
-11.293156 1.600000 1.494235 0.000000 0.000E+00 36.182142 3.703E+02 -9.960E-01 0.000E+00 1.000E+99 0.000E+00 -1.337E+27 0.000E+00
1.5446E+07 13.000000 5.645E+02 0.000000 -2.339E+46 4.453187 9.471E+02 -4.039E-01 0.000E+00 0.000E+00 0.000E+00 0.000E+00 1
solver rejected trial model
s% model_number 877
s% solver_call_number 937
nz 2037
s% num_retries 60
dt 1.6067495167816040D-04
log dt/secyer -1.1293155789113721D+01
retry 877
retry: max residual jumped -- give up in solver 877
hydro_solver_step returned ierr -1
s% model_number 877
nz 2037
s% num_retries 61
retry 877
retry: logT > hydro_mtx_max_allowed_logT 1142 877
hydro_solver_step returned ierr -1
s% model_number 877
nz 2037
s% num_retries 62
retry 877
retry: logT > hydro_mtx_max_allowed_logT 1142 877
hydro_solver_step returned ierr -1
s% model_number 877
nz 2037
s% num_retries 63
retry 877
retry: logT > hydro_mtx_max_allowed_logT 1143 877
hydro_solver_step returned ierr -1
s% model_number 877
nz 2037
s% num_retries 64
retry 877
retry: logT > hydro_mtx_max_allowed_logT 1140 877
hydro_solver_step returned ierr -1
s% model_number 877
nz 2037
s% num_retries 65
retry 877
retry: logT < hydro_mtx_min_allowed_logT 464 877
solver rejected trial model
s% model_number 877
s% solver_call_number 943
nz 2037
s% num_retries 66
dt 2.5105461199712563D-06
log dt/secyer -1.3099335763097608D+01
retry 877
retry: avg corr too large -- give up in solver 877
solver rejected trial model
s% model_number 877
s% solver_call_number 944
nz 2037
s% num_retries 67
dt 1.2552730599856282D-06
log dt/secyer -1.3400365758761589D+01
retry 877
dt 6.2763652999281408D-07
min_timestep_limit 9.9999999999999995D-07
stopping because of problems dt < min_timestep_limit
retry: avg corr too large -- give up in solver 877
dt 6.2763652999281408D-07
min_timestep_limit 9.9999999999999995D-07
stopping because of problems dt < min_timestep_limit
bin 876 14.600000 1.687E+03 2.100E+03 0.000E+00 8.125000 0 1 0.000E+00 1.000000 1.351E+54 -1.337E+27 0.000E+00
-11.293156 1.600000 1.494235 0.000000 0.000E+00 36.182142 3.703E+02 -9.960E-01 0.000E+00 1.000E+99 0.000E+00 -1.337E+27 0.000E+00
1.5446E+07 13.000000 5.645E+02 0.000000 -2.339E+46 4.453187 9.471E+02 -4.039E-01 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0
save photos/b_x876, photos/1_x876, photos/2_x876 for model 876
runtime (minutes), retries, steps 50.57 67 877
terminated evolution: cannot find acceptable model
termination code: min_timestep_limit
terminated evolution: cannot find acceptable model
termination code: min_timestep_limit
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_project
Type: application/octet-stream
Size: 428 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240323/e68823ec/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist1
Type: application/octet-stream
Size: 871 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240323/e68823ec/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist2
Type: application/octet-stream
Size: 791 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240323/e68823ec/attachment-0002.obj>
More information about the Mesa-users
mailing list