[Mesa-users] Science Fair

[Karan Jain]

Hello MESA users,

I am having a bit of trouble with my science fair project. I am not too
sure exactly what I should look for. I registered for the physics division
of my science fair competition as there was no astrophysics or astronomy
category. I am using MESA to evolve a star until the end of AGB. I have all
the data, which are stars with a variable initial metallicity and initial
mass of 2 solar masses, and I have successfully converted it all to NumPy
arrays in python. I do not know where to go from here as my project is
based around metalicity but I need to relate it to a physics concept. I was
considering using machine learning and k-means clustering to see if I could
find any patterns but I realized that would not help me too much. I was
wondering if you have any advice as to how I should use the data that I
have to relate the effects of a variable metallicity to a concept in
physics and how I could use python to display that relationship. Attached
is an inlist_defaults file that I am using for one model.

Thank You
Karan Jain
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20210201/43728d96/attachment.htm>
-------------- next part --------------
! MESA Revision 11701  
! Part of MESA-Web (c) Fields, C. E. & Timmes, F. X. 2015-2019 (http://mesa-web.asu.edu)
! Cite/offer co-authorship when appropriate.
! Questions/Comments: Carl Fields at: f i e l d s c 9 (at) msu (dot) edu
!-----------------------------------------------------------------------------------

&star_job

      create_pre_main_sequence_model = .true.
      pre_ms_relax_num_steps = 50
      set_initial_age = .true.
      initial_age = 0 ! in years
      set_initial_model_number = .true.
      initial_model_number = 0
      change_initial_net = .true.      
      new_net_name = 'basic.net'
      dump_missing_metals_into_heaviest = .false.
      new_omega_div_omega_crit = 0.0
      set_near_zams_omega_div_omega_crit_steps = 50
      pgstar_flag = .true.
      rates_cache_dir = './cache'
      rate_tables_dir = './rate_tables'
      rate_cache_suffix = 'user'
      show_log_description_at_start = .false.

/ !end of star_job namelist

&controls

  ! starting specifications

      profile_interval = 50
      history_interval = 1
      terminal_interval = 100000000
      write_header_frequency = 100000000
      photo_interval = 50000000
      max_num_profile_models = 500
      warn_when_large_rel_run_E_err = 1.0d99 

      initial_mass = 2 ! in Msun units
      initial_z = 0.03  

      ! custom stopping condition
      fe_core_infall_limit = 1d12
      when_to_stop_rtol = 1d-3
      when_to_stop_atol = 1d-3

      use_gold_tolerances = .false.
      use_eosELM = .false.

      ! a factor of ten smaller, although run_star_extras will 
      ! quit the run if dt < min_timestep = 1d-6 (the default)
      min_timestep_limit = 1d-7

      ! winds
      cool_wind_RGB_scheme = 'Reimers'
      cool_wind_AGB_scheme = 'Blocker'
      Reimers_scaling_factor = 0.1
      Blocker_scaling_factor = 0.1

      ! resolution - spatial (mass) - global
      mesh_delta_coeff = 2.0

      ! resolution - temporal (time) - global
      varcontrol_target = 1d-2
      dX_nuc_drop_min_X_limit = 1d-2
        
      ! resolution - temporal (time) - hydrogen
      delta_lg_XH_cntr_min = -6
      delta_lg_XH_cntr_max = -3
      delta_lg_XH_cntr_limit = 0.05d0
      delta_lg_XH_cntr_hard_limit = -1


      ! mixing options
      mixing_length_alpha = 2.0
      MLT_option = 'Cox'
      use_Ledoux_criterion = .true.
      
      ! overshoot
      ! nonburning regions
      overshoot_f_above_nonburn_core = 0.0001
      overshoot_f0_above_nonburn_core = 0.0005
      overshoot_f_above_nonburn_shell = 0.0001
      overshoot_f0_above_nonburn_shell = 0.0005
      overshoot_f_below_nonburn_shell = 0.0001
      overshoot_f0_below_nonburn_shell = 0.0005

      ! h burning regions
      overshoot_f_above_burn_h_core = 0.0001
      overshoot_f0_above_burn_h_core = 0.0005
      overshoot_f_above_burn_h_shell = 0.0001
      overshoot_f0_above_burn_h_shell = 0.0005
      overshoot_f_below_burn_h_shell = 0.0001
      overshoot_f0_below_burn_h_shell = 0.0005

      ! he burning regions
      overshoot_f_above_burn_he_core = 0.0001
      overshoot_f0_above_burn_he_core = 0.0005
      overshoot_f_above_burn_he_shell = 0.0001
      overshoot_f0_above_burn_he_shell = 0.0005
      overshoot_f_below_burn_he_shell = 0.0001
      overshoot_f0_below_burn_he_shell = 0.0005

      ! metal burning regions
      overshoot_f_above_burn_z_core = 0.0001
      overshoot_f0_above_burn_z_core = 0.0005
      overshoot_f_above_burn_z_shell = 0.0001
      overshoot_f0_above_burn_z_shell = 0.0005
      overshoot_f_below_burn_z_shell = 0.0001
      overshoot_f0_below_burn_z_shell = 0.0005

      ! semiconvection
      alpha_semiconvection = 0.01
      semiconvection_option = 'Langer_85 mixing; gradT = gradr'

      ! thermohaline
      thermohaline_coeff = 1.0
      thermohaline_option = 'Kippenhahn'

      ! mlt ++ 
      okay_to_reduce_gradT_excess = .true.

      ! convective pre-mixing
      do_conv_premix = .false
      conv_premix_avoid_increase = .true.
      conv_premix_fix_pgas = .true.

      am_D_mix_factor = 0.033333333d0
      D_DSI_factor = 1
      D_SH_factor  = 1
      D_SSI_factor = 1
      D_ES_factor  = 1
      D_GSF_factor = 1
      D_ST_factor  = 1

/ ! end of controls namelist