v0.14.0
param_file.petsc
Go to the documentation of this file.
1 ###############################
2 #### Essential parameters #####
3 
4 # set Griffith's energy
5 -my_gc 1.4e-4
6 
7 # arc length step
8 -arc_s 1e-2
9 
10 ## arc-length step adaptation
11 
12 # desired number of iterations
13 -adapt_step_its_d 8
14 # minimal step size
15 -adapt_step_min_s 5e-5
16 # maximal step size
17 -adapt_step_max_s 1e-2
18 
19 ## Steps and cutting
20 # number of cutting steps
21 -nb_cut_steps 100
22 # number of load steps between cutting
23 -nb_load_steps 3
24 
25 ## Approximation order and mg
26 
27 # global approximation order
28 -my_order 2
29 # increase of appox. at crack front (global_order + ref_order)
30 -my_ref_order 1
31 # h-refinemnt at crack tip
32 -my_ref 1
33 
34 # If one or true crack propagation algorithm on
35 -my_propagate_crack 1
36 
37 # Initial load step
38 -load_scale 1
39 
40 # Solvers propagation
41 
42 -propagation_pc_type lu
43 #-propagation_pc_type fieldsplit
44 #-propagation_pc_fieldsplit_type multiplicative
45 #-propagation_pc_fieldsplit_schur_precondition selfp
46 
47 
48 ######################################
49 ##### Internal stresses ##############
50 
51 -my_residual_stress_block 1
52 -mwls_dm 1
53 -mwls_number_of_base_functions 4
54 
55 ######################################
56 ##### Critical cutting parameters ####
57 
58 -fraction_level 1
59 
60 # Cut
61 #-tol_cut 1e-2
62 #-tol_cut_close 2e-1
63 
64 # Trim
65 #-tol_trim_close 2e-1
66 
67 -cut_mesh 1
68 -cut_surface_side_set 400
69 -edges_block_set 2
70 -vertex_block_set 3
71 -cut_factor 1
72 
73 #-cutting_surf_corner_factor 0.2
74 #-cutting_surf_skin_factor 0.2
75 #-cutting_surf_edge_delta 0.2
76 
77 ##############################
78 #### Sometimes change this ####
79 
80 # Arc-length
81 
82 -arc_beta 0
83 
84 # Crack front regularisation
85 
86 -griffith_E 1e1
87 -griffith_r 1.0
88 -gc_fix_threshold 0.5
89 
90 # Smoothing parameters
91 
92 -arc_alpha 1e2
93 -smoother_alpha 1e-2
94 -smoothing_stabilise 0
95 -adapt_min_smoother_alpha 1e-3
96 
97 # Approximation and order
98 
99 -my_geom_order 1
100 -material HOOKE
101 -my_add_singularity 1
102 -my_max_post_proc_ref_level 0
103 
104 ###################################
105 #### Almost never change this ####
106 
107 # Multi-grid
108 
109 -mofem_mg_verbose 1
110 -mofem_mg_coarse_order 1
111 -mofem_mg_levels 2
112 
113 -volume_length_type BARRIER_AND_CHANGE_QUALITY_SCALED_BY_VOLUME
114 
115 -ksp_type gmres
116 -ksp_atol 1e-12
117 -ksp_rtol 1e-12
118 -ksp_max_it 10
119 -pc_type lu
120 
121 -arc_length_ksp_type fgmres
122 -arc_length_ksp_max_it 40
123 -arc_length_ksp_gmres_restart 500
124 -arc_length_ksp_atol 1e-12
125 -arc_length_ksp_rtol 1e-6
126 #-arc_length_ksp_monitor
127 #-arc_length_pc_type lu
128 #-arc_length_pc_factor_mat_solver_package superlu_dist
129 
130 
131 -arc_length_pc_factor_shift_type POSITIVE_DEFINITE
132 -arc_length_pc_factor_reuse_ordering 0
133 -arc_length_pc_factor_mat_ordering_type nd
134 
135 
136 -elastic_snes_monitor
137 -elastic_snes_atol 1e-10
138 -elastic_snes_rtol 1e-10
139 -elastic_snes_converged_reason
140 -elastic_snes_linesearch_type basic
141 -elastic_snes_max_linear_solve_fail 100
142 -elastic_snes_max_fail 100
143 
144 -elastic_ksp_monitor
145 -elastic_ksp_type cg
146 -elastic_ksp_atol 1e-12
147 -elastic_ksp_rtol 1e-12
148 -elastic_ksp_stol 0
149 -elastic_ksp_max_it 100
150 #-elastic_pc_factor_shift_type POSITIVE_DEFINITE
151 #-elastic_pc_factor_reuse_ordering 0
152 #-elastic_pc_factor_mat_ordering_type nd
153 
154 
155 -elastic_pc_type lu
156 -elastic_pc_factor_mat_solver_package superlu_dist
157 -mat_superlu_dist_statprint
158 -elastic_mat_superlu_dist_statprint
159 
160 #-elastic_pc_type mg
161 -elastic_pc_mg_log
162 -elastic_pc_mg_type multiplicative
163 -elastic_pc_mg_smoothup 3
164 -elastic_pc_mg_smoothdown 3
165 -elastic_mg_coarse_ksp_type gmres
166 -elastic_mg_coarse_pc_type lu
167 #-elastic_ksp_converged_reason
168 
169 -material_snes_atol 1e-8
170 -material_snes_rtol 0
171 -material_snes_stol 0
172 -material_snes_converged_reason
173 -material_snes_max_linear_solve_fail 100
174 -material_snes_max_fail 100
175 -material_snes_max_it 100
176 
177 -material_snes_linesearch_type bt
178 -material_snes_linesearch_minlambda 1e-4
179 -material_snes_linesearch_damping 1
180 -material_snes_linesearch_max_it 3
181 
182 -material_ksp_type gmres
183 -material_ksp_atol 1e-12
184 -material_ksp_rtol 1e-12
185 -material_ksp_stol 0
186 -material_ksp_max_it 100
187 #-material_ksp_monitor
188 -material_pc_type lu
189 
190 #-propagation_snes_type ngmres
191 #-propagation_snes_ngmres_select_type difference
192 #-propagation_ngmres_monitor 1
193 #-propagation_snes_type test
194 
195 -propagation_snes_atol 1e-8
196 -propagation_snes_rtol 1e-6
197 -propagation_snes_stol 0
198 -propagation_snes_converged_reason
199 -propagation_snes_max_linear_solve_fail 100
200 -propagation_snes_max_fail 100
201 -propagation_snes_max_it 16
202 
203 -propagation_snes_linesearch_type bt
204 -propagation_snes_linesearch_minlambda 1e-3
205 -propagation_snes_linesearch_damping 1
206 -propagation_snes_linesearch_max_it 3
207 -propagation_snes_linesearch_monitor
208 -propagation_snes_monitor
209 
210 -propagation_ksp_monitor
211 -propagation_ksp_type fgmres
212 -propagation_ksp_atol 1e-9
213 -propagation_ksp_rtol 1e-6
214 -propagation_ksp_stol 0
215 -propagation_ksp_max_it 100
216 -propagation_ksp_gmres_restart 500
217 
218 -propagation_fieldsplit_0_ksp_type cg
219 -propagation_fieldsplit_0_ksp_atol 1e-16
220 -propagation_fieldsplit_0_ksp_rtol 0
221 -propagation_fieldsplit_0_ksp_max_it 5
222 -propagation_fieldsplit_0_pc_type lu
223 -propagation_fieldsplit_0_pc_factor_shift_type POSITIVE_DEFINITE
224 #-propagation_fieldsplit_0_pc_type bjacobi
225 #-propagation_fieldsplit_0_pc_type mg
226 #-propagation_fieldsplit_0_ksp_monitor
227 #-propagation_fieldsplit_0_ksp_converged_reason
228 -propagation_fieldsplit_0_pc_mg_type multiplicative
229 -propagation_fieldsplit_0_pc_mg_smoothup 4
230 -propagation_fieldsplit_0_pc_mg_smoothdown 4
231 -propagation_fieldsplit_0_mg_coarse_ksp_type cg
232 -propagation_fieldsplit_0_mg_coarse_pc_type lu
233 
234 -propagation_fieldsplit_1_ksp_type fgmres
235 -propagation_fieldsplit_1_ksp_atol 1e-16
236 -propagation_fieldsplit_1_ksp_rtol 0
237 -propagation_fieldsplit_1_ksp_max_it 5
238 #-propagation_fieldsplit_1_ksp_monitor
239 -propagation_fieldsplit_1_pc_type lu
240 
241 
242 #-pc_factor_mat_ordering_type natural
243 #-material_pc_factor_mat_ordering_type natural
244 #-elastic_pc_factor_mat_ordering_type natural
245 #-propagation_pc_factor_mat_ordering_type natural
246 #-propagation_fieldsplit_0_mg_coarse_pc_factor_mat_ordering_type natural
247 #-propagation_fieldsplit_0_pc_factor_mat_ordering_type natural
248 #-propagation_fieldsplit_1_pc_factor_mat_ordering_type natural
249 
250 #-material_mat_mumps_icntl_14 10000
251 #-elastic_mat_mumps_icntl_14 10000
252 #-propagation_mat_mumps_icntl_14 10000
253 #-propagation_fieldsplit_0_mat_mumps_icntl_14 10000
254 #-propagation_fieldsplit_1_mat_mumps_icntl_14 10000
255 #-mat_mumps_icntl_14 10000
256 #-mat_mumps_icntl_24 1
257 #-mat_mumps_icntl_7 4
258 #-mat_mumps_icntl_28 1
259 
260 
261 #-mat_mumps_icntl_14 200
262 #-mat_mumps_icntl_24 1
263 #-mat_mumps_icntl_13 1