PHOTON USE
p ; ; ; ; ;
msg Computational Grid
mgr 1 z 1 x 2 6 y 2 6 col 4
mgr 2 z 1 x 1 5 y 2 6 col 6
mgr 3 z 1 x 2 6 y 2 6 col 9
mgr 4 z 1 x 2 6 y 2 6 col 12
mgr 5 z 1 x 2 6 y 1 5 col 14
mgr 1 ou z m x 2 6 y 2 6
mgr 2 ou z m x 1 5 y 2 6
mgr 3 ou z m x 2 6 y 2 6
mgr 4 ou z m x 2 6 y 2 6
mgr 5 ou z m x 2 6 y 1 5
mgr 2 ou x 1 y 2 6
mgr 3 ou x m y 2 6
mgr 4 ou y m x 2 6
mgr 5 ou y 1 x 2 6
msg Press return for Velocity Vectors
pause
mvec 1 x 4 sh
mvec 4 x 4 sh
mvec 5 x 4 sh
red
msg Vectors on Centre Plane
msg Press return for Vectors at entry plane
pause
vec cl;vi z
mvec 1 z 1 sh
mvec 2 z 1 sh
mvec 3 z 1 sh
mvec 4 z 1 sh
mvec 5 z 1 sh
red
msg Vectors at Entry Plane
msg Press return for Vectors near exit plane
pause
vec cl;vi y;up -z
mvec 1 z 13 sh
mvec 2 z 13 sh
mvec 3 z 13 sh
mvec 4 z 13 sh
mvec 5 z 13 sh
red
msg Vectors near Exit Plane
msg Press return for Pressure Contours
pause
vec cl;vi 1 1 1;up y
mcon 1 p1 x 4 fi
.001
mcon 4 p1 x 4 fi
.001
mcon 5 p1 x 4 fi
.001
msg Pressure contours
msg Press return for LTLS contours near exit plane
pause
vec cl;con cl; vi y;up -z
mcon 1 ltls z 13 fi
.001
mcon 2 ltls z 13 fi
.001
mcon 3 ltls z 13 fi
.001
mcon 4 ltls z 13 fi
.001
mcon 5 ltls z 13 fi
.001
red
msg Press E then return to end
ENDUSE
************************************************************
Q1 created by GENERAL menu, Version 2.2, Date 17/09/96
CPVNAM=GENERAL
************************************************************
Group 1. Run Title
TEXT(Pipe Flow with Multi-Block Grid and LVEL:585
TITLE
************************************************************
Group 6. Body-Fitted coordinates
BFC=T
* Set points
XPO= 0.0000E+00;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P1 )
XPO= 5.0000E-01;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P2 )
XPO= 5.0000E-01;YPO= 5.0000E-01;ZPO= 0.0000E+00;GSET(P,P3 )
XPO= 0.0000E+00;YPO= 5.0000E-01;ZPO= 0.0000E+00;GSET(P,P4 )
XPO= 9.5710E-01;YPO= 9.5710E-01;ZPO= 0.0000E+00;GSET(P,P5 )
XPO= 9.5710E-01;YPO=-4.5710E-01;ZPO= 0.0000E+00;GSET(P,P6 )
XPO=-4.5710E-01;YPO=-4.5710E-01;ZPO= 0.0000E+00;GSET(P,P7 )
XPO=-4.5710E-01;YPO= 9.5710E-01;ZPO= 0.0000E+00;GSET(P,P8 )
XPO= 2.5000E-01;YPO= 1.2500E+00;ZPO= 0.0000E+00;GSET(P,P9 )
XPO= 1.2500E+00;YPO= 2.5000E-01;ZPO= 0.0000E+00;GSET(P,P10 )
XPO= 2.5000E-01;YPO=-7.5000E-01;ZPO= 0.0000E+00;GSET(P,P11 )
XPO=-7.5000E-01;YPO= 2.5000E-01;ZPO= 0.0000E+00;GSET(P,P12 )
* Set lines/arcs
GSET(L,L1,P1,P2,5,1.0)
GSET(L,L2,P2,P3,5,1.0)
GSET(L,L3,P3,P4,5,1.0)
GSET(L,L4,P4,P1,5,1.0)
GSET(L,L5,P4,P8,5,1.0)
GSET(L,L6,P3,P5,5,1.0)
GSET(L,L7,P2,P6,5,1.0)
GSET(L,L8,P1,P7,5,1.0)
GSET(L,A9,P8,P5,5,1.0,ARC,P9)
GSET(L,A10,P5,P6,5,1.0,ARC,P10)
GSET(L,A11,P6,P7,5,1.0,ARC,P11)
GSET(L,A12,P7,P8,5,1.0,ARC,P12)
* Set frames
GSET(F,F1,P1,-,P2,-,P3,-,P4,-)
GSET(F,F2,P7,-,P1,-,P4,-,P8,-)
GSET(F,F3,P2,-,P6,-,P5,-,P3,-)
GSET(F,F4,P4,-,P3,-,P5,-,P8,-)
GSET(F,F5,P7,-,P6,-,P2,-,P1,-)
* Match a grid mesh
GSET(D,7,7,15)
GSET(M,F1,+I+J,2,2,1,TRANS)
GSET(C,K6,F,K1,1,NX,1,NY,+,0,0,1,INC,1)
GSET(C,K11,F,K6,1,NX,1,NY,RX,-1.5710E+00,3,1,INC,1)
GSET(C,K16,F,K11,1,NX,1,NY,+,0,1,0,INC,1)
GSET(C,I1,F,I2,1,NY,1,NZ,+,0,0,0)
GSET(C,I:NX+1:,F,I:NX:,1,NY,1,NZ,+,0,0,0)
GSET(C,J1,F,J2,1,NX,1,NZ,+,0,0,0)
GSET(C,J:NY+1:,F,J:NY:,1,NX,1,NZ,+,0,0,0)
DUMPC(GRID1)
GSET(D,6,7,15)
GSET(M,F2,+I+J,1,2,1,TRANS)
GSET(C,K6,F,K1,1,NX,1,NY,+,0,0,1,INC,1)
GSET(C,K11,F,K6,1,NX,1,NY,RX,-1.5710E+00,3,1,INC,1)
GSET(C,K16,F,K11,1,NX,1,NY,+,0,1,0,INC,1)
GSET(C,I:NX+1:,F,I:NX:,1,NY,1,NZ,+,0,0,0)
GSET(C,J1,F,J2,1,NX,1,NZ,+,0,0,0)
GSET(C,J:NY+1:,F,J:NY:,1,NX,1,NZ,+,0,0,0)
DUMPC(GRID2)
GSET(D,6,7,15)
GSET(M,F3,+I+J,2,2,1,TRANS)
GSET(C,K6,F,K1,1,NX,1,NY,+,0,0,1,INC,1)
GSET(C,K11,F,K6,1,NX,1,NY,RX,-1.5710E+00,3,1,INC,1)
GSET(C,K16,F,K11,1,NX,1,NY,+,0,1,0,INC,1)
GSET(C,I1,F,I2,1,NY,1,NZ,+,0,0,0)
GSET(C,J1,F,J2,1,NX,1,NZ,+,0,0,0)
GSET(C,J:NY+1:,F,J:NY:,1,NX,1,NZ,+,0,0,0)
DUMPC(GRID3)
GSET(D,7,6,15)
GSET(M,F4,+I+J,2,2,1,TRANS)
GSET(C,K6,F,K1,1,NX,1,NY,+,0,0,1,INC,1)
GSET(C,K11,F,K6,1,NX,1,NY,RX,-1.5710E+00,3,1,INC,1)
GSET(C,K16,F,K11,1,NX,1,NY,+,0,1,0,INC,1)
GSET(C,I1,F,I2,1,NY,1,NZ,+,0,0,0)
GSET(C,I:NX+1:,F,I:NX:,1,NY,1,NZ,+,0,0,0)
GSET(C,J1,F,J2,1,NX,1,NZ,+,0,0,0)
DUMPC(GRID4)
GSET(D,7,6,15)
GSET(M,F5,+I+J,2,1,1,TRANS)
GSET(C,K6,F,K1,1,NX,1,NY,+,0,0,1,INC,1)
GSET(C,K11,F,K6,1,NX,1,NY,RX,-1.5710E+00,3,1,INC,1)
GSET(C,K16,F,K11,1,NX,1,NY,+,0,1,0,INC,1)
GSET(C,I1,F,I2,1,NY-1,1,NZ,+,0,0,0)
GSET(C,I:NX+1:,F,I:NX:,1,NY,1,NZ,+,0,0,0)
GSET(C,J:NY+1:,F,J:NY:,1,NX,1,NZ,+,0,0,0)
DUMPC(GRID5)
**********
NUMBLK=5
### Make default grid BIG, so no cells are left-handed
XSI=5; YSI=5; ZSI=5; RSET(D,CHAM)
READCO(GRID+)
MESG(Show the grid? (N/y)
READVDU(ANS,CHAR,N)
IF(:ANS:.EQ.Y) THEN
+ READCO(VIEWONLY)
+ VIEW
ENDIF
NONORT = T
* X-cyclic boundaries switched
MPATCH(1,MBL1-2,WEST , 2, 2,2,NY-1,1,NZ,1,1)
MPATCH(2,MBL2-1,EAST ,NX-1,NX-1,2,NY-1,1,NZ,1,1)
MPATCH(1,MBL1-3,EAST ,NX-1,NX-1,2,NY-1,1,NZ,1,1)
MPATCH(3,MBL3-1,WEST , 2, 2,2,NY-1,1,NZ,1,1)
MPATCH(1,MBL1-4,NORTH,2,NX-1,NY-1,NY-1,1,NZ,1,1)
MPATCH(4,MBL4-1,SOUTH,2,NX-1, 2, 2,1,NZ,1,1)
MPATCH(1,MBL1-5,SOUTH,2,NX-1, 2, 2,1,NZ,1,1)
MPATCH(5,MBL5-1,NORTH,2,NX-1,NY-1,NY-1,1,NZ,1,1)
MPATCH(2,MBL2-4,NORTH,1,NX-1,NY-1,NY-1,1,NZ,1,1)
MPATCH(4,MBL4-2,WEST ,2, 2, 2,NY ,1,NZ,1,1)
SPEDAT(SET,GCV,MBL4-2,C,WNL)
MPATCH(2,MBL2-5,SOUTH,1,NX-1,2, 2,1,NZ,1,1)
MPATCH(5,MBL5-2,WEST ,2, 2,1,NY-1,1,NZ,1,1)
SPEDAT(SET,GCV,MBL5-2,C,ESL)
MPATCH(3,MBL3-4,NORTH, 2,NX ,NY-1,NY-1,1,NZ,1,1)
MPATCH(4,MBL4-3,EAST ,NX-1,NX-1, 2,NY ,1,NZ,1,1)
SPEDAT(SET,GCV,MBL4-3,C,ESL)
MPATCH(3,MBL3-5,SOUTH, 2, NX,2, 2,1,NZ,1,1)
MPATCH(5,MBL5-3,EAST ,NX-1,NX-1,1,NY-1,1,NZ,1,1)
SPEDAT(SET,GCV,MBL5-3,C,WNL)
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,W1 )
* Stored variables list
STORE(PRPS,DEN1,UCRT,VCRT,WCRT,WDIS,WGAP,ENUT,EL1)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
TURMOD(LVEL)
************************************************************
Group 8. Terms & Devices
NEWRH1 = T ; NEWENL = T
************************************************************
Group 9. Properties
RHO1=1.189
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(PRPS) = 0.000E+00
INIADD = T
************************************************************
Group 13. Boundary & Special Sources
REAL(WIN,FAC)
WIN=5; FAC=1
MPATCH(1,IN1,LOW,2,NX-1,2,NY-1,1,1,1,1)
COVAL(IN1,P1,FIXFLU,RHO1*WIN)
COVAL(IN1,W1,ONLYMS,WIN)
MPATCH(2,IN2,LOW,1,NX-1,2,NY-1,1,1,1,1)
COVAL(IN2,P1,FIXFLU,RHO1*WIN)
COVAL(IN2,W1,ONLYMS,WIN)
COVAL(IN2,V1,ONLYMS,WIN/FAC)
MPATCH(3,IN3,LOW,2,NX,2,NY-1,1,1,1,1)
COVAL(IN3,P1,FIXFLU,RHO1*WIN)
COVAL(IN3,W1,ONLYMS,WIN)
COVAL(IN3,V1,ONLYMS,-WIN/FAC)
MPATCH(4,IN4,LOW,2,NX-1,2,NY,1,1,1,1)
COVAL(IN4,P1,FIXFLU,RHO1*WIN)
COVAL(IN4,W1,ONLYMS,WIN)
COVAL(IN4,U1,ONLYMS,WIN/FAC)
MPATCH(5,IN5,LOW,2,NX-1,1,NY-1,1,1,1,1)
COVAL(IN5,P1,FIXFLU,RHO1*WIN)
COVAL(IN5,W1,ONLYMS,WIN)
COVAL(IN5,U1,ONLYMS,-WIN/FAC)
MPATCH(1,OUT1,HIGH,2,NX-1,2,NY-1,NZ,NZ,1,1)
COVAL(OUT1,P1,1E3,0)
MPATCH(2,OUT2,HIGH,1,NX-1,2,NY-1,NZ,NZ,1,1)
COVAL(OUT2,P1,1E3,0)
MPATCH(3,OUT3,HIGH,2,NX ,2,NY-1,NZ,NZ,1,1)
COVAL(OUT3,P1,1E3,0)
MPATCH(4,OUT4,HIGH,2,NX-1,2,NY ,NZ,NZ,1,1)
COVAL(OUT4,P1,1E3,0)
MPATCH(5,OUT5,HIGH,2,NX-1,1,NY-1,NZ,NZ,1,1)
COVAL(OUT5,P1,1E3,0)
MPATCH(2,WALL1,WWALL,1,1,2,NY-1,1,NZ,1,1)
COVAL(WALL1,W1,GRND2,0)
MPATCH(3,WALL2,EWALL,NX,NX,2,NY-1,1,NZ,1,1)
COVAL(WALL2,W1,GRND2,0)
MPATCH(4,WALL3,NWALL,2,NX-1,NY,NY,1,NZ,1,1)
COVAL(WALL3,W1,GRND2,0)
MPATCH(5,WALL4,SWALL,2,NX-1,1,1,1,NZ,1,1)
COVAL(WALL4,W1,GRND2,0)
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 100; SELREF = T; RESFAC = 1.000E-03
************************************************************
Group 17. Relaxation
RELAX(U1 ,FALSDT, 3/WIN)
RELAX(V1 ,FALSDT, 3/WIN)
RELAX(W1 ,FALSDT, 3/WIN)
************************************************************
Group 19. EARTH Calls To GROUND Station
GCV=T
************************************************************
Group 22. Monitor Print-Out
IXMON = 18 ;IYMON = 3 ;IZMON = 7
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 24. Dumps For Restarts
************************************************************