DISPLAY
*********************************************************
This case concerns the prediction of air flow and temperature
distribution in and around a closed metal box, which contains
two heat sources. These liberate 5W and 7W respectively. The
box is cooled from the outside by air at 20 deg C being blown
at 1.5 m/s.
The case takes about 45 min to run on a 486/50. The conjugate
gradient solver is used to improve convergence of the
temperature equation.
The results can be plotted in PHOTON by typing 'use q1'
ENDDIS
PHOTON USE
p
GR OU X 1; GR OU X M; GR OU Y 1; GR OU Y M
gr ou z 4 y 1 10 x 5 20; gr ou z 12 y 1 10 x 5 20
gr ou x 5 y 1 10 z 4 11; gr ou x 21 y 1 10 z 4 11
gr ou z 5 y 1 9 x 6 19; gr ou z 11 y 1 9 x 6 19
gr ou x 6 y 1 9 z 5 10; gr ou x 20 y 1 9 z 5 10
GR OU X 10 Y 2 3 Z 7 8 COL 13
GR OU X 8 Y 2 3 Z 7 8 COL 13
GR OU Y 4 X 8 9 Z 7 8 COL 13
GR OU Y 2 X 8 9 Z 7 8 COL 13
GR OU X 16 Y 2 3 Z 7 8 COL 13
GR OU X 14 Y 2 3 Z 7 8 COL 13
GR OU Y 4 X 14 15 Z 7 8 COL 13
GR OU Y 2 X 14 15 Z 7 8 COL 13
GR OU X 1 Y 1 16 Z 1 15 COL 5
GR OU X 34 Y 1 16 Z 1 15 COL 12
msg The geometry - Two heat sources in a closed metal box
msg which is force-cooled
pause
vec z 8 sh
con tem1 z 8 sh;int 28
msg Vectors and temperatures on the center plane
pause
vi z
msg Viewed from the side
pause
con cl;vec cl;mag gr 3
0.94204E+03 0.15046E+04 CR
set vec ref .1
vec z 8 x 6 19 y 1 9 sh
msg Flow field inside the closed box
pause
vec cl; vi 1 1 1
set con sca pl
con tem1 x 25 fi;.001
msg Temperature rise in wake of box
ENDUSE
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(Heat Sources; Air-Cooled Box; EGWF
TITLE
Originated by JCL Jan 1994
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,33,16,15)
* Overall domain extent, RSET(D,name,XULAST,YVLAST,ZWLAST)
RSET(D,AIR,1.500E+00,3.000E-01,5.000E-01)
* Set objects: name x0 y0 z0
* dx dy dz
RSET(B,BOX1 , 2.000E-01, 0.000E+00, 1.000E-01 , $
1.500E-02, 1.000E-01, 2.500E-01)
RSET(B,BOX2 , 2.150E-01, 0.000E+00, 1.000E-01 , $
4.000E-01, 1.000E-01, 1.500E-02)
RSET(B,BOX3 , 2.150E-01, 0.000E+00, 3.350E-01 , $
4.000E-01, 1.000E-01, 1.500E-02)
RSET(B,BOX4 , 6.150E-01, 0.000E+00, 1.000E-01 , $
1.500E-02, 1.000E-01, 2.500E-01)
RSET(B,TOP , 2.000E-01, 1.000E-01, 1.000E-01 , $
4.300E-01, 5.000E-03, 2.500E-01)
RSET(B,CHP1 , 2.774E-01, 1.000E-02, 1.889E-01 , $
5.570E-02, 2.000E-02, 6.334E-02)
RSET(B,CHP2 , 4.641E-01, 1.000E-02, 1.889E-01 , $
5.570E-02, 2.000E-02, 6.334E-02)
RSET(B,FAN1 , 0.000E+00, 0.000E+00, 0.000E+00 , $
0.000E+00, 3.000E-01, 5.000E-01)
RSET(B,OUT1 , 1.500E+00, 0.000E+00, 0.000E+00 , $
0.000E+00, 3.000E-01, 5.000E-01)
* Modify default grid
RSET(X,1,4,-1.500E+00)
RSET(X,3,2,1.000E+00)
RSET(X,4,2,1.000E+00)
RSET(X,5,4,1.000E+00)
RSET(X,6,2,1.000E+00)
RSET(X,7,4,1.000E+00)
RSET(X,9,13,1.700E+00)
RSET(Y,2,2,1.000E+00)
RSET(Y,3,6,1.000E+00)
RSET(Y,5,6,1.700E+00)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(47) =BLOK ; NAME(48) =PRPS
NAME(49) =DEN1 ; NAME(50) =TEM1
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,W1 ,TEM1)
* Stored variables list
STORE(DEN1,PRPS,BLOK)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,N)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
IVARBK = -1 ; ISOLBK = 1
************************************************************
Group 8. Terms & Devices
NEWRH1 = T
NEWENL = T
ISOLX = 0 ; ISOLY = 0 ; ISOLZ = 0
************************************************************
Group 9. Properties
RHO1 = FILE
PRESS0 = 1.000E+05
TEMP0 = 2.730E+02
ENUL = FILE ; ENUT = 0.000E+00
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(BLOK) = 1.000E+00 ; FIINIT(PRPS) = 2.000E+00
FIINIT(TEM1) = 2.000E+01
CONPOR(BOX1 , -1.00,CELL ,-#2,-#2,-#1,-#3,-#2,-#6)
INIT(BOX1 ,BLOK, 0.000E+00, 2.000E+00)
INIT(BOX1 ,PRPS, 0.000E+00, 1.110E+02)
CONPOR(BOX2 , -1.00,CELL ,-#3,-#7,-#1,-#3,-#2,-#2)
INIT(BOX2 ,BLOK, 0.000E+00, 3.000E+00)
INIT(BOX2 ,PRPS, 0.000E+00, 1.110E+02)
CONPOR(BOX3 , -1.00,CELL ,-#3,-#7,-#1,-#3,-#6,-#6)
INIT(BOX3 ,BLOK, 0.000E+00, 4.000E+00)
INIT(BOX3 ,PRPS, 0.000E+00, 1.110E+02)
CONPOR(BOX4 , -1.00,CELL ,-#8,-#8,-#1,-#3,-#2,-#6)
INIT(BOX4 ,BLOK, 0.000E+00, 5.000E+00)
INIT(BOX4 ,PRPS, 0.000E+00, 1.110E+02)
CONPOR(TOP , -1.00,CELL ,-#2,-#8,-#4,-#4,-#2,-#6)
INIT(TOP ,BLOK, 0.000E+00, 6.000E+00)
INIT(TOP ,PRPS, 0.000E+00, 1.110E+02)
CONPOR(CHP1 , -1.00,CELL ,-#4,-#4,-#2,-#2,-#4,-#4)
INIT(CHP1 ,BLOK, 0.000E+00, 7.000E+00)
INIT(CHP1 ,PRPS, 0.000E+00, 1.060E+02)
CONPOR(CHP2 , -1.00,CELL ,-#6,-#6,-#2,-#2,-#4,-#4)
INIT(CHP2 ,BLOK, 0.000E+00, 8.000E+00)
INIT(CHP2 ,PRPS, 0.000E+00, 1.060E+02)
INIADD = F
************************************************************
Group 12. Convection and diffusion adjustments
************************************************************
Group 13. Boundary & Special Sources
PATCH (HOTCHP1 ,VOLUME,#4,#4,#2,#2,#4,#4,#1,#1)
COVAL (HOTCHP1 ,TEM1, FIXFLU , 7.087E+04)
PATCH (HOTCHP2 ,VOLUME,#6,#6,#2,#2,#4,#4,#1,#1)
COVAL (HOTCHP2 ,TEM1, FIXFLU , 9.922E+04)
INLET (FAN1 ,WEST ,#1,#1,#1,#5,#1,#7,1,1);
VALUE (FAN1 ,P1 , 1.783E+00)
VALUE (FAN1 ,U1 , 1.500E+00)
VALUE (FAN1 ,TEM1, 2.000E+01)
PATCH (OUT1 ,EAST ,#9,#9,#1,#5,#1,#7,1,1)
COVAL (OUT1 ,P1 , 1.000E+00, 0.000E+00)
COVAL (OUT1 ,TEM1, 0.000E+00, 2.000E+01)
#gravity
gravdir=2; rhoref=1.189
#densdiff
PATCH (BUOYANCY,PHASEM,1,33,1,16,1,15,1,1)
COVAL (BUOYANCY,V1 , FIXFLU ,DENSDIFF)
BUOYA = 0.000E+00 ; BUOYB =-9.810E+00 ; BUOYC = 0.000E+00
BUOYD = 1.189E+00
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 200
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, 1.000E-01)
RELAX(V1 ,FALSDT, 1.000E-01)
RELAX(W1 ,FALSDT, 1.000E-01)
RELAX(TEM1,FALSDT, 1.000E+04)
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
NAMGRD =HTBX
CSG3=CRGR
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
OUTPUT(BLOK,N,N,N,N,N,N)
OUTPUT(PRPS,N,N,N,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 5 ;IYMON = 5 ;IZMON = 5
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NXPRIN = 6
NYPRIN = 3
NZPRIN = 3
NPLT = 1
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
mesg(Do you want to see the geometry? (N/y)
readvdu(ans,char,n)
if(:ans:.eq.y) then
+ rset(v,0); gview(p,1,1,1); gdom(1,nx+1,1,ny+1,1,nz+1,1,0)
+ gdom(5,21,1,11,4,12,6,0); gdom(6,20,1,10,5,11,6,0)
+ gpatch(chp1,15,0); gpatch(chp2,15,0)
+ gpatch(fan1,4,0); gpatch(out1,12,0); gdraw
endif