Nodes of different colours represent the following:
Solid arrows point from a parent (sub)module to the submodule which is descended from it. Dashed arrows point from a module being used to the module or program unit using it. Where possible, edges connecting nodes are given different colours to make them easier to distinguish in large graphs.
Nodes of different colours represent the following:
Solid arrows point from a procedure to one which it calls. Dashed arrows point from an interface to procedures which implement that interface. This could include the module procedures in a generic interface or the implementation in a submodule of an interface in a parent module. Where possible, edges connecting nodes are given different colours to make them easier to distinguish in large graphs.
Nodes of different colours represent the following:
Solid arrows point from a procedure to one which it calls. Dashed arrows point from an interface to procedures which implement that interface. This could include the module procedures in a generic interface or the implementation in a submodule of an interface in a parent module. Where possible, edges connecting nodes are given different colours to make them easier to distinguish in large graphs.
SUBROUTINE DistributeHTRadGains
! SUBROUTINE INFORMATION:
! AUTHOR Rick Strand
! DATE WRITTEN February 2001
! MODIFIED April 2010 Brent Griffith, max limit to protect surface temperature calcs
! RE-ENGINEERED na
! PURPOSE OF THIS SUBROUTINE:
! To distribute the gains from the high temperature radiant heater
! as specified in the user input file. This includes distribution
! of long wavelength radiant gains to surfaces and "people" as well
! as latent, lost, and convective portions of the total gain.
! METHODOLOGY EMPLOYED:
! We must cycle through all of the radiant systems because each
! surface could feel the effect of more than one radiant system.
! Note that the energy radiated to people is assumed to affect them
! but them it is assumed to be convected to the air. This is why
! the convective portion shown below has two parts to it.
! REFERENCES:
! na
! USE STATEMENTS:
USE DataGlobals, ONLY : NumOfZones
USE DataHeatBalance, ONLY : Zone
USE DataHeatBalFanSys, ONLY : SumConvHTRadSys, SumLatentHTRadSys, &
QHTRadSysToPerson, QHTRadSysSurf, MaxRadHeatFlux
USE DataSurfaces, ONLY : Surface, TotSurfaces
USE General, ONLY : RoundSigDigits
USE DataInterfaces, ONLY : ShowContinueError, ShowWarningError, ShowSevereError, &
ShowFatalError
IMPLICIT NONE ! Enforce explicit typing of all variables in this routine
! SUBROUTINE ARGUMENT DEFINITIONS:
! na
! SUBROUTINE PARAMETER DEFINITIONS:
REAL(r64), PARAMETER :: SmallestArea = 0.001d0 ! Smallest area in meters squared (to avoid a divide by zero)
! INTERFACE BLOCK SPECIFICATIONS
! na
! DERIVED TYPE DEFINITIONS
! na
! SUBROUTINE LOCAL VARIABLE DECLARATIONS:
INTEGER :: RadSurfNum ! Counter for surfaces receiving radiation from radiant heater
INTEGER :: RadSysNum ! Counter for the radiant systems
INTEGER :: SurfNum ! Pointer to the Surface derived type
INTEGER :: ZoneNum ! Pointer to the Zone derived type
REAL(R64) :: ThisSurfIntensity ! temporary for W/m2 term for rad on a surface
! FLOW:
! Initialize arrays
SumConvHTRadSys = 0.0D0
SumLatentHTRadSys = 0.0D0
QHTRadSysSurf = 0.0D0
QHTRadSysToPerson = 0.0D0
DO RadSysNum = 1, NumOfHighTempRadSys
ZoneNum = HighTempRadSys(RadSysNum)%ZonePtr
QHTRadSysToPerson(ZoneNum) = QHTRadSource(RadSysNum) &
*HighTempRadSys(RadSysNum)%FracRadiant &
*HighTempRadSys(RadSysNum)%FracDistribPerson
SumConvHTRadSys(ZoneNum) = SumConvHTRadSys(ZoneNum) &
+QHTRadSource(RadSysNum)*HighTempRadSys(RadSysNum)%FracConvect
SumLatentHTRadSys(ZoneNum) = SumLatentHTRadSys(ZoneNum) &
+QHTRadSource(RadSysNum)*HighTempRadSys(RadSysNum)%FracLatent
DO RadSurfNum = 1, HighTempRadSys(RadSysNum)%TotSurfToDistrib
SurfNum = HighTempRadSys(RadSysNum)%SurfacePtr(RadSurfNum)
IF (Surface(SurfNum)%Area > SmallestArea) THEN
ThisSurfIntensity = (QHTRadSource(RadSysNum) &
*HighTempRadSys(RadSysNum)%FracRadiant &
*HighTempRadSys(RadSysNum)%FracDistribToSurf(RadSurfNum) &
/Surface(SurfNum)%Area )
QHTRadSysSurf(SurfNum) = QHTRadSysSurf(SurfNum) + ThisSurfIntensity
IF (ThisSurfIntensity > MaxRadHeatFlux) THEN ! CR 8074, trap for excessive intensity (throws off surface balance )
CALL ShowSevereError('DistributeHTRadGains: excessive thermal radiation heat flux intensity detected')
CALL ShowContinueError('Surface = '//TRIM(Surface(SurfNum)%Name) )
CALL ShowContinueError('Surface area = '//TRIM(RoundSigDigits(Surface(SurfNum)%Area,3))//' [m2]')
CALL ShowContinueError('Occurs in ZoneHVAC:HighTemperatureRadiant = '//Trim(HighTempRadSys(RadSysNum)%Name))
CALL ShowContinueError('Radiation intensity = '//Trim(RoundSigDigits(ThisSurfIntensity,2))//' [W/m2]')
CALL ShowContinueError('Assign a larger surface area or more surfaces in ZoneHVAC:HighTemperatureRadiant' )
CALL ShowFatalError('DistributeHTRadGains: excessive thermal radiation heat flux intensity detected')
ENDIF
ELSE ! small surface
CALL ShowSevereError('DistributeHTRadGains: surface not large enough to receive thermal radiation heat flux')
CALL ShowContinueError('Surface = '//TRIM(Surface(SurfNum)%Name) )
CALL ShowContinueError('Surface area = '//TRIM(RoundSigDigits(Surface(SurfNum)%Area,3))//' [m2]')
CALL ShowContinueError('Occurs in ZoneHVAC:HighTemperatureRadiant = '//Trim(HighTempRadSys(RadSysNum)%Name))
CALL ShowContinueError('Assign a larger surface area or more surfaces in ZoneHVAC:HighTemperatureRadiant' )
CALL ShowFatalError('DistributeHTRadGains: surface not large enough to receive thermal radiation heat flux')
ENDIF
END DO
END DO
! Here an assumption is made regarding radiant heat transfer to people.
! While the QHTRadSysToPerson array will be used by the thermal comfort
! routines, the energy transfer to people would get lost from the perspective
! of the heat balance. So, to avoid this net loss of energy which clearly
! gets added to the zones, we must account for it somehow. This assumption
! that all energy radiated to people is converted to convective energy is
! not very precise, but at least it conserves energy.
DO ZoneNum = 1, NumOfZones
SumConvHTRadSys(ZoneNum) = SumConvHTRadSys(ZoneNum) + QHTRadSysToPerson(ZoneNum)
END DO
RETURN
END SUBROUTINE DistributeHTRadGains