DISPLAY
This case concerns flow entering a sequence of slots.
\ Approaching fluid
\
\|
-
----- ----- ----- -----
| | | | | |
| | | | | |
| | | | | |
Use of the XCYCLE condition means that only one slot need be
modelled.
BFC is needed because the slots hve radiused corners.
The Q1 contains PHOTON instructions.
ENDDIS
PHOTON USE
p
gr z 1
gr z 1 x 1 10 y 1 30 col 4 ; gr z 1 x 31 m y 1 30 col 4
msg The geometry
msg Press RETURN for Velocity vectors
pause
cl
gr ou z 1;gr ou z 1 x 1 10 y 1 30; gr ou z 1 x 31 m y 1 30
vec z 1 sh
msg Velocity Vectors
msg Press RETURN for Pressure Contours
pause
vec cl;red;con p1 z 1 fi;.0005
msg Pressure Contours
msg Press RETURN for contours of Turbulrnt Viscosity
pause
con cl;red;con enut z 1 fi;.0005
msg Press E the RETURN to end
ENDUSE
MESG(Show the grid? (N/y)
READVDU(ANS,CHAR,N)
************************************************************
Group 1. Run Title
TEXT(Radiused slot entry :565
TITLE
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,40/1,50/1,1)
* Set overall domain extent:
* xulast yvlast zwlast name
XSI= 5.317E-02;YSI= 1.500E-01;ZSI= 1.000E-03;RSET(D,CHAM )
************************************************************
Group 6. Body-Fitted coordinates
BFC=T
* Set poInts
XPO= 0.0000E+00;YPO= 1.5000E-01;ZPO= 0.0000E+00;GSET(P,P4 )
XPO= 9.2360E-03;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P5 )
XPO= 9.2360E-03;YPO= 9.5000E-02;ZPO= 0.0000E+00;GSET(P,P6 )
XPO= 7.7720E-03;YPO= 9.8540E-02;ZPO= 0.0000E+00;GSET(P,P7 )
XPO= 4.2360E-03;YPO= 1.0000E-01;ZPO= 0.0000E+00;GSET(P,P8 )
XPO= 4.3940E-02;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P10 )
XPO= 4.3940E-02;YPO= 9.5000E-02;ZPO= 0.0000E+00;GSET(P,P11 )
XPO= 4.5400E-02;YPO= 9.8540E-02;ZPO= 0.0000E+00;GSET(P,P12 )
XPO= 4.8940E-02;YPO= 1.0000E-01;ZPO= 0.0000E+00;GSET(P,P13 )
XPO= 5.3170E-02;YPO= 1.5000E-01;ZPO= 0.0000E+00;GSET(P,P18 )
XPO= 5.3170E-02;YPO= 1.0000E-01;ZPO= 0.0000E+00;GSET(P,P48 )
XPO= 0.0000E+00;YPO= 1.0000E-01;ZPO= 0.0000E+00;GSET(P,P66 )
XPO= 0.0000E+00;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P75 )
XPO= 5.3170E-02;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P76 )
XPO= 6.2821E-03;YPO= 9.9516E-02;ZPO= 0.0000E-00;GSET(P,P30 )
XPO= 8.7718E-03;YPO= 9.7137E-02;ZPO= 0.0000E-00;GSET(P,P31 )
XPO= 4.6963E-02;YPO= 9.9582E-02;ZPO= 0.0000E-00;GSET(P,P32 )
XPO= 4.4518E-02;YPO= 9.7328E-02;ZPO= 0.0000E-00;GSET(P,P33 )
XPO= 7.9402E-03;YPO= 1.5000E-01;ZPO= 0.0000E-00;GSET(P,P34 )
XPO= 4.5291E-02;YPO= 1.5000E-01;ZPO= 0.0000E-00;GSET(P,P35 )
* Set lines/arcs
GSET(L,A2,P6,P7,3,1.0,ARC,P31)
GSET(L,A3,P7,P8,3,1.0,ARC,P30)
GSET(L,A5,P13,P12,3,1.0,ARC,P32)
GSET(L,A6,P12,P11,3,1.0,ARC,P33)
GSET(L,L7,P75,P5,10/1,1.0)
GSET(L,L8,P5,P10,20/1,1.0)
GSET(L,L9,P10,P76,10/1,1.0)
GSET(L,L10,P75,P66,30/1,1.0)
GSET(L,L11,P5,P6,27-15*0,1.0)
GSET(L,L12,P66,P8,7-5*0,1.0)
GSET(L,L13,P10,P11,27-15*0,1.0)
GSET(L,L14,P76,P48,30/1,1.0)
GSET(L,L15,P66,P4,20/1,1.0)
GSET(L,L17,P18,P48,20/1,1.0)
GSET(L,L18,P48,P13,7-5*0,1.0)
GSET(L,L22,P4,P34,10/1,1.0)
GSET(L,L23,P34,P35,20/1,1.0)
GSET(L,L24,P35,P18,10/1,1.0)
* Set frames
GSET(F,F5,P75,-,P5,P6,P7,P8,P66,-)
GSET(F,F7,P10,-,P76,-,P48,P13,P12,P11)
GSET(F,F8,P75,P5.P10,P76,P48,P18,P35.P34,P4,P66)
* Match a grid mesh
GSET(M,F5,+I+J,1,1,1,TRANS)
GSET(M,F7,+I+J,30/1+1,1,1,TRANS)
GSET(M,F8,+I+J,1,1,1,Lap10)
* Copy/Transfer/BLock grid pLanes
GSET(C,K2,F,K1,1,40/1,1,50/1,+,0,0,1.0000E-03,INC,1)
IF(:ANS:.EQ.Y) THEN
gclear; gview(z)
ggrid(z,1,1,nx+1,1,ny+1,4,0); ggrid(z,1,1,11,1,31,13,0)
ggrid(z,1,31,nx+1,1,31,13,0); gdraw
ENDIF
**********
NONORT = T
* X-cyclic boundaries swItched
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
* Solved variables list
STORE(ENUT,LEN1,PCOR,PRPS)
SOLVE(P1 ,U1 ,V1 )
* Stored variables list
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
TURMOD(KEMODL)
************************************************************
Group 9. Properties
RHO1 = 1.189E+00; CP1 = 1.005E+03; ENUL = 1.000E-05
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(KE ) = 1.152E-01 ;FIINIT(EP ) = 1.194E+01
CONPOR(B1 , 0.00,VOLUME,-#1,-#2,-#1,-#2,-#1,-#1)
CONPOR(B2 , 0.00,VOLUME,-#4,-#5,-#1,-#2,-#1,-#1)
INIADD = F
************************************************************
Group 13. Boundary & Special Sources
INLET (BFCIN,NORTH ,#1,#5,#3,#3,#1,#1,1,1)
VALUE (BFCIN,P1 , GRND1 ); VALUE (BFCIN,U1 , GRND1 )
VALUE (BFCIN,V1 , GRND1 ); VALUE (BFCIN,KE , 1.152E-01)
VALUE (BFCIN,EP , 1.194E+01); VALUE (BFCIN,VCRT,-2.400E+01)
VALUE (BFCIN,UCRT, 9.300E+01)
OUTLET(OUT ,SOUTH ,#3,#3,#1,#1,#1,#1,1,1)
VALUE (OUT ,P1 , 0.000E+00)
VALUE (OUT ,KE , 1.152E-01)
VALUE (OUT ,EP , 1.194E+01)
BFCA = 1.189E+00
XCYCLE = T
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 300
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1, LINRLX,1.0)
RELAX(U1 ,FALSDT, 1.E-04); RELAX(V1 ,FALSDT, 1.e-04)
RELAX(KE ,LINRLX, 2.000E-01); RELAX(EP ,LINRLX, 2.000E-01)
KELIN = 3
************************************************************
Group 19. EARTH Calls To GROUND Station
GENK = T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 20 ;IYMON = 25 ;IZMON = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
ITABL = 1 ; NXPRIN = 1; NYPRIN = 1