TALK=T;RUN(1,1)
PHOTON USE
p;;;;
view x
msg compare
msg fluid-average values: avef, tmpm, frfm, prdm, oxdm, smom
msg with corresponding single-fluid EBU values
msg mixf, tmp1, fuel, prod, oxid, smos
pause
msg x-plane mixture fraction avef according to MFM
con avef x 1 fi;0.001;pause;con off;red
msg x-plane mixture fraction mixf according to EBU
msg This is, as it should be, similar to avef
con mixf x 1 fi;0.001;pause;con off;red
msg x-plane tmpm according to MFM
con tmpm x 1 fi;0.001;pause;con off;red
msg x-plane tmp1 according to EBU
con tmp1 x 1 fi;0.001;pause;con off;red
msg x-plane frfm according to MFM
con frfm x 1 fi;0.001;pause;con off;red
msg x-plane fuel according to EBU
con fuel x 1 fi;0.001;pause;con off;red
msg x-plane prdm according to MFM
con prdm x 1 fi;0.001;pause;con off;red
msg x-plane prod according to EBU
con prod x 1 fi;0.001;pause;con off;red
msg x-plane oxdm according to MFM
con oxdm x 1 fi;0.001;pause;con off;red
msg x-plane oxid according to EBU
con oxidd x 1 fi;0.001;pause;con off;red
msg x-plane smoke concentration according to eddy-break-up 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 MFM fluid-1 concentration
con f1 x 1 fi;0.001;pause;con off;red
msg x-plane MFM fluid-11 concentration
con f11 x 1 fi;0.001;pause;con off;red
msg x-plane MFM fluid-6 concentration
con f6 x 1 fi;0.001
ENDUSE
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
display
This 2D version of the 3D combustor facilitates exploration of the
influences of the model variables, and of the boundary conditions.
Some of these are explored in the following cases, 401, etc.
enddis
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Case 400 With Enlarged EBU Constant, Viz)
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 400
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
NX = 1
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(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(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(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(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(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(V1)=1.0E-10
FIINIT(W1)=10. ;FIINIT(KE)=11.25
FIINIT(EP)=353.255707 ;FIINIT(H1)=1.1595E+06
FIINIT(PRPS)=-1. ;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, 1, 7, 7, 1, 1, 1, 1)
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, 1, 8, 8, 1, 1, 1, 1)
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, 1, 9, 10, 1, 1, 1, 1)
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, 1, 8, 8, 2, 2, 1, 1)
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, 1, 9, 9, 2, 2, 1, 1)
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, 1, 10, 10, 2, 2, 1, 1)
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 ,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 ,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 ,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 ,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 ,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 ,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 ,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 ,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,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,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 ,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,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,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(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(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(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,0.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(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 =100.
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,199,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=5893. ;EX(V1)=6.939
EX(W1)=24.77 ;EX(KE)=23.219999
EX(EP)=1.066E+04 ;EX(H1)=1.765E+06
EX(PRPS)=0.9385 ;EX(PRDM)=0.3405
EX(OXDM)=0.5908 ;EX(FRFM)=7.153E-03
EX(TMPM)=1015. ;EX(SMOS)=3.601E-06
EX(SMOM)=7.359E-06 ;EX(F1)=0.4534
EX(F2)=0.1159 ;EX(F3)=0.04251
EX(F4)=0.03242 ;EX(F5)=0.03352
EX(F6)=0.03745 ;EX(F7)=0.03901
EX(F8)=0.05495 ;EX(F9)=0.05943
EX(F10)=0.02534 ;EX(F11)=0.04501
EX(MNSQ)=0.2738 ;EX(AVEF)=0.0262
EX(MIXL)=0. ;EX(RATE)=1176.
EX(DEGF)=0. ;EX(RRAT)=0.
EX(EPKE)=1176. ;EX(HPOR)=0.
EX(NPOR)=0. ;EX(EPOR)=0.
EX(VPOR)=0. ;EX(ENUT)=5.209E-03
EX(RHO1)=2.618 ;EX(TMP1)=983.799988
EX(PROD)=0.3687 ;EX(OXID)=0.5638
EX(FUEL)=5.911E-03 ;EX(MIXF)=0.02613
************************************************************
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(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
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