mesa-users] mesa and late stages of evolution of giant stars

[Bill P. <pa...@ki...>]

On Jun 30, 2010, at 8:24 AM, Michelle Dolan wrote:

> late stages of evolution of giant stars

Hi Shelley,

It takes a rather long time for mesa to evolve a massive star from pre-ms to onset of core collapse.
Perhaps it would be useful to have a model that was saved about 1000 steps or so before collapse.
You can play with that a bit to see how things go before you do the entire evolution yourself.

Here's a model that started at 20 Msun and has evolved for 42000 timesteps.   


You can use it as a starting model by editing your inlist_massive:

         load_saved_model = .true.
         !saved_model_name = '25m_z2m2_pre_ms.mod'
         saved_model_name = '20_z2m2_si_burn.mod'


then just do
	
	./rn

you should get a terminal line saying

	load saved model 20_z2m2_si_burn.mod


and then it will start with output something like the following:

 
       step    lg_Tcntr    Teff       lg_LH     lg_Lnuc     Mass       H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr   pts  retry
   lg_dt_yr    lg_Dcntr    lg_Rsurf   lg_L3a    lg_Lneu     lg_Mdot    H_poor     He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  jacs  bckup
     age_yr    lg_Pcntr    lg_L       lg_LZ     lg_Psurf    lg_Dsurf   He_poor    C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
__________________________________________________________________________________________________________________________________________________
 
      42023   9.817492   3178.888  18.726917  26.365246  15.499273   9.327864   0.000000   0.000000   0.386320   0.353310   5.811636   1426      0
  -7.851571   9.034721   3.098539  15.152327  13.568666 -99.000000   6.171409   0.000917   0.000003   0.019955   0.383153   5.146621     10      0
 8.9557E+06  26.713394   5.159661  26.365246   2.503420  -8.754533   0.000000   0.000000   0.000000   0.000001  2.635E-01  0.189E+02    varcontrol



This plot shows the abundances



and this one shows that the infall velocity is about 1 km/sec -- the run will automatically stop when that reaches 1000 km/sec.


I hope that will be useful!

Cheers,
Bill