PHOTON USE
p;;;;
msg temperature contours
con tem1 y 1 fi;0.01
con tem1 y 5 fi;0.01
con tem1 y 10 fi;0.01
con tem1 y 15 fi;0.01
con tem1 y 20 fi;0.01
con tem1 y 25 fi;0.01
con tem1 y 30 fi;0.01
pause
msg velocity vectors at the top
vec y m
enduse
************************************************************
Group 1. Run Title
TEXT(Free Convectn And Radiatn Heat Sink
TITLE
LIBREF=608
DISPLAY
A vertical metal plate is fitted with vertical fins, which cool
it by free convection and radiation. The flow is steady. The
grid size is 20 * 30 * 30; so the computer time is not small.
PHOTON USE commands are included in the Q1 file.
ENDDIS
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
REAL(TOL,LX,LY,LZ,X0,Y0,Z0)
INTEGER(NXREG,NYREG,NZREG)
* Overall number of cells
NX=30; NY=30; NZ=20; TOL=1.E-5
RSET(M,NX,NY,NZ,TOL)
* Overall domain extent
XULAST=0.05;YVLAST=0.2;ZWLAST=0.1
RSET(D,CHAM,XULAST,YVLAST,ZWLAST)
* Set objects: name x0 y0 z0
* lx ly lz
object 1
X0=0.0; Y0=0.05; Z0=0.0455; LX=0.035; LY=0.1; LZ=0.003
RSET(B,BASE,X0,Y0,Z0,LX,LY,LZ)
object 2
X0=0.0; Y0=0.08; Z0=0.0485; LX=0.009; LY=0.04; LZ=0.002
RSET(B,HEATB,X0,Y0,Z0,LX,LY,LZ)
object 3
X0=0.0; Y0=0.08; Z0=0.0505; LX=0.0046; LY=0.04; LZ=0.008
RSET(B,HEATT,X0,Y0,Z0,LX,LY,LZ)
object 4
X0=0.0145; Y0=0.05; Z0=0.0485; LX=0.001475; LY=0.1; LZ=0.024
RSET(B,FIN1,X0,Y0,Z0,LX,LY,LZ)
object 5
X0=0.0237; Y0=0.05; Z0=0.0485; LX=0.001475; LY=0.1; LZ=0.024
RSET(B,FIN2,X0,Y0,Z0,LX,LY,LZ)
object 6
X0=0.0329; Y0=0.05; Z0=0.0485; LX=0.001475; LY=0.1; LZ=0.024
RSET(B,FIN3,X0,Y0,Z0,LX,LY,LZ)
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
#solvel
#soltem
STORE(DEN1,BLOK)
************************************************************
Group 8. Terms & Devices
CSG3=CNGR
************************************************************
Group 9. Properties
#use_props
PRESS0 = 1.0E+05; TEMP0 = 2.730E+02
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(V1 ) = 1.0E-01 ; FIINIT(BLOK) = 1.0
:fluid: =airideal ; FIINIT(DEN1) = 1.189E+00
FIINIT(TEM1) = 2.0E+01
CONPOR(BASE , -1.00,CELL ,#1,#9,#2,#4,#2,#2)
INIT(BASE ,BLOK, 0.0, 2.0); INIT(BASE ,PRPS, 0.0, 1.0E+02)
CONPOR(HEATT , -1.00,CELL ,#1,#1,#3,#3,#4,#4)
INIT(HEATT ,BLOK, 0.0, 3.0); INIT(HEATT ,PRPS, 0.0, 1.0E+02)
CONPOR(HEATB , -1.00,CELL ,#1,#2,#3,#3,#3,#3)
INIT(HEATB ,BLOK, 0.0, 4.0); INIT(HEATB ,PRPS, 0.0, 1.0E+02)
CONPOR(FIN1 , -1.00,CELL ,#4,#4,#2,#4,#3,#5)
INIT(FIN1 ,BLOK, 0.0, 5.0); INIT(FIN1 ,PRPS, 0.0, 1.0E+02)
CONPOR(FIN2 , -1.00,CELL ,#6,#6,#2,#4,#3,#5)
INIT(FIN2 ,BLOK, 0.0, 6.0); INIT(FIN2 ,PRPS, 0.0, 1.0E+02)
CONPOR(FIN3 , -1.00,CELL ,#8,#8,#2,#4,#3,#5)
INIT(FIN3 ,BLOK, 0.0, 7.0); INIT(FIN3 ,PRPS, 0.0, 1.0E+02)
INIADD = F
************************************************************
Group 13. Boundary & Special Sources
PATCH (X2 ,EAST ,nx,nx,1,ny,1,nz,1,1)
COVAL (X2 ,P1 , 1.0, 0.0)
COVAL (X2 ,U1 , 0.0, 0.0 )
COVAL (X2 ,V1 , 0.0, 0.0 )
COVAL (X2 ,W1 , 0.0, 0.0 )
COVAL (X2 ,TEM1, 0.0, 2.0E+01)
PATCH (Y1 ,SOUTH ,1,nx,1,1,1,nz,1,1)
COVAL (Y1 ,P1 , 1.0, 0.0)
COVAL (Y1 ,U1 , 0.0, 0.0 )
COVAL (Y1 ,V1 , 0.0, 0.0 )
COVAL (Y1 ,W1 , 0.0, 0.0 )
COVAL (Y1 ,TEM1, 0.0, 2.0E+01)
PATCH (Y2 ,NORTH ,1,nx,ny,ny,1,nz,1,1)
COVAL (Y2 ,P1 , 1.0, 0.0)
COVAL (Y2 ,U1 , 0.0, 0.0 )
COVAL (Y2 ,V1 , 0.0, 0.0 )
COVAL (Y2 ,W1 , 0.0, 0.0 )
COVAL (Y2 ,TEM1, 0.0, 2.0E+01)
PATCH (Z1 ,LOW ,1,nx,1,ny,1,1,1,1)
COVAL (Z1 ,P1 , 1.0, 0.0)
COVAL (Z1 ,U1 , 0.0, 0.0 )
COVAL (Z1 ,V1 , 0.0, 0.0 )
COVAL (Z1 ,W1 , 0.0, 0.0 )
COVAL (Z1 ,TEM1, 0.0, 2.0E+01)
PATCH (Z2 ,HIGH ,1,nx,1,ny,nz,nz,1,1)
COVAL (Z2 ,P1 , 1.0, 0.0)
COVAL (Z2 ,U1 , 0.0, 0.0 )
COVAL (Z2 ,V1 , 0.0, 0.0 )
COVAL (Z2 ,W1 , 0.0, 0.0 )
COVAL (Z2 ,TEM1, 0.0, 2.0E+01)
PATCH (HOTHEATT,VOLUME,#1,#1,#3,#3,#4,#4,1,1)
COVAL (HOTHEATT,TEM1, FIXFLU , 5.989E+06)
PATCH (HOTHEATB,VOLUME,#1,#2,#3,#3,#3,#3,1,1)
COVAL (HOTHEATB,TEM1, FIXFLU , 5.989E+06)
#gravity
gravdir=2; rhoref=1.189
#densdiff
PATCH (*RADBASE,LOW ,#1,#9,#2,#4,#3,#3,1,1)
COVAL (*RADBASE,TEM1, 5.670E-09, 2.0E+01)
PATCH (*RADFIN ,EAST ,#8,#8,#2,#4,#3,#5,1,1)
COVAL (*RADFIN ,TEM1, 5.670E-09, 2.0E+01)
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 100; SELREF = T; RESFAC = 1.0E-02
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(U1 ,FALSDT, 2.0E-02); RELAX(V1 ,FALSDT, 2.0E-02)
RELAX(W1 ,FALSDT, 2.0E-02); RELAX(TEM1,FALSDT, 1.0E+02)
************************************************************
Group 18. Limits on variables or increments to them
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ; IYMON = NY-1 ; IZMON = 11
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NXPRIN = 6; NYPRIN = 6; NZPRIN = 4; NPLT = 1
************************************************************
Group 24. Dumps For Restarts