thermal_conduction.f90

Go to the documentation of this file.

# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/thermal_conduction.f90"
# 1 "<built-in>"
# 1 "<command-line>"
# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/thermal_conduction.f90"
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
module thermal_conduction
  use prec
  use material
  use config
  use lattice
  use results
  use crystallite
  use source_thermal_dissipation
  use source_thermal_externalheat
 
  implicit none
  private
 
  type :: tparameters
    character(len=pStringLen), allocatable, dimension(:) :: &
      output
  end type tparameters
 
  type(tparameters),             dimension(:), allocatable :: &
    param
 
  public :: &
    thermal_conduction_init, &
    thermal_conduction_getsourceanditstangent, &
    thermal_conduction_getconductivity, &
    thermal_conduction_getspecificheat, &
    thermal_conduction_getmassdensity, &
    thermal_conduction_puttemperatureanditsrate, &
    thermal_conduction_results
 
contains
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine thermal_conduction_init
 
  integer :: ninstance,nofmyhomog,h
 
  write(6,'(/,a)') ' <<<+-  thermal_'//thermal_conduction_label//' init  -+>>>'; flush(6)
 
  ninstance = count(thermal_type == thermal_conduction_id)
  allocate(param(ninstance))
 
  do h = 1, size(config_homogenization)
    if (thermal_type(h) /= thermal_conduction_id) cycle
    associate(prm => param(thermal_typeinstance(h)),config => config_homogenization(h))
 
    prm%output = config%getStrings('(output)',defaultval=emptystringarray)
 
    nofmyhomog=count(material_homogenizationat==h)
    thermalstate(h)%sizeState = 0
    allocate(thermalstate(h)%state0   (0,nofmyhomog))
    allocate(thermalstate(h)%subState0(0,nofmyhomog))
    allocate(thermalstate(h)%state    (0,nofmyhomog))
 
    thermalmapping(h)%p => material_homogenizationmemberat
    deallocate(temperature(h)%p)
    allocate  (temperature(h)%p(nofmyhomog), source=thermal_initialt(h))
    deallocate(temperaturerate(h)%p)
    allocate  (temperaturerate(h)%p(nofmyhomog), source=0.0_preal)
 
    end associate
  enddo
 
end subroutine thermal_conduction_init
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine thermal_conduction_getsourceanditstangent(Tdot, dTdot_dT, T, ip, el)
 
  integer, intent(in) :: &
    ip, &                                                                                           !< integration point number
    el
  real(preal), intent(in) :: &
    t
  real(preal), intent(out) :: &
    tdot, dtdot_dt
  real(preal) :: &
    my_tdot, my_dtdot_dt
  integer :: &
    phase, &
    homog, &
    offset, &
    instance, &
    grain, &
    source, &
    constituent
 
  homog  = material_homogenizationat(el)
  offset = material_homogenizationmemberat(ip,el)
  instance = thermal_typeinstance(homog)
 
  tdot = 0.0_preal
  dtdot_dt = 0.0_preal
  do grain = 1, homogenization_ngrains(homog)
    phase = material_phaseat(grain,el)
    constituent = material_phasememberat(grain,ip,el)
    do source = 1, phase_nsources(phase)
      select case(phase_source(source,phase))
        case (source_thermal_dissipation_id)
          call source_thermal_dissipation_getrateanditstangent(my_tdot, my_dtdot_dt, &
                                                               crystallite_s(1:3,1:3,grain,ip,el), &
                                                               crystallite_lp(1:3,1:3,grain,ip,el), &
                                                               phase)
 
        case (source_thermal_externalheat_id)
          call source_thermal_externalheat_getrateanditstangent(my_tdot, my_dtdot_dt, &
                                                                phase, constituent)
        case default
         my_tdot = 0.0_preal
         my_dtdot_dt = 0.0_preal
 
      end select
      tdot = tdot + my_tdot
      dtdot_dt = dtdot_dt + my_dtdot_dt
    enddo
  enddo
 
  tdot = tdot/real(homogenization_ngrains(homog),preal)
  dtdot_dt = dtdot_dt/real(homogenization_ngrains(homog),preal)
 
end subroutine thermal_conduction_getsourceanditstangent
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
function thermal_conduction_getconductivity(ip,el)
 
  integer, intent(in) :: &
    ip, &                                                                                           !< integration point number
    el
  real(preal), dimension(3,3) :: &
    thermal_conduction_getconductivity
  integer :: &
    grain
 
 
  thermal_conduction_getconductivity = 0.0_preal
  do grain = 1, homogenization_ngrains(material_homogenizationat(el))
    thermal_conduction_getconductivity = thermal_conduction_getconductivity + &
     crystallite_push33toref(grain,ip,el,lattice_thermalconductivity(:,:,material_phaseat(grain,el)))
  enddo
 
  thermal_conduction_getconductivity = thermal_conduction_getconductivity &
                                     / real(homogenization_ngrains(material_homogenizationat(el)),preal)
 
end function thermal_conduction_getconductivity
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
function thermal_conduction_getspecificheat(ip,el)
 
  integer, intent(in) :: &
    ip, &                                                                                           !< integration point number
    el
  real(preal) :: &
    thermal_conduction_getspecificheat
  integer :: &
    grain
 
  thermal_conduction_getspecificheat = 0.0_preal
 
  do grain = 1, homogenization_ngrains(material_homogenizationat(el))
    thermal_conduction_getspecificheat = thermal_conduction_getspecificheat &
                                       + lattice_specificheat(material_phaseat(grain,el))
  enddo
 
  thermal_conduction_getspecificheat = thermal_conduction_getspecificheat &
                                     / real(homogenization_ngrains(material_homogenizationat(el)),preal)
 
end function thermal_conduction_getspecificheat
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
function thermal_conduction_getmassdensity(ip,el)
 
  integer, intent(in) :: &
    ip, &                                                                                           !< integration point number
    el
  real(preal) :: &
    thermal_conduction_getmassdensity
  integer :: &
    grain
 
  thermal_conduction_getmassdensity = 0.0_preal
 
 
  do grain = 1, homogenization_ngrains(material_homogenizationat(el))
    thermal_conduction_getmassdensity = thermal_conduction_getmassdensity &
                                      + lattice_massdensity(material_phaseat(grain,el))
  enddo
 
  thermal_conduction_getmassdensity = thermal_conduction_getmassdensity &
                                    / real(homogenization_ngrains(material_homogenizationat(el)),preal)
 
end function thermal_conduction_getmassdensity
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine thermal_conduction_puttemperatureanditsrate(T,Tdot,ip,el)
 
  integer, intent(in) :: &
    ip, &                                                                                           !< integration point number
    el
  real(preal),   intent(in) :: &
    t, &
    tdot
  integer :: &
    homog, &
    offset
 
  homog  = material_homogenizationat(el)
  offset = thermalmapping(homog)%p(ip,el)
  temperature(homog)%p(offset) = t
  temperaturerate(homog)%p(offset) = tdot
 
end subroutine thermal_conduction_puttemperatureanditsrate
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine thermal_conduction_results(homog,group)
 
  integer,          intent(in) :: homog
  character(len=*), intent(in) :: group
 
  integer :: o
 
  associate(prm => param(damage_typeinstance(homog)))
  outputsloop: do o = 1,size(prm%output)
    select case(trim(prm%output(o)))
      case('temperature')                                                                           ! ToDo: should be 'T'
        call results_writedataset(group,temperature(homog)%p,'T',&
                                  'temperature','K')
    end select
  enddo outputsloop
  end associate
 
end subroutine thermal_conduction_results
 
end module thermal_conduction