TALK=T;RUN(1,1)
DISPLAY
This case modeling a heat exchange by heat conductivity about
a filament. The filament is simulated by the cylinder shape with
properties of steel. It is surrounded by air. The air temperature
has constant significance 18 degree at all domain boundaries.
The filament is heated up by a electrical current which occurs
in the conductor at presence the filament votage.
The emitted heat is calculated as the filament votage divided
on thermal resistance which is function of the conductor
temperature.
Source = Voltage^2 / Resistance
where
Resistance = 0.1 * Temperature
The heating source for the HEATING object is set by In-Form as
per unit volume
(SOURCE of TEM1 at HEATING is 220^2/(0.1*TEM1) with VOLU)
The Q1 contains PHOTON USE commands
ENDDIS
PHOTON USE
p;;;
con tem1 z 1 fi;0.0001
enduse
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Heat exchange about a steel conductor )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 55
XULAST =0.1
XFRAC(1)=0.018 ;XFRAC(2)=0.036
XFRAC(3)=0.054 ;XFRAC(4)=0.072
XFRAC(5)=0.09 ;XFRAC(6)=0.108
XFRAC(7)=0.126 ;XFRAC(8)=0.144
XFRAC(9)=0.162 ;XFRAC(10)=0.18
XFRAC(11)=0.198 ;XFRAC(12)=0.216
XFRAC(13)=0.234 ;XFRAC(14)=0.252
XFRAC(15)=0.27 ;XFRAC(16)=0.288
XFRAC(17)=0.306 ;XFRAC(18)=0.324
XFRAC(19)=0.342 ;XFRAC(20)=0.36
XFRAC(21)=0.378 ;XFRAC(22)=0.396
XFRAC(23)=0.414 ;XFRAC(24)=0.432
XFRAC(25)=0.45 ;XFRAC(26)=0.466667
XFRAC(27)=0.483333 ;XFRAC(28)=0.5
XFRAC(29)=0.516667 ;XFRAC(30)=0.533333
XFRAC(31)=0.55 ;XFRAC(32)=0.56875
XFRAC(33)=0.5875 ;XFRAC(34)=0.60625
XFRAC(35)=0.625 ;XFRAC(36)=0.64375
XFRAC(37)=0.6625 ;XFRAC(38)=0.68125
XFRAC(39)=0.7 ;XFRAC(40)=0.71875
XFRAC(41)=0.7375 ;XFRAC(42)=0.75625
XFRAC(43)=0.775 ;XFRAC(44)=0.79375
XFRAC(45)=0.8125 ;XFRAC(46)=0.83125
XFRAC(47)=0.85 ;XFRAC(48)=0.86875
XFRAC(49)=0.8875 ;XFRAC(50)=0.90625
XFRAC(51)=0.925 ;XFRAC(52)=0.94375
XFRAC(53)=0.9625 ;XFRAC(54)=0.98125
XFRAC(55)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 55
YVLAST =0.1
YFRAC(1)=0.018 ;YFRAC(2)=0.036
YFRAC(3)=0.054 ;YFRAC(4)=0.072
YFRAC(5)=0.09 ;YFRAC(6)=0.108
YFRAC(7)=0.126 ;YFRAC(8)=0.144
YFRAC(9)=0.162 ;YFRAC(10)=0.18
YFRAC(11)=0.198 ;YFRAC(12)=0.216
YFRAC(13)=0.234 ;YFRAC(14)=0.252
YFRAC(15)=0.27 ;YFRAC(16)=0.288
YFRAC(17)=0.306 ;YFRAC(18)=0.324
YFRAC(19)=0.342 ;YFRAC(20)=0.36
YFRAC(21)=0.378 ;YFRAC(22)=0.396
YFRAC(23)=0.414 ;YFRAC(24)=0.432
YFRAC(25)=0.45 ;YFRAC(26)=0.466667
YFRAC(27)=0.483333 ;YFRAC(28)=0.5
YFRAC(29)=0.516667 ;YFRAC(30)=0.533333
YFRAC(31)=0.55 ;YFRAC(32)=0.56875
YFRAC(33)=0.5875 ;YFRAC(34)=0.60625
YFRAC(35)=0.625 ;YFRAC(36)=0.64375
YFRAC(37)=0.6625 ;YFRAC(38)=0.68125
YFRAC(39)=0.7 ;YFRAC(40)=0.71875
YFRAC(41)=0.7375 ;YFRAC(42)=0.75625
YFRAC(43)=0.775 ;YFRAC(44)=0.79375
YFRAC(45)=0.8125 ;YFRAC(46)=0.83125
YFRAC(47)=0.85 ;YFRAC(48)=0.86875
YFRAC(49)=0.8875 ;YFRAC(50)=0.90625
YFRAC(51)=0.925 ;YFRAC(52)=0.94375
YFRAC(53)=0.9625 ;YFRAC(54)=0.98125
YFRAC(55)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST =0.1
ZFRAC(1)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(149)=PRPS ;NAME(150)=TEM1
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(PRPS,Y,N,N,N,N,Y)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
PRPS = 149
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
TERMS(TEM1,N,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 =1.189 ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =273. ;PRESS0 =1.0E+05
DVO1DT =3.41E-03 ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =1.544E-05 ;ENUT =0.
PRNDTL(TEM1)=-0.0258
PRT(TEM1)=1.
CP1 =1005. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(PRPS)=-1. ;FIINIT(TEM1)=1.0E-10
PATCH(^OB1 ,INIVAL, 1, 0, 0, 0, 0, 0, 1, 1)
INIT(^OB1 ,PRPS,0. ,111. )
INIT(^OB1 ,TEM1,0. ,18. )
INIADD = F
FSWEEP = 1
NAMFI =CHAM
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources
PATCH(^OD1 ,CELL , 1, 18, 0, 0, 0, 0, 1, 1)
COVAL(^OD1 ,TEM1,In-Form:source - see Grp 19)
PATCH(OB2 ,WWALL , 1, 1, 1, 55, 1, 1, 1, 1)
COVAL(OB2 ,TEM1,1. ,18. )
PATCH(OB3 ,EWALL , 55, 55, 1, 55, 1, 1, 1, 1)
COVAL(OB3 ,TEM1,1. ,18. )
PATCH(OB4 ,SWALL , 1, 55, 1, 1, 1, 1, 1, 1)
COVAL(OB4 ,TEM1,1. ,18. )
PATCH(OB5 ,NWALL , 1, 55, 55, 55, 1, 1, 1, 1)
COVAL(OB5 ,TEM1,1. ,18. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 500 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
LITER(TEM1)=500
ENDIT(TEM1)=1.0E-20
************************************************************
Group 17. Relaxation
RELAX(PRPS,LINRLX,1.)
RELAX(TEM1,FALSDT,1.0E+09)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
VARMAX(TEM1)=3000. ;VARMIN(TEM1)=-204.75
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
ISG62 = 1
SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,0)
SPEDAT(SET,SOURCE,TEM1!HEATING,C,=220^2/(0.1*TEM1)!VOLU)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,OBJNAM,^OB1,C,HEATING)
SPEDAT(SET,OBJTYP,^OB1,C,BLOCKAGE)
SPEDAT(SET,HEATING,MATERIAL,R,111.)
SPEDAT(SET,OBJNAM,^OD1,C,HEATING)
SPEDAT(SET,OBJTYP,^OD1,C,BLOCKAGE)
SPEDAT(SET,OBJNAM,!OB2,C,W_WALL)
SPEDAT(SET,OBJTYP,!OB2,C,PLATE)
SPEDAT(SET,OBJNAM,!OB3,C,E_WALL)
SPEDAT(SET,OBJTYP,!OB3,C,PLATE)
SPEDAT(SET,OBJNAM,!OB4,C,S_WALL)
SPEDAT(SET,OBJTYP,!OB4,C,PLATE)
SPEDAT(SET,OBJNAM,!OB5,C,N_WALL)
SPEDAT(SET,OBJTYP,!OB5,C,PLATE)
SPEDAT(SET,FACETDAT,NUMOBJ,I,5)
SPEDAT(SET,MATERIAL,111,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(PRPS)=1.582 ;EX(TEM1)=20.93
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
OUTPUT(PRPS,Y,N,Y,N,N,N)
OUTPUT(TEM1,Y,N,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 41 ;IYMON = 30 ;IZMON = 1
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NXPRIN = 1 ;IXPRF = 25 ;IXPRL = 32
NYPRIN = 1 ;IYPRF = 25 ;IYPRL = 32
IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP