This test case checks Riemann HLLC solver can hold an envelope model in hydrostatic equilibrium.
This test case has 1 parts. Click to see a larger version of a plot.
- Part 1 (
inlist_finish) loads a model containing the envelope from a 12 msun model with its core excised. The model is then evolved for 1e5 s, and the
run_star_extras.f90then checks that the peak mach number is less than 0.1. This is an example of a null test; if the Riemann HLLC solver is working, then the profiles should remain static and the peak mach number should be small:
pgstar commands used for the plots that ake the movie above:
&pgstar file_white_on_black_flag = .true. ! white_on_black flags -- true means white foreground color on black background !file_device = 'png' ! png !file_extension = 'png' file_device = 'vcps' ! postscript file_extension = 'ps' pgstar_interval = 2 Profile_Panels1_win_flag = .true. Profile_Panels1_win_width = 12 Profile_Panels1_win_aspect_ratio = 1.0 Profile_Panels1_xaxis_name = 'logR_cm' Profile_Panels1_xmin = -101d0 Profile_Panels1_xmax = -101d0 Profile_Panels1_num_panels = 3 Profile_Panels1_yaxis_name(1) = 'logRho' Profile_Panels1_ymin(1) = -101d0 Profile_Panels1_ymax(1) = -101d0 Profile_Panels1_ymargin(1) = 0.1 Profile_Panels1_other_yaxis_name(1) = 'entropy' Profile_Panels1_other_ymin(1) = -101d0 Profile_Panels1_other_ymax(1) = -101d0 Profile_Panels1_other_ymargin(1) = 0.1 Profile_Panels1_yaxis_name(2) = 'logT' Profile_Panels1_ymin(2) = -101d0 Profile_Panels1_ymax(2) = -101d0 Profile_Panels1_ymargin(2) = 0.1 Profile_Panels1_yaxis_name(3) = 'v_div_csound' Profile_Panels1_ymin(3) = -0.06 Profile_Panels1_ymax(3) = 0.06 Profile_Panels1_ymargin(3) = 0.1 Profile_Panels1_file_flag = .true. Profile_Panels1_file_dir = 'pgstar_out' Profile_Panels1_file_prefix = 'profile1_' Profile_Panels1_file_interval = 5 ! output when mod(model_number,file_interval)==0 Profile_Panels1_file_width = -1 ! (inches) negative means use same value as for window Profile_Panels1_file_aspect_ratio = -1 ! negative means use same value as for window / ! end of pgstar namelist
Last-Updated: 14Jun2021 (MESA 5be9e57) by fxt.