|
32 | 32 | #define pulse |
33 | 33 | P_amp = 10.E+6 |
34 | 34 | P_len = 45 #length of the impulse |
35 | | -theta = -math.pi/2 #direction of propagation |
| 35 | +theta = -math.pi/2 #direction of propagation |
36 | 36 |
|
37 | 37 | #non-dim |
38 | 38 |
|
|
44 | 44 | stress_char = rho_char*c_char*c_char/gammag |
45 | 45 |
|
46 | 46 | #non-dim the properties |
47 | | -rhog_n = rhog/rho_char |
48 | | -c_g_n = c_g/c_char |
| 47 | +rhog_n = rhog/rho_char |
| 48 | +c_g_n = c_g/c_char |
49 | 49 | rhol_n = rhol/rho_char |
50 | 50 | c_l_n = c_l/c_char |
51 | 51 | Bg_n = Bg/stress_char |
|
73 | 73 | time_end = 2.5 |
74 | 74 | cfl = 0.5 |
75 | 75 |
|
76 | | -dt = cfl * dx/c_l |
| 76 | +dt = cfl * dx/c_l |
77 | 77 | Nt = int(time_end/dt) |
78 | 78 | Nframes = 50000 |
79 | 79 | tstart = 0 |
80 | 80 | tstop = Nt |
81 | 81 | tsave = int(Nt/Nframes) |
82 | 82 |
|
83 | 83 | #interface profile |
84 | | -interface_amp = 0.5 |
| 84 | +interface_amp = 0.5 |
85 | 85 |
|
86 | 86 | # Configuring case dictionary |
87 | 87 | print(json.dumps({ |
|
108 | 108 | 'model_eqns' : 2, |
109 | 109 | 'alt_soundspeed' : 'F', |
110 | 110 | 'num_fluids' : 2, |
111 | | - 'adv_alphan' : 'T', |
112 | 111 | 'mpp_lim' : 'T', |
113 | 112 | 'mixture_err' : 'T', |
114 | 113 | 'time_stepper' : 3, |
115 | 114 | 'weno_order' : 5, |
116 | 115 | 'weno_eps' : 1.E-16, |
117 | | - 'weno_Re_flux' : 'F', |
| 116 | + 'weno_Re_flux' : 'F', |
118 | 117 | 'weno_avg' : 'F', |
119 | 118 | 'mapped_weno' : 'T', |
120 | 119 | 'null_weights' : 'F', |
|
134 | 133 | 'prim_vars_wrt' :'T', |
135 | 134 | 'parallel_io' :'T', |
136 | 135 | # ========================================================================== |
137 | | - |
| 136 | + |
138 | 137 | # Monopole setting ========================================================= |
139 | | - 'Monopole' : 'T', # creating an acoustic wave |
140 | | - 'num_mono' : 1, # place in the middle and expand |
141 | | - 'Mono(1)%pulse' : 3, # sine wave |
142 | | - 'Mono(1)%npulse' : 1, # 1 pulse |
143 | | - 'Mono(1)%mag' : 10.0*patmos_n, # magnitude |
144 | | - 'Mono(1)%length' : 1*dlengx, # impulse length |
145 | | - 'Mono(1)%loc(1)' : dlengx/2, # x_center of the domain |
146 | | - 'Mono(1)%loc(2)' : 5.0*dlengx, # upper boundary of the domain |
147 | | - 'Mono(1)%dir' : -math.pi/2, # direction: -pi/2 |
| 138 | + 'acoustic_source' : 'T', # creating an acoustic wave |
| 139 | + 'num_source' : 1, # place in the middle and expand |
| 140 | + 'acoustic(1)%pulse' : 3, # sine wave |
| 141 | + 'acoustic(1)%npulse' : 1, # 1 pulse |
| 142 | + 'acoustic(1)%mag' : 10.0*patmos_n, # magnitude |
| 143 | + 'acoustic(1)%length' : 1*dlengx, # impulse length |
| 144 | + 'acoustic(1)%loc(1)' : dlengx/2, # x_center of the domain |
| 145 | + 'acoustic(1)%loc(2)' : 5.0*dlengx, # upper boundary of the domain |
| 146 | + 'acoustic(1)%dir' : -math.pi/2, # direction: -pi/2 |
148 | 147 | #============================================================================== |
149 | | - |
| 148 | + |
150 | 149 | # Patch 1: Background ====================================================== |
151 | 150 | 'patch_icpp(1)%geometry' : 3, |
152 | 151 | 'patch_icpp(1)%x_centroid' : dlengx/2, |
|
186 | 185 | 'fluid_pp(1)%pi_inf' : gammal*Bl_n/(gammal-1.E+00), |
187 | 186 | 'fluid_pp(2)%gamma' : 1.E+00/(gammag-1.E+00), |
188 | 187 | 'fluid_pp(2)%pi_inf' : gammag*Bg_n/(gammag-1.E+00), |
189 | | - # ========================================================================== |
190 | | - |
| 188 | + # ========================================================================== |
| 189 | + |
191 | 190 | })) |
192 | 191 |
|
193 | 192 | # ============================================================================== |
0 commit comments