[Mesa-users] Increased irradiation during planetary evolution
Marit Mol Lous
marit.mollous at gmail.com
Wed Aug 7 08:05:30 EDT 2019
Dear MESA users,
Using make_planets in test_suite I am trying to evolve irradiated planets
based on migration. With the subroutine that I wrote for this in
run_star_extras this effectively means that the model receives an increase
in flux from 6e4 erg/s to 6e8 erg/s over a timescale 1e5 yrs. I've noticed
that when I let this increase in radiation take place at a later time the
hydro_code fails, which I think is because the model has cooled to a lower
temperature and the sudden increase becomes relatively too high. But the
transition from a converging to a non-converging code seems very sudden to
me, in my case specifically the code seems to work fine if the increased
irradiation starts at 1.4e6 years, but crashes when it does at 1.5e6 yrs.
I've included three plots that show the received irradiation, logarithmic
timesteps and effective temperature during the evolution. For the
converging models there is a slight decrease in timestep during the highest
increase in radiation, but for the non-converging models this timestep
cascades. When I use *report_ierr = .true. *this turns out to be due to a
negative luminosity and the terminal shows the following:
*reduce dt because of lgTeff 195*
*hydro_call_number, s% dt, dt, dt/secyer, log dt/yr, tim 110
3.0092625467543285D+03 3.0092625467543285D+03 9.5356113976263702D-05
-4.0206514559370339D+00 1.0122168632241566D+14*
* model iter var zone corr norm max corr
xscale dx new-prev new prev
mass loc log dt/yr lgE lgT lgRho*
* 196 1 lnPgas 1 0.58351089E-03 0.20147429E-01
0.11314896E+02 0.22796607E+00 0.22796607E+00 0.11542862E+02
0.11314896E+02 0.92445190E-03 -0.40206515E+01 0.10889323E+02
0.29629335E+01 -0.56044180E+01*
* 196 1 coeff 1.0000 slope 0.000E+00 f 0.000E+00 avg resid
0.123E-04 max resid dv_dt 1 0.901E-01 avg corr 0.584E-03
max corr lnPgas 1 0.201E-01 lg dt/yr -4.02 avg+max corr+resid*
* 196 2 lnPgas 1 0.63288810E-02 -0.22002018E+00
0.11314896E+02 -0.24895055E+01 -0.24895055E+01 0.88253905E+01
0.11314896E+02 0.92445190E-03 -0.40206515E+01 0.10919772E+02
0.29921793E+01 -0.55346595E+01*
* get_surf_PT: L_surf <= 0 196
-8.5587266038618357D+28*
* L u kap logT logRho 1
-8.5587266038618357D+28 NaN 1.5045286385974722D-03
2.7343843484821400D+00 -6.3580430235840470D+00*
* L u kap logT logRho 2
-8.5583147776500534D+28 NaN 1.5045790089649306D-03
2.7343858068970479D+00 -6.3580317608739918D+00*
* L u kap logT logRho 3
-8.5579029640878850D+28 NaN 1.5046671590873324D-03
2.7343883591139426D+00 -6.3580120515545628D+00*
* L u kap logT logRho 4
-8.5568734855225617D+28 NaN 1.5047930921615925D-03
2.7343920051135533D+00 -6.3579838963149147D+00*
* L u kap logT logRho 5
-8.5558440860097758D+28 NaN 1.5049820013973831D-03
2.7343974740666934D+00 -6.3579416655149323D+00*
* L u kap logT logRho 6
-8.5537855241217367D+28 NaN 1.5053598544494447D-03
2.7344084117997891D+00 -6.3578572113232799D+00*
* L u kap logT logRho 7
-8.5496693487497867D+28 NaN 1.5062416912031706D-03
2.7344339323138827D+00 -6.3576601899745278D+00*
* L u kap logT logRho 8
-8.5393844408057417D+28 NaN 1.5075018906701839D-03
2.7344703877444689D+00 -6.3573788239340558D+00*
* L u kap logT logRho 9
-8.5291074287981111D+28 NaN 1.5093931518066987D-03
2.7345250662783829D+00 -6.3569569805611765D+00*
* L u kap logT logRho 10
-8.5085770725569926D+28 NaN 1.5131791290736154D-03
2.7346344061310695D+00 -6.3561140330404697D+00*
* logT =
2.7343843484821400D+00*
* logRho =
-6.3580430235840470D+00*
* x =
7.4000000000000021D-01*
* z =
1.9999999999999796D-02*
* abar =
1.2480945159693049D+00*
* zbar =
1.0858475417883164D+00*
* s% use_compression_outer_BC F*
* s% use_T_Paczynski_outer_BC F*
* s% use_T_black_body_outer_BC F*
* s% use_fixed_L_for_BB_outer_BC F*
* s% tau_for_L_BB
-1.0000000000000000D+00*
I am new to MESA and since there are many functions in model_planet that I
don't understand I was wondering if there is anything that I can enable or
disable to get this to converge without losing credibility, or if I am
maybe deliberately using something that is causing the crash. Or maybe I am
asking to much of the code and it has a fair reason to crash. I am also
curious why the hydro_code crashes for some of these models while it seems
to work fine if the irradiation starts to increase only slightly earlier,
and if this is a sign that something is wrong with these working models as
well.
For now I am using surface radiation only, no interior heating. I am using
the most recent mesa-r11701. I've also attached my inlist files (which are
barely different from the make_planet default except for the use of
*use_other_energy
= .true.* and evolving to 5e9 yrs) and run_star_extra. Hopefully I've
provided enough information and if someone could give insight into why the
code crashes so suddenly when it receives the radiation at a slightly later
time or point me towards something I am overlooking I'd be grateful.
Kind regards and thanks in advance,
Marit Mol Lous
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190807/78b3bdf4/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: F_rec.png
Type: image/png
Size: 40811 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190807/78b3bdf4/attachment-0003.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: T_eff.png
Type: image/png
Size: 39827 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190807/78b3bdf4/attachment-0004.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: log_dt.png
Type: image/png
Size: 30383 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190807/78b3bdf4/attachment-0005.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_create
Type: application/octet-stream
Size: 1449 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190807/78b3bdf4/attachment-0006.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_create_1.0_MJ_10.0_ME_2.0_RJ
Type: application/octet-stream
Size: 3815 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190807/78b3bdf4/attachment-0007.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_core_1.0_MJ_10.0_ME_2.0_RJ
Type: application/octet-stream
Size: 2291 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190807/78b3bdf4/attachment-0008.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_evolve
Type: application/octet-stream
Size: 1435 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190807/78b3bdf4/attachment-0009.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_core
Type: application/octet-stream
Size: 1624 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190807/78b3bdf4/attachment-0010.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_evolve_1.0_MJ_10.0_ME_2.0_RJ
Type: application/octet-stream
Size: 4512 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190807/78b3bdf4/attachment-0011.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: run_star_extras.f
Type: text/x-fortran
Size: 11461 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190807/78b3bdf4/attachment-0001.bin>
More information about the Mesa-users
mailing list