[Mesa-users] Problems With High Irradiation

[Francis Timmes]

hi tuhin,

the 2017 labs used version 9793. is this the version you are running?
if not, my suggestion is to get the lab running (which includes solutions) 
in version 9793 and then port it to, say, the current public version 10398.

fxt





> On Nov 19, 2018, at 10:50 AM, tuhin.ghosh via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> 
> 
> Hi everyone,
> 
> I am trying to make highly irradiated sub-Neptune size planetary models. For this I am 
> closely following the Summer School 2017 instructions by 
> Leslie A. Rogers. But when I tried to run my models  for core mass 4 M_earth and 20 % 
> envelope at very close distance (0.01 A.U.) from a sun 
> twin star, I always ran into some trouble.
> 
> I used this formula flux_dayside = Lsun/(4*pi*d**2)
> 
> For :
> 
> years_for_initial_dt = 1.d-6
> relax_irradiation_max_yrs_dt = 1d-1
> irrad_col_depth =  1 !  g cm^-2
> set_to_this_tau_factor = 100
> 
> I get the following at a certain point
> 
>         relax irradiation, model, step, fraction done, flux  1052    51      0.510E+00   
> 0.696E+10
> hydro_newton_step failed to accept
>       retry log10(dt/yr), log10(dt), retry_factor     1  1053       -1.602060        
> 5.897052        0.500000
> hydro_newton_step failed to accept
>       retry log10(dt/yr), log10(dt), retry_factor     2  1053       -1.903090        
> 5.596022        0.500000
> have reached retry limit so now backup
>       retry log10(dt/yr), log10(dt), retry_factor     1  1052       -2.204120        
> 5.294992        0.500000
>       retry log10(dt/yr), log10(dt), retry_factor     2  1052       -2.505150        
> 4.993962        0.500000
> have reached retry limit so now backup
>       retry log10(dt/yr), log10(dt), retry_factor     1  1052       -3.408240        
> 4.090872        0.500000
>       retry log10(dt/yr), log10(dt), retry_factor     2  1052       -3.709270        
> 3.789842        0.500000
> have reached retry limit so now backup
> hydro_newton_step failed to accept
>       retry log10(dt/yr), log10(dt), retry_factor     1  1052       -5.214420        
> 2.284692        0.500000
> hydro_newton_step failed to accept
>       retry log10(dt/yr), log10(dt), retry_factor     2  1052       -5.515450        
> 1.983662        0.500000
> have reached retry limit so now backup
> hydro_newton_step failed to accept
>       retry log10(dt/yr), log10(dt), retry_factor     1  1052       -7.622660       
> -0.123548        0.500000
>     relax irradiation, model, step, frac, flux, wait for dt  1052    51      0.510E+00   
> 0.696E+10
> 
> 
> The terminal kept showing me this untill I manually stop the process after many models 
> :
> 
> relax irradiation, model, step, frac, flux, wait for dt  3552    51      0.510E+00   
> 0.696E+10
> 
> 
> Now when I used irrad_col_depth =  10 (or 100) I get similar output but this time it is 
> stuck at 11% (or 1%).
> 
> However when I reduce tau factor drastically (irrad_col_depth =  0.01 !  g cm^-2) the 
> relax irradiation process gets done but then while gets 
> stuck at evolution :
> 
> _______________________________________________________________________________________
> _______________________________________________________
> ____
> 
>       step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass       H_rich     
> H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones 
> retry
>   lg_dt_yr    lg_Dcntr    lg_R       lg_L3a    lg_Lneu     lg_Mdot    He_core    
> He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters 
> bckup
>     age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf   C_core     
> C_cntr     Ne_cntr    Z_cntr     Z_avg     v_div_cs     
> dt_limit
> _______________________________________________________________________________________
> _______________________________________________________
> ____
> 
>       2700   3.875057   1753.615 -25.432549 -25.432549   0.000015   0.000015   
> 0.720000   0.001512   0.250000   0.720000 -19.888005    673      
> 6
>  -4.817735  -0.606211  -0.910594 -99.000000 -26.830359 -99.000000   0.000000   
> 0.250000   0.014041   0.030000   0.250000  21.117045      4      
> 2
> 1.0000E+06  11.179237  -3.891969 -99.000000   4.574482  -6.906418   0.000000   
> 0.005162   0.003149  3.000E-02  3.000E-02  0.146E-01    
> varcontrol
> 
> 
> As you can see the lg_dt_yr is very low and and it oscillates around and the process 
> stops due to reaching maximum model number. I have tested 
> this 
> by increasing the maximum model number significantly, but no effect.
> 
> 
> So I want to ask that what am I doing wrong. I also want to ask what should be the 
> optimal values of these parameters (years_for_initial_dt, 
> relax_irradiation_max_yrs_dt, irrad_col_depth,set_to_this_tau_factor) for highly 
> irradiated situation and in general how they relate to 
> successful compilation of the code and what should I keep in mind while changing these 
> parameters.
> 
> I know I am asking for a lot but any help would be greatly appreciated.
> 
> 
> Yours sincerely,
> Tuhin Ghosh
> <inlist_relaxirrad_4.0_ME_0.2_9_1.0_Lsun_0.01_AU_><planet_removeLc_4.0_ME_0.2_9.mod>_______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
>