TALK=T;RUN(1,1)
PHOTON USE
p
gr ou y 1
MSG Velocity vectors
vec y 1 sh
msg
msg Press return to plot pressure contours
pause
cont p1 y 1 fil;.01
msg
msg Press return to plot streamlines
pause
cont cl; red
stream sh
y 1
y
posit
3
0.32683E+03 0.18303E+04 CR
0.34606E+03 0.16414E+04 CR
0.33965E+03 0.14590E+04 CR
0.32683E+03 0.12636E+04 CR
0.32042E+03 0.10878E+04 CR
0.94845E+03 0.95098E+03 t
exit
msg Type e to End
ENDUSE
DISPLAY
This case represents the flow between two turbine blades in a
cascade. The blade profile information was supplied in an
arbitrary coordinate system and the coding here transforms
this into two surfaces set at the appropriate pitch and stagger
angle. The non-orthogonal sheared grid has a portion upstream
of the blades inclined at the gas inlet angle and a downstream
portion inclined at the estimated outlet angle. Cyclic boundary
conditions are used in these regions.
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(FLOW THROUGH A REACTION TURBINE: B525)
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 9
XULAST =1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 1
YVLAST =1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 28
ZWLAST =1.
************************************************************
Group 6. Body-Fitted Coordinates
BFC = T ;NONORT = T
NCRT = -1
RSTGEO = F ;SAVGEO = F
UUP = F ;VUP = F ;WUP = F
NGEOM =CHAM
NAMXYZ =CHAM
ANGMIN =20.
DOMAIN( 1, 10, 1, 2, 6, 23)
* Set fixed sub-domain
FIXDOM( 1, 0, 0, 0, 0, 0, 0)
FIXDOM( 2, 0, 0, 0, 0, 0, 0)
FIXDOM( 3, 0, 0, 0, 0, 0, 0)
FIXDOM( 4, 0, 0, 0, 0, 0, 0)
FIXDOM( 5, 0, 0, 0, 0, 0, 0)
FIXDOM( 6, 0, 0, 0, 0, 0, 0)
FIXDOM( 7, 0, 0, 0, 0, 0, 0)
FIXDOM( 8, 0, 0, 0, 0, 0, 0)
FIXDOM( 9, 0, 0, 0, 0, 0, 0)
FIXDOM( 10, 0, 0, 0, 0, 0, 0)
LIJ = F ;LJK = T ;LIK = F
* X-cyclic boundaries switched
XCYIZ( 1, T ) ; XCYIZ( 2, T ) ; XCYIZ( 3, T )
XCYIZ( 4, T ) ; XCYIZ( 5, T ) ; XCYIZ( 24, T )
XCYIZ( 25, T ) ; XCYIZ( 26, T ) ; XCYIZ( 27, T )
XCYIZ( 28, T )
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(7)=W1 ;NAME(147)=RHO1
NAME(148)=WCRT ;NAME(149)=VCRT
NAME(150)=UCRT
* 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(W1,Y,Y,N,N,N,Y)
SOLUTN(RHO1,Y,N,N,N,N,Y)
SOLUTN(WCRT,Y,N,N,N,N,N)
SOLUTN(VCRT,Y,N,N,N,N,N)
SOLUTN(UCRT,Y,N,N,N,N,N)
DEN1 = 147
************************************************************
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(W1,Y,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWRH1 = T
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 = GRND3 ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =1.0E+05
DVO1DT =0. ;DRH1DP = GRND3
RHO1A =2.682695E-04 ;RHO1B =0.714286 ;RHO1C =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =1.0E-10 ;ENUT =0.
PRNDTL(U1)=1. ;PRNDTL(W1)=1.
PRT(U1)=1. ;PRT(W1)=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(W1)=190. ;FIINIT(RHO1)=1.35
FIINIT(WCRT)=1.0E-10 ;FIINIT(VCRT)=1.0E-10
FIINIT(UCRT)=1.0E-10
No PATCHes yet used for this Group
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(INLET ,LOW , 1, 9, 1, 1, 1, 1, 1, 1)
COVAL(INLET ,P1 , FIXFLU ,256.5 )
COVAL(INLET ,U1 ,0. ,-99.533348 )
COVAL(INLET ,W1 ,0. ,190. )
PATCH(OUTLET ,HIGH , 1, 9, 1, 1, 28, 28, 1, 1)
COVAL(OUTLET ,P1 ,1000. ,0. )
COVAL(OUTLET ,U1 ,0. ,0. )
COVAL(OUTLET ,W1 ,0. ,0. )
XCYCLE = T
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 500 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =0.1
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(U1)=10
LITER(W1)=10
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(W1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,0.3)
RELAX(U1,FALSDT,5.0E-04)
RELAX(W1,FALSDT,5.0E-04)
RELAX(RHO1,LINRLX,0.3)
RELAX(WCRT,LINRLX,1.)
RELAX(VCRT,LINRLX,1.)
RELAX(UCRT,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(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=-1.0E+10
VARMAX(WCRT)=1.0E+10 ;VARMIN(WCRT)=-1.0E+10
VARMAX(VCRT)=1.0E+10 ;VARMIN(VCRT)=-1.0E+10
VARMAX(UCRT)=1.0E+10 ;VARMIN(UCRT)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,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(W1,Y,N,Y,Y,Y,Y)
OUTPUT(RHO1,N,N,N,N,N,N)
OUTPUT(WCRT,Y,N,Y,N,N,N)
OUTPUT(VCRT,Y,N,Y,N,N,N)
OUTPUT(UCRT,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 3 ;IYMON = 1 ;IZMON = 10
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NXPRIN = 3 ;IXPRF = 1 ;IXPRL = 10000
NZPRIN = -1 ;IZPRF = 4 ;IZPRL = 24
XZPR = F ;YZPR = F
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
PATCH(XZ ,CONTUR, 1, 9, 1, 1, 1, 28, 1, 1)
PLOT(XZ ,P1 ,0. ,20. )
PLOT(XZ ,W1 ,0. ,20. )
PATCH(PRESSURE,PROFIL, 9, 9, 1, 1, 1, 28, 1, 1)
PLOT(PRESSURE,P1 ,0. ,0. )
PATCH(SUCTION ,PROFIL, 1, 1, 1, 1, 1, 28, 1, 1)
PLOT(SUCTION ,P1 ,0. ,0. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP