[Mesa-users] Increased irradiation during planetary evolution
Marit Mol Lous
marit.mollous at gmail.com
Thu Aug 8 05:47:54 EDT 2019
One more thing: I can get the evolution to complete using *delta_lgTeff_limit
= 1.d0 *in the &controls of inlist_evolve, but that is still suspicious as
it requires timesteps much smaller than the other models, see the figure
attached.
Op wo 7 aug. 2019 om 14:05 schreef Marit Mol Lous <marit.mollous at gmail.com>:
> 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/20190808/851d9e86/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: log_dt_lgTeff_limit.png
Type: image/png
Size: 30877 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190808/851d9e86/attachment-0001.png>
More information about the Mesa-users
mailing list