fix HLLC ADC parity, precision, docs, examples

Author: ChrisZYJ

docs/documentation/case.md

@@ -550,7 +550,7 @@ It is recommended to set `weno_eps` to $10^{-6}$ for WENO-JS, and to $10^{-40}$
 
 - `riemann_solver` specifies the choice of the Riemann solver that is used in simulation by an integer from 1 through 4.
 `riemann_solver = 1`, `2`, and `3` correspond to HLL, HLLC, and Exact Riemann solver, respectively (\cite Toro09).
-`riemann_solver = 4` is only for MHD simulations. It resolves 5 of the full seven-wave structure of the MHD equations (\cite Miyoshi05).
+`riemann_solver = 4` is the HLLD solver for MHD or hypoelasticity simulations. For MHD it resolves 5 of the full seven-wave structure of the MHD equations (\cite Miyoshi05).
 
 - `hll_u_interface`: Selects between two HLL discretizations of volume fraction advection (`riemann_solver = 1`):
   - **Default** (``'F'``): \f$\partial_t \alpha_k + u\,\partial_x \alpha_k = 0\f$

examples/2D_axisym_hypo_hlld/case.py

@@ -1,18 +1,15 @@
 #!/usr/bin/env python3
 import json
-import math
 
 # 2D axisymmetric hypoelastic case with HLLD solver.
 # Solid sphere in liquid hit by focused acoustic pulse.
-# Adapted from HLL_new/circle_1_interface for cylindrical coordinates.
-# x = axial, y = radial (axis of symmetry at y = 0).
 
 config = {
     "run_time_info": "T",
     # Computational Domain
-    "x_domain%beg": 0,
+    "x_domain%beg": 0.0,
     "x_domain%end": 1.0,
-    "y_domain%beg": 0,
+    "y_domain%beg": 0.0,
     "y_domain%end": 1.0,
     "cyl_coord": "T",
     "m": 49,
@@ -35,7 +32,6 @@
     "null_weights": "F",
     "mp_weno": "F",
     "riemann_solver": 4,
-    # "hypo_hll_interface_rhs": "T",
     "wave_speeds": 1,
     "avg_state": 2,
     "bc_x%beg": -2,

examples/2D_axisym_hypo_hlld_riemann/case.py

@@ -1,14 +1,12 @@
 #!/usr/bin/env python3
 import json
-import math
 
 config = {
-    # Logistics
     "run_time_info": "T",
-    # Computational Domain Parameters
-    "x_domain%beg": 0,
+    # Computational Domain
+    "x_domain%beg": 0.0,
     "x_domain%end": 1.0,
-    "y_domain%beg": 0,
+    "y_domain%beg": 0.0,
     "y_domain%end": 1.0,
     "m": 49,
     "n": 49,
@@ -17,46 +15,35 @@
     "t_step_start": 0,
     "t_step_stop": 250,
     "t_step_save": 50,
-    # Simulation Algorithm Parameters
+    # Simulation Algorithm
     "num_patches": 2,
     "model_eqns": 2,
     "alt_soundspeed": "F",
     "num_fluids": 2,
-    # "mpp_lim": "F",
     "mpp_lim": "T",
     "mixture_err": "F",
     "time_stepper": 1,
     "weno_order": 1,
     "weno_eps": 1.0e-20,
-    # "mapped_weno": "T",
     "null_weights": "F",
     "mp_weno": "F",
-    # "riemann_solver": 2,
     "riemann_solver": 4,
-    # "riemann_hypo_ADC": 'T',
-    # "ADC_kappa": 1.0,
     "wave_speeds": 1,
     "avg_state": 2,
-    # "bc_x%beg": -6,
-    # "bc_x%end": -6,
-    # "bc_y%beg": -6,
-    # "bc_y%end": -6,
     "bc_x%beg": -3,
     "bc_x%end": -3,
     "bc_y%beg": -3,
     "bc_y%end": -3,
-    # !!! BC=-6 doesn't work with 1st order upwind !!!
-    # Formatted Database Files Structure Parameters
+    # Output
     "format": 1,
-    # "format": 2,
     "precision": 2,
     "prim_vars_wrt": "T",
     "rho_wrt": "T",
     "parallel_io": "T",
     # Hypoelasticity
     "hypoelasticity": "T",
     "fd_order": 4,
-    # Patch 1 Liquid
+    # Patch 1: Liquid background
     "patch_icpp(1)%geometry": 3,
     "patch_icpp(1)%x_centroid": 0.5,
     "patch_icpp(1)%y_centroid": 0.5,
@@ -70,7 +57,7 @@
     "patch_icpp(1)%alpha(1)": 1.0 - 1e-8,
     "patch_icpp(1)%alpha_rho(2)": 1000 * 1e-8,
     "patch_icpp(1)%alpha(2)": 1e-8,
-    # Patch 1 Solid
+    # Patch 2: Solid circle
     "patch_icpp(2)%alter_patch(1)": "T",
     "patch_icpp(2)%geometry": 2,
     "patch_icpp(2)%x_centroid": 0.6,
@@ -84,9 +71,6 @@
     "patch_icpp(2)%alpha(1)": 1e-8,
     "patch_icpp(2)%alpha_rho(2)": 1000 * (1.0 - 1e-8),
     "patch_icpp(2)%alpha(2)": 1.0 - 1e-8,
-    # "patch_icpp(2)%smoothen": "T",
-    # "patch_icpp(2)%smooth_patch_id": 1,
-    # "patch_icpp(2)%smooth_coeff": 4.0e00,
     # Acoustic source
     "acoustic_source": "T",
     "num_source": 1,

examples/2D_hypo_hlld/case.py

@@ -1,16 +1,17 @@
 #!/usr/bin/env python3
 import json
-import math
+
+# 3D hypoelastic case with HLLD solver.
+# Solid sphere in liquid hit by focused acoustic pulse.
 
 config = {
-    # Logistics
     "run_time_info": "T",
-    # Computational Domain Parameters
-    "x_domain%beg": 0,
+    # Computational Domain
+    "x_domain%beg": 0.0,
     "x_domain%end": 1.0,
-    "y_domain%beg": 0,
+    "y_domain%beg": 0.0,
     "y_domain%end": 1.0,
-    "z_domain%beg": 0,
+    "z_domain%beg": 0.0,
     "z_domain%end": 1.0,
     "m": 24,
     "n": 24,