PHOTON USE
p
up y
VIEW 1 1 1
red
msg pressure contours at ground level
con p1 y 1 fi;0.1
set prop off
set por off
msg the buildings
surf prps y 198 col 4
surf prps x 198 col 5
surf prps z 198 col 6
gr ou x 1
gr ou y 1
gr z 1
msg a smoke plume
surf smok y 0.005 col 7
surf smok x 0.005 col 8
surf smok z 0.005 col 9
msg press return to view from the top
pause
con off;gr off
view y
msg press return to view from the side
pause
con off
view x
msg press return to view from the end
pause
con off
view z
enduse
Group 1. Run Title
TEXT(Flow Around A Group Of Buildings
TITLE
DISPLAY
The flow is steady. Smoke from a chimney is tracked.
A PHOTON USE file is supplied.
ENDDIS
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,15,10,48)
* Overall domain extent, RSET(D,name,XULAST,YVLAST,ZWLAST)
RSET(D,CHAM,1.5E+01,1.0E+01,4.8E+01)
* Set objects: name x0 y0 z0
* dx dy dz
RSET(B,X1, 2, 0, 0 ,0.0, 1.E+01, 4.8E+01)
RSET(B,X2, 3., 0.0, 0.0 ,0.0, 1.E+01, 4.8E+01)
RSET(B,X3, 4., 0.0, 0.0 ,0.0, 1.E+01, 4.8E+01)
RSET(B,X4, 6., 0.0, 0.0 ,0.0, 1.E+01, 4.8E+01)
RSET(B,X5, 8., 0.0, 0.0 ,0.0, 1.E+01, 4.8E+01)
RSET(B,X6, 1.E+01, 0.0, 0.0 ,0.0, 1.E+01, 4.8E+01)
RSET(B,X7, 1.1E+01, 0.0, 0.0 ,0.0, 1.E+01, 4.8E+01)
RSET(B,X8, 1.3E+01, 0.0, 0.0 ,0.0, 1.E+01, 4.8E+01)
RSET(B,Y1, 0.0, 3., 0.0 ,1.5E+01, 0.0, 4.8E+01)
RSET(B,Y2, 0.0, 4., 0.0 ,1.5E+01, 0.0, 4.8E+01)
RSET(B,Y3, 0.0, 5., 0.0 ,1.5E+01, 0.0, 4.8E+01)
RSET(B,Y4, 0.0, 8., 0.0 ,1.5E+01, 0.0, 4.8E+01)
RSET(B,Z1, 0.0, 0.0, 4. ,1.5E+01, 1.E+01, 0.0)
RSET(B,Z2, 0.0, 0.0, 6. ,1.5E+01, 1.E+01, 0.0)
RSET(B,Z3, 0.0, 0.0, 8. ,1.5E+01, 1.E+01, 0.0)
RSET(B,Z4, 0.0, 0.0, 1.2E+01 ,1.5E+01, 1.E+01, 0.0)
RSET(B,Z5, 0.0, 0.0, 1.6E+01 ,1.5E+01, 1.E+01, 0.0)
RSET(B,Z6, 0.0, 0.0, 1.8E+01 ,1.5E+01, 1.E+01, 0.0)
RSET(B,Z7, 0.0, 0.0, 2.E+01 , 1.5E+01, 1.E+01, 0.0)
RSET(B,Z8, 0.0, 0.0, 2.4E+01 ,1.5E+01, 1.E+01, 0.0)
RSET(B,Z9, 0.0, 0.0, 2.8E+01 ,1.5E+01, 1.E+01, 0.0)
RSET(B,Z10 , 0.0, 0.0, 3.E+01 ,1.5E+01, 1.E+01, 0.0)
RSET(B,Z11 , 0.0, 0.0, 3.2E+01 ,1.5E+01, 1.E+01, 0.0)
RSET(B,Z12 , 0.0, 0.0, 3.6E+01 , 1.5E+01, 1.E+01, 0.0)
RSET(B,Z13 , 0.0, 0.0, 4.E+01 , 1.5E+01, 1.E+01, 0.0)
RSET(B,Z14 , 0.0, 0.0, 4.2E+01 , 1.5E+01, 1.E+01, 0.0)
RSET(B,Z15 , 0.0, 0.0, 4.4E+01 , 1.5E+01, 1.E+01, 0.0)
RSET(B,GROUND , 0.0, 0.0, 0.0 , 1.5E+01, 0.0, 4.8E+01)
RSET(B,SKY , 0.0, 1.E+01, 0.0 , 1.5E+01, 0.0, 4.8E+01)
RSET(B,LOWIN1 , 0.0, 0.0, 0.0 , 1.5E+01, 1.E+01, 0.0)
RSET(B,WESTIN2 , 0.0, 0.0, 0.0 , 0.0, 1.E+01, 4.8E+01)
RSET(B,HIGHOT1 , 0.0, 0.0, 4.8E+01 , 1.5E+01, 1.E+01, 0.0)
RSET(B,EASTOT2 , 1.5E+01, 0.0, 0.0 , 0.0, 1.E+01, 4.8E+01)
RSET(B,SUBD-Z1 , 0.0, 0.0, 0.0 , 1.5E+01, 1.E+01, 1.2E+01)
RSET(B,SUBD-Z2 , 0.0, 0.0, 1.2E+01 , 1.5E+01, 1.E+01, 1.2E+01)
RSET(B,SUBD-Z3 , 0.0, 0.0, 2.4E+01 , 1.5E+01, 1.E+01, 1.2E+01)
RSET(B,SUBD-Z4 , 0.0, 0.0, 3.6E+01 , 1.5E+01, 1.E+01, 1.2E+01)
* Modify default grid
RSET(X,1,2,1.0)
RSET(X,2,1,1.0)
RSET(X,3,1,1.0)
RSET(X,4,2,1.0)
RSET(X,5,2,1.0)
RSET(X,6,2,1.0)
RSET(X,7,1,1.0)
RSET(X,8,2,1.0)
RSET(X,9,2,1.0)
RSET(Y,1,3,1.0)
RSET(Y,2,1,1.0)
RSET(Y,3,1,1.0)
RSET(Y,4,3,1.0)
RSET(Y,5,2,1.0)
RSET(Z,1,4,1.0)
RSET(Z,2,2,1.0)
RSET(Z,3,2,1.0)
RSET(Z,4,4,1.0)
RSET(Z,5,4,1.0)
RSET(Z,6,2,1.0)
RSET(Z,7,2,1.0)
RSET(Z,8,4,1.0)
RSET(Z,9,4,1.0)
RSET(Z,10,2,1.0)
RSET(Z,11,2,1.0)
RSET(Z,12,4,1.0)
RSET(Z,13,4,1.0)
RSET(Z,14,2,1.0)
RSET(Z,15,2,1.0)
RSET(Z,16,4,1.0)
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,W1 ,smok)
* Stored variables list
STORE(VPOR)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,N)
SOLUTN(SMOK,Y,Y,Y,N,N,N)
Group 8. Terms & Devices
CSG3=CNGR
Group 9. Properties
RHO1 = 1.0
PRESS0 = 1.E+05
ENUL = 1.E-05 ;ENUT = 0.0
Group 11.Initialise Var/Porosity Fields
FIINIT(U1 ) = 5. ;FIINIT(V1 ) = 0.0
FIINIT(W1 ) = 5. ;FIINIT(VPOR) = 1.0
CONPOR(S1B1 , 0.00,CELL ,#3,#4,#1,#2,#2,#2)
CONPOR(S1B2 , 0.00,CELL ,#2,#3,#1,#3,#3,#3)
CONPOR(S1B3 , 0.00,CELL ,#6,#7,#1,#2,#2,#2)
CONPOR(S1B4 , 0.00,CELL ,#7,#8,#1,#3,#3,#3)
CONPOR(S2B1 , 0.00,CELL ,#5,#5,#1,#4,#5,#5)
CONPOR(S2B2 , 0.00,CELL ,#4,#6,#1,#2,#6,#6)
CONPOR(S3B1 , 0.00,CELL ,#2,#3,#1,#1,#8,#8)
CONPOR(S3B2 , 0.00,CELL ,#7,#8,#1,#1,#8,#8)
CONPOR(S34B1 , 0.00,CELL ,#2,#2,#1,#1,#9,#11)
CONPOR(S34B2 , 0.00,CELL ,#8,#8,#1,#1,#9,#11)
CONPOR(S34B3 , 0.00,CELL ,#5,#5,#1,#4,#10,#10)
CONPOR(S4B1 , 0.00,CELL ,#2,#3,#1,#1,#12,#12)
CONPOR(S4B2 , 0.00,CELL ,#7,#8,#1,#1,#12,#12)
INIADD = F
Group 13. Boundary & Special Sources
PATCH (GROUND ,SOUTH ,#1,#9,#1,#1,#1,#16,1,1)
PATCH (SKY ,NORTH ,#1,#9,#5,#5,#1,#16,1,1)
COVAL (SKY ,P1 , 1.0, 0.0)
INLET (LOWIN1 ,LOW ,#1,#9,#1,#5,#1,#1,1,1)
VALUE (LOWIN1 ,P1 , 7.070)
VALUE (LOWIN1 ,U1 , 5.)
VALUE (LOWIN1 ,W1 , 5.)
INLET (WESTIN2 ,WEST ,#1,#1,#1,#5,#1,#16,1,1)
VALUE (WESTIN2 ,P1 , 7.070)
VALUE (WESTIN2 ,U1 , 5.)
VALUE (WESTIN2 ,W1 , 5.)
PATCH (HIGHOT1 ,HIGH ,#1,#9,#1,#5,#16,#16,1,1)
COVAL (HIGHOT1 ,P1 , 1.0, 0.0)
PATCH (EASTOT2 ,EAST ,#9,#9,#1,#5,#1,#16,1,1)
COVAL (EASTOT2 ,P1 , 1.0, 0.0)
patch(chimney,cell,1,1,ny/2,ny/2,1,1,1,1)
coval(chimney,smok,fixval,1.0)
Group 15. Terminate Sweeps
LSWEEP = 100; SELREF = T; RESFAC = 1.E-02
Group 16. Terminate Iterations
LITER(P1)=100
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.0); RELAX(U1 ,FALSDT, 1.0)
RELAX(V1 ,FALSDT, 1.0); RELAX(W1 ,FALSDT, 1.0)
Group 21. Print-out of Variables
OUTPUT(VPOR,N,N,N,Y,N,N)
Group 22. Monitor Print-Out
IXMON = 8 ;IYMON = 6 ;IZMON = 10
TSTSWP=-1
Group 23.Field Print-Out & Plot Control
NZPRIN = 1
NPLT = 1
No PATCHes used for this Group
Group 24. Dumps For Restarts