Skip to content

Commit a3bffbd

Browse files
mrodrig6mrodrig6
committed
added examples to be made tests after verifying
1 parent 5f3adcb commit a3bffbd

File tree

6 files changed

+503
-9
lines changed

6 files changed

+503
-9
lines changed

‎examples/1D_hyper_impact_strong/case.py‎

Lines changed: 102 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,102 @@
1+
#!/usr/bin/env python3
2+
import math
3+
import json
4+
5+
#Numerical setup
6+
Nx = 201
7+
dx = 1./(1.*(Nx+1))
8+
9+
Tend = 41E-06
10+
Nt = 4000
11+
mydt = Tend/(1.*Nt)
12+
13+
# Configuring case dictionary
14+
print(json.dumps({
15+
# Logistics ================================================
16+
'run_time_info' : 'T',
17+
# ==========================================================
18+
19+
# Computational Domain Parameters ==========================
20+
'x_domain%beg' : 0.E+00,
21+
'x_domain%end' : 1.E+00,
22+
# 'y_domain%beg' : 0.E+00,
23+
# 'y_domain%end' : 0.002,
24+
'm' : Nx,
25+
'n' : 0,
26+
'p' : 0,
27+
'dt' : mydt,
28+
't_step_start' : 0,
29+
't_step_stop' : int(Nt),
30+
't_step_save' : int(Nt/200),
31+
# ==========================================================
32+
33+
# Simulation Algorithm Parameters ==========================
34+
'num_patches' : 2,
35+
'model_eqns' : 2,
36+
'alt_soundspeed' : 'F',
37+
'num_fluids' : 1,
38+
'mpp_lim' : 'F',
39+
'mixture_err' : 'F',
40+
'time_stepper' : 3,
41+
'weno_order' : 5,
42+
'weno_eps' : 1.E-16,
43+
'weno_Re_flux' : 'F',
44+
'weno_avg' : 'F',
45+
'mapped_weno' : 'F',
46+
'null_weights' : 'F',
47+
'mp_weno' : 'T',
48+
'riemann_solver' : 2,
49+
'wave_speeds' : 1,
50+
'avg_state' : 2,
51+
'bc_x%beg' : -3,
52+
'bc_x%end' : -3,
53+
#'bc_y%beg' : -3,
54+
#'bc_y%end' : -3,
55+
# ==========================================================
56+
57+
# Turning on Hypoelasticity ================================
58+
'hyperelasticity' : 'T',
59+
# ==========================================================
60+
61+
# Formatted Database Files Structure Parameters ============
62+
'format' : 1,
63+
'precision' : 2,
64+
'prim_vars_wrt' :'T',
65+
'parallel_io' :'F',
66+
# ==========================================================
67+
68+
# Patch 1 L ================================================
69+
'patch_icpp(1)%geometry' : 1,
70+
'patch_icpp(1)%x_centroid' : 0.25,
71+
# 'patch_icpp(1)%y_centroid' : 0.001,
72+
'patch_icpp(1)%length_x' : 0.5,
73+
# 'patch_icpp(1)%length_y' : 0.002,
74+
'patch_icpp(1)%vel(1)' : 1000,
75+
# 'patch_icpp(1)%vel(2)' : 100*0,
76+
'patch_icpp(1)%pres' : 1.E+5,
77+
'patch_icpp(1)%alpha_rho(1)' : 1000,
78+
'patch_icpp(1)%alpha(1)' : 1.,
79+
'patch_icpp(1)%tau_e(1)' : 0.0,
80+
# ==========================================================
81+
82+
# Patch 2 R ================================================
83+
'patch_icpp(2)%geometry' : 1,
84+
'patch_icpp(2)%x_centroid' : 0.75,
85+
# 'patch_icpp(2)%y_centroid' : 0.001,
86+
'patch_icpp(2)%length_x' : 0.5,
87+
# 'patch_icpp(2)%length_y' : 0.002,
88+
'patch_icpp(2)%vel(1)' : 1000,
89+
# 'patch_icpp(2)%vel(2)' : -100*0,
90+
'patch_icpp(2)%pres' : 1.E+05,
91+
'patch_icpp(2)%alpha_rho(1)' : 1000,
92+
'patch_icpp(2)%alpha(1)' : 1.,
93+
'patch_icpp(2)%tau_e(1)' : 0.0,
94+
# ==========================================================
95+
96+
# Fluids Physical Parameters ===============================
97+
'fluid_pp(1)%gamma' : 1.E+00/(4.4E+00-1.E+00),
98+
'fluid_pp(1)%pi_inf' : 4.4E+00*6.E+08/(4.4E+00 - 1.E+00),
99+
'fluid_pp(1)%G' : 1E+010,
100+
# ==========================================================
101+
}))
102+
# ==============================================================================

‎examples/1D_hyper_impact_weak/case.py‎

Lines changed: 102 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,102 @@
1+
#!/usr/bin/env python3
2+
import math
3+
import json
4+
5+
#Numerical setup
6+
Nx = 201
7+
dx = 1./(1.*(Nx+1))
8+
9+
Tend = 64E-06
10+
Nt = 4000
11+
mydt = Tend/(1.E+00*Nt)
12+
#print(mydt)
13+
# Configuring case dictionary
14+
print(json.dumps({
15+
# Logistics ================================================
16+
'run_time_info' : 'T',
17+
# ==========================================================
18+
19+
# Computational Domain Parameters ==========================
20+
'x_domain%beg' : 0.E+00,
21+
'x_domain%end' : 1.E+00,
22+
'y_domain%beg' : 0.E+00,
23+
'y_domain%end' : 0.002,
24+
'm' : Nx,
25+
'n' : 0,
26+
'p' : 0,
27+
'dt' : mydt,
28+
't_step_start' : 0,
29+
't_step_stop' : int(Nt),
30+
't_step_save' : int(Nt/200),
31+
# ==========================================================
32+
33+
# Simulation Algorithm Parameters ==========================
34+
'num_patches' : 2,
35+
'model_eqns' : 3,
36+
'alt_soundspeed' : 'F',
37+
'num_fluids' : 1,
38+
'mpp_lim' : 'F',
39+
'mixture_err' : 'F',
40+
'time_stepper' : 3,
41+
'weno_order' : 5,
42+
'weno_eps' : 1.E-16,
43+
'weno_Re_flux' : 'F',
44+
'weno_avg' : 'F',
45+
'mapped_weno' : 'F',
46+
'null_weights' : 'F',
47+
'mp_weno' : 'T',
48+
'riemann_solver' : 2,
49+
'wave_speeds' : 1,
50+
'avg_state' : 2,
51+
'bc_x%beg' : -3,
52+
'bc_x%end' : -3,
53+
'bc_y%beg' : -3,
54+
'bc_y%end' : -3,
55+
# ==========================================================
56+
57+
# Turning on Hyperelasticity ================================
58+
'hyperelasticity' : 'T',
59+
# ==========================================================
60+
61+
# Formatted Database Files Structure Parameters ============
62+
'format' : 1,
63+
'precision' : 2,
64+
'prim_vars_wrt' :'T',
65+
'parallel_io' :'F',
66+
# ==========================================================
67+
68+
# Patch 1 L ================================================
69+
'patch_icpp(1)%geometry' : 1,
70+
'patch_icpp(1)%x_centroid' : 0.25,
71+
'patch_icpp(1)%y_centroid' : 0.001,
72+
'patch_icpp(1)%length_x' : 0.5,
73+
'patch_icpp(1)%length_y' : 0.002,
74+
'patch_icpp(1)%vel(1)' : 10,
75+
'patch_icpp(1)%vel(2)' : 0,
76+
'patch_icpp(1)%pres' : 1.E+5,
77+
'patch_icpp(1)%alpha_rho(1)' : 1000,
78+
'patch_icpp(1)%alpha(1)' : 1.,
79+
'patch_icpp(1)%tau_e(1)' : 0.0,
80+
# ==========================================================
81+
82+
# Patch 2 R ================================================
83+
'patch_icpp(2)%geometry' : 1,
84+
'patch_icpp(2)%x_centroid' : 0.75,
85+
'patch_icpp(2)%y_centroid' : 0.001,
86+
'patch_icpp(2)%length_x' : 0.5,
87+
'patch_icpp(2)%length_y' : 0.002,
88+
'patch_icpp(2)%vel(1)' : -10, #10,
89+
'patch_icpp(2)%vel(2)' : 0,
90+
'patch_icpp(2)%pres' : 1.E+05,
91+
'patch_icpp(2)%alpha_rho(1)' : 1000,
92+
'patch_icpp(2)%alpha(1)' : 1.,
93+
'patch_icpp(2)%tau_e(1)' : 0.0,
94+
# ==========================================================
95+
96+
# Fluids Physical Parameters ===============================
97+
'fluid_pp(1)%gamma' : 1.E+00/(4.4E+00-1.E+00),
98+
'fluid_pp(1)%pi_inf' : 4.4E+00*6.E+08/(4.4E+00 - 1.E+00),
99+
'fluid_pp(1)%G' : 1.E+00, #.E+010,
100+
# ==========================================================
101+
}))
102+
# ==============================================================================

‎examples/1D_hyper_impact_weak/original.py‎

Lines changed: 102 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,102 @@
1+
#!/usr/bin/env python3
2+
import math
3+
import json
4+
5+
#Numerical setup
6+
Nx = 201
7+
dx = 1./(1.*(Nx+1))
8+
9+
Tend = 64E-06
10+
Nt = 4000
11+
mydt = Tend/(1.E+00*Nt)
12+
#print(mydt)
13+
# Configuring case dictionary
14+
print(json.dumps({
15+
# Logistics ================================================
16+
'run_time_info' : 'T',
17+
# ==========================================================
18+
19+
# Computational Domain Parameters ==========================
20+
'x_domain%beg' : 0.E+00,
21+
'x_domain%end' : 1.E+00,
22+
# 'y_domain%beg' : 0.E+00,
23+
# 'y_domain%end' : 0.002,
24+
'm' : Nx,
25+
'n' : 0,
26+
'p' : 0,
27+
'dt' : mydt,
28+
't_step_start' : 0,
29+
't_step_stop' : int(Nt),
30+
't_step_save' : int(Nt/200),
31+
# ==========================================================
32+
33+
# Simulation Algorithm Parameters ==========================
34+
'num_patches' : 2,
35+
'model_eqns' : 3,
36+
'alt_soundspeed' : 'F',
37+
'num_fluids' : 1,
38+
'mpp_lim' : 'F',
39+
'mixture_err' : 'F',
40+
'time_stepper' : 3,
41+
'weno_order' : 5,
42+
'weno_eps' : 1.E-16,
43+
'weno_Re_flux' : 'F',
44+
'weno_avg' : 'F',
45+
'mapped_weno' : 'F',
46+
'null_weights' : 'F',
47+
'mp_weno' : 'T',
48+
'riemann_solver' : 2,
49+
'wave_speeds' : 1,
50+
'avg_state' : 2,
51+
'bc_x%beg' : -3,
52+
'bc_x%end' : -3,
53+
#'bc_y%beg' : -3,
54+
#'bc_y%end' : -3,
55+
# ==========================================================
56+
57+
# Turning on Hyperelasticity ================================
58+
'hyperelasticity' : 'T',
59+
# ==========================================================
60+
61+
# Formatted Database Files Structure Parameters ============
62+
'format' : 1,
63+
'precision' : 2,
64+
'prim_vars_wrt' :'T',
65+
'parallel_io' :'F',
66+
# ==========================================================
67+
68+
# Patch 1 L ================================================
69+
'patch_icpp(1)%geometry' : 1,
70+
'patch_icpp(1)%x_centroid' : 0.25,
71+
# 'patch_icpp(1)%y_centroid' : 0.001,
72+
'patch_icpp(1)%length_x' : 0.5,
73+
# 'patch_icpp(1)%length_y' : 0.002,
74+
'patch_icpp(1)%vel(1)' : 10,
75+
# 'patch_icpp(1)%vel(2)' : 100*0,
76+
'patch_icpp(1)%pres' : 1.E+5,
77+
'patch_icpp(1)%alpha_rho(1)' : 1000,
78+
'patch_icpp(1)%alpha(1)' : 1.,
79+
'patch_icpp(1)%tau_e(1)' : 0.0,
80+
# ==========================================================
81+
82+
# Patch 2 R ================================================
83+
'patch_icpp(2)%geometry' : 1,
84+
'patch_icpp(2)%x_centroid' : 0.75,
85+
# 'patch_icpp(2)%y_centroid' : 0.001,
86+
'patch_icpp(2)%length_x' : 0.5,
87+
# 'patch_icpp(2)%length_y' : 0.002,
88+
'patch_icpp(2)%vel(1)' : 10, #-10,
89+
# 'patch_icpp(2)%vel(2)' : -100*0,
90+
'patch_icpp(2)%pres' : 1.E+05,
91+
'patch_icpp(2)%alpha_rho(1)' : 1000,
92+
'patch_icpp(2)%alpha(1)' : 1.,
93+
'patch_icpp(2)%tau_e(1)' : 0.0,
94+
# ==========================================================
95+
96+
# Fluids Physical Parameters ===============================
97+
'fluid_pp(1)%gamma' : 1.E+00/(4.4E+00-1.E+00),
98+
'fluid_pp(1)%pi_inf' : 4.4E+00*6.E+08/(4.4E+00 - 1.E+00),
99+
'fluid_pp(1)%G' : 1.E+00, #.E+010,
100+
# ==========================================================
101+
}))
102+
# ==============================================================================

0 commit comments

Comments
 (0)