TALK=T;RUN(1,1)
photon use
p;;;;;
use patgeo; msg Velocity vectors. Press RETURN for pressure
vec z 1 sh;pause
msg Pressure contours. Press RETURN for streamlines
con p1 z 1 fi;0.001
pause; msg streamlines
stre 2d z 1 ;int 30
enduse
DISPLAY
This Q1, which sets up a steady 2D flow through a passage created
by blockages, has been set up by means of the general input menu.
Enter the menu (Type M, return, R) to see the geometry and
settings, and to change them if desired.
On re-entering command mode, another means of interpreting the
settings can be activated by typing: L(20). This is always
available.
PHOTON USE commands have been included in the Q1.
ENDDIS
DISPLAY
In this case, the false-time-step rlaxation factors set by
the General Menu for case 293 are first multiplied by 1.E-3.
Then LSWEEP is raised to 200.
This proves to be not enough to secure convergence.
ENDDIS
DISPLAY
Now EXPERT is called in, by L(30), to choose the optimum
relaxation. Convergence is obtained after 135 sweeps.
The relaxation factors chosen by the menu (in case 293) prove
to have been nearly optimal.
ENDDIS
DISPLAY
Library case 30 loaded, to activate EXPERT
expsol f .... expert is not used in the solver
expdtf t .... expert is used to adjust false time step
indtf 2 ... adjust the false time step of all velocities
inres 3 ... monitor the residuals of variable 3
ifrsts 20 ... first adjustment sweep
ifrequ 10 ... adjustments will be made every ifrequ sweeps
facdec 0.5 .. decreases will be by a factor of facdec
facinc 2.0 .. increases will be by a factor of facinc
exprin t .... print adjusted false time steps in RESULT file
expend ..... end of the inputs to expert via Q1
xprtend
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Simple Labyrinth; Expert Relaxation )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 294
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 30
XULAST =2.5
XFRAC(1)=0.033333 ;XFRAC(2)=0.066667
XFRAC(3)=0.1 ;XFRAC(4)=0.133333
XFRAC(5)=0.166667 ;XFRAC(6)=0.2
XFRAC(7)=0.233333 ;XFRAC(8)=0.266667
XFRAC(9)=0.3 ;XFRAC(10)=0.332
XFRAC(11)=0.364 ;XFRAC(12)=0.396
XFRAC(13)=0.428 ;XFRAC(14)=0.46
XFRAC(15)=0.486667 ;XFRAC(16)=0.513333
XFRAC(17)=0.54 ;XFRAC(18)=0.572
XFRAC(19)=0.604 ;XFRAC(20)=0.636
XFRAC(21)=0.668 ;XFRAC(22)=0.7
XFRAC(23)=0.7375 ;XFRAC(24)=0.775
XFRAC(25)=0.8125 ;XFRAC(26)=0.85
XFRAC(27)=0.8875 ;XFRAC(28)=0.925
XFRAC(29)=0.9625 ;XFRAC(30)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 25
YVLAST =1.5
YFRAC(1)=0.038462 ;YFRAC(2)=0.076923
YFRAC(3)=0.115385 ;YFRAC(4)=0.153846
YFRAC(5)=0.192308 ;YFRAC(6)=0.230769
YFRAC(7)=0.269231 ;YFRAC(8)=0.307692
YFRAC(9)=0.346154 ;YFRAC(10)=0.384615
YFRAC(11)=0.423077 ;YFRAC(12)=0.461538
YFRAC(13)=0.5 ;YFRAC(14)=0.541667
YFRAC(15)=0.583333 ;YFRAC(16)=0.625
YFRAC(17)=0.666667 ;YFRAC(18)=0.708333
YFRAC(19)=0.75 ;YFRAC(20)=0.791667
YFRAC(21)=0.833333 ;YFRAC(22)=0.875
YFRAC(23)=0.916667 ;YFRAC(24)=0.958333
YFRAC(25)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST =1.
ZFRAC(1)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(49)=DEN1
NAME(149)=PRPS
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(P1,Y,Y,Y,N,N,N)
SOLUTN(U1,Y,Y,N,N,N,Y)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(DEN1,Y,N,N,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,N)
DEN1 = 49
PRPS = 149
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
TERMS(P1,Y,Y,Y,N,Y,Y)
TERMS(U1,Y,Y,Y,Y,Y,Y)
TERMS(V1,Y,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 =1.161 ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =1.589E-05 ;ENUT =3.178E-03
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRT(U1)=1. ;PRT(V1)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10
FIINIT(V1)=1.0E-10 ;FIINIT(DEN1)=1.0E-10
FIINIT(PRPS)=-1.
PATCH(B1 ,INIVAL, 1, 9, 1, 17, 1, 1, 1, 1)
INIT(B1 ,PRPS,0. ,199. )
PATCH(B2 ,INIVAL, 15, 17, 14, 25, 1, 1, 1, 1)
INIT(B2 ,PRPS,0. ,199. )
PATCH(B3 ,INIVAL, 23, 30, 1, 17, 1, 1, 1, 1)
INIT(B3 ,PRPS,0. ,199. )
INIADD = F
FSWEEP = 1
NAMFI =CHAM
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources
PATCH(KESOURCE,PHASEM, 1, 30, 1, 25, 1, 1, 1, 1)
PATCH(IN1 ,WEST , 1, 1, 18, 25, 1, 1, 1, 1)
COVAL(IN1 ,P1 , FIXFLU ,5.805 )
COVAL(IN1 ,U1 ,0. ,5. )
COVAL(IN1 ,V1 ,0. ,0. )
PATCH(OUT1 ,EAST , 30, 30, 18, 25, 1, 1, 1, 1)
COVAL(OUT1 ,P1 ,1. ,0. )
PATCH(TWALL ,NWALL , 1, 30, 25, 25, 1, 1, 1, 1)
COVAL(TWALL ,U1 , GRND2 ,0. )
PATCH(BWALL ,SWALL , 1, 30, 1, 1, 1, 1, 1, 1)
COVAL(BWALL ,U1 , GRND2 ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 300 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-02
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(U1)=10
LITER(V1)=10
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(U1,FALSDT,6.0E-03)
RELAX(V1,FALSDT,6.0E-03)
RELAX(DEN1,LINRLX,1.)
RELAX(PRPS,LINRLX,1.)
OVRRLX =0.
EXPERT = T ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
GENK = T
PARSOL = F
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,MATERIAL,199,L,T)
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
OUTPUT(P1,Y,N,Y,Y,Y,Y)
OUTPUT(U1,Y,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(DEN1,Y,N,Y,N,N,N)
OUTPUT(PRPS,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 19 ;IYMON = 23 ;IZMON = 1
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP