thermal_adiabatic.f90

Go to the documentation of this file.

# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/thermal_adiabatic.f90"
# 1 "<built-in>"
# 1 "<command-line>"
# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/thermal_adiabatic.f90"
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
module thermal_adiabatic
  use prec
  use config
  use numerics
  use material
  use results
  use source_thermal_dissipation
  use source_thermal_externalheat
  use crystallite
  use lattice
 
  implicit none
  private
 
  type :: tparameters
    character(len=pStringLen), allocatable, dimension(:) :: &
      output
  end type tparameters
  
  type(tparameters),             dimension(:), allocatable :: &
    param
 
  public :: &
    thermal_adiabatic_init, &
    thermal_adiabatic_updatestate, &
    thermal_adiabatic_getsourceanditstangent, &
    thermal_adiabatic_getspecificheat, &
    thermal_adiabatic_getmassdensity, &
    thermal_adiabatic_results
 
contains
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine thermal_adiabatic_init
 
  integer :: maxninstance,h,nofmyhomog   
 
  write(6,'(/,a)') ' <<<+-  thermal_'//thermal_adiabatic_label//' init  -+>>>'; flush(6)
  
  maxninstance = count(thermal_type == thermal_adiabatic_id)
  if (maxninstance == 0) return
  
  allocate(param(maxninstance))
  
  do h = 1, size(thermal_type)
    if (thermal_type(h) /= thermal_adiabatic_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 = 1
    allocate(thermalstate(h)%state0   (1,nofmyhomog), source=thermal_initialt(h))
    allocate(thermalstate(h)%subState0(1,nofmyhomog), source=thermal_initialt(h))
    allocate(thermalstate(h)%state    (1,nofmyhomog), source=thermal_initialt(h))
 
    thermalmapping(h)%p => material_homogenizationmemberat
    deallocate(temperature(h)%p)
    temperature(h)%p => thermalstate(h)%state(1,:)
    deallocate(temperaturerate(h)%p)
    allocate  (temperaturerate(h)%p(nofmyhomog), source=0.0_preal)
  
    end associate
  enddo
 
end subroutine thermal_adiabatic_init
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
function thermal_adiabatic_updatestate(subdt, ip, el)
  
  integer,     intent(in) :: &
    ip, &                                                                                           !< integration point number
    el
  real(preal), intent(in) :: &
    subdt
 
  logical, dimension(2)  :: &
    thermal_adiabatic_updatestate
  integer :: &
    homog, &
    offset
  real(preal) :: &
    t, tdot, dtdot_dt  
 
  homog  = material_homogenizationat(el)
  offset = material_homogenizationmemberat(ip,el)
  
  t = thermalstate(homog)%subState0(1,offset)
  call thermal_adiabatic_getsourceanditstangent(tdot, dtdot_dt, t, ip, el)
  t = t + subdt*tdot/(thermal_adiabatic_getspecificheat(ip,el)*thermal_adiabatic_getmassdensity(ip,el))
  
  thermal_adiabatic_updatestate = [     abs(t - thermalstate(homog)%state(1,offset)) &
                                     <= err_thermal_tolabs &
                                   .or. abs(t - thermalstate(homog)%state(1,offset)) &
                                     <= err_thermal_tolrel*abs(thermalstate(homog)%state(1,offset)), &
                                   .true.]
 
  temperature(homog)%p(thermalmapping(homog)%p(ip,el)) = t  
  temperaturerate(homog)%p(thermalmapping(homog)%p(ip,el)) = &
    (thermalstate(homog)%state(1,offset) - thermalstate(homog)%subState0(1,offset))/(subdt+tiny(0.0_preal))
  
end function thermal_adiabatic_updatestate
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine thermal_adiabatic_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, &
    instance, &
    grain, &
    source, &
    constituent
    
  homog  = material_homogenizationat(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_adiabatic_getsourceanditstangent
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
function thermal_adiabatic_getspecificheat(ip,el)
 
  integer, intent(in) :: &
    ip, &                                                                                           !< integration point number
    el
  
  real(preal) :: &
    thermal_adiabatic_getspecificheat
  integer :: &
    grain
   
  thermal_adiabatic_getspecificheat = 0.0_preal
  
  do grain = 1, homogenization_ngrains(material_homogenizationat(el))
    thermal_adiabatic_getspecificheat = thermal_adiabatic_getspecificheat & 
                                      + lattice_specificheat(material_phaseat(grain,el))
  enddo
 
  thermal_adiabatic_getspecificheat = thermal_adiabatic_getspecificheat &
                                    / real(homogenization_ngrains(material_homogenizationat(el)),preal)
  
end function thermal_adiabatic_getspecificheat
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
function thermal_adiabatic_getmassdensity(ip,el)
    
  integer, intent(in) :: &
    ip, &                                                                                           !< integration point number
    el
  real(preal) :: &
    thermal_adiabatic_getmassdensity
  integer :: &
    grain
   
  thermal_adiabatic_getmassdensity = 0.0_preal
 
  do grain = 1, homogenization_ngrains(material_homogenizationat(el))
    thermal_adiabatic_getmassdensity = thermal_adiabatic_getmassdensity &
                                     + lattice_massdensity(material_phaseat(grain,el))
  enddo
 
  thermal_adiabatic_getmassdensity = thermal_adiabatic_getmassdensity &
                                   / real(homogenization_ngrains(material_homogenizationat(el)),preal)
 
end function thermal_adiabatic_getmassdensity
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine thermal_adiabatic_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_adiabatic_results
 
end module thermal_adiabatic