source_thermal_externalheat.f90

Go to the documentation of this file.

# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/source_thermal_externalheat.f90"
# 1 "<built-in>"
# 1 "<command-line>"
# 1 "/home/damask_user/GitLabCI_Pipeline_4301/DAMASK/src/source_thermal_externalheat.f90"
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
module source_thermal_externalheat
  use prec
  use debug
  use discretization
  use material
  use config
 
  implicit none
  private
 
  integer,           dimension(:),   allocatable :: &
    source_thermal_externalheat_offset, &                                                           !< which source is my current thermal dissipation mechanism?
    source_thermal_externalheat_instance
 
  type :: tparameters                                                                              
    real(preal), dimension(:), allocatable :: &
      time, &
      heat_rate
    integer :: &
     nintervals
  end type tparameters
 
  type(tparameters), dimension(:), allocatable  :: param
 
 
  public :: &
    source_thermal_externalheat_init, &
    source_thermal_externalheat_dotstate, &
    source_thermal_externalheat_getrateanditstangent
 
contains
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine source_thermal_externalheat_init
 
  integer :: ninstance,sourceoffset,nipcmyphase,p
 
  write(6,'(/,a)') ' <<<+-  source_'//source_thermal_externalheat_label//' init  -+>>>'; flush(6)
 
  ninstance = count(phase_source == source_thermal_externalheat_id)
  if (iand(debug_level(debug_constitutive),debug_levelbasic) /= 0) &
    write(6,'(a16,1x,i5,/)') '# instances:',ninstance
 
  allocate(source_thermal_externalheat_offset(size(config_phase)), source=0)
  allocate(source_thermal_externalheat_instance(size(config_phase)), source=0)
  allocate(param(ninstance))
 
  do p = 1, size(config_phase)
    source_thermal_externalheat_instance(p) = count(phase_source(:,1:p) == source_thermal_externalheat_id)
    do sourceoffset = 1, phase_nsources(p)
      if (phase_source(sourceoffset,p) == source_thermal_externalheat_id) then
        source_thermal_externalheat_offset(p) = sourceoffset
        exit
      endif
    enddo
 
    if (all(phase_source(:,p) /= source_thermal_externalheat_id)) cycle
    associate(prm => param(source_thermal_externalheat_instance(p)), &
              config => config_phase(p))
 
    prm%time       = config%getFloats('externalheat_time')
    prm%nIntervals = size(prm%time) - 1
 
    prm%heat_rate = config%getFloats('externalheat_rate',requiredsize = size(prm%time))
 
    nipcmyphase = count(material_phaseat==p) * discretization_nip
    call material_allocatesourcestate(p,sourceoffset,nipcmyphase,1,1,0)
 
    end associate
  enddo
 
end subroutine source_thermal_externalheat_init
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine source_thermal_externalheat_dotstate(phase, of)
 
  integer, intent(in) :: &
    phase, &
    of
 
  integer :: &
    sourceoffset
 
  sourceoffset = source_thermal_externalheat_offset(phase)
 
  sourcestate(phase)%p(sourceoffset)%dotState(1,of) = 1.0_preal                                     ! state is current time
 
end subroutine source_thermal_externalheat_dotstate
 
 
!--------------------------------------------------------------------------------------------------
!--------------------------------------------------------------------------------------------------
subroutine source_thermal_externalheat_getrateanditstangent(TDot, dTDot_dT, phase, of)
 
  integer, intent(in) :: &
    phase, &
    of
  real(preal),  intent(out) :: &
    tdot, &
    dtdot_dt
 
  integer :: &
    sourceoffset, interval
  real(preal) :: &
    frac_time
 
  sourceoffset = source_thermal_externalheat_offset(phase)
 
  associate(prm => param(source_thermal_externalheat_instance(phase)))
  do interval = 1, prm%nIntervals                                                                   ! scan through all rate segments
    frac_time = (sourcestate(phase)%p(sourceoffset)%state(1,of) - prm%time(interval)) &
              / (prm%time(interval+1) - prm%time(interval))                                         ! fractional time within segment
    if (     (frac_time <  0.0_preal .and. interval == 1) &
        .or. (frac_time >= 1.0_preal .and. interval == prm%nIntervals) &
        .or. (frac_time >= 0.0_preal .and. frac_time < 1.0_preal) ) &
      tdot = prm%heat_rate(interval  ) * (1.0_preal - frac_time) + &
             prm%heat_rate(interval+1) * frac_time                                                  ! interpolate heat rate between segment boundaries...
                                                                                                    ! ...or extrapolate if outside of bounds
  enddo
  dtdot_dt = 0.0
  end associate
 
end subroutine source_thermal_externalheat_getrateanditstangent
 
end module source_thermal_externalheat