HLL[C] & Gamma Methods Testing

Author: henryleberre

examples/1D_inert_shocktube/case.py

@@ -28,7 +28,7 @@
 L    = 0.10
 Nx   = 400
 dx   = L / Nx
-dt   = 5e-9
+dt   = 20e-8
 Tend = 40e-6
 
 NT         = int(Tend / dt)
@@ -66,9 +66,9 @@
     'weno_avg'                     : 'F',
     'mapped_weno'                  : 'T',
     'mp_weno'                      : 'T',
-    'riemann_solver'               : 1,
-    'wave_speeds'                  : 1,
-    'avg_state'                    : 2,
+    'riemann_solver'               : 2,
+    'wave_speeds'                  : 2,
+    'avg_state'                    : 1,
     'bc_x%beg'                     :-2,
     'bc_x%end'                     :-3,
     # ==========================================================================
@@ -106,7 +106,7 @@
     # ==========================================================================
 
     # Fluids Physical Parameters ===============================================
-    'fluid_pp(1)%gamma'            : 1.0E+00/(4.4E+00-1.0E+00),
+    'fluid_pp(1)%gamma'            : 1.0E+00/(1.55E+00-1.0E+00),
     'fluid_pp(1)%pi_inf'           : 0,
     # ==========================================================================
 

examples/1D_reactive_shocktube/case.py

@@ -30,9 +30,9 @@
 u_r = -487.34
 
 L  = 0.12
-Nx = 800 * args.scale
+Nx = 400 * args.scale
 dx = L/Nx
-dt = dx/abs(u_r)*0.01
+dt = dx/abs(u_r)*0.02
 Tend=230e-6
 
 NT=int(Tend/dt)
@@ -55,7 +55,7 @@
     't_step_stop'                  : NT,
     't_step_save'                  : NS,
     't_step_print'                 : NS,
-    'parallel_io'                  : 'T' if args.mfc.get("mpi", True) else 'F',
+    'parallel_io'                  : 'F' if args.mfc.get("mpi", True) else 'F',
 
     # Simulation Algorithm Parameters ==========================================
     'model_eqns'                   : 2,
@@ -69,7 +69,7 @@
     'weno_avg'                     : 'F',
     'mapped_weno'                  : 'T',
     'mp_weno'                      : 'T',
-    'riemann_solver'               : 1,
+    'riemann_solver'               : 2,
     'wave_speeds'                  : 1,
     'avg_state'                    : 2,
     'bc_x%beg'                     :-2,