debugging hyper schemes

.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_1material/case.py

@@ -0,0 +1,107 @@
+#!/usr/bin/env python3
+
+# 1D -> Hyperelasticity -> 1 Fluid(s)
+
+import math, json
+
+print(json.dumps({
+    'run_time_info': 'T',
+    'm': 299,
+    'n': 0,
+    'p': 0,
+    'dt': 1e-07,
+    't_step_start': 0,
+    't_step_stop': 50,
+    't_step_save': 25,
+    'num_patches': 3,
+    'model_eqns': 3,
+    'alt_soundspeed': 'F',
+    'num_fluids': 1,
+    'mpp_lim': 'F',
+    'mixture_err': 'F',
+    'time_stepper': 3,
+    'weno_order': 5,
+    'weno_eps': 1e-16,
+    'mapped_weno': 'F',
+    'null_weights': 'F',
+    'mp_weno': 'F',
+    'riemann_solver': 2,
+    'wave_speeds': 1,
+    'avg_state': 2,
+    'format': 1,
+    'precision': 2,
+    'prim_vars_wrt': 'F',
+    'parallel_io': 'F',
+    'patch_icpp(1)%pres': 1000000.0,
+    'patch_icpp(1)%alpha_rho(1)': 1000.0,
+    'patch_icpp(1)%alpha(1)': 1.0,
+    'patch_icpp(2)%pres': 100000.0,
+    'patch_icpp(2)%alpha_rho(1)': 1000.0,
+    'patch_icpp(2)%alpha(1)': 1.0,
+    'patch_icpp(3)%pres': 500000.0,
+    'patch_icpp(3)%alpha_rho(1)': 1000.0,
+    'patch_icpp(3)%alpha(1)': 1.0,
+    'fluid_pp(1)%gamma': 0.3,
+    'fluid_pp(1)%pi_inf': 780000.0,
+    'Ca': 0.9769178386380458,
+    'Web': 13.927835051546392,
+    'Re_inv': 0.009954269975623245,
+    'pref': 101325.0,
+    'rhoref': 1000.0,
+    'bubble_model': 3,
+    'polytropic': 'T',
+    'polydisperse': 'F',
+    'thermal': 3,
+    'R0ref': 1e-05,
+    'patch_icpp(1)%r0': 1,
+    'patch_icpp(1)%v0': 0,
+    'patch_icpp(2)%r0': 1,
+    'patch_icpp(2)%v0': 0,
+    'patch_icpp(3)%r0': 1,
+    'patch_icpp(3)%v0': 0,
+    'qbmm': 'F',
+    'dist_type': 2,
+    'poly_sigma': 0.3,
+    'R0_type': 1,
+    'sigR': 0.1,
+    'sigV': 0.1,
+    'rhoRV': 0.0,
+    'x_domain%beg': 0.0,
+    'x_domain%end': 1.0,
+    'bc_x%beg': -3,
+    'bc_x%end': -3,
+    'patch_icpp(1)%geometry': 1,
+    'patch_icpp(1)%x_centroid': 0.05,
+    'patch_icpp(1)%length_x': 0.1,
+    'patch_icpp(2)%x_centroid': 0.45,
+    'patch_icpp(2)%length_x': 0.7,
+    'patch_icpp(3)%x_centroid': 0.9,
+    'patch_icpp(3)%length_x': 0.2,
+    'patch_icpp(1)%vel(1)': 0.0,
+    'patch_icpp(2)%geometry': 1,
+    'patch_icpp(2)%vel(1)': 0.0,
+    'patch_icpp(3)%geometry': 1,
+    'patch_icpp(3)%vel(1)': 0.0,
+    'hyperelasticity': 'T',
+    'hyper_model': 1,
+    'fd_order': 4,
+    'patch_icpp(1)%tau_e(1)': 0.0,
+    'patch_icpp(2)%tau_e(1)': 0.0,
+    'patch_icpp(3)%tau_e(1)': 0.0,
+    'fluid_pp(1)%G': 1000000000.0,
+    'parallel_io'  : 'T', 
+    'cons_vars_wrt'   : 'T',
+    'prim_vars_wrt': 'T', 
+    'alpha_rho_wrt(1)': 'T',
+    'rho_wrt'      : 'T', 
+    'mom_wrt(1)'      : 'T',
+    'vel_wrt(1)'   : 'T', 
+    'E_wrt'           : 'T',
+    'pres_wrt'     : 'T',
+    'alpha_wrt(1)'    : 'T',
+    'gamma_wrt'    : 'T',
+    'heat_ratio_wrt'  : 'T',
+    'pi_inf_wrt'   : 'T', 
+    'pres_inf_wrt'    : 'T',
+    'c_wrt'        : 'T',
+}))

examples/1D_hyper_2materials/case.py

@@ -0,0 +1,110 @@
+#!/usr/bin/env python3
+
+# 1D -> Hyperelasticity -> 2 Fluid(s)
+
+import math, json
+import argparse
+
+print(json.dumps({
+    'run_time_info': 'T',
+    'm': 299,
+    'n': 0,
+    'p': 0,
+    'dt': 1e-07,
+    't_step_start': 0,
+    't_step_stop': 50,
+    't_step_save': 25,
+    'num_patches': 3,
+    'model_eqns': 3,
+    'alt_soundspeed': 'F',
+    'num_fluids': 2,
+    'mpp_lim': 'F',
+    'mixture_err': 'F',
+    'time_stepper': 3,
+    'weno_order': 5,
+    'weno_eps': 1e-16,
+    'mapped_weno': 'F',
+    'null_weights': 'F',
+    'mp_weno': 'F',
+    'riemann_solver': 2,
+    'wave_speeds': 1,
+    'avg_state': 2,
+    'format': 1,
+    'precision': 2,
+    'prim_vars_wrt': 'F',
+    'parallel_io': 'F',
+    'patch_icpp(1)%pres': 1000000.0,
+    'patch_icpp(1)%alpha_rho(1)': 900.0,
+    'patch_icpp(1)%alpha(1)': 0.9,
+    'patch_icpp(2)%pres': 100000.0,
+    'patch_icpp(2)%alpha_rho(1)': 100,
+    'patch_icpp(2)%alpha(1)': 0.1,
+    'patch_icpp(3)%pres': 500000.0,
+    'patch_icpp(3)%alpha_rho(1)': 900,
+    'patch_icpp(3)%alpha(1)': 0.9,
+    'fluid_pp(1)%gamma': 0.3,
+    'fluid_pp(1)%pi_inf': 780000.0,
+#    'Ca': 0.9769178386380458,
+#    'Web': 13.927835051546392,
+#    'Re_inv': 0.009954269975623245,
+    'pref': 101325.0,
+    'rhoref': 1000.0,
+#    'bubble_model': 3,
+#    'polytropic': 'T',
+#    'polydisperse': 'F',
+#    'thermal': 3,
+    'R0ref': 1e-05,
+    'patch_icpp(1)%r0': 1,
+    'patch_icpp(1)%v0': 0,
+    'patch_icpp(2)%r0': 1,
+    'patch_icpp(2)%v0': 0,
+    'patch_icpp(3)%r0': 1,
+    'patch_icpp(3)%v0': 0,
+#    'qbmm': 'F',
+#    'dist_type': 2,
+#    'poly_sigma': 0.3,
+#    'R0_type': 1,
+#    'sigR': 0.1,
+#    'sigV': 0.1,
+#    'rhoRV': 0.0,
+    'x_domain%beg': 0.0,
+    'x_domain%end': 1.0,
+    'bc_x%beg': -3,
+    'bc_x%end': -3,
+    'patch_icpp(1)%geometry': 1,
+    'patch_icpp(1)%x_centroid': 0.05,
+    'patch_icpp(1)%length_x': 0.1,
+    'patch_icpp(2)%x_centroid': 0.45,
+    'patch_icpp(2)%length_x': 0.7,
+    'patch_icpp(3)%x_centroid': 0.9,
+    'patch_icpp(3)%length_x': 0.2,
+    'patch_icpp(1)%vel(1)': 0.0,
+    'patch_icpp(2)%geometry': 1,
+    'patch_icpp(2)%vel(1)': 0.0,
+    'patch_icpp(3)%geometry': 1,
+    'patch_icpp(3)%vel(1)': 0.0,
+    'hyperelasticity': 'T',
+    'hyper_model': 1,
+    'fd_order': 4,
+    'patch_icpp(1)%tau_e(1)': 0.0,
+    'patch_icpp(2)%tau_e(1)': 0.0,
+    'patch_icpp(3)%tau_e(1)': 0.0,
+    'fluid_pp(1)%G': 1000000000.0,
+    'fluid_pp(2)%gamma': 0.3,
+    'fluid_pp(2)%pi_inf': 780000.0,
+    'patch_icpp(1)%alpha_rho(2)': 100,
+    'patch_icpp(1)%alpha(2)': 0.1,
+    'patch_icpp(2)%alpha_rho(2)': 900,
+    'patch_icpp(2)%alpha(2)': 0.9,
+    'patch_icpp(3)%alpha_rho(2)': 100,
+    'patch_icpp(3)%alpha(2)': 0.1,
+    'fluid_pp(2)%G': 0.0, #50000.0
+    'parallel_io'  : 'T', 'cons_vars_wrt'   : 'T',
+    'prim_vars_wrt': 'T', 'alpha_rho_wrt(1)': 'T',
+    'rho_wrt'      : 'T', 'mom_wrt(1)'      : 'T',
+    'vel_wrt(1)'   : 'T', 'E_wrt'           : 'T',
+    'pres_wrt'     : 'T', 'alpha_wrt(1)'    : 'T',
+    'gamma_wrt'    : 'T', 'heat_ratio_wrt'  : 'T',
+    'pi_inf_wrt'   : 'T', 'pres_inf_wrt'    : 'T',
+    'c_wrt'        : 'T',
+ }))   

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,
+        }
+    )
+)