hse_riemann

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.f90 then 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.