TALK=T;RUN(1,1)
PHOTON USE
p;;;;
view z
msg exit smoke concentration according to single-fluid model
con smos z m fi;0.001
pause;con off;red
msg exit smoke concentration according to multi-fluid model
con smom z m fi;0.001;pause;con off;red;view x
msg x-plane smoke concentration according to single-fluid model
con smos x 1 fi;0.001;pause;con off;red
msg x-plane smoke concentration according to multi-fluid model
con smom x 1 fi;0.001;pause;con off;red
msg x-plane mixture fraction avef according to multi-fluid model
con avef x 1 fi;0.001;pause;con off;red
msg x-plane mixture fraction mixf according to single-fluid model
con mixf x 1 fi;0.001;pause;con off;red
msg x-plane fluid-1 concentration according to multi-fluid model
con f1 x 1 fi;0.001;pause;con off;red
msg x-plane fluid-11 concentration according to multi-fluid model
con f11 x 1 fi;0.001;pause;con off;red
msg x-plane fluid-6 concentration according to multi-fluid model
con f6 x 1 fi;0.001
ENDUSE
DISPLAY
This simulation shows how the multi-fluid model of turbulence
predicts a smoke-production pattern in a gas-turbine combustor
which is different from, and probably more realistic than, that
predicted by a conventional single-fluid model.
It is explained and discussed extensively in a report which is
published on CHAM's website www.cham.co.uk (Click on short-cuts,
then MFM then 1998b).
The geometry is that of library case 492, which is loaded first.
Note the $ in #$492; it ensures that the above PHOTON USE
commands remain in the Q1 file, and therefore take priority over
those which appear in case 493.
The injected fuel-air composition, which is stoichiometric in
case 492, is made twice stoichiometric here.
To reduce computation time, it is best to answer N to the question
about radiation.
ENDDIS
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
DISPLAY
Notes:
(1) four macros are provided in order reduce the repetiveness of
MFM q1s.
(2) they are named mfm, mfm1, mfm2 and mfm3, and are normally
called in that order at the top of the q1 file, which then
appears as:
#mfm
settings of nflr, nflf, etc
#mfm1
settings of
#mfm2
#mfm3
(3) the character variable mfm is declared and set in the
always-loaded core-library macro 014.htm, which, in its turn,
declares and sets the character variables mfm1, mfm2, mfm3.
(4) the values are:
mfm = $L004
mfm1 = $L003
mfm2 = $L002
mfm3 = $L001
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(MFM; 3D Cmbstr; Conmix=5. )
************************************************************
************************************************************
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(2)=6.5E-03
YFRAC(3)=1. ;YFRAC(4)=5.0E-03
YFRAC(5)=1. ;YFRAC(6)=8.0E-03
YFRAC(7)=3. ;YFRAC(8)=6.5E-03
************************************************************
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(115)=PRPS
NAME(116)=PRDM ;NAME(117)=OXDM
NAME(118)=FRFM ;NAME(119)=TMPM
NAME(120)=SMOS ;NAME(121)=SMOM
NAME(122)=F1 ;NAME(123)=F2
NAME(124)=F3 ;NAME(125)=F4
NAME(126)=F5 ;NAME(127)=F6
NAME(128)=F7 ;NAME(129)=F8
NAME(130)=F9 ;NAME(131)=F10
NAME(132)=F11 ;NAME(133)=MNSQ
NAME(134)=AVEF ;NAME(135)=MIXL
NAME(136)=RATE ;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,Y)
SOLUTN(PRDM,Y,N,N,N,N,Y)
SOLUTN(OXDM,Y,N,N,N,N,Y)
SOLUTN(FRFM,Y,N,N,N,N,Y)
SOLUTN(TMPM,Y,N,N,N,N,Y)
SOLUTN(SMOS,Y,Y,N,N,N,Y)
SOLUTN(SMOM,Y,Y,N,N,N,Y)
SOLUTN(F1,Y,Y,N,N,N,Y)
SOLUTN(F2,Y,Y,N,N,N,Y)
SOLUTN(F3,Y,Y,N,N,N,Y)
SOLUTN(F4,Y,Y,N,N,N,Y)
SOLUTN(F5,Y,Y,N,N,N,Y)
SOLUTN(F6,Y,Y,N,N,N,Y)
SOLUTN(F7,Y,Y,N,N,N,Y)
SOLUTN(F8,Y,Y,N,N,N,Y)
SOLUTN(F9,Y,Y,N,N,N,Y)
SOLUTN(F10,Y,Y,N,N,N,Y)
SOLUTN(F11,Y,Y,N,N,N,Y)
SOLUTN(MNSQ,Y,N,N,N,N,Y)
SOLUTN(AVEF,Y,N,N,N,N,Y)
SOLUTN(MIXL,Y,N,N,N,N,N)
SOLUTN(RATE,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,N,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 = 115
************************************************************
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(SMOS,N,Y,Y,Y,Y,Y)
TERMS(SMOM,N,Y,Y,Y,N,Y)
TERMS(F1,N,Y,Y,Y,Y,Y)
TERMS(F2,N,Y,Y,Y,N,Y)
TERMS(F3,N,Y,Y,Y,Y,Y)
TERMS(F4,N,Y,Y,Y,N,Y)
TERMS(F5,N,Y,Y,Y,Y,Y)
TERMS(F6,N,Y,Y,Y,N,Y)
TERMS(F7,N,Y,Y,Y,Y,Y)
TERMS(F8,N,Y,Y,Y,N,Y)
TERMS(F9,N,Y,Y,Y,Y,Y)
TERMS(F10,N,Y,Y,Y,N,Y)
TERMS(F11,N,Y,Y,Y,Y,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 = GRND7 ;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.054825
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(SMOS)=1. ;PRNDTL(SMOM)=1.
PRNDTL(F1)=1. ;PRNDTL(F2)=1.
PRNDTL(F3)=1. ;PRNDTL(F4)=1.
PRNDTL(F5)=1. ;PRNDTL(F6)=1.
PRNDTL(F7)=1. ;PRNDTL(F8)=1.
PRNDTL(F9)=1. ;PRNDTL(F10)=1.
PRNDTL(F11)=1. ;PRNDTL(MIXF)=1.
PRT(U1)=1. ;PRT(V1)=1.
PRT(W1)=1. ;PRT(KE)=1.
PRT(EP)=1.314 ;PRT(H1)=1.
PRT(SMOS)=1. ;PRT(SMOM)=1.
PRT(F1)=1. ;PRT(F2)=1.
PRT(F3)=1. ;PRT(F4)=1.
PRT(F5)=1. ;PRT(F6)=1.
PRT(F7)=1. ;PRT(F8)=1.
PRT(F9)=1. ;PRT(F10)=1.
PRT(F11)=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.0E-10
FIINIT(PRDM)=1.0E-10 ;FIINIT(OXDM)=1.0E-10
FIINIT(FRFM)=1.0E-10 ;FIINIT(TMPM)=1.0E-10
FIINIT(SMOS)=0. ;FIINIT(SMOM)=0.
FIINIT(F1)=0.090909 ;FIINIT(F2)=0.090909
FIINIT(F3)=0.090909 ;FIINIT(F4)=0.090909
FIINIT(F5)=0.090909 ;FIINIT(F6)=0.090909
FIINIT(F7)=0.090909 ;FIINIT(F8)=0.090909
FIINIT(F9)=0.090909 ;FIINIT(F10)=0.090909
FIINIT(F11)=0.090909 ;FIINIT(MNSQ)=1.0E-10
FIINIT(AVEF)=1.0E-10 ;FIINIT(MIXL)=1.0E-10
FIINIT(RATE)=31.400507 ;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.109649
FIINIT(MIXF)=0.109649
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 ,180.300003 )
COVAL(FOIN ,U1 ,0. ,0. )
COVAL(FOIN ,V1 ,0. ,50. )
COVAL(FOIN ,W1 ,0. ,0. )
COVAL(FOIN ,KE ,0. ,0.0625 )
COVAL(FOIN ,EP ,0. ,3.672038 )
COVAL(FOIN ,H1 ,0. ,4.001728E+06 )
COVAL(FOIN ,F11 ,0. ,1. )
COVAL(FOIN ,FUEL,0. ,0.058005 )
COVAL(FOIN ,MIXF,0. ,0.109649 )
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 ,F1 ,0. ,1. )
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 ,F1 ,0. ,1. )
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 ,F1 ,0. ,1. )
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 ,SMOS,0. , SAME )
COVAL(OUTLET ,SMOM,0. , SAME )
COVAL(OUTLET ,F1 ,0. , SAME )
COVAL(OUTLET ,F2 ,0. , SAME )
COVAL(OUTLET ,F3 ,0. , SAME )
COVAL(OUTLET ,F4 ,0. , SAME )
COVAL(OUTLET ,F5 ,0. , SAME )
COVAL(OUTLET ,F6 ,0. , SAME )
COVAL(OUTLET ,F7 ,0. , SAME )
COVAL(OUTLET ,F8 ,0. , SAME )
COVAL(OUTLET ,F9 ,0. , SAME )
COVAL(OUTLET ,F10 ,0. , SAME )
COVAL(OUTLET ,F11 ,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)
PATCH(SMOKE ,PHASEM, 1, 6, 1, 10, 1, 13, 1, 1)
COVAL(SMOKE ,SMOS, FIXFLU , GRND )
COVAL(SMOKE ,SMOM, FIXFLU , GRND )
XCYCLE = F
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
SARAH =1.
SELREF = T
RESFAC =1.0E-03
************************************************************
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(SMOS)=20
LITER(SMOM)=20 ;LITER(F1)=20
LITER(F2)=20 ;LITER(F3)=20
LITER(F4)=20 ;LITER(F5)=20
LITER(F6)=20 ;LITER(F7)=20
LITER(F8)=20 ;LITER(F9)=20
LITER(F10)=20 ;LITER(F11)=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(SMOS)=1.0E-03
ENDIT(SMOM)=1.0E-03 ;ENDIT(F1)=1.0E-03
ENDIT(F2)=1.0E-03 ;ENDIT(F3)=1.0E-03
ENDIT(F4)=1.0E-03 ;ENDIT(F5)=1.0E-03
ENDIT(F6)=1.0E-03 ;ENDIT(F7)=1.0E-03
ENDIT(F8)=1.0E-03 ;ENDIT(F9)=1.0E-03
ENDIT(F10)=1.0E-03 ;ENDIT(F11)=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(PRDM,LINRLX,1.)
RELAX(OXDM,LINRLX,1.)
RELAX(FRFM,LINRLX,1.)
RELAX(TMPM,LINRLX,1.)
RELAX(SMOS,FALSDT,1.0E+09)
RELAX(SMOM,FALSDT,1.0E+09)
RELAX(F1,LINRLX,1.)
RELAX(F2,LINRLX,1.)
RELAX(F3,LINRLX,1.)
RELAX(F4,LINRLX,1.)
RELAX(F5,LINRLX,1.)
RELAX(F6,LINRLX,1.)
RELAX(F7,LINRLX,1.)
RELAX(F8,LINRLX,1.)
RELAX(F9,LINRLX,1.)
RELAX(F10,LINRLX,1.)
RELAX(F11,LINRLX,1.)
RELAX(MNSQ,LINRLX,1.)
RELAX(AVEF,LINRLX,1.)
RELAX(MIXL,LINRLX,1.)
RELAX(RATE,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,LINRLX,1.)
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(PRDM)=1.0E+10 ;VARMIN(PRDM)=-1.0E+10
VARMAX(OXDM)=1.0E+10 ;VARMIN(OXDM)=-1.0E+10
VARMAX(FRFM)=1.0E+10 ;VARMIN(FRFM)=-1.0E+10
VARMAX(TMPM)=1.0E+10 ;VARMIN(TMPM)=-1.0E+10
VARMAX(SMOS)=1.0E+10 ;VARMIN(SMOS)=-1.0E+10
VARMAX(SMOM)=1.0E+10 ;VARMIN(SMOM)=-1.0E+10
VARMAX(F1)=1. ;VARMIN(F1)=0.
VARMAX(F2)=1. ;VARMIN(F2)=0.
VARMAX(F3)=1. ;VARMIN(F3)=0.
VARMAX(F4)=1. ;VARMIN(F4)=0.
VARMAX(F5)=1. ;VARMIN(F5)=0.
VARMAX(F6)=1. ;VARMIN(F6)=0.
VARMAX(F7)=1. ;VARMIN(F7)=0.
VARMAX(F8)=1. ;VARMIN(F8)=0.
VARMAX(F9)=1. ;VARMIN(F9)=0.
VARMAX(F10)=1. ;VARMIN(F10)=0.
VARMAX(F11)=1. ;VARMIN(F11)=0.
VARMAX(MNSQ)=1. ;VARMIN(MNSQ)=0.
VARMAX(AVEF)=1. ;VARMIN(AVEF)=0.
VARMAX(MIXL)=1.0E+10 ;VARMIN(MIXL)=-1.0E+10
VARMAX(RATE)=1.0E+10 ;VARMIN(RATE)=-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,STORED,MIXL,C,=LEN1!ZSLFIN)
SPEDAT(SET,STORED,RATE,C,=EPKE!ZSLFIN)
SPEDAT(SET,MFM,FFUIN,R,0.109649)
SPEDAT(SET,MFM,TAIRIN,R,773.)
SPEDAT(SET,MFM,TFUIN,R,773.)
SPEDAT(SET,MFM,TSTOIC,R,2106.90918)
SPEDAT(SET,MFM,FACST,R,0.054825)
SPEDAT(SET,MFM,GASCON,R,8314.)
SPEDAT(SET,MFM,SPVLMX,R,100.)
SPEDAT(SET,MFM,SPVLMN,R,1.0E-02)
SPEDAT(SET,MFM,SMOCON,R,1.)
SPEDAT(SET,MFM,SMOEXP,R,7.)
SPEDAT(SET,MFM,MFMMOD,C,SCRS)
SPEDAT(SET,MFM,NFLUIDS,I,11)
SPEDAT(SET,MFM,NFLR,I,1)
SPEDAT(SET,MFM,NFLF,I,11)
SPEDAT(SET,MFM,CONREA,R,0.)
SPEDAT(SET,MFM,CONMIX,R,5.)
SPEDAT(SET,MFM,VISCON,R,0.)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,MATERIAL,0,L,T)
SPEDAT(SET,MATERIAL,199,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=204.5 ;EX(U1)=1.01
EX(V1)=5.748 ;EX(W1)=7.893
EX(KE)=52.830002 ;EX(EP)=1.995E+04
EX(H1)=1.93E+06 ;EX(PRPS)=0.9385
EX(PRDM)=0.5259 ;EX(OXDM)=0.4056
EX(FRFM)=7.007E-03 ;EX(TMPM)=1266.
EX(SMOS)=3.38E-06 ;EX(SMOM)=7.983E-06
EX(F1)=0.07249 ;EX(F2)=0.1613
EX(F3)=0.2186 ;EX(F4)=0.1438
EX(F5)=0.07713 ;EX(F6)=0.04748
EX(F7)=0.03901 ;EX(F8)=0.05254
EX(F9)=0.05453 ;EX(F10)=0.02652
EX(F11)=0.04501 ;EX(MNSQ)=0.3369
EX(AVEF)=0.03621 ;EX(MIXL)=0.
EX(RATE)=290.5 ;EX(DEGF)=0.
EX(RRAT)=0. ;EX(EPKE)=290.5
EX(HPOR)=0. ;EX(NPOR)=0.
EX(EPOR)=0. ;EX(VPOR)=0.
EX(ENUT)=0.01362 ;EX(RHO1)=2.198
EX(TMP1)=1121. ;EX(PROD)=0.4996
EX(OXID)=0.4337 ;EX(FUEL)=5.078E-03
EX(MIXF)=0.03247
************************************************************
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(PRDM,Y,N,Y,N,N,N)
OUTPUT(OXDM,Y,N,Y,N,N,N)
OUTPUT(FRFM,Y,N,Y,N,N,N)
OUTPUT(TMPM,Y,N,Y,N,N,N)
OUTPUT(SMOS,Y,N,Y,Y,Y,Y)
OUTPUT(SMOM,Y,N,Y,Y,Y,Y)
OUTPUT(F1,Y,N,Y,Y,Y,Y)
OUTPUT(F2,Y,N,Y,Y,Y,Y)
OUTPUT(F3,Y,N,Y,Y,Y,Y)
OUTPUT(F4,Y,N,Y,Y,Y,Y)
OUTPUT(F5,Y,N,Y,Y,Y,Y)
OUTPUT(F6,Y,N,Y,Y,Y,Y)
OUTPUT(F7,Y,N,Y,Y,Y,Y)
OUTPUT(F8,Y,N,Y,Y,Y,Y)
OUTPUT(F9,Y,N,Y,Y,Y,Y)
OUTPUT(F10,Y,N,Y,Y,Y,Y)
OUTPUT(F11,Y,N,Y,Y,Y,Y)
OUTPUT(MNSQ,Y,N,Y,N,N,N)
OUTPUT(AVEF,Y,N,Y,N,N,N)
OUTPUT(MIXL,Y,N,Y,N,N,N)
OUTPUT(RATE,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,N,N,N,N,N,N)
OUTPUT(NPOR,N,N,N,N,N,N)
OUTPUT(EPOR,N,N,N,N,N,N)
OUTPUT(VPOR,N,N,N,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 = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 500 ;NUMCLS = 5
NXPRIN = -1 ;IXPRF = 3 ;IXPRL = 3
NYPRIN = 1 ;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 = 1 ;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. )
PATCH(PROF1 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
PLOT(PROF1 ,F1 ,0. ,0. )
PATCH(PROF2 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
PLOT(PROF2 ,F2 ,0. ,0. )
PATCH(PROF3 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
PLOT(PROF3 ,F3 ,0. ,0. )
PATCH(PROF4 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
PLOT(PROF4 ,F4 ,0. ,0. )
PATCH(PROF5 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
PLOT(PROF5 ,F5 ,0. ,0. )
PATCH(PROF6 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
PLOT(PROF6 ,F6 ,0. ,0. )
PATCH(PROF7 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
PLOT(PROF7 ,F7 ,0. ,0. )
PATCH(PROF8 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
PLOT(PROF8 ,F8 ,0. ,0. )
PATCH(PROF9 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
PLOT(PROF9 ,F9 ,0. ,0. )
PATCH(PROF10 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
PLOT(PROF10 ,F10 ,0. ,0. )
PATCH(PROF11 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
PLOT(PROF11 ,F11 ,0. ,0. )
PATCH(MIDDLE ,PROFIL, 1, 1, 5, 5, 1, 13, 1, 1)
PLOT(MIDDLE ,MNSQ,0. ,0. )
PLOT(MIDDLE ,AVEF,0. ,0. )
PLOT(MIDDLE ,ENUT,0. ,0. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP