-
Notifications
You must be signed in to change notification settings - Fork 145
Expand file tree
/
Copy pathcase.py
More file actions
107 lines (106 loc) · 3.9 KB
/
case.py
File metadata and controls
107 lines (106 loc) · 3.9 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
#!/usr/bin/env python3
import json
# Configuring case dictionary
print(
json.dumps(
{
# Logistics
"run_time_info": "T",
# Computational Domain Parameters
"x_domain%beg": 0.0e00,
"x_domain%end": 4.0e-03 / 1.0e-03,
"y_domain%beg": 0.0e00,
"y_domain%end": 4.0e-03 / 1.0e-03,
"z_domain%beg": 0.0e00,
"z_domain%end": 4.0e-03 / 1.0e-03,
"stretch_x": "T",
"a_x": 4.0e00,
"x_stretch%beg": -1.5e-03 / 1.0e-03,
"x_stretch%end": 1.5e-03 / 1.0e-03,
"stretch_y": "T",
"a_y": 4.0e00,
"y_stretch%beg": -1.5e-03 / 1.0e-03,
"y_stretch%end": 1.5e-03 / 1.0e-03,
"stretch_z": "T",
"a_z": 4.0e00,
"z_stretch%beg": -1.5e-03 / 1.0e-03,
"z_stretch%end": 1.5e-03 / 1.0e-03,
"cyl_coord": "F",
"m": 200,
"n": 200,
"p": 200,
"dt": 0.2e-09 / 1.0e-03,
"t_step_start": 0,
"t_step_stop": 30,
"t_step_save": 30,
# Simulation Algorithm Parameters
"num_patches": 2,
"model_eqns": 2,
"alt_soundspeed": "F",
"num_fluids": 2,
"mpp_lim": "T",
"mixture_err": "T",
"time_stepper": 3,
"weno_order": 5,
"weno_eps": 1.0e-16,
"weno_Re_flux": "F",
"weno_avg": "F",
"avg_state": 2,
"mapped_weno": "T",
"null_weights": "F",
"mp_weno": "F",
"riemann_solver": 2,
"wave_speeds": 1,
"bc_x%beg": -2,
"bc_x%end": -6,
"bc_y%beg": -2,
"bc_y%end": -6,
"bc_z%beg": -2,
"bc_z%end": -6,
# Formatted Database Files Structure Parameters
"format": 1,
"precision": 2,
"prim_vars_wrt": "T",
"parallel_io": "T",
# Patch 1: High pressured water
"patch_icpp(1)%geometry": 9,
"patch_icpp(1)%x_centroid": 80.0e-03 / 1.0e-03,
"patch_icpp(1)%y_centroid": 80.0e-03 / 1.0e-03,
"patch_icpp(1)%z_centroid": 80.0e-03 / 1.0e-03,
"patch_icpp(1)%length_x": 160.0e-03 / 1.0e-03,
"patch_icpp(1)%length_y": 160.0e-03 / 1.0e-03,
"patch_icpp(1)%length_z": 160.0e-03 / 1.0e-03,
"patch_icpp(1)%vel(1)": 0.0e00,
"patch_icpp(1)%vel(2)": 0.0e00,
"patch_icpp(1)%vel(3)": 0.0e00,
"patch_icpp(1)%pres": 1.0e05,
"patch_icpp(1)%alpha_rho(1)": 1000.0e00,
"patch_icpp(1)%alpha_rho(2)": 0.1e00,
"patch_icpp(1)%alpha(1)": 0.9e00,
"patch_icpp(1)%alpha(2)": 0.1e00,
# Patch 3: Air bubble
"patch_icpp(2)%geometry": 8,
"patch_icpp(2)%smoothen": "T",
"patch_icpp(2)%smooth_patch_id": 1,
"patch_icpp(2)%smooth_coeff": 0.5e00,
"patch_icpp(2)%x_centroid": 0.0e00,
"patch_icpp(2)%y_centroid": 0.0e00,
"patch_icpp(2)%z_centroid": 0.0e00,
"patch_icpp(2)%radius": 1.0e-03 / 1.0e-03,
"patch_icpp(2)%alter_patch(1)": "T",
"patch_icpp(2)%vel(1)": 0.0e00,
"patch_icpp(2)%vel(2)": 0.0e00,
"patch_icpp(2)%vel(3)": 0.0e00,
"patch_icpp(2)%pres": 1.0e03,
"patch_icpp(2)%alpha_rho(1)": 100.0e00,
"patch_icpp(2)%alpha_rho(2)": 0.9e00,
"patch_icpp(2)%alpha(1)": 0.1e00,
"patch_icpp(2)%alpha(2)": 0.9e00,
# 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(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00),
"fluid_pp(2)%pi_inf": 0.0e00,
}
)
)