This subroutine is used for nonuniform convective heat flux, characterized by distributed load labels of the form FxNUy, cf *FILM. The load label can be up to 20 characters long. In particular, y can be used to distinguish different nonuniform film patterns. The header and variable description is as follows:
subroutine film(h,sink,temp,kstep,kinc,time,noel,npt,
& coords,jltyp,field,nfield,loadtype,node,area,vold,mi,
& ipkon,kon,lakon,iponoel,inoel,ielprop,prop,ielmat,
& shcon,nshcon,rhcon,nrhcon,ntmat_,cocon,ncocon)
!
! user subroutine film
!
!
! INPUT:
!
! sink most recent sink temperature
! temp current temperature value
! kstep step number
! kinc increment number
! time(1) current step time
! time(2) current total time
! noel element number
! npt integration point number
! coords(1..3) global coordinates of the integration point
! jltyp loading face kode:
! 11 = face 1
! 12 = face 2
! 13 = face 3
! 14 = face 4
! 15 = face 5
! 16 = face 6
! field currently not used
! nfield currently not used (value = 1)
! loadtype load type label
! node network node (only for forced convection)
! area area covered by the integration point
! vold(0..4,1..nk) solution field in all nodes;
! for structural nodes:
! 0: temperature
! 1: displacement in global x-direction
! 2: displacement in global y-direction
! 3: displacement in global z-direction
! 4: static pressure
! for network nodes
! 0: total temperature (at end nodes)
! = static temperature for liquids
! 1: mass flow (at middle nodes)
! 2: total pressure (at end nodes)
! = static pressure for liquids
! 3: static temperature (at end nodes; only for gas)
! mi(1) max # of integration points per element (max
! over all elements)
! mi(2) max degree of freedom per node (max over all
! nodes) in fields like v(0:mi(2))...
! ipkon(i) points to the location in field kon preceding
! the topology of element i
! kon(*) contains the topology of all elements. The
! topology of element i starts at kon(ipkon(i)+1)
! and continues until all nodes are covered. The
! number of nodes depends on the element label
! lakon(i) contains the label of element i
! iponoel(i) the network elements to which node i belongs
! are stored in inoel(1,iponoel(i)),
! inoel(1,inoel(2,iponoel(i)))...... until
! inoel(2,inoel(2,inoel(2......)=0
! inoel(1..2,*) field containing the network elements
! ielprop(i) points to the location in field prop preceding
! the properties of element i
! prop(*) contains the properties of all network elements. The
! properties of element i start at prop(ielprop(i)+1)
! and continues until all properties are covered. The
! appropriate amount of properties depends on the
! element label. The kind of properties, their
! number and their order corresponds
! to the description in the user's manual,
! cf. the sections "Fluid Section Types"
! ielmat(i) contains the material number for element i
! shcon(0,j,i) temperature at temperature point j of material i
! shcon(1,j,i) specific heat at constant pressure at the
! temperature point j of material i
! shcon(2,j,i) dynamic viscosity at the temperature point j of
! material i
! shcon(3,1,i) specific gas constant of material i
! nshcon(i) number of temperature data points for the specific
! heat of material i
! rhcon(0,j,i) temperature at density temperature point j of
! material i
! rhcon(1,j,i) density at the density temperature point j of
! material i
! nrhcon(i) number of temperature data points for the density
! of material i
! ntmat_ maximum number of temperature data points for
! any material property for any material
! ncocon(1,i) number of conductivity constants for material i
! ncocon(2,i) number of temperature data points for the
! conductivity coefficients of material i
! cocon(0,j,i) temperature at conductivity temperature point
! j of material i
! cocon(k,j,i) conductivity coefficient k at conductivity
! temperature point j of material i
!
! OUTPUT:
!
! h(1) magnitude of the film coefficient
! h(2) not used; please do NOT assign any value
! sink (updated) sink temperature (need not be
! defined for forced convection)
! ntmat_ maximum number of temperature data points for
! any material property for any material
! ncocon(1,i) number of conductivity constants for material i
! ncocon(2,i) number of temperature data points for the
! conductivity coefficients of material i
! cocon(0,j,i) temperature at conductivity temperature point
! j of material i
! cocon(k,j,i) conductivity coefficient k at conductivity
! temperature point j of material i
!
! OUTPUT:
!
! h(1) magnitude of the film coefficient
! h(2) not used; please do NOT assign any value
! sink (updated) sink temperature (need not be
! defined for forced convection)
! heatnod extra heat flow going to the network node
! (zero if not specified)
! heatfac extra heat flow going to the structural face
! (zero if not specified)