|
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