Skip to content

Commit 8f96e30

Browse files
haocheyhaochey
authored
Ibm+stl pr (#26)
* IBM+STL Upstream infrastructure * Add models * Still need to fix checker and docs * Fix magic number and add examples * Fix format and spell * fix spell * fix trailing white space * fix some docstring in m_patches
1 parent b41386d commit 8f96e30

File tree

14 files changed

+42264
-24
lines changed

14 files changed

+42264
-24
lines changed

‎examples/2D_IBM_STL_wedge/20deg_wedge.stl‎

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,9 @@
1+
solid 20deg_wedge
2+
facet normal 0.000000e+00 -0.000000e+00 1.000000e+00
3+
outer loop
4+
vertex 1.511113e-02 0.000000e+00 0.000000e+00
5+
vertex 1.511113e-02 1.100000e-02 0.000000e+00
6+
vertex 0.000000e+00 5.500000e-03 0.000000e+00
7+
endloop
8+
endfacet
9+
endsolid 20deg_wedge

‎examples/2D_IBM_STL_wedge/case.py‎

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

0 commit comments

Comments
 (0)