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 Type e to End
ENDUSE
DISPLAY
The flow considered is supersonic flow through a cascade of
wedges with inlet Mach number 3.0 and completely supersonic
flow. A leading-edge shock reflects off the pressure surface
and should be exactly cancelled at the upstream corner giving
a uniform parallel flow through the two surfaces. The flow
then expands off the downstream corner and exits through the
blade row where two compression waves are formed at the trailing
edge. Cyclic boundary conditions are applied upstream and
downstream of the cascade. The geometry is as follows:
wall
//////////////////////////////////////////////////////
------------------------------------------------------
---------> zero
---> ---------> ---> pressure
____________________
.'////////////////////`.
---> .'| wall -`. --->
.'| wall wall -`.
___________.'| -`._______
/////////// ///////
^ wall wall
x|
|--->
z
For simplicity, the flow is treated as isentropic. However,
shock theory indicates that there is a significant entropy
change across the shocks for the given approach Mach number
and wedge angle. Therefore, in future work the isentropic
treatment will be replaced with one which allows for entropy
changes across shock fronts.
The exit boundary condition is one of fixed pressure according
to the post-expansion pressure calculated from gas-dynamic
theory; this neglects the presence of trailing-edge shocks.
Strictly, the flow is hyperbolic and so the exit boundary
condition should be modified accordingly.
The system of units used are the same as those used in case 522.
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(SUPERSONIC FLOW THRU WEDGE CASCADE: B523)
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 14
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 10
XULAST =1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 1
YVLAST =1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 28
ZWLAST =5.6
************************************************************
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, 11, 1, 2, 1, 29)
* 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( 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 =0.
DVO1DT =0. ;DRH1DP = GRND3
RHO1A =1. ;RHO1B =0.714286 ;RHO1C =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =0. ;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)=0.027224 ;FIINIT(U1)=1.0E-10
FIINIT(W1)=2.12132 ;FIINIT(RHO1)=0.076226
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, 10, 1, 1, 1, 1, 1, 1)
COVAL(INLET ,P1 , FIXFLU ,0.1617 )
COVAL(INLET ,U1 ,0. ,0. )
COVAL(INLET ,W1 ,0. ,2.12132 )
PATCH(OUTLET ,HIGH , 1, 10, 1, 1, 28, 28, 1, 1)
COVAL(OUTLET ,P1 ,5.0E+04 ,0.0377 )
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 = 100 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-02
************************************************************
Group 16. Terminate Iterations
LITER(P1)=15 ;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.8)
RELAX(U1,FALSDT,0.5)
RELAX(W1,FALSDT,0.5)
RELAX(RHO1,LINRLX,1.)
RELAX(WCRT,LINRLX,1.)
RELAX(VCRT,LINRLX,1.)
RELAX(UCRT,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1. ;VARMIN(P1)=2.722368E-04
VARMAX(U1)=50. ;VARMIN(U1)=-50.
VARMAX(W1)=50. ;VARMIN(W1)=-50.
VARMAX(RHO1)=1. ;VARMIN(RHO1)=7.62263E-03
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,Y,N,Y,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 = 2 ;IYMON = 1 ;IZMON = 9
NPRMON = 100 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100 ;NUMCLS = 5
NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
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(PLOT1 ,PROFIL, 5, 5, 1, 1, 1, 28, 1, 1)
PLOT(PLOT1 ,P1 ,0. ,0. )
PATCH(CASCADE ,CONTUR, 1, 10, 1, 1, 1, 28, 1, 1)
PLOT(CASCADE ,P1 ,0. ,20. )
PLOT(CASCADE ,W1 ,0. ,20. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP