Skip to content

Commit 095bedf

Browse files
Close #258 (#260)
Co-authored-by: Henry LE BERRE <[email protected]>
1 parent 37137c3 commit 095bedf

File tree

1 file changed

+110
-0
lines changed

1 file changed

+110
-0
lines changed

examples/2D_riemann_test/case.py

Lines changed: 110 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,110 @@
1+
#!/usr/bin/env python3
2+
3+
import json
4+
import math
5+
6+
# Configuring case dictionary
7+
print(json.dumps({
8+
# Logistics ================================================
9+
'run_time_info' : 'T',
10+
# ==========================================================
11+
12+
# Computational Domain Parameters ==========================
13+
'x_domain%beg' : 0.E+00,
14+
'x_domain%end' : 1.E+00,
15+
'y_domain%beg' : 0.E+00,
16+
'y_domain%end' : 1.E+00,
17+
'm' : 499,
18+
'n' : 499,
19+
'p' : 0,
20+
'dt' : 8e-05,
21+
't_step_start' : 0,
22+
't_step_stop' : 10000,
23+
't_step_save' : 100,
24+
# ==========================================================
25+
26+
# Simulation Algorithm Parameters ==========================
27+
'num_patches' : 4,
28+
'model_eqns' : 2,
29+
'alt_soundspeed' : 'F',
30+
'num_fluids' : 1,
31+
'adv_alphan' : 'T',
32+
'mpp_lim' : 'F',
33+
'mixture_err' : 'T',
34+
'time_stepper' : 3,
35+
'mp_weno' : 'F',
36+
'recon_type' : 1,
37+
'weno_order' : 5,
38+
'weno_eps' : 1e-16,
39+
#'muscl_order' : 2,
40+
#'muscl_lim' : 1,
41+
'riemann_solver' : 2,
42+
'wave_speeds' : 1,
43+
'avg_state' : 2,
44+
'bc_x%beg' : -3,
45+
'bc_x%end' : -3,
46+
'bc_y%beg' : -3,
47+
'bc_y%end' : -3,
48+
# ==========================================================
49+
50+
# Formatted Database Files Structure Parameters ============
51+
'format' : 1,
52+
'precision' : 2,
53+
'prim_vars_wrt' :'T',
54+
'parallel_io' :'T',
55+
# ==========================================================
56+
57+
# Patch 1: Base ============================================
58+
'patch_icpp(1)%geometry' : 3,
59+
'patch_icpp(1)%x_centroid' : 0.4,
60+
'patch_icpp(1)%y_centroid' : 0.4,
61+
'patch_icpp(1)%length_x' : 0.8,
62+
'patch_icpp(1)%length_y' : 0.8,
63+
'patch_icpp(1)%vel(1)' : 4/math.sqrt(11),
64+
'patch_icpp(1)%vel(2)' : 4/math.sqrt(11),
65+
'patch_icpp(1)%pres' : 9/310,
66+
'patch_icpp(1)%alpha_rho(1)' : 77/558,
67+
'patch_icpp(1)%alpha(1)' : 1,
68+
# ==========================================================
69+
# Patch 1: Base ============================================
70+
'patch_icpp(2)%geometry' : 3,
71+
'patch_icpp(2)%x_centroid' : 0.4,
72+
'patch_icpp(2)%y_centroid' : 0.9,
73+
'patch_icpp(2)%length_x' : 0.8,
74+
'patch_icpp(2)%length_y' : 0.2,
75+
'patch_icpp(2)%vel(1)' : 4/math.sqrt(11),
76+
'patch_icpp(2)%vel(2)' : 0,
77+
'patch_icpp(2)%pres' : 0.3,
78+
'patch_icpp(2)%alpha_rho(1)' : 33/62,
79+
'patch_icpp(2)%alpha(1)' : 1,
80+
# ==========================================================
81+
# Patch 1: Base ============================================
82+
'patch_icpp(3)%geometry' : 3,
83+
'patch_icpp(3)%x_centroid' : 0.9E+00,
84+
'patch_icpp(3)%y_centroid' : 0.4E+00,
85+
'patch_icpp(3)%length_x' : 0.2,
86+
'patch_icpp(3)%length_y' : 0.8,
87+
'patch_icpp(3)%vel(1)' : 0,
88+
'patch_icpp(3)%vel(2)' : 4/math.sqrt(11),
89+
'patch_icpp(3)%pres' : 0.3,
90+
'patch_icpp(3)%alpha_rho(1)' : 33/62,
91+
'patch_icpp(3)%alpha(1)' : 1,
92+
# ==========================================================
93+
# Patch 1: Base ============================================
94+
'patch_icpp(4)%geometry' : 3,
95+
'patch_icpp(4)%x_centroid' : 0.9E+00,
96+
'patch_icpp(4)%y_centroid' : 0.9E+00,
97+
'patch_icpp(4)%length_x' : 0.2,
98+
'patch_icpp(4)%length_y' : 0.2,
99+
'patch_icpp(4)%vel(1)' : 0,
100+
'patch_icpp(4)%vel(2)' : 0,
101+
'patch_icpp(4)%pres' : 1.5,
102+
'patch_icpp(4)%alpha_rho(1)' : 1.5,
103+
'patch_icpp(4)%alpha(1)' : 1.,
104+
# ==========================================================
105+
# Fluids Physical Parameters ===============================
106+
'fluid_pp(1)%gamma' : 1.E+00/(1.4E+00-1.E+00),
107+
'fluid_pp(1)%pi_inf' : 0.E+00,
108+
# ==========================================================
109+
}))
110+
# ==============================================================================

0 commit comments

Comments
 (0)