Skip to content

Commit d373ebf

Browse files
haocheyhaochey
committed
Test tolerance and more examples
1 parent 2f891f6 commit d373ebf

File tree

16 files changed

+42463
-77
lines changed

16 files changed

+42463
-77
lines changed

‎examples/2D_ibm_stl_MFCCharacter/Character_IBM.STL‎

Lines changed: 1017 additions & 0 deletions
Large diffs are not rendered by default.

‎examples/2D_ibm_stl_MFCCharacter/case.py‎

Lines changed: 95 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,95 @@
1+
#!/usr/bin/env python3
2+
3+
import json
4+
import math
5+
6+
# Dynamic Viscosity
7+
Mu1 = 0.0000184
8+
#Mu2 = 0.01
9+
rho1 = 1.19 #0.2199
10+
gam_a = 1.4
11+
# Patch Design
12+
D = 0.1
13+
14+
# Configuring case dictionary
15+
print(json.dumps({
16+
# Logistics ================================================================
17+
'run_time_info' : 'T',
18+
# ==========================================================================
19+
20+
# Computational Domain Parameters ==========================================
21+
'x_domain%beg' : -1*D,
22+
'x_domain%end' : 1*D,
23+
'y_domain%beg' : -0.3*D,
24+
'y_domain%end' : 0.3*D,
25+
'm' : 499,
26+
'n' : 149,
27+
'p' : 0,
28+
'dt' : 1.0E-6,
29+
't_step_start' : 0,
30+
't_step_stop' : 1000,
31+
't_step_save' : 10,
32+
# ==========================================================================
33+
34+
# Simulation Algorithm Parameters ==========================================
35+
'num_patches' : 1,
36+
'model_eqns' : 2,
37+
'alt_soundspeed' : 'F',
38+
'num_fluids' : 1,
39+
'time_stepper' : 3,
40+
'weno_order' : 5,
41+
'weno_eps' : 1.E-16,
42+
'weno_Re_flux' : 'F',
43+
'weno_avg' : 'F',
44+
'avg_state' : 2,
45+
'mapped_weno' : 'T',
46+
'null_weights' : 'F',
47+
'mp_weno' : 'F',
48+
'riemann_solver' : 2,
49+
'wave_speeds' : 1,
50+
'viscous' : 'T',
51+
'bc_x%beg' : -3,
52+
'bc_x%end' : -3,
53+
'bc_y%beg' : -3,
54+
'bc_y%end' : -3,
55+
'ib' : 'T',
56+
'num_ibs' : 1,
57+
# ==========================================================================
58+
59+
# Formatted Database Files Structure Parameters ============================
60+
# Export primitive variables in double precision with parallel
61+
# I/O to minimize I/O computational time during large simulations
62+
'format' : 1,
63+
'precision' : 2,
64+
'prim_vars_wrt' :'T',
65+
'parallel_io' :'T',
66+
# ==========================================================================
67+
68+
# Patch: Middle ============================================================
69+
'patch_icpp(1)%geometry' : 3,
70+
'patch_icpp(1)%x_centroid' : 0,
71+
'patch_icpp(1)%y_centroid' : 0,
72+
'patch_icpp(1)%length_x' : 1000*D,
73+
'patch_icpp(1)%length_y' : 1000*D,
74+
'patch_icpp(1)%vel(1)' : 10, #527.2E+00,
75+
'patch_icpp(1)%vel(2)' : 0.0E+00,
76+
'patch_icpp(1)%pres' : 100000, #10918.2549,
77+
'patch_icpp(1)%alpha_rho(1)' : (1.0)*rho1,
78+
'patch_icpp(1)%alpha(1)' : 1.0,
79+
# ==========================================================================
80+
'patch_ib(1)%geometry' : 5,
81+
'patch_ib(1)%model%filepath' : 'Character_IBM.stl',
82+
'patch_ib(1)%model%translate(1)' : -0.0469828735,
83+
'patch_ib(1)%model%translate(2)' : -0.0094892600,
84+
'patch_ib(1)%model%spc' : 100,
85+
'patch_ib(1)%model%threshold' : 0.95,
86+
'patch_ib(1)%slip' : 'F',
87+
# # ========================================================================
88+
89+
# Fluids Physical Parameters ===============================================
90+
'fluid_pp(1)%gamma' : 1.E+00/(gam_a-1.E+00),
91+
'fluid_pp(1)%pi_inf' : 0,
92+
'fluid_pp(1)%Re(1)' : 1000, #7535533.2,
93+
# ==========================================================================
94+
}))
95+

‎examples/2D_ibm_stl_test/case.py‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -9,7 +9,7 @@
99
rho1 = 1.19 #0.2199
1010
gam_a = 1.4
1111
# Patch Design
12-
D = 0.1
12+
D = 2
1313

1414
# Configuring case dictionary
1515
print(json.dumps({
@@ -42,7 +42,7 @@
4242
'weno_Re_flux' : 'F',
4343
'weno_avg' : 'F',
4444
'avg_state' : 2,
45-
'mapped_weno' : 'F',
45+
'mapped_weno' : 'T',
4646
'null_weights' : 'F',
4747
'mp_weno' : 'F',
4848
'riemann_solver' : 2,
@@ -82,7 +82,7 @@
8282
'patch_ib(1)%model%translate(1)' : -0.05,
8383
'patch_ib(1)%model%translate(2)' : -0.05,
8484
'patch_ib(1)%model%spc' : 200,
85-
'patch_ib(1)%model%threshold' : 0.99,
85+
'patch_ib(1)%model%threshold' : 0.95,
8686
'patch_ib(1)%slip' : 'F',
8787
# # ========================================================================
8888

‎examples/2D_ibm_stl_wedge/Wedge2D_IBM.STL‎

Lines changed: 16 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,16 @@
1+
solid Wedge2D_IBM
2+
facet normal -0.000000e+00 -0.000000e+00 1.000000e+00
3+
outer loop
4+
vertex 3.673940e-19 3.639702e-02 0.000000e+00
5+
vertex 1.000000e-01 0.000000e+00 0.000000e+00
6+
vertex 0.000000e+00 3.839702e-02 0.000000e+00
7+
endloop
8+
endfacet
9+
facet normal -0.000000e+00 0.000000e+00 1.000000e+00
10+
outer loop
11+
vertex 0.000000e+00 3.839702e-02 0.000000e+00
12+
vertex 1.000000e-01 0.000000e+00 0.000000e+00
13+
vertex 1.000000e-01 7.479405e-02 0.000000e+00
14+
endloop
15+
endfacet
16+
endsolid Wedge2D_IBM

‎examples/2D_ibm_stl_wedge/case.py‎

Lines changed: 95 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,95 @@
1+
#!/usr/bin/env python3
2+
3+
import json
4+
import math
5+
6+
# Dynamic Viscosity
7+
Mu1 = 0.0000184
8+
#Mu2 = 0.01
9+
rho1 = 1.19 #0.2199
10+
gam_a = 1.4
11+
# Patch Design
12+
D = 0.1
13+
14+
# Configuring case dictionary
15+
print(json.dumps({
16+
# Logistics ================================================================
17+
'run_time_info' : 'T',
18+
# ==========================================================================
19+
20+
# Computational Domain Parameters ==========================================
21+
'x_domain%beg' : -3*D,
22+
'x_domain%end' : 3*D,
23+
'y_domain%beg' : -1.5*D,
24+
'y_domain%end' : 1.5*D,
25+
'm' : 399,
26+
'n' : 199,
27+
'p' : 0,
28+
'dt' : 1.0E-6,
29+
't_step_start' : 0,
30+
't_step_stop' : 1000,
31+
't_step_save' : 10,
32+
# ==========================================================================
33+
34+
# Simulation Algorithm Parameters ==========================================
35+
'num_patches' : 1,
36+
'model_eqns' : 2,
37+
'alt_soundspeed' : 'F',
38+
'num_fluids' : 1,
39+
'time_stepper' : 3,
40+
'weno_order' : 5,
41+
'weno_eps' : 1.E-16,
42+
'weno_Re_flux' : 'F',
43+
'weno_avg' : 'F',
44+
'avg_state' : 2,
45+
'mapped_weno' : 'T',
46+
'null_weights' : 'F',
47+
'mp_weno' : 'F',
48+
'riemann_solver' : 2,
49+
'wave_speeds' : 1,
50+
'viscous' : 'T',
51+
'bc_x%beg' : -3,
52+
'bc_x%end' : -3,
53+
'bc_y%beg' : -3,
54+
'bc_y%end' : -3,
55+
'ib' : 'T',
56+
'num_ibs' : 1,
57+
# ==========================================================================
58+
59+
# Formatted Database Files Structure Parameters ============================
60+
# Export primitive variables in double precision with parallel
61+
# I/O to minimize I/O computational time during large simulations
62+
'format' : 1,
63+
'precision' : 2,
64+
'prim_vars_wrt' :'T',
65+
'parallel_io' :'T',
66+
# ==========================================================================
67+
68+
# Patch: Middle ============================================================
69+
'patch_icpp(1)%geometry' : 3,
70+
'patch_icpp(1)%x_centroid' : 0,
71+
'patch_icpp(1)%y_centroid' : 0,
72+
'patch_icpp(1)%length_x' : 1000*D,
73+
'patch_icpp(1)%length_y' : 1000*D,
74+
'patch_icpp(1)%vel(1)' : 527.2E+00,
75+
'patch_icpp(1)%vel(2)' : 0.0E+00,
76+
'patch_icpp(1)%pres' : 10918.2549,
77+
'patch_icpp(1)%alpha_rho(1)' : (1.0)*rho1,
78+
'patch_icpp(1)%alpha(1)' : 1.0,
79+
# ==========================================================================
80+
'patch_ib(1)%geometry' : 5,
81+
'patch_ib(1)%model%filepath' : 'Wedge2D_IBM.stl',
82+
'patch_ib(1)%model%translate(1)' : -0.0500000000,
83+
'patch_ib(1)%model%translate(2)' : -0.0373970250,
84+
'patch_ib(1)%model%spc' : 200,
85+
'patch_ib(1)%model%threshold' : 0.90,
86+
'patch_ib(1)%slip' : 'F',
87+
# # ========================================================================
88+
89+
# Fluids Physical Parameters ===============================================
90+
'fluid_pp(1)%gamma' : 1.E+00/(gam_a-1.E+00),
91+
'fluid_pp(1)%pi_inf' : 0,
92+
'fluid_pp(1)%Re(1)' : 7535533.2,
93+
# ==========================================================================
94+
}))
95+

0 commit comments

Comments
 (0)