Skip to content

Commit 929a820

Browse files
mcarcanamcarcana
committed
debugging hyperelasticity
1 parent 2d0d572 commit 929a820

File tree

7 files changed

+825
-53
lines changed

7 files changed

+825
-53
lines changed

‎examples/1D_hyper_1material/case.py‎

Lines changed: 107 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,107 @@
1+
#!/usr/bin/env python3
2+
3+
# 1D -> Hyperelasticity -> 1 Fluid(s)
4+
5+
import math, json
6+
7+
print(json.dumps({
8+
'run_time_info': 'T',
9+
'm': 299,
10+
'n': 0,
11+
'p': 0,
12+
'dt': 1e-07,
13+
't_step_start': 0,
14+
't_step_stop': 50,
15+
't_step_save': 25,
16+
'num_patches': 3,
17+
'model_eqns': 3,
18+
'alt_soundspeed': 'F',
19+
'num_fluids': 1,
20+
'mpp_lim': 'F',
21+
'mixture_err': 'F',
22+
'time_stepper': 3,
23+
'weno_order': 5,
24+
'weno_eps': 1e-16,
25+
'mapped_weno': 'F',
26+
'null_weights': 'F',
27+
'mp_weno': 'F',
28+
'riemann_solver': 2,
29+
'wave_speeds': 1,
30+
'avg_state': 2,
31+
'format': 1,
32+
'precision': 2,
33+
'prim_vars_wrt': 'F',
34+
'parallel_io': 'F',
35+
'patch_icpp(1)%pres': 1000000.0,
36+
'patch_icpp(1)%alpha_rho(1)': 1000.0,
37+
'patch_icpp(1)%alpha(1)': 1.0,
38+
'patch_icpp(2)%pres': 100000.0,
39+
'patch_icpp(2)%alpha_rho(1)': 1000.0,
40+
'patch_icpp(2)%alpha(1)': 1.0,
41+
'patch_icpp(3)%pres': 500000.0,
42+
'patch_icpp(3)%alpha_rho(1)': 1000.0,
43+
'patch_icpp(3)%alpha(1)': 1.0,
44+
'fluid_pp(1)%gamma': 0.3,
45+
'fluid_pp(1)%pi_inf': 780000.0,
46+
'Ca': 0.9769178386380458,
47+
'Web': 13.927835051546392,
48+
'Re_inv': 0.009954269975623245,
49+
'pref': 101325.0,
50+
'rhoref': 1000.0,
51+
'bubble_model': 3,
52+
'polytropic': 'T',
53+
'polydisperse': 'F',
54+
'thermal': 3,
55+
'R0ref': 1e-05,
56+
'patch_icpp(1)%r0': 1,
57+
'patch_icpp(1)%v0': 0,
58+
'patch_icpp(2)%r0': 1,
59+
'patch_icpp(2)%v0': 0,
60+
'patch_icpp(3)%r0': 1,
61+
'patch_icpp(3)%v0': 0,
62+
'qbmm': 'F',
63+
'dist_type': 2,
64+
'poly_sigma': 0.3,
65+
'R0_type': 1,
66+
'sigR': 0.1,
67+
'sigV': 0.1,
68+
'rhoRV': 0.0,
69+
'x_domain%beg': 0.0,
70+
'x_domain%end': 1.0,
71+
'bc_x%beg': -3,
72+
'bc_x%end': -3,
73+
'patch_icpp(1)%geometry': 1,
74+
'patch_icpp(1)%x_centroid': 0.05,
75+
'patch_icpp(1)%length_x': 0.1,
76+
'patch_icpp(2)%x_centroid': 0.45,
77+
'patch_icpp(2)%length_x': 0.7,
78+
'patch_icpp(3)%x_centroid': 0.9,
79+
'patch_icpp(3)%length_x': 0.2,
80+
'patch_icpp(1)%vel(1)': 0.0,
81+
'patch_icpp(2)%geometry': 1,
82+
'patch_icpp(2)%vel(1)': 0.0,
83+
'patch_icpp(3)%geometry': 1,
84+
'patch_icpp(3)%vel(1)': 0.0,
85+
'hyperelasticity': 'T',
86+
'hyper_model': 1,
87+
'fd_order': 4,
88+
'patch_icpp(1)%tau_e(1)': 0.0,
89+
'patch_icpp(2)%tau_e(1)': 0.0,
90+
'patch_icpp(3)%tau_e(1)': 0.0,
91+
'fluid_pp(1)%G': 1000000000.0,
92+
'parallel_io' : 'T',
93+
'cons_vars_wrt' : 'T',
94+
'prim_vars_wrt': 'T',
95+
'alpha_rho_wrt(1)': 'T',
96+
'rho_wrt' : 'T',
97+
'mom_wrt(1)' : 'T',
98+
'vel_wrt(1)' : 'T',
99+
'E_wrt' : 'T',
100+
'pres_wrt' : 'T',
101+
'alpha_wrt(1)' : 'T',
102+
'gamma_wrt' : 'T',
103+
'heat_ratio_wrt' : 'T',
104+
'pi_inf_wrt' : 'T',
105+
'pres_inf_wrt' : 'T',
106+
'c_wrt' : 'T',
107+
}))

‎examples/1D_hyper_2materials/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+
# 1D -> Hyperelasticity -> 2 Fluid(s)
4+
5+
import math, json
6+
import argparse
7+
8+
print(json.dumps({
9+
'run_time_info': 'T',
10+
'm': 299,
11+
'n': 0,
12+
'p': 0,
13+
'dt': 1e-07,
14+
't_step_start': 0,
15+
't_step_stop': 50,
16+
't_step_save': 25,
17+
'num_patches': 3,
18+
'model_eqns': 3,
19+
'alt_soundspeed': 'F',
20+
'num_fluids': 2,
21+
'mpp_lim': 'F',
22+
'mixture_err': 'F',
23+
'time_stepper': 3,
24+
'weno_order': 5,
25+
'weno_eps': 1e-16,
26+
'mapped_weno': 'F',
27+
'null_weights': 'F',
28+
'mp_weno': 'F',
29+
'riemann_solver': 2,
30+
'wave_speeds': 1,
31+
'avg_state': 2,
32+
'format': 1,
33+
'precision': 2,
34+
'prim_vars_wrt': 'F',
35+
'parallel_io': 'F',
36+
'patch_icpp(1)%pres': 1000000.0,
37+
'patch_icpp(1)%alpha_rho(1)': 900.0,
38+
'patch_icpp(1)%alpha(1)': 0.9,
39+
'patch_icpp(2)%pres': 100000.0,
40+
'patch_icpp(2)%alpha_rho(1)': 100,
41+
'patch_icpp(2)%alpha(1)': 0.1,
42+
'patch_icpp(3)%pres': 500000.0,
43+
'patch_icpp(3)%alpha_rho(1)': 900,
44+
'patch_icpp(3)%alpha(1)': 0.9,
45+
'fluid_pp(1)%gamma': 0.3,
46+
'fluid_pp(1)%pi_inf': 780000.0,
47+
# 'Ca': 0.9769178386380458,
48+
# 'Web': 13.927835051546392,
49+
# 'Re_inv': 0.009954269975623245,
50+
'pref': 101325.0,
51+
'rhoref': 1000.0,
52+
# 'bubble_model': 3,
53+
# 'polytropic': 'T',
54+
# 'polydisperse': 'F',
55+
# 'thermal': 3,
56+
'R0ref': 1e-05,
57+
'patch_icpp(1)%r0': 1,
58+
'patch_icpp(1)%v0': 0,
59+
'patch_icpp(2)%r0': 1,
60+
'patch_icpp(2)%v0': 0,
61+
'patch_icpp(3)%r0': 1,
62+
'patch_icpp(3)%v0': 0,
63+
# 'qbmm': 'F',
64+
# 'dist_type': 2,
65+
# 'poly_sigma': 0.3,
66+
# 'R0_type': 1,
67+
# 'sigR': 0.1,
68+
# 'sigV': 0.1,
69+
# 'rhoRV': 0.0,
70+
'x_domain%beg': 0.0,
71+
'x_domain%end': 1.0,
72+
'bc_x%beg': -3,
73+
'bc_x%end': -3,
74+
'patch_icpp(1)%geometry': 1,
75+
'patch_icpp(1)%x_centroid': 0.05,
76+
'patch_icpp(1)%length_x': 0.1,
77+
'patch_icpp(2)%x_centroid': 0.45,
78+
'patch_icpp(2)%length_x': 0.7,
79+
'patch_icpp(3)%x_centroid': 0.9,
80+
'patch_icpp(3)%length_x': 0.2,
81+
'patch_icpp(1)%vel(1)': 0.0,
82+
'patch_icpp(2)%geometry': 1,
83+
'patch_icpp(2)%vel(1)': 0.0,
84+
'patch_icpp(3)%geometry': 1,
85+
'patch_icpp(3)%vel(1)': 0.0,
86+
'hyperelasticity': 'T',
87+
'hyper_model': 1,
88+
'fd_order': 4,
89+
'patch_icpp(1)%tau_e(1)': 0.0,
90+
'patch_icpp(2)%tau_e(1)': 0.0,
91+
'patch_icpp(3)%tau_e(1)': 0.0,
92+
'fluid_pp(1)%G': 1000000000.0,
93+
'fluid_pp(2)%gamma': 0.3,
94+
'fluid_pp(2)%pi_inf': 780000.0,
95+
'patch_icpp(1)%alpha_rho(2)': 100,
96+
'patch_icpp(1)%alpha(2)': 0.1,
97+
'patch_icpp(2)%alpha_rho(2)': 900,
98+
'patch_icpp(2)%alpha(2)': 0.9,
99+
'patch_icpp(3)%alpha_rho(2)': 100,
100+
'patch_icpp(3)%alpha(2)': 0.1,
101+
'fluid_pp(2)%G': 0.0, #50000.0
102+
'parallel_io' : 'T', 'cons_vars_wrt' : 'T',
103+
'prim_vars_wrt': 'T', 'alpha_rho_wrt(1)': 'T',
104+
'rho_wrt' : 'T', 'mom_wrt(1)' : 'T',
105+
'vel_wrt(1)' : 'T', 'E_wrt' : 'T',
106+
'pres_wrt' : 'T', 'alpha_wrt(1)' : 'T',
107+
'gamma_wrt' : 'T', 'heat_ratio_wrt' : 'T',
108+
'pi_inf_wrt' : 'T', 'pres_inf_wrt' : 'T',
109+
'c_wrt' : 'T',
110+
}))

‎examples/2D_hyper_1material/case.py‎

Lines changed: 125 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,125 @@
1+
#!/usr/bin/env python3
2+
3+
# 2D -> Hyperelasticity -> 1 Fluid(s)
4+
5+
import math, json
6+
7+
print(json.dumps({
8+
'run_time_info': 'T',
9+
'm': 49,
10+
'n': 39,
11+
'p': 0,
12+
'dt': 1e-08, #1e-06,
13+
# 'cfl_adap_dt': 'T',
14+
# 'cfl_target': 0.1,
15+
# 'n_start': 0,
16+
# 't_stop': 1e-06,
17+
# 't_save': 1e-07,
18+
't_step_start': 0,
19+
't_step_stop': 50,
20+
't_step_save': 1,
21+
'num_patches': 3,
22+
'model_eqns': 3,
23+
'alt_soundspeed': 'F',
24+
'num_fluids': 1,
25+
'mpp_lim': 'F',
26+
'mixture_err': 'F',
27+
'time_stepper': 3,
28+
'weno_order': 5,
29+
'weno_eps': 1e-16,
30+
'mapped_weno': 'F',
31+
'null_weights': 'F',
32+
'mp_weno': 'F',
33+
'riemann_solver': 2,
34+
'wave_speeds': 1,
35+
'avg_state': 2,
36+
'format': 1,
37+
'precision': 2,
38+
'prim_vars_wrt': 'F',
39+
'parallel_io': 'F',
40+
'patch_icpp(1)%pres': 1000000.0,
41+
'patch_icpp(1)%alpha_rho(1)': 1000.0,
42+
'patch_icpp(1)%alpha(1)': 1.0,
43+
'patch_icpp(2)%pres': 100000.0,
44+
'patch_icpp(2)%alpha_rho(1)': 1000.0,
45+
'patch_icpp(2)%alpha(1)': 1.0,
46+
'patch_icpp(3)%pres': 500000.0,
47+
'patch_icpp(3)%alpha_rho(1)': 1000.0,
48+
'patch_icpp(3)%alpha(1)': 1.0,
49+
'fluid_pp(1)%gamma': 0.3,
50+
'fluid_pp(1)%pi_inf': 780000.0,
51+
# 'Ca': 0.9769178386380458,
52+
# 'Web': 13.927835051546392,
53+
# 'Re_inv': 0.009954269975623245,
54+
'pref': 101325.0,
55+
'rhoref': 1000.0,
56+
# 'bubble_model': 3,
57+
# 'polytropic': 'T',
58+
# 'polydisperse': 'F',
59+
# 'thermal': 3,
60+
'R0ref': 1e-05,
61+
'patch_icpp(1)%r0': 1,
62+
'patch_icpp(1)%v0': 0,
63+
'patch_icpp(2)%r0': 1,
64+
'patch_icpp(2)%v0': 0,
65+
'patch_icpp(3)%r0': 1,
66+
'patch_icpp(3)%v0': 0,
67+
# 'qbmm': 'F',
68+
# 'dist_type': 2,
69+
# 'poly_sigma': 0.3,
70+
# 'R0_type': 1,
71+
# 'sigR': 0.1,
72+
# 'sigV': 0.1,
73+
# 'rhoRV': 0.0,
74+
'x_domain%beg': 0.0,
75+
'x_domain%end': 1.0,
76+
'y_domain%beg': 0.0,
77+
'y_domain%end': 1.0,
78+
'bc_x%beg': -3,
79+
'bc_x%end': -3,
80+
'bc_y%beg': -3,
81+
'bc_y%end': -3,
82+
'patch_icpp(1)%geometry': 3,
83+
'patch_icpp(1)%y_centroid': 0.05,
84+
'patch_icpp(1)%length_y': 0.1,
85+
'patch_icpp(2)%y_centroid': 0.45,
86+
'patch_icpp(2)%length_y': 0.7,
87+
'patch_icpp(3)%y_centroid': 0.9,
88+
'patch_icpp(3)%length_y': 0.2,
89+
'patch_icpp(1)%x_centroid': 0.5,
90+
'patch_icpp(1)%length_x': 1,
91+
'patch_icpp(1)%vel(1)': 0.0,
92+
'patch_icpp(1)%vel(2)': 0.0,
93+
'patch_icpp(2)%geometry': 3,
94+
'patch_icpp(2)%x_centroid': 0.5,
95+
'patch_icpp(2)%length_x': 1,
96+
'patch_icpp(2)%vel(1)': 0.0,
97+
'patch_icpp(2)%vel(2)': 0.0,
98+
'patch_icpp(3)%geometry': 3,
99+
'patch_icpp(3)%x_centroid': 0.5,
100+
'patch_icpp(3)%length_x': 1,
101+
'patch_icpp(3)%vel(1)': 0.0,
102+
'patch_icpp(3)%vel(2)': 0.0,
103+
'hyperelasticity': 'T',
104+
'hyper_model': 1,
105+
'fd_order': 4,
106+
'patch_icpp(1)%tau_e(1)': 0.0,
107+
'patch_icpp(2)%tau_e(1)': 0.0,
108+
'patch_icpp(3)%tau_e(1)': 0.0,
109+
'fluid_pp(1)%G': 1000000000.0,
110+
'patch_icpp(1)%tau_e(2)': 0.0,
111+
'patch_icpp(1)%tau_e(3)': 0.0,
112+
'patch_icpp(2)%tau_e(2)': 0.0,
113+
'patch_icpp(2)%tau_e(3)': 0.0,
114+
'patch_icpp(3)%tau_e(2)': 0.0,
115+
'patch_icpp(3)%tau_e(3)': 0.0,
116+
'parallel_io' : 'T', 'cons_vars_wrt' : 'T',
117+
'prim_vars_wrt': 'T', 'alpha_rho_wrt(1)': 'T',
118+
'rho_wrt' : 'T', 'mom_wrt(1)' : 'T',
119+
'vel_wrt(1)' : 'T', 'E_wrt' : 'T',
120+
'pres_wrt' : 'T', 'alpha_wrt(1)' : 'T',
121+
'gamma_wrt' : 'T', 'heat_ratio_wrt' : 'T',
122+
'pi_inf_wrt' : 'T', 'pres_inf_wrt' : 'T',
123+
'c_wrt' : 'T',
124+
}))
125+

0 commit comments

Comments
 (0)