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
PLANTBEGIN
STORE(ABSR)
The above command provides the 3D-storage for aborption
coefficient.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
RG(1)=0.2;RG(2)=0.3
The use of RG's, as above, in PLANTed expression is
especially convenient for parametric studies.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
ABSR=RG(1)*FUEL+RG(2)*PROD
The above statement makes the absorbtion coefficient to
be the function of the fuel and combustion products
mass fractions. It is done for the whole domain at the
end of each iz-slab.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
PRNDTL(RADX)=-GRND
LAMPR(RADX)=1/(:SCAT:+ABSR)
PRNDTL(RADZ)=-GRND
LAMPR(RADZ)=1/(:SCAT:+ABSR)
The above two statements followed by their pointers
instruct PLANT to make the diffusivities of radiation
fluxes in X- and Z-directions reciprocal to the sum of
scattering and absorption coefficients.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
PRNDTL(RADY)=-GRND
LAMPR(RADY)=1/(:SCAT:+ABSR)
IF(CARTES)
LAMPR(RADY)=RV2D**2/(1.+RV2D*(:SCAT:+ABSR))
IF(.NOT.CARTES)
PLANT makes the provision for the calculations of
diffusivity of lateral radiation flux. The first
statement of the above two is applied for cartesian
coordinates, while the second will be active for
cylindrical polar ones.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
PLANTEND
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Idealised Gas-Turbine Combustion Chamber)
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 492
************************************************************
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(135)=PRPS
NAME(136)=ABSR ;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(ABSR,Y,N,N,N,N,Y)
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 = 135
************************************************************
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(ABSR)=1.0E-10 ;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 = 40 ;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(ABSR,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(ABSR)=1.0E+10 ;VARMIN(ABSR)=-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
NAMSAT =MOSG
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)
RG( 1) =0.2
RG( 2) =0.3
************************************************************
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(ABSR,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