param_file.petsc

Go to the documentation of this file.

###############################
#### Essential parameters #####
# set Griffith's energy
# Pay attention to the size
-my_gc 0.23
# arc length step
# Impact by the size
-arc_s 5
## arc-length step adaptation
# desired number of iterations
-adapt_step_its_d 6
# minimal step size
-adapt_step_min_s 1
# maximal step size
-adapt_step_max_s 20
## Steps and cutting
# number of cutting steps
-nb_cut_steps 500
# number of load steps between cutting
-nb_load_steps 3
## Approximation order and mg
# global approximation order
-my_order 1
# increase of appox. at crack front (global_order + ref_order)
-my_ref_order 2
# If one or true crack propagation algorithm on
-my_propagate_crack 1
# Initial load step
-load_scale 1
# Contact parameters
-my_contact_order 2              # default value 2
-my_lambda_contact_order 1       # default value 1 (not to be changed)
-my_cn_value 1.e4                # default value 1.0
-my_r_value 1.0                  # default value 1.0 (not to be changed)
-my_cracked_body_block_set 1001  # Block ID of the fractured solid
# Logging parameters
-log_no_color 1
# Solvers propagation
-mofem_mg_verbose 1
-mofem_mg_coarse_order 1
-mofem_mg_levels 3
-arc_length_ksp_type fgmres
-arc_length_ksp_max_it 20
-arc_length_ksp_gmres_restart 500
-arc_length_ksp_atol 1e-10
-arc_length_ksp_rtol 0
-elastic_pc_type lu
-propagation_pc_type lu
-propagation_fieldsplit_0_ksp_type cg
-propagation_fieldsplit_0_ksp_atol 1e-16
-propagation_fieldsplit_0_ksp_rtol 0
-propagation_fieldsplit_0_ksp_max_it 5
-propagation_fieldsplit_0_pc_type lu
-propagation_fieldsplit_0_pc_mg_type multiplicative
-propagation_fieldsplit_0_pc_mg_smoothup 1
-propagation_fieldsplit_0_pc_mg_smoothdown 1
-propagation_fieldsplit_0_mg_coarse_ksp_type cg
-propagation_fieldsplit_0_mg_coarse_pc_type lu
-propagation_fieldsplit_1_ksp_type fgmres
-propagation_fieldsplit_1_ksp_atol 1e-16
-propagation_fieldsplit_1_ksp_rtol 0
-propagation_fieldsplit_1_ksp_max_it 3
-propagation_fieldsplit_1_pc_type lu
######################################
##### Critical cutting parameters ####
-fraction_level 2
# Cut
#-tol_cut 4e-2
#-tol_cut_close 5e-1
# Trim
#-tol_trim 4e-1
#-tol_trim_close 4e-2
-cut_mesh 1
-cut_surface_side_set 200
-edges_block_set 900
-vertex_block_set 300
-body_skin 3000
-crack_side_set 4000
-front_side_set 4001
-cut_factor 0
-cutting_end_node_edge_delta_factor 0.2
-cutting_surf_skin_factor 0.5
-cutting_snap_to_fixed_edge_rtol 0.1
-ref_before_cut 1
-ref_before_trim 1
-griffith_force_push 0  # True/False
##############################
#### Sometimes change this ####
# Arc-length
-arc_beta 0.0
# Crack front regularisation
-griffith_E 1
-griffith_r 1.01
-gc_fix_threshold 0.5
# Smoothing parameters
-arc_alpha 1
-smoother_alpha 1
-smoothing_stabilise 1
-adapt_desired_alpha 1e-1
-adapt_min_smoother_alpha 1e-6
# Approximation and order
-my_ref 0
-my_geom_order 1
-material HOOKE
-my_add_singularity 0
-my_max_post_proc_ref_level 0
-post_proc_level 4
###################################
#### Almost never change this  ####
-volume_length_type BARRIER_AND_CHANGE_QUALITY_SCALED_BY_VOLUME
-ksp_type gmres
-ksp_atol 1e-12
-ksp_rtol 1e-12
-ksp_max_it 10
-pc_type lu
-elastic_snes_monitor
-elastic_snes_atol 1e-10
-elastic_snes_rtol 1e-10
-elastic_snes_converged_reason
-elastic_snes_linesearch_type basic
-elastic_snes_max_linear_solve_fail 100
-elastic_snes_max_fail 100
-elastic_ksp_monitor
-elastic_ksp_type gmres
-elastic_ksp_atol 1e-12
-elastic_ksp_rtol 1e-12
-elastic_ksp_stol 0
-elastic_ksp_max_it 100
-elastic_pc_mg_log
-elastic_pc_mg_type multiplicative
-elastic_pc_mg_smoothup 1
-elastic_pc_mg_smoothdown 1
-elastic_mg_coarse_ksp_type gmres
-elastic_mg_coarse_pc_type lu
-material_snes_atol 1e-8
-material_snes_rtol 0
-material_snes_stol 0
-material_snes_converged_reason
-material_snes_max_linear_solve_fail 100
-material_snes_max_fail 100
-material_snes_max_it 100
-material_snes_linesearch_type bt
-material_snes_linesearch_minlambda 1e-6
-material_snes_linesearch_damping 1
-material_snes_linesearch_max_it 6
-material_ksp_type gmres
-material_ksp_atol 1e-12
-material_ksp_rtol 1e-12
-material_ksp_max_it 100
-material_pc_type lu
-propagation_snes_atol 1e-4
-propagation_snes_rtol 1e-4
-propagation_snes_stol 0
-propagation_snes_converged_reason
-propagation_snes_max_linear_solve_fail 100
-propagation_snes_max_fail 100
-propagation_snes_max_it 20
-propagation_snes_linesearch_type l2
-propagation_snes_linesearch_minlambda 1e-3
-propagation_snes_linesearch_damping 1
-propagation_snes_linesearch_max_it 1
-propagation_snes_linesearch_monitor
-propagation_snes_monitor
-propagation_ksp_monitor
-propagation_ksp_type fgmres
-propagation_ksp_atol 1e-10
-propagation_ksp_rtol 1e-6
-propagation_ksp_max_it 100
-propagation_ksp_gmres_restart 500
-pc_factor_mat_ordering_type natural
-material_pc_factor_mat_ordering_type natural
-elastic_pc_factor_mat_ordering_type natural
-propagation_pc_factor_mat_ordering_type natural
-propagation_fieldsplit_0_mg_coarse_pc_factor_mat_ordering_type natural
-propagation_fieldsplit_0_pc_factor_mat_ordering_type natural
-propagation_fieldsplit_1_pc_factor_mat_ordering_type natural