v0.10.0
droplet_adaptivity_analysis.sh
Go to the documentation of this file.
1 # analysis of automatic adapptivity performed on droplet application */
2 
3 
4 rm -rf droplet_adaptivity_pressure
5 rm -rf droplet_reference_pressure
6 rm -rf droplet_h_refinment_pressure
7 
8 # # remove all the previous time dependant results.
9 
10 mkdir droplet_adaptivity_pressure
11 mkdir droplet_reference_pressure
12 mkdir droplet_h_refinment_pressure
13 
14 # COMPUTE THE AUTOMATIC ADAPTIVITY SCHEME
15 
16 # rm -rf droplet_pressure
17 #
18 #
19 ~/build_moFEM1/usr/tools/mofem_part -my_file ./coarse_droplet.cub -my_nparts 8
20 #
21 # # 4 MHz
22 mpirun -np 8 ./fe_field_split \
23 -my_file out.h5m \
24 -my_is_partitioned true \
25 -duffy true \
26 -adaptivity true \
27 -error_level 1 \
28 -my_order 6 \
29 -save_postproc_mesh true \
30 -wave_number 6283.1 \
31 -my_max_post_proc_ref_level 0 \
32 -add_incident_wave true \
33 -wave_direction 1,0,0 \
34 -wave_oscilation_direction -0.3756,-0.9268,0.000000001 \
35 -monochromatic_wave true \
36 -signal_length 0.001 \
37 -radiation_field true \
38 -frequency 4000000 \
39 -signal_duration 0.00000025 \
40 -velocity 1500 \
41 -complex_wave_number 79.8387 \
42 -transmission_coefficient 0.1491 \
43 -rayleigh_wave false \
44 -density 998 \
45 -reynolds_stress false \
46 -radiation_force false \
47 -material_coefficient1 -1348.8 \
48 -space_data acoustic_data.txt \
49 -nb_of_time_step 14 \
50 -material_coefficient2 \
51 -fieldsplit_1_ksp_type fgmres \
52 -fieldsplit_1_pc_type lu \
53 -fieldsplit_1_pc_factor_mat_solver_package superlu_dist \
54 -fieldsplit_1_ksp_max_it 25 \
55 -fieldsplit_0_ksp_type fgmres \
56 -fieldsplit_0_pc_type lu \
57 -fieldsplit_0_pc_factor_mat_solver_package superlu_dist \
58 -fieldsplit_0_ksp_max_it 25 \
59 -pc_fieldsplit_type schur \
60 -ksp_type fgmres \
61 -ksp_atol 1e-8 \
62 -ksp_rtol 1e-8 \
63 -ksp_max_it 800 \
64 -ksp_monitor \
65 -ksp_converged_reason \
66 -droplet_radius 0.00210 \
67 -attenuation 0.067077277922418 | tee log_droplet_adaptivity_4MHZ
68 
69 
70 
71 mbconvert fe_solution_mesh_post_proc.h5m fe_solution_mesh_post_proc.vtk
72 mbconvert fe_solution_mesh_post_proc.vtk fe_solution_mesh_post_proc_4MHZ.vtk
73 mv fe_solution_mesh_post_proc_4MHZ.vtk droplet_adaptivity_pressure
74 mv log_droplet_adaptivity_4MHZ droplet_adaptivity_pressure
75 
76 # -complex_wave_number 2768 \
77 
78 #
79 mpirun -np 8 ./fe_field_split \
80 -my_file out.h5m \
81 -my_is_partitioned true \
82 -ksp_final_residual \
83 -ksp_monitor \
84 -ksp_converged_reason \
85 -pc_type lu \
86 -ksp_type fgmres \
87 -pc_factor_mat_solver_package superlu_dist \
88 -duffy true \
89 -adaptivity true \
90 -error_level 1 \
91 -my_order 6 \
92 -save_postproc_mesh true \
93 -wave_number 15707.96326794897 \
94 -my_max_post_proc_ref_level 0 \
95 -add_incident_wave true \
96 -wave_direction 1,0,0 \
97 -wave_oscilation_direction -0.3756,-0.9268,3.2051e-09 \
98 -monochromatic_wave true \
99 -signal_length 0.0004 \
100 -signal_duration 0.0000001 \
101 -frequency 10000000 \
102 -velocity 1500 \
103 -complex_wave_number 199.5968 \
104 -transmission_coefficient 0.1491 \
105 -rayleigh_wave false \
106 -density 998 \
107 -droplet_radius 0.00210 \
108 -attenuation 0.265146221848709 \
109 -material_coefficient1 -1348.8 \
110 -ksp_atol 1e-8 \
111 -ksp_rtol 1e-8 \
112 -ksp_max_it 2000 \
113 -space_data acoustic_data.txt \
114 -nb_of_time_step 13 \
115 -material_coefficient2 | tee log_droplet_adaptivity_10MHZ
116 
117 mbconvert fe_solution_mesh_post_proc.h5m fe_solution_mesh_post_proc.vtk
118 mbconvert fe_solution_mesh_post_proc.vtk fe_solution_mesh_post_proc_10MHZ.vtk
119 mv fe_solution_mesh_post_proc_10MHZ.vtk droplet_adaptivity_pressure
120 mv log_droplet_adaptivity_10MHZ droplet_adaptivity_pressure
121 
122 
123 
124 mpirun -np 8 ./fe_field_split \
125 -my_file out.h5m \
126 -my_is_partitioned true \
127 -ksp_final_residual \
128 -ksp_monitor \
129 -ksp_converged_reason \
130 -pc_type lu \
131 -ksp_type fgmres \
132 -pc_factor_mat_solver_package superlu_dist \
133 -duffy true \
134 -adaptivity true \
135 -error_level 1 \
136 -my_order 6 \
137 -save_postproc_mesh true \
138 -wave_number 31415.92653589793 \
139 -my_max_post_proc_ref_level 0 \
140 -add_incident_wave true \
141 -wave_direction 1,0,0 \
142 -wave_oscilation_direction -0.3756,-0.9268,3.2051e-09 \
143 -monochromatic_wave true \
144 -signal_length 0.0002 \
145 -signal_duration 0.00000005 \
146 -frequency 20000000 \
147 -velocity 1500 \
148 -complex_wave_number 399.1935 \
149 -transmission_coefficient 0.1491 \
150 -rayleigh_wave false \
151 -density 998 \
152 -droplet_radius 0.00210 \
153 -attenuation 0.749946765900858 \
154 -material_coefficient1 -1348.8 \
155 -ksp_atol 1e-8 \
156 -ksp_rtol 1e-8 \
157 -ksp_max_it 2000 \
158 -space_data acoustic_data.txt \
159 -nb_of_time_step 13 \
160 -material_coefficient2 | tee log_droplet_adaptivity_20MHZ
161 
162 mbconvert fe_solution_mesh_post_proc.h5m fe_solution_mesh_post_proc.vtk
163 mbconvert fe_solution_mesh_post_proc.vtk fe_solution_mesh_post_proc_20MHZ.vtk
164 mv fe_solution_mesh_post_proc_20MHZ.vtk droplet_adaptivity_pressure
165 mv log_droplet_adaptivity_20MHZ droplet_adaptivity_pressure
166 
167 #
168 # ./do_vtk1.sh pressure_real_time_step_*.h5m
169 # mbconvert reynolds_stress.h5m reynolds_stress.vtk
170 # mkdir droplet_pressure
171 # mv pressure_real_time_step_*.vtk droplet_pressure
172 # mv reynolds_stress.vtk droplet_pressure
173 
174 
175 # -complex_wave_number 199.5968 \
176 
177 
178 ##################################################################################################################################
179 
180 ##################################################################################################################################
181 
182 ##################################################################################################################################
183 
184 ##################################################################################################################################
185 
186 ##################################################################################################################################
187 
188 
189 
190 # COMPUTE THE hhhhhhhhh ADAPTIVITY SCHEME
191 
192 # rm -rf droplet_pressure
193 #
194 #
195 ~/build_moFEM1/usr/tools/mofem_part -my_file ./dense_droplet.cub -my_nparts 8
196 #
197 # # 4 MHz
198 mpirun -np 8 ./fe_field_split \
199 -my_file out.h5m \
200 -my_is_partitioned true \
201 -duffy true \
202 -adaptivity false \
203 -error_level 1 \
204 -my_order 2 \
205 -save_postproc_mesh true \
206 -wave_number 6283.1 \
207 -my_max_post_proc_ref_level 0 \
208 -add_incident_wave true \
209 -wave_direction 1,0,0 \
210 -wave_oscilation_direction -0.3756,-0.9268,0.000000001 \
211 -monochromatic_wave true \
212 -signal_length 0.001 \
213 -radiation_field true \
214 -frequency 4000000 \
215 -signal_duration 0.00000025 \
216 -velocity 1500 \
217 -complex_wave_number 79.8387 \
218 -transmission_coefficient 0.1491 \
219 -rayleigh_wave false \
220 -density 998 \
221 -reynolds_stress false \
222 -radiation_force false \
223 -material_coefficient1 -1348.8 \
224 -ksp_atol 1e-8 \
225 -ksp_rtol 1e-8 \
226 -ksp_max_it 2000 \
227 -ksp_monitor \
228 -space_data acoustic_data.txt \
229 -nb_of_time_step 14 \
230 -material_coefficient2 \
231 -ksp_converged_reason \
232 -ksp_monitor_true_residual \
233 -pc_type lu \
234 -ksp_type fgmres \
235 -pc_factor_mat_solver_package superlu_dist \
236 -droplet_radius 0.00210 \
237 -attenuation 0.067077277922418 | tee log_droplet_h_refinment_4MHZ
238 
239 
240 
241 mbconvert fe_solution_mesh_post_proc.h5m fe_solution_mesh_post_proc.vtk
242 mbconvert fe_solution_mesh_post_proc.vtk fe_solution_mesh_post_proc_4MHZ.vtk
243 mv fe_solution_mesh_post_proc_4MHZ.vtk droplet_h_refinment_pressure
244 mv log_droplet_h_refinment_4MHZ droplet_h_refinment_pressure
245 
246 # -complex_wave_number 2768 \
247 
248 #
249 mpirun -np 8 ./fe_field_split \
250 -my_file out.h5m \
251 -my_is_partitioned true \
252 -ksp_final_residual \
253 -ksp_monitor \
254 -ksp_converged_reason \
255 -pc_type lu \
256 -ksp_type fgmres \
257 -pc_factor_mat_solver_package superlu_dist \
258 -duffy true \
259 -adaptivity false \
260 -error_level 1 \
261 -my_order 2 \
262 -save_postproc_mesh true \
263 -wave_number 15707.96326794897 \
264 -my_max_post_proc_ref_level 0 \
265 -add_incident_wave true \
266 -wave_direction 1,0,0 \
267 -wave_oscilation_direction -0.3756,-0.9268,3.2051e-09 \
268 -monochromatic_wave true \
269 -signal_length 0.0004 \
270 -signal_duration 0.0000001 \
271 -frequency 10000000 \
272 -velocity 1500 \
273 -complex_wave_number 199.5968 \
274 -transmission_coefficient 0.1491 \
275 -rayleigh_wave false \
276 -density 998 \
277 -droplet_radius 0.00210 \
278 -attenuation 0.265146221848709 \
279 -material_coefficient1 -1348.8 \
280 -ksp_atol 1e-8 \
281 -ksp_rtol 1e-8 \
282 -ksp_max_it 2000 \
283 -space_data acoustic_data.txt \
284 -nb_of_time_step 13 \
285 -material_coefficient2 | tee log_droplet_h_refinment_10MHZ
286 
287 mbconvert fe_solution_mesh_post_proc.h5m fe_solution_mesh_post_proc.vtk
288 mbconvert fe_solution_mesh_post_proc.vtk fe_solution_mesh_post_proc_10MHZ.vtk
289 mv fe_solution_mesh_post_proc_10MHZ.vtk droplet_h_refinment_pressure
290 mv log_droplet_h_refinment_10MHZ droplet_h_refinment_pressure
291 
292 
293 
294 mpirun -np 8 ./fe_field_split \
295 -my_file out.h5m \
296 -my_is_partitioned true \
297 -ksp_final_residual \
298 -ksp_monitor \
299 -ksp_converged_reason \
300 -pc_type lu \
301 -ksp_type fgmres \
302 -pc_factor_mat_solver_package superlu_dist \
303 -duffy true \
304 -adaptivity false \
305 -error_level 1 \
306 -my_order 2 \
307 -save_postproc_mesh true \
308 -wave_number 31415.92653589793 \
309 -my_max_post_proc_ref_level 0 \
310 -add_incident_wave true \
311 -wave_direction 1,0,0 \
312 -wave_oscilation_direction -0.3756,-0.9268,3.2051e-09 \
313 -monochromatic_wave true \
314 -signal_length 0.0002 \
315 -signal_duration 0.00000005 \
316 -frequency 20000000 \
317 -velocity 1500 \
318 -complex_wave_number 399.1935 \
319 -transmission_coefficient 0.1491 \
320 -rayleigh_wave false \
321 -density 998 \
322 -droplet_radius 0.00210 \
323 -attenuation 0.749946765900858 \
324 -material_coefficient1 -1348.8 \
325 -ksp_atol 1e-8 \
326 -ksp_rtol 1e-8 \
327 -ksp_max_it 2000 \
328 -space_data acoustic_data.txt \
329 -nb_of_time_step 13 \
330 -material_coefficient2 | tee log_droplet_h_refinment_20MHZ
331 
332 mbconvert fe_solution_mesh_post_proc.h5m fe_solution_mesh_post_proc.vtk
333 mbconvert fe_solution_mesh_post_proc.vtk fe_solution_mesh_post_proc_20MHZ.vtk
334 mv fe_solution_mesh_post_proc_20MHZ.vtk droplet_h_refinment_pressure
335 mv log_droplet_h_refinment_20MHZ droplet_h_refinment_pressure
336 
337 ##################################################################################################################################
338 
339 ##################################################################################################################################
340 
341 ##################################################################################################################################
342 
343 ##################################################################################################################################
344 
345 ##################################################################################################################################
346 
347 ########################################## Reference solution - Pure Uniform p adaptivity ############################################
348 
349 
350 #
351 # ./do_vtk1.sh pressure_real_time_step_*.h5m
352 # mbconvert reynolds_stress.h5m reynolds_stress.vtk
353 # mkdir droplet_pressure
354 # mv pressure_real_time_step_*.vtk droplet_pressure
355 # mv reynolds_stress.vtk droplet_pressure
356 
357 
358 # -complex_wave_number 199.5968 \