IGR User Guide test cases

docs/documentation/case.md

@@ -637,6 +637,7 @@ To restart the simulation from $k$-th time step, see @ref running "Restarting Ca
 | `cf_wrt`                | Logical | Write color function field |
 | `chem_wrt_T`            | Logical | Write temperature field for chemistry output |
 | `fft_wrt`               | Logical | Enable FFT output |
+| `double_mach`           | Logical | Enable Double Mach Reflection case |
 | `sim_data`              | Logical | Write interface and energy data files (post_process) |
 | `integral_wrt`          | Logical | Write integral data |
 | `num_integrals`         | Integer | Number of integral regions |

examples/2D_IGR_double_mach/README.md

@@ -0,0 +1,3 @@
+# Double Mach reflection case referenced in IGR User Guide paper 
+
+Mach number set to 10 and shock angle at 60 degrees with igr_order = 3 and alf_factor = 2

examples/2D_IGR_double_mach/case.py

@@ -0,0 +1,83 @@
+import json
+import math
+
+# Numerical setup
+Nx = 799
+Ny = int((Nx + 1) / 4) - 1
+dx = 4.0 / (1.0 * (Nx + 1))
+
+c = 1.0
+cfl = 0.1
+mydt = cfl * dx / (11.0 * c)
+Nt = int(0.2 / (1 * mydt))
+
+# Configuring case dictionary
+print(
+    json.dumps(
+        {
+            # Logistics
+            "run_time_info": "T",
+            # Computational Domain Parameters
+            "x_domain%beg": 0,
+            "x_domain%end": 4,
+            "y_domain%beg": 0,
+            "y_domain%end": 1,
+            "m": Nx,
+            "n": Ny,
+            "p": 0,
+            "dt": mydt,
+            "t_step_start": 0,
+            "t_step_stop": Nt,
+            "t_step_save": int(Nt / 100),
+            # Simulation Algorithm Parameters
+            "num_patches": 1,
+            "model_eqns": 2,
+            "alt_soundspeed": "F",
+            "num_fluids": 1,
+            "mpp_lim": "F",
+            "mixture_err": "F",
+            "time_stepper": 1,
+            "weno_order": 3,
+            "weno_eps": 1.0e-16,
+            "mapped_weno": "T",
+            "null_weights": "F",
+            "mp_weno": "F",
+            "riemann_solver": 2,
+            "wave_speeds": 1,
+            "avg_state": 2,
+            "bc_x%beg": -17,
+            "bc_x%end": -3,
+            "bc_y%beg": -2,
+            "bc_y%end": -3,
+            "igr": "T",
+            "igr_order": 3,
+            "igr_iter_solver": 1,
+            "num_igr_iters": 5,
+            "num_igr_warm_start_iters": 150,
+            "alf_factor": 10,
+            "double_mach": "T",
+            # Formatted Database Files Structure Parameters
+            "format": 1,
+            "precision": 2,
+            "prim_vars_wrt": "T",
+            "parallel_io": "T",
+            "omega_wrt(3)": "F",
+            "fd_order": 2,
+            # Patch 1
+            "patch_icpp(1)%geometry": 3,
+            "patch_icpp(1)%x_centroid": 2.0,
+            "patch_icpp(1)%y_centroid": 0.5,
+            "patch_icpp(1)%length_x": 4.0,
+            "patch_icpp(1)%length_y": 1.0,
+            "patch_icpp(1)%vel(1)": 0.0,
+            "patch_icpp(1)%vel(2)": 0.0,
+            "patch_icpp(1)%pres": 1.0,
+            "patch_icpp(1)%alpha_rho(1)": 1.4,
+            "patch_icpp(1)%hcid": 285,
+            "patch_icpp(1)%alpha(1)": 1.0,
+            # Fluids Physical Parameters
+            "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00),
+            "fluid_pp(1)%pi_inf": 0.0,
+        }
+    )
+)

examples/2D_IGR_isentropicvortex/README.md

@@ -0,0 +1 @@
+# Isentropic Vortex with modified initial condition and IGR enabled  used in IGR User Guide paper 

examples/2D_IGR_isentropicvortex/case.py

@@ -0,0 +1,97 @@
+import json
+import math
+
+# Parameters
+epsilon = "5d0"
+alpha = "1d0"
+gamma = "1.4"
+
+
+# Numerical setup
+Nx = 149
+dx = 20.0 / (1.0 * (Nx + 1))
+
+c = 1.4**0.5
+C = 0.3
+mydt = C * dx / c
+Nt = int(20 / (0.1 * mydt))
+
+# Configuring case dictionary
+print(
+    json.dumps(
+        {
+            # Logistics
+            "run_time_info": "T",
+            # Computational Domain Parameters
+            "x_domain%beg": -10,
+            "x_domain%end": 10,
+            "y_domain%beg": -10,
+            "y_domain%end": 10,
+            "stretch_x": "F",
+            "stretch_y": "F",
+            "loops_x": 2,
+            "loops_y": 2,
+            "a_x": 1.03,
+            "a_y": 1.03,
+            "x_a": -1.5,
+            "y_a": -1.5,
+            "x_b": 1.5,
+            "y_b": 1.5,
+            "m": Nx,
+            "n": Nx,
+            "p": 0,
+            "dt": mydt,
+            "t_step_start": 0,
+            "t_step_stop": int(30 * Nt),
+            "t_step_save": int(Nt),
+            # Simulation Algorithm Parameters
+            "num_patches": 1,
+            "model_eqns": 2,
+            "alt_soundspeed": "F",
+            "num_fluids": 1,
+            "mpp_lim": "F",
+            "mixture_err": "F",
+            "time_stepper": 3,
+            "weno_order": 5,
+            "weno_eps": 1.0e-16,
+            "mapped_weno": "F",
+            "null_weights": "F",
+            "mp_weno": "F",
+            "riemann_solver": 5,
+            "wave_speeds": 1,
+            "avg_state": 2,
+            "bc_x%beg": -1,
+            "bc_x%end": -1,
+            "bc_y%beg": -1,
+            "bc_y%end": -1,
+            "igr": "T",
+            "igr_order": 5,
+            "igr_iter_solver": 1,
+            "num_igr_iters": 5,
+            "num_igr_warm_start_iters": 150,
+            "alf_factor": 2,
+            # Formatted Database Files Structure Parameters
+            "format": 1,
+            "precision": 2,
+            "prim_vars_wrt": "T",
+            "parallel_io": "T",
+            "omega_wrt(3)": "T",
+            "fd_order": 2,
+            # Patch 1
+            "patch_icpp(1)%geometry": 3,
+            "patch_icpp(1)%x_centroid": 0,
+            "patch_icpp(1)%y_centroid": 0,
+            "patch_icpp(1)%length_x": 20.0,
+            "patch_icpp(1)%length_y": 20.0,
+            "patch_icpp(1)%vel(1)": 0.0,
+            "patch_icpp(1)%vel(2)": 0.0,
+            "patch_icpp(1)%pres": 0.0,
+            "patch_icpp(1)%alpha_rho(1)": 0.0,
+            "patch_icpp(1)%hcid": 283,
+            "patch_icpp(1)%alpha(1)": 1.0,
+            # Fluids Physical Parameters
+            "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00),
+            "fluid_pp(1)%pi_inf": 0.0,
+        }
+    )
+)

examples/2D_IGR_riemann_test/README.md

@@ -0,0 +1 @@
+# 2D Riemann test  with IGR and entropy wave fluctuations referenced in IGR User Guide paper 

examples/2D_IGR_riemann_test/case.py

@@ -0,0 +1,69 @@
+#!/usr/bin/env python3
+import json
+import math
+
+# Configuring case dictionary
+print(
+    json.dumps(
+        {
+            # Logistics
+            "run_time_info": "T",
+            # Computational Domain Parameters
+            "x_domain%beg": 0.0,
+            "x_domain%end": 1.0,
+            "y_domain%beg": 0.0,
+            "y_domain%end": 1.0,
+            "m": 499,
+            "n": 499,
+            "p": 0,
+            "dt": 8e-05,
+            "t_step_start": 0,
+            "t_step_stop": 10000,
+            "t_step_save": 100,
+            # Simulation Algorithm Parameters
+            "num_patches": 1,
+            "model_eqns": 2,
+            "alt_soundspeed": "F",
+            "num_fluids": 1,
+            "mpp_lim": "F",
+            "mixture_err": "T",
+            "time_stepper": 3,
+            "mp_weno": "F",
+            "weno_order": 5,
+            "weno_eps": 1e-16,
+            "riemann_solver": 2,
+            "wave_speeds": 1,
+            "avg_state": 2,
+            "igr": "T",
+            "igr_order": 5,
+            "igr_iter_solver": 1,
+            "num_igr_iters": 5,
+            "num_igr_warm_start_iters": 150,
+            "alf_factor": 2,
+            "bc_x%beg": -3,
+            "bc_x%end": -3,
+            "bc_y%beg": -3,
+            "bc_y%end": -3,
+            # Formatted Database Files Structure Parameters
+            "format": 1,
+            "precision": 2,
+            "prim_vars_wrt": "T",
+            "parallel_io": "T",
+            # Patch 1: Base
+            "patch_icpp(1)%geometry": 3,
+            "patch_icpp(1)%x_centroid": 0.5,
+            "patch_icpp(1)%y_centroid": 0.5,
+            "patch_icpp(1)%length_x": 1.0,
+            "patch_icpp(1)%length_y": 1.0,
+            "patch_icpp(1)%vel(1)": 0.0,
+            "patch_icpp(1)%vel(2)": 0.0,
+            "patch_icpp(1)%pres": 1.0,
+            "patch_icpp(1)%alpha_rho(1)": 1.0,
+            "patch_icpp(1)%alpha(1)": 1.0,
+            "patch_icpp(1)%hcid": 284,
+            # Fluids Physical Parameters
+            "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00),
+            "fluid_pp(1)%pi_inf": 0.0e00,
+        }
+    )
+)

examples/3D_IGR_33jet/README.md

@@ -0,0 +1 @@
+# 33-jet Mach 10 case, in a configuration inspired by SpaceX Super Heavy, used in ACM Gordon Bell prize submission and referenced in IGR User Guide

examples/3D_IGR_TaylorGreenVortex/README.md

@@ -0,0 +1,5 @@
+# Taylor Green Vortex used in IGR User Guide paper 
+
+Ma = 1.25 and Re = 10,000 
+alf_factor is changed to 2
+fft_wrt is enabled to get turbulent energy cascade  

examples/3D_IGR_TaylorGreenVortex/case.py

@@ -2,20 +2,20 @@
 import json
 import math
 
-N = 99
+N = 127
 
-Re = 1600
+Re = 10000
 L = 1
 P0 = 101325
 rho0 = 1
 C0 = math.sqrt(1.4 * P0)
-V0 = 0.1 * C0
+V0 = 1.25 * C0
 mu = V0 * L / Re
 
-cfl = 0.5
+cfl = 0.3
 dx = 2 * math.pi * L / (N + 1)
 
-dt = cfl * dx / (C0)
+dt = cfl * dx / (C0 + V0)
 
 tC = L / V0
 tEnd = 20 * tC
@@ -60,7 +60,7 @@
             "igr_iter_solver": 1,
             "num_igr_iters": 3,
             "num_igr_warm_start_iters": 3,
-            "alf_factor": 10,
+            "alf_factor": 2,
             "viscous": "T",
             # Formatted Database Files Structure Parameters
             "format": 1,
@@ -70,6 +70,7 @@
             "omega_wrt(2)": "T",
             "omega_wrt(3)": "T",
             "qm_wrt": "T",
+            "fft_wrt": "T",
             "fd_order": 4,
             "parallel_io": "T",
             # Patch 1: Background (AIR - 2)