|
35 | 35 |
|
36 | 36 | # Configuring case dictionary |
37 | 37 | print(json.dumps({ |
38 | | - # Logistics ================================================ |
39 | | - 'run_time_info' : 'F', |
| 38 | + # Logistics ================================================ |
| 39 | + 'run_time_info' : 'F', |
40 | 40 | # ========================================================== |
41 | | - |
| 41 | + |
42 | 42 | # Computational Domain Parameters ========================== |
43 | | - 'x_domain%beg' : 0.E+00, |
44 | | - 'x_domain%end' : 1.E+00, |
| 43 | + 'x_domain%beg' : 0.E+00, |
| 44 | + 'x_domain%end' : 1.E+00, |
45 | 45 | 'y_domain%beg' : 0, |
46 | 46 | 'y_domain%end' : 0.5, |
47 | 47 | 'z_domain%beg' : 0., |
48 | | - 'z_domain%end' : 0.5, |
49 | | - 'm' : Nx, |
50 | | - 'n' : Ny, |
51 | | - 'p' : Nz, |
52 | | - 'dt' : mydt, |
53 | | - 't_step_start' : 0, |
| 48 | + 'z_domain%end' : 0.5, |
| 49 | + 'm' : Nx, |
| 50 | + 'n' : Ny, |
| 51 | + 'p' : Nz, |
| 52 | + 'dt' : mydt, |
| 53 | + 't_step_start' : 0, |
54 | 54 | 't_step_stop' : int(20*(95*size + 5)), |
55 | 55 | 't_step_save' : int(20*(95*size + 5)), |
56 | | - # ========================================================== |
57 | | - |
| 56 | + # ========================================================== |
| 57 | + |
58 | 58 | # Simulation Algorithm Parameters ========================== |
59 | | - 'num_patches' : 2, |
60 | | - 'model_eqns' : 2, |
61 | | - 'alt_soundspeed' : 'F', |
62 | | - 'num_fluids' : 1, |
63 | | - 'mpp_lim' : 'F', |
64 | | - 'mixture_err' : 'F', |
65 | | - 'time_stepper' : 3, |
66 | | - 'weno_order' : 5, |
| 59 | + 'num_patches' : 2, |
| 60 | + 'model_eqns' : 2, |
| 61 | + 'alt_soundspeed' : 'F', |
| 62 | + 'num_fluids' : 1, |
| 63 | + 'mpp_lim' : 'F', |
| 64 | + 'mixture_err' : 'F', |
| 65 | + 'time_stepper' : 3, |
| 66 | + 'weno_order' : 5, |
67 | 67 | 'weno_eps' : 1.E-16, |
68 | | - 'weno_Re_flux' : 'F', |
| 68 | + 'weno_Re_flux' : 'F', |
69 | 69 | 'weno_avg' : 'F', |
70 | | - 'mapped_weno' : 'F', |
71 | | - 'null_weights' : 'F', |
72 | | - 'mp_weno' : 'F', |
73 | | - 'riemann_solver' : 1, |
74 | | - 'wave_speeds' : 1, |
75 | | - 'avg_state' : 2, |
76 | | - 'bc_x%beg' : -3, |
| 70 | + 'mapped_weno' : 'F', |
| 71 | + 'null_weights' : 'F', |
| 72 | + 'mp_weno' : 'F', |
| 73 | + 'riemann_solver' : 1, |
| 74 | + 'wave_speeds' : 1, |
| 75 | + 'avg_state' : 2, |
| 76 | + 'bc_x%beg' : -3, |
77 | 77 | 'bc_x%end' : -3, |
78 | 78 | 'bc_y%beg' : -3, |
79 | 79 | 'bc_y%end' : -3, |
80 | 80 | 'bc_z%beg' : -3, |
81 | | - 'bc_z%end' : -3, |
| 81 | + 'bc_z%end' : -3, |
82 | 82 | # ========================================================== |
83 | 83 |
|
84 | 84 | # Turning on Hypoelasticity ================================ |
85 | | - 'hypoelasticity' : 'T', |
| 85 | + 'hypoelasticity' : 'T', |
| 86 | + 'fd_order' : 4, |
86 | 87 | # ========================================================== |
87 | | - |
| 88 | + |
88 | 89 | # Formatted Database Files Structure Parameters ============ |
89 | | - 'format' : 1, |
90 | | - 'precision' : 2, |
91 | | - 'prim_vars_wrt' :'T', |
92 | | - 'parallel_io' :'T', |
93 | | - # ========================================================== |
94 | | - |
95 | | - # Patch 1 L ================================================ |
96 | | - 'patch_icpp(1)%geometry' : 9, |
97 | | - 'patch_icpp(1)%x_centroid' : 0.25, |
| 90 | + 'format' : 1, |
| 91 | + 'precision' : 2, |
| 92 | + 'prim_vars_wrt' :'T', |
| 93 | + 'parallel_io' :'T', |
| 94 | + # ========================================================== |
| 95 | + |
| 96 | + # Patch 1 L ================================================ |
| 97 | + 'patch_icpp(1)%geometry' : 9, |
| 98 | + 'patch_icpp(1)%x_centroid' : 0.25, |
98 | 99 | 'patch_icpp(1)%y_centroid' : 0.25, |
99 | | - 'patch_icpp(1)%z_centroid' : 0.25, |
| 100 | + 'patch_icpp(1)%z_centroid' : 0.25, |
100 | 101 | 'patch_icpp(1)%length_x' : 0.5, |
101 | 102 | 'patch_icpp(1)%length_y' : 0.5, |
102 | | - 'patch_icpp(1)%length_z' : 0.5, |
103 | | - 'patch_icpp(1)%vel(1)' : 0.0, |
| 103 | + 'patch_icpp(1)%length_z' : 0.5, |
| 104 | + 'patch_icpp(1)%vel(1)' : 0.0, |
104 | 105 | 'patch_icpp(1)%vel(2)' : 0, |
105 | | - 'patch_icpp(1)%vel(3)' : 0, |
106 | | - 'patch_icpp(1)%pres' : 1.E+8, |
107 | | - 'patch_icpp(1)%alpha_rho(1)' : 1000, |
108 | | - 'patch_icpp(1)%alpha(1)' : 1., |
109 | | - 'patch_icpp(1)%tau_e(1)' : 0.0, |
| 106 | + 'patch_icpp(1)%vel(3)' : 0, |
| 107 | + 'patch_icpp(1)%pres' : 1.E+8, |
| 108 | + 'patch_icpp(1)%alpha_rho(1)' : 1000, |
| 109 | + 'patch_icpp(1)%alpha(1)' : 1., |
| 110 | + 'patch_icpp(1)%tau_e(1)' : 0.0, |
110 | 111 | # ========================================================== |
111 | 112 |
|
112 | 113 | # Patch 2 R ================================================ |
113 | | - 'patch_icpp(2)%geometry' : 9, |
| 114 | + 'patch_icpp(2)%geometry' : 9, |
114 | 115 | 'patch_icpp(2)%x_centroid' : 0.75, |
115 | 116 | 'patch_icpp(2)%y_centroid' : 0.25, |
116 | | - 'patch_icpp(2)%z_centroid' : 0.25, |
117 | | - 'patch_icpp(2)%length_x' : 0.5, |
| 117 | + 'patch_icpp(2)%z_centroid' : 0.25, |
| 118 | + 'patch_icpp(2)%length_x' : 0.5, |
118 | 119 | 'patch_icpp(2)%length_y' : 0.5, |
119 | | - 'patch_icpp(2)%length_z' : 0.5, |
120 | | - 'patch_icpp(2)%vel(1)' : 0, |
| 120 | + 'patch_icpp(2)%length_z' : 0.5, |
| 121 | + 'patch_icpp(2)%vel(1)' : 0, |
121 | 122 | 'patch_icpp(2)%vel(2)' : 0, |
122 | | - 'patch_icpp(2)%vel(3)' : 0, |
123 | | - 'patch_icpp(2)%pres' : 1.E+05, |
124 | | - 'patch_icpp(2)%alpha_rho(1)' : 1000, |
125 | | - 'patch_icpp(2)%alpha(1)' : 1., |
126 | | - 'patch_icpp(2)%tau_e(1)' : 0.0, |
| 123 | + 'patch_icpp(2)%vel(3)' : 0, |
| 124 | + 'patch_icpp(2)%pres' : 1.E+05, |
| 125 | + 'patch_icpp(2)%alpha_rho(1)' : 1000, |
| 126 | + 'patch_icpp(2)%alpha(1)' : 1., |
| 127 | + 'patch_icpp(2)%tau_e(1)' : 0.0, |
127 | 128 | # ========================================================== |
128 | 129 |
|
129 | 130 | # Fluids Physical Parameters =============================== |
130 | | - 'fluid_pp(1)%gamma' : 1.E+00/(4.4E+00-1.E+00), |
131 | | - 'fluid_pp(1)%pi_inf' : 4.4E+00*6.E+08/(4.4E+00 - 1.E+00), |
132 | | - 'fluid_pp(1)%G' : 10E+09, |
| 131 | + 'fluid_pp(1)%gamma' : 1.E+00/(4.4E+00-1.E+00), |
| 132 | + 'fluid_pp(1)%pi_inf' : 4.4E+00*6.E+08/(4.4E+00 - 1.E+00), |
| 133 | + 'fluid_pp(1)%G' : 10E+09, |
133 | 134 | # ========================================================== |
134 | 135 | })) |
135 | 136 | # ============================================================================== |
0 commit comments