1.4M_ms_op_mono

The test checks the functionality of OP mono opacities. The test vehicle is a 1.4 \({\rm M}_\odot\) solar metallicity model.

This test case has three parts. Click to see a larger view of a plot.

  • Part 1 (inlist_1.4M_ms_initial_model) creates the pre-main-sequence model and evolves the model until the central hydrogen mass fractions falls below 0.5. During this evolution the OP mono opacities are off.

../_images/initial_model.svg

The upper plot shows the fraction of the opacity that is from OP mono (yellow curve) and the opacity (blue curve) as a function of the temperature profile - center on right, surface on left. The lower plot shows the opacity derivatives.

  • Part 2 (inlist_1.4M_ms_op_mono) continues the evolution of inlist_1.4M_ms_initial_model for 10 timesteps using the OP mono opacities, as shown by the fraction of the opacity that is from OP mono (yellow curve, upper plot)

../_images/op_mono.svg

  • Part 3 (inlist_1.4M_ms_op_mono) continues the evolution of inlist_1.4M_ms_initial_model for 10 timesteps using an alternative formulation of the OP mono opacities, as shown by the smoother opacity derivatives in the lower plot. See $MESA_DIR/kap/public/kap_lib.f.90 for the differences between the regular and alternative OP mono formulations.

../_images/op_mono_alt.svg

pgstar commands used for the plots above:

&pgstar

  pgstar_interval = 10

 ! device

  file_white_on_black_flag = .true.
  file_device = 'vcps'
  file_extension = 'ps'

 ! two profile panels
   Profile_Panels1_win_flag = .true.

   Profile_Panels1_win_width = 15
   Profile_Panels1_win_aspect_ratio = 0.75

   Profile_Panels1_xleft = 0.15
   Profile_Panels1_xright = 0.85
   Profile_Panels1_ybot = 0.15
   Profile_Panels1_ytop = 0.85
   Profile_Panels1_txt_scale = 1.0
   Profile_Panels1_title = '1.4_ms_op_mono'

   Profile_Panels1_num_panels = 2

   Profile_Panels1_xaxis_name = 'logT'
   Profile_Panels1_xaxis_reversed = .false.
   Profile_Panels1_xmin = -101d0
   Profile_Panels1_xmax = -101d0

   Profile_Panels1_yaxis_name(1) = 'kap_frac_op_mono'
   Profile_Panels1_other_yaxis_name(1) = 'log_opacity'

   Profile_Panels1_yaxis_name(2) = 'dlnkap_dlnd'
   Profile_Panels1_other_yaxis_name(2) = 'dlnkap_dlnT'

   Profile_Panels1_file_flag = .true.
   Profile_Panels1_file_dir = 'pgstar_out'
   Profile_Panels1_file_prefix = 'profile_panels1_'
   Profile_Panels1_file_interval = 10
   Profile_Panels1_file_width = 15
   Profile_Panels1_file_aspect_ratio = -1

/ ! end of pgstar namelist

Last-Updated: 27May2021 (MESA ebecc10) by fxt