Skip to content

Commit 8383415

Browse files
wilfonbawilfonba
committed
new example
1 parent e4b8976 commit 8383415

File tree

1 file changed

+117
-0
lines changed

1 file changed

+117
-0
lines changed

‎examples/2D_triple_point/case.py‎

Lines changed: 117 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,117 @@
1+
#!/usr/bin/env python3
2+
import math
3+
import json
4+
5+
eps = 1e-8
6+
Nx = 699
7+
Ny = 299
8+
9+
print(json.dumps({
10+
# Logistics ================================================
11+
'run_time_info' : 'F',
12+
# ==========================================================
13+
14+
# Computational Domain Parameters ==========================
15+
'x_domain%beg' : 0,
16+
'x_domain%end' : 7,
17+
'y_domain%beg' : 0,
18+
'y_domain%end' : 3,
19+
'm' : int(Nx),
20+
'n' : int(Ny),
21+
'p' : 0,
22+
'cfl_adap_dt' : 'T',
23+
'cfl_target' : 0.8,
24+
'n_start' : 0,
25+
't_stop' : 4.0,
26+
't_save' : 0.04,
27+
# ==========================================================
28+
29+
# Simulation Algorithm Parameters ==========================
30+
'num_patches' : 3,
31+
'model_eqns' : 2,
32+
'alt_soundspeed' : 'F',
33+
'num_fluids' : 2,
34+
'mpp_lim' : 'T',
35+
'mixture_err' : 'T',
36+
'time_stepper' : 3,
37+
'weno_order' : 7,
38+
'weno_eps' : 1.E-16,
39+
'weno_Re_flux' : 'F',
40+
'weno_avg' : 'F',
41+
'mapped_weno' : 'T',
42+
'null_weights' : 'F',
43+
'mp_weno' : 'F',
44+
'riemann_solver' : 2,
45+
'wave_speeds' : 1,
46+
'avg_state' : 2,
47+
'bc_x%beg' : -3,#11,
48+
'bc_x%end' : -3,#12
49+
'bc_y%beg' : -3,
50+
'bc_y%end' : -3,
51+
# ==========================================================
52+
53+
# Formatted Database Files Structure Parameters ============
54+
'format' : 1,
55+
'precision' : 2,
56+
'prim_vars_wrt' :'T',
57+
'schlieren_wrt' :'T',
58+
'fd_order' : 4,
59+
'schlieren_alpha(1)' : 0.5,
60+
'schlieren_alpha(2)' : 0.5,
61+
'parallel_io' :'T',
62+
# ==========================================================
63+
64+
# Patch 1: Left state ====================================
65+
'patch_icpp(1)%geometry' : 3,
66+
'patch_icpp(1)%x_centroid' : 0.5,
67+
'patch_icpp(1)%y_centroid' : 3,
68+
'patch_icpp(1)%length_x' : 1,
69+
'patch_icpp(1)%length_y' : 6,
70+
'patch_icpp(1)%vel(1)' : 0.,
71+
'patch_icpp(1)%vel(2)' : 0.,
72+
'patch_icpp(1)%pres' : 1.,
73+
'patch_icpp(1)%alpha_rho(1)' : (1 - eps)*1.,
74+
'patch_icpp(1)%alpha_rho(2)' : eps,
75+
'patch_icpp(1)%alpha(1)' : 1 - eps,
76+
'patch_icpp(1)%alpha(2)' : eps,
77+
# ==========================================================
78+
79+
# Patch 2: Top right state ================================
80+
'patch_icpp(2)%geometry' : 3,
81+
'patch_icpp(2)%alter_patch(1)' : 'T',
82+
'patch_icpp(2)%x_centroid' : 4,
83+
'patch_icpp(2)%y_centroid' : 2.25,
84+
'patch_icpp(2)%length_x' : 6,
85+
'patch_icpp(2)%length_y' : 1.5,
86+
'patch_icpp(2)%vel(1)' : 0.,
87+
'patch_icpp(2)%vel(2)' : 0.,
88+
'patch_icpp(2)%pres' : 0.1,
89+
'patch_icpp(2)%alpha_rho(1)' : (1-eps)*0.125,
90+
'patch_icpp(2)%alpha_rho(2)' : eps,
91+
'patch_icpp(2)%alpha(1)' : 1 - eps,
92+
'patch_icpp(2)%alpha(2)' : eps,
93+
# ==========================================================
94+
95+
# Patch 3: Bottom right state =============================
96+
'patch_icpp(3)%geometry' : 3,
97+
'patch_icpp(3)%alter_patch(1)' : 'T',
98+
'patch_icpp(3)%x_centroid' : 4,
99+
'patch_icpp(3)%y_centroid' : 0.75,
100+
'patch_icpp(3)%length_x' : 6,
101+
'patch_icpp(3)%length_y' : 1.5,
102+
'patch_icpp(3)%vel(1)' : 0.,
103+
'patch_icpp(3)%vel(2)' : 0.,
104+
'patch_icpp(3)%pres' : 0.1,
105+
'patch_icpp(3)%alpha_rho(1)' : eps,
106+
'patch_icpp(3)%alpha_rho(2)' : (1 - eps)*1.,
107+
'patch_icpp(3)%alpha(1)' : eps,# 0.95
108+
'patch_icpp(3)%alpha(2)' : 1 - eps,#0.05,
109+
# ==========================================================
110+
111+
# Fluids Physical Parameters ===============================
112+
'fluid_pp(1)%gamma' : 1./(1.5 - 1.),
113+
'fluid_pp(1)%pi_inf' : 0,
114+
'fluid_pp(2)%gamma' : 1./(1.4 - 1.),
115+
'fluid_pp(2)%pi_inf' : 0.,
116+
# ==========================================================
117+
}))

0 commit comments

Comments
 (0)