TALK=T;RUN(1,1)
PHOTON USE
ext;;;;
gr ou z 1; use patgeo
msg boundary condition patches. Press RETURN
pause
msg temperature contours. Press RETURN
con tmp1 y m fi;0.001;con tmp1 x m fi;0.001;con tmp1 z m fi;0.001
gr ou y m;gr ou z m
msg press RETURN for view x
pause; con off; view x; con tmp1 x m fi;0.001
con tmp1 y m fi;0.001; gr ou y m; gr ou x m; gr ou z 1
msg This is view x. Press RETURN for view y
pause; con off; view y; con tmp1 y m fi;0.001
msg This is view y. Press RETURN for view z
pause; con off; view z; con tmp1 z m fi;0.001; gr ou z m
msg This is view z. Press RETURN for view x and velocity vectors
pause; con off; view x
msg velocity vectors. Type menu for menu and further possibilities
msg Press e to end. Otherwise enter photon-readable commands
vec x m sh; gr ou y m; gr ou x m; gr ou z 1
ENDUSE
DISPLAY
READQ1_BEGIN
The following text is provided as an example of what, by use of
the readq1 command, can be transmitted to EARTH for writing near
the top of the RESULT file. Such text should leave columns 1 and 2
blank; and it should not extend beyond the 68th column.
****************************************************************
* This library case dates from the earliest days of PHOENICS, *
* when Professor WU Chung-Hua ('turbomachinery Wu', who had *
* returned to China from the USA) visited CHAM in 1982. *
* *
* The configuration of the combustion chamber, and its being *
* supplied with premixed fuel vapour and air, was proposed by *
* Professor Wu's accompanying assistant. *
****************************************************************
READQ1_END
The shape of the combustion chamber is as shown.
Pre-mixed fuel vapour and air enter near the axis on left.
Secondary and dilution air enter through holes in outer wall.
A 36-degree sector is simulated.
secondary ox. inlets dilution inlet
____________ 1,2 _______________ ________
blocked /
region /
/
/ outlet
|
|_____ _
fuel-ox. inlet| Symmetry axis
-- - -- -| -- - -- - -- - -- - -- - -- - -- - -- -
The flow is turbulent; the Simple Chemical Reaction Scheme is
used; and the reaction-rate is physically controlled by means of
the Eddy-Breakup Model.
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Idealised Gas-Turbine Combustion Chamber)
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
NX = 6
XULAST =0.628319
Method of pairs used for grid setting.
XFRAC(1)=-6. ;XFRAC(2)=0.166666
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 10
YVLAST =1.
Method of pairs used for grid setting.
YFRAC(1)=-5. ;YFRAC(3)=1.
YFRAC(5)=1. ;YFRAC(7)=3.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 13
ZWLAST =1.
Method of pairs used for grid setting.
ZFRAC(1)=-2. ;ZFRAC(3)=1.
ZFRAC(5)=1. ;ZFRAC(7)=1.
ZFRAC(9)=1. ;ZFRAC(11)=7.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(7)=W1
NAME(12)=KE ;NAME(13)=EP
NAME(14)=H1 ;NAME(136)=PRPS
NAME(137)=DEGF ;NAME(138)=RRAT
NAME(139)=EPKE ;NAME(140)=HPOR
NAME(141)=NPOR ;NAME(142)=EPOR
NAME(143)=VPOR ;NAME(144)=ENUT
NAME(145)=RHO1 ;NAME(146)=TMP1
NAME(147)=PROD ;NAME(148)=OXID
NAME(149)=FUEL ;NAME(150)=MIXF
* 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(W1,Y,Y,N,N,N,Y)
SOLUTN(KE,Y,Y,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
SOLUTN(H1,Y,Y,N,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,N)
SOLUTN(DEGF,Y,N,N,N,N,N)
SOLUTN(RRAT,Y,N,N,N,N,N)
SOLUTN(EPKE,Y,N,N,N,N,N)
SOLUTN(HPOR,Y,N,N,N,N,N)
SOLUTN(NPOR,Y,N,N,N,N,N)
SOLUTN(EPOR,Y,N,N,N,N,N)
SOLUTN(VPOR,Y,N,N,N,N,N)
SOLUTN(ENUT,Y,N,N,N,N,Y)
SOLUTN(RHO1,Y,N,N,N,N,Y)
SOLUTN(TMP1,Y,N,N,N,N,Y)
SOLUTN(PROD,Y,N,N,N,N,Y)
SOLUTN(OXID,Y,N,N,N,N,Y)
SOLUTN(FUEL,Y,Y,N,N,N,Y)
SOLUTN(MIXF,Y,Y,N,N,N,Y)
DEN1 = 145
VIST = 144
EPOR = 142 ;HPOR = 140 ;NPOR = 141 ;VPOR = 143
TEMP1 = 146
PRPS = 136
************************************************************
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)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(KE,N,Y,Y,Y,Y,N)
TERMS(EP,N,Y,Y,Y,Y,N)
TERMS(H1,N,Y,Y,N,Y,N)
TERMS(FUEL,N,Y,Y,Y,N,Y)
TERMS(MIXF,N,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWRH1 = T
NEWENT = T
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
DENPCO = T
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 = GRND6 ;TMP1 = GRND8 ;EL1 = GRND4
TSURR =0. ;TEMP0 =0. ;PRESS0 =8.0E+05
DVO1DT =0. ;DRH1DP =0.
RHO1A =16. ;RHO1B =29. ;RHO1C =28.
TMP1A =0. ;TMP1B =0. ;TMP1C =0.
TMP2A =0.054825 ;TMP2B =4.9E+07 ;TMP2C =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0. ;EL1B =0. ;EL1C =0.
ENUL =9.983361E-06 ;ENUT = GRND3
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(W1)=1. ;PRNDTL(KE)=1.
PRNDTL(EP)=1. ;PRNDTL(H1)=1.
PRNDTL(FUEL)=1. ;PRNDTL(MIXF)=1.
PRT(U1)=1. ;PRT(V1)=1.
PRT(W1)=1. ;PRT(KE)=1.
PRT(EP)=1.314 ;PRT(H1)=1.
PRT(FUEL)=1. ;PRT(MIXF)=1.
CP1 = GRND10 ;CP2 =1.
CP1A =1500. ;CP1B =1500. ;CP1C =1500. ;CP1D =0.
************************************************************
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(W1)=10.
FIINIT(KE)=11.25 ;FIINIT(EP)=353.255707
FIINIT(H1)=1.1595E+06 ;FIINIT(PRPS)=-1.
FIINIT(DEGF)=1.0E-10 ;FIINIT(RRAT)=1.0E-10
FIINIT(EPKE)=1.0E-10 ;FIINIT(HPOR)=1.
FIINIT(NPOR)=1. ;FIINIT(EPOR)=1.
FIINIT(VPOR)=1. ;FIINIT(ENUT)=1.0E-10
FIINIT(RHO1)=3.606 ;FIINIT(TMP1)=773.
FIINIT(PROD)=1.0E-10 ;FIINIT(OXID)=1.0E-10
FIINIT(FUEL)=0.054825 ;FIINIT(MIXF)=0.054825
PATCH(CMP1 ,INIVAL, 1, 6, 1, 2, 1, 3, 1, 1)
INIT(CMP1 ,PRPS,0. ,199. )
PATCH(CMP2 ,INIVAL, 1, 5, 7, 7, 1, 1, 1, 1)
INIT(CMP2 ,EPOR,0. ,0.709 )
PATCH(CMP3 ,INIVAL, 1, 6, 7, 7, 1, 1, 1, 1)
INIT(CMP3 ,VPOR,0. ,0.709 )
PATCH(CMP4 ,INIVAL, 1, 6, 7, 7, 1, 1, 1, 1)
INIT(CMP4 ,NPOR,0. ,0.418 )
PATCH(CMP5 ,INIVAL, 1, 6, 7, 7, 1, 1, 1, 1)
INIT(CMP5 ,HPOR,0. ,1. )
PATCH(CMP6 ,INIVAL, 1, 5, 8, 8, 1, 1, 1, 1)
INIT(CMP6 ,EPOR,0. ,0.185 )
PATCH(CMP7 ,INIVAL, 1, 6, 8, 8, 1, 1, 1, 1)
INIT(CMP7 ,VPOR,0. ,0.185 )
PATCH(CMP8 ,INIVAL, 1, 6, 8, 8, 1, 1, 1, 1)
INIT(CMP8 ,V1 ,0. ,0. )
INIT(CMP8 ,NPOR,0. ,0. )
PATCH(CMP9 ,INIVAL, 1, 6, 8, 8, 1, 1, 1, 1)
INIT(CMP9 ,HPOR,0. ,0.877 )
PATCH(CMP10 ,INIVAL, 1, 5, 9, 10, 1, 1, 1, 1)
INIT(CMP10 ,U1 ,0. ,0. )
INIT(CMP10 ,EPOR,0. ,0. )
PATCH(CMP11 ,INIVAL, 1, 6, 9, 10, 1, 1, 1, 1)
INIT(CMP11 ,PRPS,0. ,199. )
PATCH(CMP12 ,INIVAL, 1, 6, 9, 9, 1, 1, 1, 1)
INIT(CMP12 ,V1 ,0. ,0. )
INIT(CMP12 ,NPOR,0. ,0. )
PATCH(CMP13 ,INIVAL, 1, 6, 9, 10, 1, 1, 1, 1)
INIT(CMP13 ,W1 ,0. ,0. )
INIT(CMP13 ,HPOR,0. ,0. )
PATCH(CMP14 ,INIVAL, 1, 5, 8, 8, 2, 2, 1, 1)
INIT(CMP14 ,EPOR,0. ,0.997 )
PATCH(CMP15 ,INIVAL, 1, 6, 8, 8, 2, 2, 1, 1)
INIT(CMP15 ,VPOR,0. ,0.997 )
PATCH(CMP16 ,INIVAL, 1, 6, 8, 8, 2, 2, 1, 1)
INIT(CMP16 ,NPOR,0. ,0.946 )
PATCH(CMP17 ,INIVAL, 1, 6, 8, 8, 2, 2, 1, 1)
INIT(CMP17 ,HPOR,0. ,1. )
PATCH(CMP18 ,INIVAL, 1, 5, 9, 9, 2, 2, 1, 1)
INIT(CMP18 ,EPOR,0. ,0.706 )
PATCH(CMP19 ,INIVAL, 1, 6, 9, 9, 2, 2, 1, 1)
INIT(CMP19 ,VPOR,0. ,0.706 )
PATCH(CMP20 ,INIVAL, 1, 6, 9, 9, 2, 2, 1, 1)
INIT(CMP20 ,NPOR,0. ,0.473 )
PATCH(CMP21 ,INIVAL, 1, 6, 9, 9, 2, 2, 1, 1)
INIT(CMP21 ,HPOR,0. ,1. )
PATCH(CMP22 ,INIVAL, 1, 5, 10, 10, 2, 2, 1, 1)
INIT(CMP22 ,EPOR,0. ,0.236 )
PATCH(CMP23 ,INIVAL, 1, 6, 10, 10, 2, 2, 1, 1)
INIT(CMP23 ,VPOR,0. ,0.236 )
PATCH(CMP24 ,INIVAL, 1, 6, 10, 10, 2, 2, 1, 1)
PATCH(CMP25 ,INIVAL, 1, 6, 10, 10, 2, 2, 1, 1)
INIT(CMP25 ,HPOR,0. ,1. )
PATCH(WALL4 ,INIVAL, 1, 6, 2, 2, 1, 3, 1, 1)
INIT(WALL4 ,V1 ,0. ,0. )
INIT(WALL4 ,NPOR,0. ,0. )
PATCH(WALL6 ,INIVAL, 1, 6, 1, 2, 3, 3, 1, 1)
INIT(WALL6 ,W1 ,0. ,0. )
INIT(WALL6 ,HPOR,0. ,0. )
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, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(KESOURCE,KE , GRND4 , GRND4 )
COVAL(KESOURCE,EP , GRND4 , GRND4 )
PATCH(FOIN ,SOUTH , 1, 6, 3, 3, 3, 3, 1, 1)
COVAL(FOIN ,P1 , FIXFLU ,360.600006 )
COVAL(FOIN ,U1 ,0. ,0. )
COVAL(FOIN ,V1 ,0. ,100. )
COVAL(FOIN ,W1 ,0. ,0. )
COVAL(FOIN ,KE ,0. ,0.25 )
COVAL(FOIN ,EP ,0. ,58.752605 )
COVAL(FOIN ,H1 ,0. ,3.845904E+06 )
COVAL(FOIN ,FUEL,0. ,0.054825 )
COVAL(FOIN ,MIXF,0. ,0.054825 )
PATCH(SOIN1 ,CELL , 1, 1, 10, 10, 7, 7, 1, 1)
COVAL(SOIN1 ,P1 , FIXFLU ,8.157E-03 )
COVAL(SOIN1 ,U1 ,0. ,0. )
COVAL(SOIN1 ,V1 ,0. ,-40. )
COVAL(SOIN1 ,W1 ,0. ,0. )
COVAL(SOIN1 ,KE ,0. ,0.123 )
COVAL(SOIN1 ,EP ,0. ,8.429416 )
COVAL(SOIN1 ,H1 ,0. ,1.1595E+06 )
COVAL(SOIN1 ,FUEL,0. ,0. )
COVAL(SOIN1 ,MIXF,0. ,0. )
PATCH(SOIN2 ,CELL , 4, 4, 10, 10, 7, 7, 1, 1)
COVAL(SOIN2 ,P1 , FIXFLU ,8.157E-03 )
COVAL(SOIN2 ,U1 ,0. ,0. )
COVAL(SOIN2 ,V1 ,0. ,-40. )
COVAL(SOIN2 ,W1 ,0. ,0. )
COVAL(SOIN2 ,KE ,0. ,0.123 )
COVAL(SOIN2 ,EP ,0. ,8.429416 )
COVAL(SOIN2 ,H1 ,0. ,1.1595E+06 )
COVAL(SOIN2 ,FUEL,0. ,0. )
COVAL(SOIN2 ,MIXF,0. ,0. )
PATCH(DILUIN ,CELL , 1, 1, 10, 10, 10, 10, 1, 1)
COVAL(DILUIN ,P1 , FIXFLU ,0.0145 )
COVAL(DILUIN ,U1 ,0. ,0. )
COVAL(DILUIN ,V1 ,0. ,-40. )
COVAL(DILUIN ,W1 ,0. ,0. )
COVAL(DILUIN ,KE ,0. ,0.123 )
COVAL(DILUIN ,EP ,0. ,6.299674 )
COVAL(DILUIN ,H1 ,0. ,1.1595E+06 )
COVAL(DILUIN ,FUEL,0. ,0. )
COVAL(DILUIN ,MIXF,0. ,0. )
PATCH(OUTLET ,HIGH , 1, 6, 1, 10, 13, 13, 1, 1)
COVAL(OUTLET ,P1 ,1000. ,0. )
COVAL(OUTLET ,U1 ,0. ,0. )
COVAL(OUTLET ,V1 ,0. ,0. )
COVAL(OUTLET ,W1 ,0. ,0. )
COVAL(OUTLET ,KE ,0. , SAME )
COVAL(OUTLET ,EP ,0. , SAME )
COVAL(OUTLET ,H1 ,0. , SAME )
COVAL(OUTLET ,FUEL,0. , SAME )
COVAL(OUTLET ,MIXF,0. , SAME )
PATCH(WALL1 ,NWALL , 1, 6, 10, 10, 3, 6, 1, 1)
COVAL(WALL1 ,U1 , GRND2 ,0. )
COVAL(WALL1 ,W1 , GRND2 ,0. )
COVAL(WALL1 ,KE , GRND2 , GRND2 )
COVAL(WALL1 ,EP , GRND2 , GRND2 )
PATCH(WALL2 ,NWALL , 1, 6, 10, 10, 8, 9, 1, 1)
COVAL(WALL2 ,U1 , GRND2 ,0. )
COVAL(WALL2 ,W1 , GRND2 ,0. )
COVAL(WALL2 ,KE , GRND2 , GRND2 )
COVAL(WALL2 ,EP , GRND2 , GRND2 )
PATCH(WALL3 ,NWALL , 1, 6, 10, 10, 11, 13, 1, 1)
COVAL(WALL3 ,U1 , GRND2 ,0. )
COVAL(WALL3 ,W1 , GRND2 ,0. )
COVAL(WALL3 ,KE , GRND2 , GRND2 )
COVAL(WALL3 ,EP , GRND2 , GRND2 )
PATCH(WALL4-NW,NWALL , 1, 6, 2, 2, 1, 3, 1, 1)
COVAL(WALL4-NW,U1 , GRND2 ,0. )
COVAL(WALL4-NW,W1 , GRND2 ,0. )
COVAL(WALL4-NW,KE , GRND2 , GRND2 )
COVAL(WALL4-NW,EP , GRND2 , GRND2 )
PATCH(WALL4-SW,SWALL , 1, 6, 3, 3, 1, 3, 1, 1)
COVAL(WALL4-SW,U1 , GRND2 ,0. )
COVAL(WALL4-SW,W1 , GRND2 ,0. )
COVAL(WALL4-SW,KE , GRND2 , GRND2 )
COVAL(WALL4-SW,EP , GRND2 , GRND2 )
PATCH(WALL5 ,LWALL , 1, 6, 3, 6, 1, 1, 1, 1)
COVAL(WALL5 ,U1 , GRND2 ,0. )
COVAL(WALL5 ,V1 , GRND2 ,0. )
COVAL(WALL5 ,KE , GRND2 , GRND2 )
COVAL(WALL5 ,EP , GRND2 , GRND2 )
PATCH(WALL6-HW,HWALL , 1, 6, 1, 2, 3, 3, 1, 1)
COVAL(WALL6-HW,U1 , GRND2 ,0. )
COVAL(WALL6-HW,V1 , GRND2 ,0. )
COVAL(WALL6-HW,KE , GRND2 , GRND2 )
COVAL(WALL6-HW,EP , GRND2 , GRND2 )
PATCH(WALL6-LW,LWALL , 1, 6, 1, 2, 4, 4, 1, 1)
COVAL(WALL6-LW,U1 , GRND2 ,0. )
COVAL(WALL6-LW,V1 , GRND2 ,0. )
COVAL(WALL6-LW,KE , GRND2 , GRND2 )
COVAL(WALL6-LW,EP , GRND2 , GRND2 )
PATCH(WHOLE ,PHASEM, 1, 6, 1, 10, 1, 13, 1, 1)
COVAL(WHOLE ,FUEL,In-Form:source - see Grp 19)
XCYCLE = F
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
SARAH =1.
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(U1)=20
LITER(V1)=20 ;LITER(W1)=20
LITER(KE)=20 ;LITER(EP)=20
LITER(H1)=20 ;LITER(FUEL)=20
LITER(MIXF)=20
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03
ENDIT(H1)=1.0E-03 ;ENDIT(FUEL)=1.0E-03
ENDIT(MIXF)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(U1,FALSDT,0.02)
RELAX(V1,FALSDT,0.02)
RELAX(W1,FALSDT,0.02)
RELAX(KE,LINRLX,0.5)
RELAX(EP,LINRLX,0.5)
RELAX(H1,FALSDT,10.)
RELAX(PRPS,LINRLX,1.)
RELAX(DEGF,LINRLX,1.)
RELAX(RRAT,LINRLX,1.)
RELAX(EPKE,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
RELAX(RHO1,LINRLX,1.)
RELAX(TMP1,LINRLX,1.)
RELAX(PROD,LINRLX,1.)
RELAX(OXID,LINRLX,1.)
RELAX(FUEL,FALSDT,1.0E+09)
RELAX(MIXF,FALSDT,1.0E+09)
KELIN = 3
OVRRLX =0.
EXPERT = F ;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(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
VARMAX(DEGF)=1.0E+10 ;VARMIN(DEGF)=-1.0E+10
VARMAX(RRAT)=1.0E+10 ;VARMIN(RRAT)=-1.0E+10
VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
VARMAX(HPOR)=1.0E+10 ;VARMIN(HPOR)=-1.0E+10
VARMAX(NPOR)=1.0E+10 ;VARMIN(NPOR)=-1.0E+10
VARMAX(EPOR)=1.0E+10 ;VARMIN(EPOR)=-1.0E+10
VARMAX(VPOR)=1.0E+10 ;VARMIN(VPOR)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=-1.0E+10
VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=773.
VARMAX(PROD)=1.0E+10 ;VARMIN(PROD)=0.
VARMAX(OXID)=1.0E+10 ;VARMIN(OXID)=0.
VARMAX(FUEL)=1.0E+10 ;VARMIN(FUEL)=0.
VARMAX(MIXF)=1.0E+10 ;VARMIN(MIXF)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
GENK = T
READQ1 = T
PARSOL = F
ISG62 = 1
CHSOA =0.054825 ;CHSOB =1.
SPEDAT(SET,PROPERTY,TMP1,C,=MAX(773.&(H1-4.9E+07*FUEL)/(1500.*FUE$)
SPEDAT(SET,PROPERTY,TMP1,C,L+1500.*OXID+1500.*PROD))!IMAT<100)
SPEDAT(SET,PROPERTY,RHO1,C,=(P1+8.0E+05)*1.0/(FUEL/16.+OXID/29.+P$)
SPEDAT(SET,PROPERTY,RHO1,C,ROD/28.)/(TMP1*8313.4))
SPEDAT(SET,STORED,PROD,C,=(MIXF-FUEL)*(1+17.24)!ZSLFIN)
SPEDAT(SET,STORED,OXID,C,=1-FUEL-PROD!ZSLFIN)
SPEDAT(SET,SOURCE,FUEL!WHOLE,C,=COVAL(1.*EPKE&0.0))
SPEDAT(SET,STORED,RRAT,C,=1.*EPKE*FUEL!ZSLFIN)
SPEDAT(SET,LONGNAME,TMP1,C,absolute_temperature_of_the_gas_Kelvin)
SPEDAT(SET,STORED,DEGF,C,=(9./5.)*(TMP1-273)!ZSLFIN)
SPEDAT(SET,LONGNAME,DEGF,C,degrees_Fahrenheit)
SPEDAT(SET,LONGNAME,RRAT,C,rate_of_consumption_of_fuel_kg/m^3sec)
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(W1,Y,N,Y,Y,Y,Y)
OUTPUT(KE,Y,N,Y,Y,Y,Y)
OUTPUT(EP,Y,N,Y,Y,Y,Y)
OUTPUT(H1,Y,N,Y,Y,Y,Y)
OUTPUT(PRPS,Y,N,Y,N,N,N)
OUTPUT(DEGF,Y,N,Y,N,N,N)
OUTPUT(RRAT,Y,N,Y,N,N,N)
OUTPUT(EPKE,Y,N,Y,N,N,N)
OUTPUT(HPOR,Y,N,Y,N,N,N)
OUTPUT(NPOR,Y,N,Y,N,N,N)
OUTPUT(EPOR,Y,N,Y,N,N,N)
OUTPUT(VPOR,Y,N,Y,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,N,N)
OUTPUT(RHO1,Y,N,Y,N,N,N)
OUTPUT(TMP1,Y,N,Y,N,N,N)
OUTPUT(PROD,Y,N,Y,N,N,N)
OUTPUT(OXID,Y,N,Y,N,N,N)
OUTPUT(FUEL,Y,N,Y,Y,Y,Y)
OUTPUT(MIXF,Y,N,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 3 ;IYMON = 5 ;IZMON = 5
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -5
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100 ;NUMCLS = 5
NXPRIN = -1 ;IXPRF = 3 ;IXPRL = 3
NYPRIN = 2 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = 2 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = T
IPLTF = 2 ;IPLTL = 100 ;NPLT = 5
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 40
PATCH(IZ4 ,CONTUR, 1, 6, 1, 10, 4, 4, 1, 1)
PLOT(IZ4 ,TMP1,1. ,10. )
PATCH(XSECIN1 ,CONTUR, 1, 1, 1, 10, 1, 13, 1, 1)
PLOT(XSECIN1 ,W1 ,0. ,10. )
PLOT(XSECIN1 ,TMP1,0. ,10. )
PLOT(XSECIN1 ,FUEL,0. ,10. )
PLOT(XSECIN1 ,MIXF,0. ,10. )
PATCH(XSECIN2 ,CONTUR, 4, 4, 1, 10, 1, 13, 1, 1)
PLOT(XSECIN2 ,W1 ,0. ,10. )
PLOT(XSECIN2 ,TMP1,0. ,10. )
PLOT(XSECIN2 ,FUEL,0. ,10. )
PLOT(XSECIN2 ,MIXF,0. ,10. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP