[mesa-users] Convergence of models in the 11 - 14 Msun range
Manos Chatzopoulos
manolis at flash.uchicago.edu
Thu Jun 26 16:51:06 EDT 2014
Dear Bill et al.,
Attached is the inlist I am using, with the latest version of MESA
(6596) compiled with the latest version of MESASDK
for MacOSX. Running on 4 threads. Thank you for looking into this, in
the meantime I will try to play around with
the overshooting parameter as you mentioned.
Manos
> Hi Manos,
>
> Please send inlists and all details necessary to reproduce a sample case.
> I'd like to know how far the evolution has gone when convergence fails.
> In the past, this mass range has ignited Ne20 off-center after C/O burning.
> It then takes a very large number of timesteps for the burning front to advance inward to the center.
> The advanced is limited by not being able to deplete the Ne20 over a significant mass range at once.
> With convection that works to help the burning go relatively quickly.
> For inward burning, there is no help from mixing --- unless you crank up overshooting below the burn.
> Without the help of overshooting to drive things inward faster, it has to creep inward bit by bit.
>
> -Bill
>
>
>
>
>
>
> On Jun 26, 2014, at 1:38 PM, Manos Chatzopoulos wrote:
>
>> Dear all,
>>
>> Has anyone managed to have a converged MESA model all the way from
>> pre-MS to pre-supernova (core-collapse)
>> in this rather interesting mass range of 11 - 14 Msun? Stars around this
>> mass range are thought to make electron-capture
>> SN progenitors with O/Ne/Mg cores, but when I try running them they
>> always seem to fail to converge past
>> core C/O burning. Is this relevant to the core becoming significantly
>> degenerate? Any insight on this would be
>> appreciated.
>>
>> Best,
>>
>> Manos
>>
>> ------------------------------------------------------------------------------
>> Open source business process management suite built on Java and Eclipse
>> Turn processes into business applications with Bonita BPM Community Edition
>> Quickly connect people, data, and systems into organized workflows
>> Winner of BOSSIE, CODIE, OW2 and Gartner awards
>> http://p.sf.net/sfu/Bonitasoft
>> _______________________________________________
>> mesa-users mailing list
>> mesa-users at lists.sourceforge.net
>> https://lists.sourceforge.net/lists/listinfo/mesa-users
-------------- next part --------------
! inlist_rotation
&star_job
!! BASIC MODEL PARAMETERS -- I/O
!------------------------------------------------------------!
!mesa_data_dir='../../mesa_4298/data'
show_log_description_at_start = .false.
create_pre_main_sequence_model = .true.
show_log_description_at_start = .false.
save_model_number = 284
save_model_filename = 'ZAMS.mod'
!load_saved_model = .true.
!saved_model_name = 'pre_bounce.mod'
saved_model_name = 'ZAMS.mod'
!! NUCLEAR REACTION NETWORK
!------------------------------------------------------------!
change_initial_net = .true.
new_net_name = 'approx21.net'
!new_net_name = 'pp_extras.net'
!new_net_name = 'approx21.net'
!new_net_name = 'extended.net'
!new_net_name = 'basic.net'
!new_net_name = 'cno_extras.net'
!auto_extend_net = .true.
eos_file_prefix = 'mesa'
kappa_file_prefix = 'gs98'
!! STELLAR ROTATION PARAMETERS
!------------------------------------------------------------!
change_v_flag = .true.
new_v_flag = .true.
!new_v_flag = .false.
!change_rotation_flag = .true.
change_rotation_flag = .false.
!new_rotation_flag = .true.
new_rotation_flag = .false.
!set_initial_omega_div_omega_crit = .true. ! only applies at start of run, not for restarts.
!new_omega_div_omega_crit = 0.3 ! omega as fraction of current omega_crit at surface
!set_initial_surface_rotation_v = .true.
!set_surface_rotation_v = .true.
!new_surface_rotation_v = 200 ! km/sec 200 for massive
!! CONVERGENCE PARAMETERS BASED ON PHOTOSPHERE OPACITY
!------------------------------------------------------------!
!set_tau_factor= .true.
!set_to_this_tau_factor= 1.d01! 7.5d04
!relax_tau_factor = .true.
!relax_to_this_tau_factor = 1.d01
!! SPECIFIC REACTION RATE MODIFICATION FOR CONVERGENCE
!------------------------------------------------------------!
!set_rate_c12ag = 'Buchmann'
!num_special_rate_factors = 1
!reaction_for_special_factor(:) = 'r_c12_ag_o16'
!special_rate_factor(1) = 1.1
!! CREATION OF HIGHER MASS MODELS COMMAND
!------------------------------------------------------------!
!relax_mass_scale=.true.
!dlgm_per_step=1d-2
!new_mass= 140.
pgstar_flag = .true.
/ ! end of star_job namelist
&controls
!! BASIC STELLAR PROPERTIES - MASS AND METALLICITY
!------------------------------------------------------------!
initial_mass = 13
initial_z = 0.02d0
!! MAGNETIC PROPERTIES
!------------------------------------------------------------!
!set to 0 for non-magnetic
D_ST_factor = 0
am_nu_ST_factor = 0
am_D_mix_factor = 0.0333333333333333d0
am_gradmu_factor = 0.05d0 ! 0.05d0
!! MASS LOSS PARAMETERS
!------------------------------------------------------------!
!RGB_wind_scheme = 'Vink'
!AGB_wind_scheme = 'Vink'
!RGB_to_AGB_wind_switch = 1d-4
!vink_wind_eta = 0.8
RGB_wind_scheme = 'Dutch'
AGB_wind_scheme = 'Dutch'
RGB_to_AGB_wind_switch = 1d-4
Dutch_wind_eta = 0.8
!RGB_wind_scheme = 'de Jager'
!AGB_wind_scheme = 'de Jager'
!RGB_to_AGB_wind_switch = 1d-4
!de_Jager_wind_eta = 0.8 ! 13.5
!de_Jager_wind_eta = 1.5 ! 12.7
!de_Jager_wind_eta = 1.7 ! 10.9
!de_Jager_wind_eta = 1.6 ! 12.5
! tune eta and smooth to match 15M mass and J(1.5) at he depletion
!de_Jager_wind_eta = 1.65
!RGB_wind_scheme = 'Reimers'
!AGB_wind_scheme = 'Blocker'
!RGB_to_AGB_wind_switch = 1d-4
!Reimers_wind_eta = 0.5d0
!Blocker_wind_eta = 0.08d0
!! GRID REFINEMENT FOR CONVERGENCE
!------------------------------------------------------------!
!mesh_delta_coeff = 1.3 ! for massive stars advanced burning
mesh_delta_coeff = 0.4
max_allowed_nz = 100000
max_iter_for_resid_tol1 = 4
max_iter_for_resid_tol2 = 7
max_iter_for_resid_tol3 = 10
!edit during Si burning.
!mesh_delta_coeff = 4
!max_iter_for_resid_tol1 = 3
!max_iter_for_resid_tol2 = -1
!max_iter_for_resid_tol3 = -1
log_center_density_limit = 12
log_center_temp_limit = 10.3d0 ! max T about 2d10
!xa_central_lower_limit_species(1) = 'si28' ! isotope name as defined in chem_def
!xa_central_lower_limit(1) = 1.d-04
!! OUTPUT DATA PARAMETERS
!------------------------------------------------------------!
photostep = 1
profile_interval = 1
history_interval = 1
terminal_cnt = 2
write_header_frequency = 5
max_num_profile_models = 1000
max_model_number = 63000
!! MIXING AND CONVECTION PARAMETERS
!------------------------------------------------------------!
use_Ledoux_criterion = .false.
!use_Ledoux_criterion = .true.
!nsmooth_gradL_composition_term = 1 ! comment out for Schw
!use_gradmu_alt_for_mlt = .true.
!use_gradmu_X_for_mlt = .false.
!use_gradmu_alt_for_mlt = .false.
!use_gradmu_X_for_mlt = .true.
!gradL_comp_term_logRho_switch = 3
!alpha_semiconvection = 0.1! 4d-2 ! comment out for Schw
!thermo_haline_coeff = 2 ! comment out for Schw
dX_mix_dist_limit = 1d-4
omega_function_weight = 20
mdot_omega_power = 0.43
mixing_length_alpha = 1.0 ! 1.0 for M > 30
clip_D_limit = 10 ! zero mixing diffusion coeffs that are smaller than this
min_overshoot_q = 1d-3
T_mix_limit = 1d4
!set_min_D_mix = .true.
!min_D_mix = 1d8 ! D_mix will be at least this large
!min_center_Ye_for_min_D_mix = 0.87 !0.47 ! min_D_mix is only used when center Ye >= this
!i.e., when center Ye drops below this, min_D_mix is 0
!!!!New stuff below on T limits for mixing at late stages of massive stellar evol.
! high center T limit to avoid negative mass fractions
!sig_min_factor_for_high_Tcenter = 0.01
! inactive when >= 1d0
! if Tcenter >= Tcenter_min_for_sig_min_factor_full_on,
! then okay to reduce sig by as much as this factor
! as needed to prevent causing negative abundances
Tcenter_min_for_sig_min_factor_full_on = 2d9
! if Tcenter >= this, factor = sig_min_factor_for_neg_abundances
! this should be > Tcenter_max_for_sig_min_factor_full_off.
Tcenter_max_for_sig_min_factor_full_off = 1d9
! if Tcenter <= this, factor = 1, so has no effect
! this should be < Tcenter_min_for_sig_min_factor_full_on.
! for T > full_off and < full_on, factor changes linearly with Tcenter
use_Type2_opacities = .true.
Zbase = 0.02d0 ! must set this in the main inlist
delta_Ye_highT_limit = 1d-3
okay_to_reduce_gradT_excess = .true.
gradT_excess_f1 = 1d-4
gradT_excess_f2 = 1d-2
gradT_excess_age_fraction = 0.9d0
gradT_excess_lambda1 = 1.0
gradT_excess_beta1 = 0.35
gradT_excess_lambda2 = 0.5
gradT_excess_beta2 = 0.25
gradT_excess_dlambda = 0.1
gradT_excess_dbeta = 0.1
!! MESH PROPERTIES AND OTHER MISCELLANEOUS PARAMETERS
!------------------------------------------------------------!
!rate_factor_for_fe56ec_to_cr56 = 1d-4
varcontrol_target = 1d-3
xa_scale = 1d-5
newton_itermin = 2
tol_correction_norm_alt = 1d-3
tol_max_correction_alt = 1d0
mass_change_full_on_dt = 1d7 ! (seconds)
mass_change_full_off_dt = 1d6 ! (seconds)
min_timestep_limit = 1d-15 ! (seconds)
/ ! end of controls namelist
!! PGSTAR CONTROLS START HERE
&pgstar
! 1. Profile plots
TRho_Profile_win_flag = .true.
!Summary_Profile_win_flag = .true.
Abundance_win_flag = .true.
Power_win_flag = .false.
!Mixing_Ds_win_flag = .true.
!Dynamo_win_flag = .true.
!Mode_Prop_win_flag = .true
! 1A. Options for Profile plots
Abundance_win_width = 7.5
Abundance_win_aspect_ratio = 0.9 ! aspect_ratio = height/width
TRho_Profile_win_width = 7.5
TRho_Profile_win_aspect_ratio = 0.9 ! aspect_ratio = height/width
TRho_Profile_xleft = 0.15
TRho_Profile_xright = 0.85
TRho_Profile_ybot = 0.15
TRho_Profile_ytop = 0.85
TRho_Profile_txt_scale = 1.0
TRho_Profile_title = 'TRho_Profile'
TRho_switch_to_Column_Depth = .false.
! if true, replace logRho for xaxis by log column depth (g/cm^2)
show_TRho_Profile_legend = .false.
TRho_Profile_legend_coord = 0.05
TRho_Profile_legend_fjust = 0.0
TRho_Profile_legend_disp1 = -2.0
TRho_Profile_legend_del_disp = -1.3
show_TRho_Profile_text_info = .false.
TRho_Profile_text_info_xfac = 0.77 ! controls x location
TRho_Profile_text_info_dxfac = 0.02 ! controls x spacing to value from text
TRho_Profile_text_info_yfac = 0.6 ! controls y location of 1st line
TRho_Profile_text_info_dyfac = -0.04 ! controls line spacing
show_TRho_Profile_mass_locs = .true.
show_TRho_Profile_kap_regions = .false.
show_TRho_Profile_eos_regions = .true.
show_TRho_Profile_degeneracy_line = .true.
show_TRho_Profile_Pgas_Prad_line = .true.
show_TRho_Profile_burn_lines = .true.
show_TRho_Profile_burn_labels = .true.
show_TRho_Profile_logQ_limit = .false.
show_TRho_Profile_annotation1 = .false.
show_TRho_Profile_annotation2 = .false.
show_TRho_Profile_annotation3 = .false.
! axis limits
TRho_Profile_xmin = 0.0
TRho_Profile_xmax = 9.0
TRho_Profile_ymin = 7.0
TRho_Profile_ymax = 10.5
! file output
TRho_Profile_file_flag = .false.
TRho_Profile_file_dir = 'TRho_png'
TRho_Profile_file_prefix = 'trho_profile_'
TRho_Profile_file_cnt = 5 ! output when mod(model_number,TRho_Profile_file_cnt)==0
TRho_Profile_file_width = -1 ! negative means use same value as for window
TRho_Profile_file_aspect_ratio = -1 ! negative means use same value as for window
Power_win_width = 6
Power_win_aspect_ratio = 0.75 ! aspect_ratio = height/width
Power_xleft = 0.15
Power_xright = 0.85
Power_ybot = 0.15
Power_ytop = 0.85
Power_txt_scale = 1.0
Power_title = 'Power'
Power_xaxis_name = 'mass'
Power_xaxis_reversed = .false.
Power_legend_max_cnt = 16
Power_legend_txt_scale_factor = 0.7 ! relative to other text
! power xaxis limits -- to override system default selections
Power_xmin = -101d0 ! only used if /= -101d0
Power_xmax = -101d0 ! only used if /= -101d0
! power yaxis limits -- to override system default selections
Power_ymin = -101d0 ! only used if /= -101d0
Power_ymax = -101d0 ! only used if /= -101d0
! file output
Power_file_flag = .true.
Power_file_dir = 'power_png'
Power_file_prefix = 'power_'
Power_file_cnt = 1 ! output when mod(model_number,Power_file_cnt)==0
Power_file_width = -1 ! (inches) negative means use same value as for window
Power_file_aspect_ratio = -1 ! negative means use same value as for window
! 2. History plots
Kipp_win_flag = .true.
Kipp_file_flag = .false.
Kipp_file_dir = 'Kipp_png'
Kipp_file_prefix = 'Kipp_profile_'
Kipp_file_cnt = 5
Kipp_file_width = -1
Kipp_file_aspect_ratio = -1
!TRho_win_flag = .true.
!HR_win_flag = .true.
/ ! end of pgstar namelist
More information about the Mesa-users
mailing list