Boundary Condition Paches

examples/2D_jet/README.md

@@ -0,0 +1,4 @@
+# Gas Jet (2D)
+
+## Final Condition
+<img src='final_condition.png' height='MAX_HEIGHT'/>

examples/2D_jet/case.py

@@ -0,0 +1,116 @@
+#!/usr/bin/env python3
+import math
+import json
+
+pA = 101325
+rhoA = 1.29
+gam = 1.4
+c_l = math.sqrt(1.4 * pA / rhoA)
+
+pS = 2.458 * pA
+velS = 0.701 * c_l
+rhoS = 1.862 * rhoA
+
+leng = 1e-2
+Ny = 250
+Nx = Ny
+dx = leng / Nx
+
+time_end = 2 * leng / velS
+cfl = 0.8
+
+dt = cfl * dx / c_l
+Nt = int(time_end / dt)
+
+eps = 1e-5
+
+# Configuring case dictionary
+print(
+    json.dumps(
+        {
+            # Logistics
+            "run_time_info": "T",
+            # Computational Domain Parameters
+            "x_domain%beg": -leng / 2.0,
+            "x_domain%end": 3 * leng / 2,
+            "y_domain%beg": -leng,
+            "y_domain%end": leng,
+            "m": int(Nx),
+            "n": int(Ny),
+            "p": 0,
+            "dt": dt,
+            "cfl_adap_dt": "T",
+            "t_stop": time_end,
+            "t_save": time_end / 50,
+            "n_start": 0,
+            "cfl_target": 0.8,
+            # Simulation Algorithm Parameters
+            "num_patches": 2,
+            "model_eqns": 2,
+            "alt_soundspeed": "F",
+            "num_fluids": 2,
+            "mpp_lim": "F",
+            "mixture_err": "T",
+            "time_stepper": 3,
+            "weno_order": 5,
+            "weno_eps": 1.0e-16,
+            "weno_Re_flux": "F",
+            "weno_avg": "F",
+            "mapped_weno": "T",
+            "null_weights": "F",
+            "mp_weno": "F",
+            "riemann_solver": 2,
+            "wave_speeds": 1,
+            "avg_state": 2,
+            "elliptic_smoothing": "T",
+            "elliptic_smoothing_iters": 50,
+            "bc_x%beg": -2,
+            "bc_x%end": -3,
+            "bc_y%beg": -3,
+            "bc_y%end": -3,
+            "num_bc_patches": 1,
+            "patch_bc(1)%dir": 1,
+            "patch_bc(1)%loc": -1,
+            "patch_bc(1)%geometry": 1,
+            "patch_bc(1)%type": -17,
+            "patch_bc(1)%centroid(2)": 0.0,
+            "patch_bc(1)%length(2)": leng / 4,
+            # Formatted Database Files Structure Parameters
+            "format": 1,
+            "precision": 2,
+            "prim_vars_wrt": "T",
+            "parallel_io": "T",
+            # Patch 1: Background
+            "patch_icpp(1)%geometry": 3,
+            "patch_icpp(1)%x_centroid": 0.0,
+            "patch_icpp(1)%y_centroid": 0.0,
+            "patch_icpp(1)%length_x": 10 * leng,
+            "patch_icpp(1)%length_y": 10 * leng,
+            "patch_icpp(1)%vel(1)": 0.0e00,
+            "patch_icpp(1)%vel(2)": 0.0e00,
+            "patch_icpp(1)%pres": pA,
+            "patch_icpp(1)%alpha_rho(1)": rhoA,
+            "patch_icpp(1)%alpha(1)": 1.0 - eps,
+            "patch_icpp(1)%alpha_rho(2)": eps,
+            "patch_icpp(1)%alpha(2)": eps,
+            "patch_icpp(2)%geometry": 3,
+            "patch_icpp(2)%alter_patch(1)": "T",
+            "patch_icpp(2)%x_centroid": -leng / 2,
+            "patch_icpp(2)%y_centroid": 0.0,
+            "patch_icpp(2)%length_x": leng / 4,
+            "patch_icpp(2)%length_y": leng / 3,
+            "patch_icpp(2)%vel(1)": velS,
+            "patch_icpp(2)%vel(2)": 0,
+            "patch_icpp(2)%pres": pS,
+            "patch_icpp(2)%alpha_rho(1)": eps,
+            "patch_icpp(2)%alpha(1)": eps,
+            "patch_icpp(2)%alpha_rho(2)": (1 - eps) * rhoS,
+            "patch_icpp(2)%alpha(2)": 1 - eps,
+            # Fluids Physical Parameters
+            "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00),
+            "fluid_pp(1)%pi_inf": 0.0,
+            "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00),
+            "fluid_pp(2)%pi_inf": 0.0,
+        }
+    )
+)

examples/2D_shockbubble/case.py

@@ -35,8 +35,8 @@
             "p": 0,
             "dt": dt,
             "t_step_start": 0,
-            "t_step_stop": Nt,
-            "t_step_save": int(Nt / 20.0),
+            "t_step_stop": 1000,
+            "t_step_save": 10,
             # Simulation Algorithm Parameters
             "num_patches": 3,
             "model_eqns": 2,