better numerical parameters for thermal solver/PETSc

works for realistic values for thermal problems
parameters copied from spectral_damage
This commit is contained in:
Martin Diehl 2021-05-05 17:08:31 +02:00
parent 18dfc9e54d
commit c8cdd7e622
4 changed files with 17 additions and 2 deletions

@ -1 +1 @@
Subproject commit ab64793bb04c506d815ebc850672ed0f2d013e67
Subproject commit 95f7faea920dd6956884e4a55f72e5d5b1ffcdc8

View File

@ -0,0 +1,7 @@
references:
- www.engineeringtoolbox.com/thermal-conductivity-metals-d_858.html
- www.engineeringtoolbox.com/specific-heat-metals-d_152.html
c_p: 0.91e3
K_11: 236.0
K_22: 236.0
K_33: 236.0

View File

@ -0,0 +1,7 @@
references:
- www.engineeringtoolbox.com/thermal-conductivity-metals-d_858.html
- www.engineeringtoolbox.com/specific-heat-metals-d_152.html
c_p: 0.49e3
K_11: 54.0
K_22: 54.0
K_33: 54.0

View File

@ -90,7 +90,8 @@ subroutine grid_thermal_spectral_init(T_0)
!--------------------------------------------------------------------------------------------------
! set default and user defined options for PETSc
call PetscOptionsInsertString(PETSC_NULL_OPTIONS,'-thermal_snes_type ngmres',ierr)
call PetscOptionsInsertString(PETSC_NULL_OPTIONS,'-thermal_snes_type newtonls -thermal_snes_mf &
&-thermal_snes_ksp_ew -thermal_ksp_type fgmres',ierr)
CHKERRQ(ierr)
call PetscOptionsInsertString(PETSC_NULL_OPTIONS,num_grid%get_asString('petsc_options',defaultVal=''),ierr)
CHKERRQ(ierr)