Hyperelastic RMT code with updated modular precision

.github/workflows/lint-source.yml

@@ -42,10 +42,15 @@ jobs:
       run: |
         ! grep -iR -e '\.\.\.' -e '\-\-\-' -e '===' ./src/*
 
+    - name: Looking for all caps comments in source
+      run: |
+        ! grep -RE '^\s+[!]\s+[A-Z]{5}' ./src/*
+
     - name: Looking for junk comments in examples
       run: |
         ! grep -R '# ===' ./benchmarks **/*.py
         ! grep -R '# ===' ./examples/**/*.py
         ! grep -R '==='   ./benchmarks/**/*.py
         ! grep -R '==='   ./examples/**/*.py
 
+

docs/documentation/case.md

@@ -508,6 +508,7 @@ To restart the simulation from $k$-th time step, see [Restarting Cases](running.
 | `vel_wrt(i)`         | Logical | Add the $i$-direction velocity to the database	  |
 | `E_wrt`              | Logical | Add the total energy to the database	 |
 | `pres_wrt`           | Logical | Add the pressure to the database	|
+| `tau_wrt`            | Logical | Add the elastic stresses to the database |
 | `alpha_wrt(i)`       | Logical | Add the volume fraction of fluid $i$ to the database	|
 | `gamma_wrt`          | Logical | Add the specific heat ratio function to the database	|
 | `heat_ratio_wrt`     | Logical | Add the specific heat ratio to the database	|

examples/1D_hyper_impact_strong/case.py

@@ -0,0 +1,92 @@
+#!/usr/bin/env python3
+import math
+import json
+
+# Numerical setup
+Nx = 201
+dx = 1.0 / (1.0 * (Nx + 1))
+
+Tend = 41e-06
+Nt = 4000
+mydt = Tend / (1.0 * Nt)
+
+# Configuring case dictionary
+print(
+    json.dumps(
+        {
+            # Logistics
+            "run_time_info": "T",
+            # Computational Domain Parameters
+            "x_domain%beg": 0.0e00,
+            "x_domain%end": 1.0e00,
+            #                    'y_domain%beg'                 : 0.E+00,
+            #                    'y_domain%end'                 : 0.002,
+            "m": Nx,
+            "n": 0,
+            "p": 0,
+            "dt": mydt,
+            "t_step_start": 0,
+            "t_step_stop": int(Nt),
+            "t_step_save": int(Nt / 200),
+            # Simulation Algorithm Parameters
+            "num_patches": 2,
+            "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,
+            "weno_Re_flux": "F",
+            "weno_avg": "F",
+            "mapped_weno": "F",
+            "null_weights": "F",
+            "mp_weno": "T",
+            "riemann_solver": 2,
+            "wave_speeds": 1,
+            "avg_state": 2,
+            "bc_x%beg": -3,
+            "bc_x%end": -3,
+            #'bc_y%beg'                     : -3,
+            #'bc_y%end'                     : -3,
+            # Turning on Hypoelasticity
+            "hyperelasticity": "T",
+            "hyper_model": 1,
+            "fd_order": 4,
+            # Formatted Database Files Structure Parameters
+            "format": 1,
+            "precision": 2,
+            "prim_vars_wrt": "T",
+            "parallel_io": "F",
+            # Patch 1 L
+            "patch_icpp(1)%geometry": 1,
+            "patch_icpp(1)%x_centroid": 0.25,
+            #                   'patch_icpp(1)%y_centroid'     : 0.001,
+            "patch_icpp(1)%length_x": 0.5,
+            #                   'patch_icpp(1)%length_y'       : 0.002,
+            "patch_icpp(1)%vel(1)": 1000,
+            #                   'patch_icpp(1)%vel(2)'         : 100*0,
+            "patch_icpp(1)%pres": 1.0e5,
+            "patch_icpp(1)%alpha_rho(1)": 1000,
+            "patch_icpp(1)%alpha(1)": 1.0,
+            "patch_icpp(1)%tau_e(1)": 0.0,
+            # Patch 2 R
+            "patch_icpp(2)%geometry": 1,
+            "patch_icpp(2)%x_centroid": 0.75,
+            #                    'patch_icpp(2)%y_centroid'     : 0.001,
+            "patch_icpp(2)%length_x": 0.5,
+            #                    'patch_icpp(2)%length_y'       : 0.002,
+            "patch_icpp(2)%vel(1)": 1000,
+            #                    'patch_icpp(2)%vel(2)'         : -100*0,
+            "patch_icpp(2)%pres": 1.0e05,
+            "patch_icpp(2)%alpha_rho(1)": 1000,
+            "patch_icpp(2)%alpha(1)": 1.0,
+            "patch_icpp(2)%tau_e(1)": 0.0,
+            # Fluids Physical Parameters
+            "fluid_pp(1)%gamma": 1.0e00 / (4.4e00 - 1.0e00),
+            "fluid_pp(1)%pi_inf": 4.4e00 * 6.0e08 / (4.4e00 - 1.0e00),
+            "fluid_pp(1)%G": 1e010,
+        }
+    )
+)

examples/1D_hyper_impact_weak/case.py

@@ -0,0 +1,82 @@
+#!/usr/bin/env python3
+import math
+import json
+
+# Numerical setup
+Nx = 201
+dx = 1.0 / (1.0 * (Nx + 1))
+
+Tend = 64e-06
+Nt = 4000
+mydt = Tend / (1.0e00 * Nt)
+# print(mydt)
+# Configuring case dictionary
+print(
+    json.dumps(
+        {
+            # Logistics
+            "run_time_info": "T",
+            # Computational Domain Parameters
+            "x_domain%beg": 0.0e00,
+            "x_domain%end": 1.0e00,
+            "m": Nx,
+            "n": 0,
+            "p": 0,
+            "dt": mydt,
+            "t_step_start": 0,
+            "t_step_stop": int(Nt),
+            "t_step_save": int(Nt / 200),
+            # Simulation Algorithm Parameters
+            "num_patches": 2,
+            "model_eqns": 3,
+            "alt_soundspeed": "F",
+            "num_fluids": 1,
+            "mpp_lim": "F",
+            "mixture_err": "F",
+            "time_stepper": 3,
+            "weno_order": 5,
+            "weno_eps": 1.0e-16,
+            "weno_Re_flux": "F",
+            "weno_avg": "F",
+            "mapped_weno": "F",
+            "null_weights": "F",
+            "mp_weno": "T",
+            "riemann_solver": 2,
+            "wave_speeds": 1,
+            "avg_state": 2,
+            "bc_x%beg": -3,
+            "bc_x%end": -3,
+            # Turning on Hyperelasticity
+            "hyperelasticity": "T",
+            "hyper_model": 1,
+            "fd_order": 4,
+            # Formatted Database Files Structure Parameters
+            "format": 1,
+            "precision": 2,
+            "prim_vars_wrt": "T",
+            "parallel_io": "F",
+            # Patch 1 L
+            "patch_icpp(1)%geometry": 1,
+            "patch_icpp(1)%x_centroid": 0.25,
+            "patch_icpp(1)%length_x": 0.5,
+            "patch_icpp(1)%vel(1)": 10,
+            "patch_icpp(1)%pres": 1.0e5,
+            "patch_icpp(1)%alpha_rho(1)": 1000,
+            "patch_icpp(1)%alpha(1)": 1.0,
+            "patch_icpp(1)%tau_e(1)": 0.0,
+            # Patch 2 R
+            "patch_icpp(2)%geometry": 1,
+            "patch_icpp(2)%x_centroid": 0.75,
+            "patch_icpp(2)%length_x": 0.5,
+            "patch_icpp(2)%vel(1)": -10,  # 10,
+            "patch_icpp(2)%pres": 1.0e05,
+            "patch_icpp(2)%alpha_rho(1)": 1000,
+            "patch_icpp(2)%alpha(1)": 1.0,
+            "patch_icpp(2)%tau_e(1)": 0.0,
+            # Fluids Physical Parameters
+            "fluid_pp(1)%gamma": 1.0e00 / (4.4e00 - 1.0e00),
+            "fluid_pp(1)%pi_inf": 4.4e00 * 6.0e08 / (4.4e00 - 1.0e00),
+            "fluid_pp(1)%G": 1.0e00,  # .E+010,
+        }
+    )
+)

examples/1D_hyper_impact_weak/original.py

@@ -0,0 +1,90 @@
+#!/usr/bin/env python3
+import math
+import json
+
+# Numerical setup
+Nx = 201
+dx = 1.0 / (1.0 * (Nx + 1))
+
+Tend = 64e-06
+Nt = 4000
+mydt = Tend / (1.0e00 * Nt)
+# print(mydt)
+# Configuring case dictionary
+print(
+    json.dumps(
+        {
+            # Logistics
+            "run_time_info": "T",
+            # Computational Domain Parameters
+            "x_domain%beg": 0.0e00,
+            "x_domain%end": 1.0e00,
+            #                    'y_domain%beg'                 : 0.E+00,
+            #                    'y_domain%end'                 : 0.002,
+            "m": Nx,
+            "n": 0,
+            "p": 0,
+            "dt": mydt,
+            "t_step_start": 0,
+            "t_step_stop": int(Nt),
+            "t_step_save": int(Nt / 200),
+            # Simulation Algorithm Parameters
+            "num_patches": 2,
+            "model_eqns": 3,
+            "alt_soundspeed": "F",
+            "num_fluids": 1,
+            "mpp_lim": "F",
+            "mixture_err": "F",
+            "time_stepper": 3,
+            "weno_order": 5,
+            "weno_eps": 1.0e-16,
+            "weno_Re_flux": "F",
+            "weno_avg": "F",
+            "mapped_weno": "F",
+            "null_weights": "F",
+            "mp_weno": "T",
+            "riemann_solver": 2,
+            "wave_speeds": 1,
+            "avg_state": 2,
+            "bc_x%beg": -3,
+            "bc_x%end": -3,
+            #'bc_y%beg'                     : -3,
+            #'bc_y%end'                     : -3,
+            # Turning on Hyperelasticity
+            "hyperelasticity": "T",
+            # Formatted Database Files Structure Parameters
+            "format": 1,
+            "precision": 2,
+            "prim_vars_wrt": "T",
+            "parallel_io": "F",
+            # Patch 1 L
+            "patch_icpp(1)%geometry": 1,
+            "patch_icpp(1)%x_centroid": 0.25,
+            #                   'patch_icpp(1)%y_centroid'     : 0.001,
+            "patch_icpp(1)%length_x": 0.5,
+            #                   'patch_icpp(1)%length_y'       : 0.002,
+            "patch_icpp(1)%vel(1)": 10,
+            #                   'patch_icpp(1)%vel(2)'         : 100*0,
+            "patch_icpp(1)%pres": 1.0e5,
+            "patch_icpp(1)%alpha_rho(1)": 1000,
+            "patch_icpp(1)%alpha(1)": 1.0,
+            "patch_icpp(1)%tau_e(1)": 0.0,
+            # Patch 2 R
+            "patch_icpp(2)%geometry": 1,
+            "patch_icpp(2)%x_centroid": 0.75,
+            #                    'patch_icpp(2)%y_centroid'     : 0.001,
+            "patch_icpp(2)%length_x": 0.5,
+            #                    'patch_icpp(2)%length_y'       : 0.002,
+            "patch_icpp(2)%vel(1)": 10,  # -10,
+            #                    'patch_icpp(2)%vel(2)'         : -100*0,
+            "patch_icpp(2)%pres": 1.0e05,
+            "patch_icpp(2)%alpha_rho(1)": 1000,
+            "patch_icpp(2)%alpha(1)": 1.0,
+            "patch_icpp(2)%tau_e(1)": 0.0,
+            # Fluids Physical Parameters
+            "fluid_pp(1)%gamma": 1.0e00 / (4.4e00 - 1.0e00),
+            "fluid_pp(1)%pi_inf": 4.4e00 * 6.0e08 / (4.4e00 - 1.0e00),
+            "fluid_pp(1)%G": 1.0e00,  # .E+010,
+        }
+    )
+)

examples/1D_hypo_impact_strong/case.py

@@ -0,0 +1,91 @@
+#!/usr/bin/env python3
+import math
+import json
+
+# Numerical setup
+Nx = 201
+dx = 1.0 / (1.0 * (Nx + 1))
+
+Tend = 41e-06
+Nt = 4000
+mydt = Tend / (1.0 * Nt)
+
+# Configuring case dictionary
+print(
+    json.dumps(
+        {
+            # Logistics
+            "run_time_info": "T",
+            # Computational Domain Parameters
+            "x_domain%beg": 0.0e00,
+            "x_domain%end": 1.0e00,
+            #                    'y_domain%beg'                 : 0.E+00,
+            #                    'y_domain%end'                 : 0.002,
+            "m": Nx,
+            "n": 0,
+            "p": 0,
+            "dt": mydt,
+            "t_step_start": 0,
+            "t_step_stop": int(Nt),
+            "t_step_save": int(Nt / 200),
+            # Simulation Algorithm Parameters
+            "num_patches": 2,
+            "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,
+            "weno_Re_flux": "F",
+            "weno_avg": "F",
+            "mapped_weno": "F",
+            "null_weights": "F",
+            "mp_weno": "T",
+            "riemann_solver": 2,
+            "wave_speeds": 1,
+            "avg_state": 2,
+            "bc_x%beg": -3,
+            "bc_x%end": -3,
+            #'bc_y%beg'                     : -3,
+            #'bc_y%end'                     : -3,
+            # Turning on Hypoelasticity
+            "hypoelasticity": "T",
+            "fd_order": 4,
+            # Formatted Database Files Structure Parameters
+            "format": 1,
+            "precision": 2,
+            "prim_vars_wrt": "T",
+            "parallel_io": "F",
+            # Patch 1 L
+            "patch_icpp(1)%geometry": 1,
+            "patch_icpp(1)%x_centroid": 0.25,
+            #                   'patch_icpp(1)%y_centroid'     : 0.001,
+            "patch_icpp(1)%length_x": 0.5,
+            #                   'patch_icpp(1)%length_y'       : 0.002,
+            "patch_icpp(1)%vel(1)": 1000,
+            #                   'patch_icpp(1)%vel(2)'         : 100*0,
+            "patch_icpp(1)%pres": 1.0e5,
+            "patch_icpp(1)%alpha_rho(1)": 1000,
+            "patch_icpp(1)%alpha(1)": 1.0,
+            "patch_icpp(1)%tau_e(1)": 0.0,
+            # Patch 2 R
+            "patch_icpp(2)%geometry": 1,
+            "patch_icpp(2)%x_centroid": 0.75,
+            #                    'patch_icpp(2)%y_centroid'     : 0.001,
+            "patch_icpp(2)%length_x": 0.5,
+            #                    'patch_icpp(2)%length_y'       : 0.002,
+            "patch_icpp(2)%vel(1)": 1000,
+            #                    'patch_icpp(2)%vel(2)'         : -100*0,
+            "patch_icpp(2)%pres": 1.0e05,
+            "patch_icpp(2)%alpha_rho(1)": 1000,
+            "patch_icpp(2)%alpha(1)": 1.0,
+            "patch_icpp(2)%tau_e(1)": 0.0,
+            # Fluids Physical Parameters
+            "fluid_pp(1)%gamma": 1.0e00 / (4.4e00 - 1.0e00),
+            "fluid_pp(1)%pi_inf": 4.4e00 * 6.0e08 / (4.4e00 - 1.0e00),
+            "fluid_pp(1)%G": 1e010,
+        }
+    )
+)