talk=f;run(1,1)
DISPLAY
Library case H105:
EXPERIMENTAL BOX
This is an example of HotBox set-up using BFC grid. It is
a box with two slanted plates. Forced flow is generated
by the two fans, one blowing-in and one extracting air.
Two PCBs generate heat. Turbulence is modelled using LVEL
turbulence model.
Please note that in order to run this case in-core the
F-array dimension in your GROUND.F file must be set to
1.000.000.
(Use PHOTON to display results. Enter USE Q1 in PHOTON
command mode, or create a file called U, which contains
the one line: USE Q1 before typing RUNPHO.)
ENDDIS
************************************************************
PHOTON USE
p;;;;
use patgeo
msg Geometry of the problem
msg
use hotspot
msg Press return to plot temperature contours
pause
cont tem1 x 12 fil;.001
msg
msg Press return to plot velocity vectors
pause
vec x 6 sh
msg
msg Type E to end
ENDUSE
************************************************************
Group 1. Run Title
TEXT(Experimental Box
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,12,28,18)
* Set overall domain extent:
* xulast yvlast zwlast name
XSI= 1.300E-01;YSI= 6.400E-01;ZSI= 4.000E-01;RSET(D,CHAM )
************************************************************
Group 6. Body-Fitted coordinates
BFC=T
* Copy/Transfer/Block grid planes
GSET(C,I6,F,I1,1,28,1,18,+,5.0000E-02,0,0,INC,1)
GSET(C,I11,F,I6,1,28,1,18,+,5.0000E-02,0,0,INC,1)
GSET(C,I13,F,I11,1,28,1,18,+,3.0000E-02,0,0,INC,1)
GSET(C,K5,F,K1,1,12,1,28,+,0,0,7.0000E-02,INC,1)
GSET(C,K19,F,K5,1,12,1,28,+,0,0,3.3000E-01,INC,1)
GSET(C,J29,F,J19,1,12,1,18,+,0,1.3000E-01,0,INC,1)
GSET(C,J11,F,J1,1,12,1,18,+,0,1.3000E-01,0,INC,1)
GSET(C,J19,F,J11,1,12,1,18,+,0,3.8000E-01,0,INC,1)
GSET(C,J29,F,J19,1,12,1,18,+,0,1.3000E-01,0,INC,1)
GSET(C,J1,F,J11,1,12,18,18,+,0,-2.5000E-02,0,INC,1)
GSET(T,K19,F,K9,1,12,1,10,1.0)
GSET(C,J29,F,J19,1,12,18,18,+,0,2.5000E-02,0,INC,1)
GSET(T,K19,F,K9,1,12,19,28,1.0)
**********
NONORT = T
* X-cyclic boundaries switched
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME( 43) =EPOR ; NAME( 44) =BLOK
NAME( 45) =PRPS ; NAME( 46) =DEN1
NAME( 47) =TEM1 ; NAME( 48) =WCRT
NAME( 49) =VCRT ; NAME( 50) =UCRT
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,W1 ,TEM1)
* Stored variables list
STORE(UCRT,VCRT,WCRT,DEN1,PRPS,BLOK,EPOR,WDIS)
SOLUTN(P1 ,Y,Y,Y,N,N,N)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
IVARBK = -1 ;ISOLBK = 1
TURMOD(LVEL)
************************************************************
Group 8. Terms & Devices
NEWRH1 = T; NEWENL = T; NEWENT = T
ISOLX = 0 ;ISOLY = 0 ;ISOLZ = 0
************************************************************
Group 9. Properties
RHO1 = FILE
RHO1A = 1.161E+00 ;RHO1B = 0.000E+00 ;RHO1C = 0.000E+00
PRESS0 = 1.000E+05; TEMP0 = 2.730E+02
CP1 = FILE; ENUL = FILE
PRNDTL(LTLS) = GRND2 ;PRNDTL(TEM1) = CONDFILE
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(WDIS) = 1.300E-02 ;FIINIT(EPOR) = 1.000E+00
FIINIT(BLOK) = 1.000E+00 ;FIINIT(PRPS) = 2.000E+00
FIINIT(TEM1) = 2.000E+01
CONPOR(PCB1 , 0.00,EAST ,-#2,-#2,-#2,-#2,-#2,-#4)
CONPOR(FRONT , -1.00,CELL ,-#1,-#3,-#2,-#2,-#1,-#1)
INIT(FRONT ,BLOK, 0.000E+00, 2.000E+00)
INIT(FRONT ,PRPS, 0.000E+00, 1.090E+02)
CONPOR(SIDE , -1.00,CELL ,-#1,-#1,-#2,-#2,-#2,-#4)
INIT(SIDE ,BLOK, 0.000E+00, 3.000E+00)
INIT(SIDE ,PRPS, 0.000E+00, 1.000E+02)
INIADD = F
************************************************************
Group 13. Boundary & Special Sources
INLET (BFCFNIN ,LOW ,#1,#3,#3,#3,#1,#1,1,1)
VALUE (BFCFNIN ,P1 , GRND1 )
VALUE (BFCFNIN ,W1 , GRND1 )
VALUE (BFCFNIN ,TEM1, 2.000E+01)
VALUE (BFCFNIN ,WCRT, 1.500E+00)
INLET (BFCFNOUT,LOW ,#1,#3,#1,#1,#1,#1,1,1)
VALUE (BFCFNOUT,P1 , GRND1 )
VALUE (BFCFNOUT,W1 , GRND1 )
VALUE (BFCFNOUT,WCRT,-1.500E+00)
PATCH (HOTPCB1 ,EAST ,#2,#2,#2,#2,#2,#4,1,1)
COVAL (HOTPCB1 ,TEM1, FIXFLU , 3.000E+02)
PATCH (PCB2 ,EWALL ,#3,#3,#2,#2,#2,#4,1,1)
COVAL (PCB2 ,V1 , LOGLAW , 0.000E+00)
COVAL (PCB2 ,W1 , LOGLAW , 0.000E+00)
COVAL (PCB2 ,LTLS, 1.000E+00, 0.000E+00)
PATCH (HOTPCB2 ,EAST ,#3,#3,#2,#2,#2,#4,1,1)
COVAL (HOTPCB2 ,TEM1, FIXFLU , 3.000E+02)
PATCH (RELIEF ,CELL ,6,6,14,14,9,9,1,1)
COVAL (RELIEF ,P1 , 1.000E+00, 0.000E+00)
COVAL (RELIEF ,TEM1, 0.000E+00, SAME )
PATCH (BUOYANCY,PHASEM,#1,#NREGX,#1,#NREGY,#1,#NREGZ,#1,#NREGT)
COVAL (BUOYANCY,U1 , FIXFLU , DENSDIFF )
COVAL (BUOYANCY,V1 , FIXFLU , DENSDIFF )
COVAL (BUOYANCY,W1 , FIXFLU , DENSDIFF )
BUOYA = 0.000E+00 ; BUOYB =-9.810E+00 ; BUOYC = 0.000E+00
BUOYD = 1.189E+00
BFCA = 1.189E+00
EGWF = T
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 50; SELREF = T; RESFAC = 1.000E-02
************************************************************
Group 16. Terminate Iterations
LITER (TEM1) = 15
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000E+00); RELAX(U1 ,FALSDT, 4.267E-03)
RELAX(V1 ,FALSDT, 4.267E-03); RELAX(W1 ,FALSDT, 4.267E-03)
RELAX(TEM1,FALSDT, 4.267E+01)
************************************************************
Group 19. EARTH Calls To GROUND Station
NAMGRD =HTBX
************************************************************
Group 21. Print-out of Variables
OUTPUT(BLOK,N,N,N,N,N,N); OUTPUT(PRPS,N,N,N,N,N,N)
OUTPUT(WDIS,Y,N,N,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 6 ;IYMON = 10 ;IZMON = 6
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NXPRIN = 2; NYPRIN = 5; NZPRIN = 3
NPLT = 1; ITABL = 3
STOP