TALK=T;RUN(1,1)
PHOTON USE
p
up z
con temp x 1 y 1 m z 7 14 fil;.001
con temp x 1 y 1 m z 15 m sh;int 10
vec x 1 y 1 m z 15 m sh
vec x 1 y 1 m z 7 14
set vec ref 4.
vec x 1 y 1 m z 1 6 sh
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Velocity vectors and
msg temperature contours
msg
msg
msg Press
msg to continue
pause;vec off;con off;red
con temp x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Temperature
msg
msg
msg Press
msg to continue
pause;con off;red
con smok x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Smoke
msg
msg
msg Press
msg to continue
pause;con off;red
con YO2 x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Oxygen
msg
msg
msg Press
msg to continue
pause;con off;red
con YN2 x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Nitrogen
msg
msg
msg Press
msg to continue
pause;con off;red
con YH2 x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Hydrogen
msg
msg
msg Press
msg to continue
pause;con off;red
con YH2O x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Water vapour
msg
msg
msg Press
msg to continue
pause;con off;red
con YCO2 x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Carbon dioxide
msg
msg
msg Press
msg to continue
pause;con off;red
con YCO x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Carbon monoxide
msg
msg
msg Press
msg to continue
enduse
PLANTBEGIN
NAMSAT=MOSG
RHO1=GRND
DEN1=PRESS0/(RMIX*TEMP+tiny)
DEN1=AMIN1(VARMAX(140),AMAX1(0.0,DEN1,VARMIN(140)))
**
Carbon mass transfer related sources:
------------------------------------
PATCH(CARGAS,VOLUME,1,NX,1,NY,1,NZ,1,1)
(1) Transfer of mass leading to increase of gas flow rate:
- VPOR is volume fraction of lump coal
VAL=RG(14)*(1.-VPOR)*(RG(15)-FCL)
COVAL(CARGAS,P1,FIXFLU,GRND)
(2) Transfer of carbon leading to increase of mixture
fraction at the same rate:
- CO=1. signifies that mass tarnsfer brings in
material which is 100% carbon
COVAL(CARGAS,FCL,ONLYMS,1.)
COVAL(CARGAS,SMOK,ONLYMS,0.0)
(3) Transfer of enthalpy and heat leading to increase of
gas enthalpy at the same rate:
- Interphase gas temperature is assumed as TEMP.
- HSUB = HCOCO2*YCO * HH2*YH2
VAL=RG(10)*TEMP+RG(18)+HSUB
COVAL(CARGAS,H1,ONLYMS,GRND)
PATCH(NOX,PHASEM,1,NX,1,NY,1,NZ,1,1)
VAL=10.*EXP(-3000./TEMP)*YN2*(YO2/TEMP)**0.5
COVAL(NOX,SMOK,FIXFLU,GRND)
** Cell-wise composition parameters
--------------------------------
FLIM=RG(8)/(RG(8)+RG(12)*RG(1)/RG(7)+ RG(11)*RG(1$
)/(2*RG(3)))
GO=RG(8)*(1-FCL)
GC=RG(12)*FCL
GH=RG(11)*FCL
GOPA=GC*RG(1)/(2*RG(7))/(1-GO+GC*RG(1)/(2*RG(7))+TINY)
GOFU=(GH*RG(1)/(2*RG(3))+GC*RG(1)/RG(7))/ (1$
.-GO+GH*RG(1)/(2*RG(3))+GC*RG(1)/RG(7)+TINY)
FRAC=(GO-GOPA)/(GOFU-GOPA+TINY)
** For all regions
---------------
YN2=RG(13)*FCL+RG(19)*(1.-FCL)
** Region 1
--------
YH2O=RG(11)*FCL*RG(5)/RG(3)
IF(FCL.LE.FLIM)
YCO2=RG(12)*FCL*RG(4)/RG(7)
IF(FCL.LE.FLIM)
YO2 =RG(8)*(1-FCL)-RG(12)*FCL*RG(1)/RG(7)- RG(11)*FC$
L*RG(1)/(2.*RG(3))
IF(FCL.LE.FLIM)
YCO=0.0
IF(FCL.LE.FLIM)
YH2=0.0
IF(FCL.LE.FLIM)
HSUB=0.0
IF(FCL.LE.FLIM)
RMIX=RG(9)*(YO2/RG(1)+YH2O/RG(5)+YCO2/RG(4)+ YN2/R$
G(2))
IF(FCL.LE.FLIM)
** Region 2
--------
YH2O=RG(11)*FCL*RG(5)/RG(3)*FRAC*(1-GOFU)/(1-GO+TINY)
IF(FCL.GT.FLIM.AND.FRAC.GE.0.)
YCO2=RG(12)*FCL*RG(4)/RG(7)*FRAC*(1-GOFU)/(1-GO+TINY)
IF(FCL.GT.FLIM.AND.FRAC.GE.0.)
YO2=0.0
IF(FCL.GT.FLIM.AND.FRAC.GE.0.)
YCO=RG(12)*FCL*RG(6)/RG(7)*(1-FRAC)* (1$
-GOPA)/(1-GO+TINY)
IF(FCL.GT.FLIM.AND.FRAC.GE.0.)
YH2=RG(11)*FCL*(1-FRAC)*(1-GOPA)/(1-GO+TINY)
IF(FCL.GT.FLIM.AND.FRAC.GE.0.)
HSUB=YCO*RG(16)+YH2*RG(17)
IF(FCL.GT.FLIM.AND.FRAC.GE.0.)
RMIX=RG(9)*(YH2O/RG(5)+YCO/RG(6)+YCO2/RG(4)+ YH2/R$
G(3)+YN2/RG(2))
IF(FCL.GT.FLIM.AND.FRAC.GE.0.)
** Region 3
--------
YH2O=0.0
IF(FCL.GT.FLIM.AND.FRAC.LT.0.)
YCO2=0.0
IF(FCL.GT.FLIM.AND.FRAC.LT.0.)
YO2=0.0
IF(FCL.GT.FLIM.AND.FRAC.LT.0.)
YCO=RG(8)*(1-FCL)*2*RG(6)/RG(1)
IF(FCL.GT.FLIM.AND.FRAC.LT.0.)
YH2=RG(11)*FCL
IF(FCL.GT.FLIM.AND.FRAC.LT.0.)
HSUB=YCO*RG(16)+YH2*RG(17)
IF(FCL.GT.FLIM.AND.FRAC.LT.0.)
RMIX=RG(9)*(YCO/RG(6)+YH2/RG(3)+YN2/RG(2))
IF(FCL.GT.FLIM.AND.FRAC.LT.0.)
** Calculation of absolute gas temperature
--------------------------------------
TEMP=(H1-HSUB)/RG(10)
TEMP=AMIN1(VARMAX(147),AMAX1(100.,TEMP,VARMIN(147)))
store(ysum)
YSUM=YN2+YO2+YCO+YCO2+YH2O+YH2
PLANTEND
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Coke Burning In A Grate Stoker )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 1
XULAST =1.
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 20
YVLAST =1.
YFRAC(1)=0.05 ;YFRAC(2)=0.1
YFRAC(3)=0.15 ;YFRAC(4)=0.2
YFRAC(5)=0.25 ;YFRAC(6)=0.3
YFRAC(7)=0.35 ;YFRAC(8)=0.4
YFRAC(9)=0.45 ;YFRAC(10)=0.5
YFRAC(11)=0.55 ;YFRAC(12)=0.6
YFRAC(13)=0.65 ;YFRAC(14)=0.7
YFRAC(15)=0.75 ;YFRAC(16)=0.8
YFRAC(17)=0.85 ;YFRAC(18)=0.9
YFRAC(19)=0.95 ;YFRAC(20)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 30
ZWLAST =1.5
ZFRAC(1)=0.033333 ;ZFRAC(2)=0.066667
ZFRAC(3)=0.1 ;ZFRAC(4)=0.133333
ZFRAC(5)=0.166667 ;ZFRAC(6)=0.2
ZFRAC(7)=0.233333 ;ZFRAC(8)=0.266667
ZFRAC(9)=0.3 ;ZFRAC(10)=0.333333
ZFRAC(11)=0.366667 ;ZFRAC(12)=0.4
ZFRAC(13)=0.433333 ;ZFRAC(14)=0.466667
ZFRAC(15)=0.5 ;ZFRAC(16)=0.533333
ZFRAC(17)=0.566667 ;ZFRAC(18)=0.6
ZFRAC(19)=0.633333 ;ZFRAC(20)=0.666667
ZFRAC(21)=0.7 ;ZFRAC(22)=0.733333
ZFRAC(23)=0.766667 ;ZFRAC(24)=0.8
ZFRAC(25)=0.833333 ;ZFRAC(26)=0.866667
ZFRAC(27)=0.9 ;ZFRAC(28)=0.933333
ZFRAC(29)=0.966667 ;ZFRAC(30)=1.
************************************************************
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(126)=PRPS ;NAME(127)=EPKE
NAME(128)=YSUM ;NAME(129)=HPOR
NAME(130)=NPOR ;NAME(131)=VPOR
NAME(132)=GOPA ;NAME(133)=GOFU
NAME(134)=GH ;NAME(135)=GC
NAME(136)=GO ;NAME(137)=FRAC
NAME(138)=FLIM ;NAME(139)=DEN1
NAME(140)=YH2O ;NAME(141)=YCO2
NAME(142)=YCO ;NAME(143)=YO2
NAME(144)=YH2 ;NAME(145)=YN2
NAME(146)=TEMP ;NAME(147)=HSUB
NAME(148)=RMIX ;NAME(149)=SMOK
NAME(150)=FCL
* 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,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,Y,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,N)
SOLUTN(EPKE,Y,N,N,N,N,Y)
SOLUTN(YSUM,Y,N,N,N,N,Y)
SOLUTN(HPOR,Y,N,N,N,N,Y)
SOLUTN(NPOR,Y,N,N,N,N,Y)
SOLUTN(VPOR,Y,N,N,N,N,Y)
SOLUTN(GOPA,Y,N,N,N,N,Y)
SOLUTN(GOFU,Y,N,N,N,N,Y)
SOLUTN(GH,Y,N,N,N,N,Y)
SOLUTN(GC,Y,N,N,N,N,Y)
SOLUTN(GO,Y,N,N,N,N,Y)
SOLUTN(FRAC,Y,N,N,N,N,Y)
SOLUTN(FLIM,Y,N,N,N,N,Y)
SOLUTN(DEN1,Y,N,N,N,N,Y)
SOLUTN(YH2O,Y,N,N,N,N,Y)
SOLUTN(YCO2,Y,N,N,N,N,Y)
SOLUTN(YCO,Y,N,N,N,N,Y)
SOLUTN(YO2,Y,N,N,N,N,Y)
SOLUTN(YH2,Y,N,N,N,N,Y)
SOLUTN(YN2,Y,N,N,N,N,Y)
SOLUTN(TEMP,Y,N,N,N,N,Y)
SOLUTN(HSUB,Y,N,N,N,N,Y)
SOLUTN(RMIX,Y,N,N,N,N,Y)
SOLUTN(SMOK,Y,Y,Y,N,N,Y)
SOLUTN(FCL,Y,Y,Y,N,N,Y)
DEN1 = 139
EPOR = 0 ;HPOR = 129 ;NPOR = 130 ;VPOR = 131
PRPS = 126
************************************************************
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,N,Y,Y,Y)
TERMS(SMOK,N,Y,N,Y,N,Y)
TERMS(FCL,N,Y,N,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWRH1 = T
NEWENT = 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 = GRND ;TMP1 =0. ;EL1 = GRND4
TSURR =0. ;TEMP0 =0. ;PRESS0 =1.0E+05
DVO1DT =0. ;DRH1DP =0.
RHO1A =0. ;RHO1B =0. ;RHO1C =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0. ;EL1B =0. ;EL1C =0.
ENUL =0. ;ENUT = GRND3
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PRNDTL(V1)=1. ;PRNDTL(W1)=1.
PRNDTL(KE)=1. ;PRNDTL(EP)=1.
PRT(V1)=1. ;PRT(W1)=1.
PRT(KE)=1. ;PRT(EP)=1.314
CP1 =1005. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.3E-04 ;FIINIT(V1)=1.0E-10
FIINIT(W1)=1. ;FIINIT(KE)=0.012093
FIINIT(EP)=9.469614E-03 ;FIINIT(H1)=3.85E+05
FIINIT(PRPS)=-1. ;FIINIT(EPKE)=1.
FIINIT(YSUM)=1.0E-10 ;FIINIT(HPOR)=1.
FIINIT(NPOR)=1. ;FIINIT(VPOR)=1.
FIINIT(GOPA)=1.0E-10 ;FIINIT(GOFU)=1.0E-10
FIINIT(GH)=1.0E-10 ;FIINIT(GC)=1.0E-10
FIINIT(GO)=1.0E-10 ;FIINIT(FRAC)=1.0E-10
FIINIT(FLIM)=1.0E-10 ;FIINIT(DEN1)=1.0E-10
FIINIT(YH2O)=1.0E-10 ;FIINIT(YCO2)=1.0E-10
FIINIT(YCO)=1.0E-10 ;FIINIT(YO2)=1.0E-10
FIINIT(YH2)=1.0E-10 ;FIINIT(YN2)=1.0E-10
FIINIT(TEMP)=1.0E-10 ;FIINIT(HSUB)=1.0E-10
FIINIT(RMIX)=1.0E-10 ;FIINIT(SMOK)=0.
FIINIT(FCL)=0.154804
PATCH(OB2 ,INIVAL, 1, 1, 1, 20, 7, 14, 1, 1)
INIT(OB2 ,PRPS,0. ,-1. )
INIT(OB2 ,HPOR,0. ,0.75 )
INIT(OB2 ,NPOR,0. ,0.75 )
INIT(OB2 ,VPOR,0. ,0.75 )
INIT(OB2 ,SMOK,0. ,0. )
INIT(OB2 ,FCL ,0. ,0. )
PATCH(^OB3 ,INIVAL, 3, 0, 0, 0, 0, 0, 1, 1)
INIT(^OB3 ,PRPS,0. ,198. )
PATCH(^OB4 ,INIVAL, 4, 0, 0, 0, 0, 0, 1, 1)
INIT(^OB4 ,PRPS,0. ,198. )
PATCH(OB5 ,INIVAL, 1, 1, 1, 1, 6, 6, 1, 1)
INIT(OB5 ,HPOR,0. ,0. )
PATCH(OB6 ,INIVAL, 1, 1, 4, 5, 6, 6, 1, 1)
INIT(OB6 ,HPOR,0. ,0. )
PATCH(OB7 ,INIVAL, 1, 1, 10, 10, 6, 6, 1, 1)
INIT(OB7 ,HPOR,0. ,0. )
PATCH(OB8 ,INIVAL, 1, 1, 11, 11, 6, 6, 1, 1)
INIT(OB8 ,HPOR,0. ,0. )
PATCH(OB9 ,INIVAL, 1, 1, 14, 15, 6, 6, 1, 1)
INIT(OB9 ,HPOR,0. ,0. )
PATCH(OBA ,INIVAL, 1, 1, 20, 20, 6, 6, 1, 1)
INIT(OBA ,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, 1, 1, 1, 20, 1, 30, 1, 1)
COVAL(KESOURCE,KE , GRND4 , GRND4 )
COVAL(KESOURCE,EP , GRND4 , GRND4 )
PATCH(CARGAS ,VOLUME, 1, 1, 1, 20, 1, 30, 1, 1)
COVAL(CARGAS ,P1 , FIXFLU , GRND )
COVAL(CARGAS ,H1 ,0. , GRND )
COVAL(CARGAS ,SMOK,0. ,0. )
COVAL(CARGAS ,FCL ,0. ,1. )
PATCH(NOX ,PHASEM, 1, 1, 1, 20, 1, 30, 1, 1)
COVAL(NOX ,SMOK, FIXFLU , GRND )
PATCH(OV2 ,PHASEM, 1, 1, 1, 20, 7, 14, 1, 1)
COVAL(OV2 ,V1 ,100. ,0. )
COVAL(OV2 ,W1 ,100. ,0. )
PATCH(OB5-L ,HWALL , 1, 1, 1, 1, 6, 6, 1, 1)
COVAL(OB5-L ,V1 , GRND2 ,0. )
COVAL(OB5-L ,W1 , FIXVAL ,0. )
COVAL(OB5-L ,KE , GRND2 , GRND2 )
COVAL(OB5-L ,EP , GRND2 , GRND2 )
PATCH(OB5-H ,LWALL , 1, 1, 1, 1, 7, 7, 1, 1)
COVAL(OB5-H ,V1 , GRND2 ,0. )
COVAL(OB5-H ,KE , GRND2 , GRND2 )
COVAL(OB5-H ,EP , GRND2 , GRND2 )
PATCH(OB6-L ,HWALL , 1, 1, 4, 5, 6, 6, 1, 1)
COVAL(OB6-L ,V1 , GRND2 ,0. )
COVAL(OB6-L ,W1 , FIXVAL ,0. )
COVAL(OB6-L ,KE , GRND2 , GRND2 )
COVAL(OB6-L ,EP , GRND2 , GRND2 )
PATCH(OB6-H ,LWALL , 1, 1, 4, 5, 7, 7, 1, 1)
COVAL(OB6-H ,V1 , GRND2 ,0. )
COVAL(OB6-H ,KE , GRND2 , GRND2 )
COVAL(OB6-H ,EP , GRND2 , GRND2 )
PATCH(OB7-L ,HWALL , 1, 1, 10, 10, 6, 6, 1, 1)
COVAL(OB7-L ,V1 , GRND2 ,0. )
COVAL(OB7-L ,W1 , FIXVAL ,0. )
COVAL(OB7-L ,KE , GRND2 , GRND2 )
COVAL(OB7-L ,EP , GRND2 , GRND2 )
PATCH(OB7-H ,LWALL , 1, 1, 10, 10, 7, 7, 1, 1)
COVAL(OB7-H ,V1 , GRND2 ,0. )
COVAL(OB7-H ,KE , GRND2 , GRND2 )
COVAL(OB7-H ,EP , GRND2 , GRND2 )
PATCH(OB8-L ,HWALL , 1, 1, 11, 11, 6, 6, 1, 1)
COVAL(OB8-L ,V1 , GRND2 ,0. )
COVAL(OB8-L ,W1 , FIXVAL ,0. )
COVAL(OB8-L ,KE , GRND2 , GRND2 )
COVAL(OB8-L ,EP , GRND2 , GRND2 )
PATCH(OB8-H ,LWALL , 1, 1, 11, 11, 7, 7, 1, 1)
COVAL(OB8-H ,V1 , GRND2 ,0. )
COVAL(OB8-H ,KE , GRND2 , GRND2 )
COVAL(OB8-H ,EP , GRND2 , GRND2 )
PATCH(OB9-L ,HWALL , 1, 1, 14, 15, 6, 6, 1, 1)
COVAL(OB9-L ,V1 , GRND2 ,0. )
COVAL(OB9-L ,W1 , FIXVAL ,0. )
COVAL(OB9-L ,KE , GRND2 , GRND2 )
COVAL(OB9-L ,EP , GRND2 , GRND2 )
PATCH(OB9-H ,LWALL , 1, 1, 14, 15, 7, 7, 1, 1)
COVAL(OB9-H ,V1 , GRND2 ,0. )
COVAL(OB9-H ,KE , GRND2 , GRND2 )
COVAL(OB9-H ,EP , GRND2 , GRND2 )
PATCH(OBA-L ,HWALL , 1, 1, 20, 20, 6, 6, 1, 1)
COVAL(OBA-L ,V1 , GRND2 ,0. )
COVAL(OBA-L ,W1 , FIXVAL ,0. )
COVAL(OBA-L ,KE , GRND2 , GRND2 )
COVAL(OBA-L ,EP , GRND2 , GRND2 )
PATCH(OBA-H ,LWALL , 1, 1, 20, 20, 7, 7, 1, 1)
COVAL(OBA-H ,V1 , GRND2 ,0. )
COVAL(OBA-H ,KE , GRND2 , GRND2 )
COVAL(OBA-H ,EP , GRND2 , GRND2 )
PATCH(OBB ,WEST , 1, 1, 1, 20, 1, 6, 1, 1)
COVAL(OBB ,P1 , FIXFLU ,2.199388 )
COVAL(OBB ,V1 ,0. ,0. )
COVAL(OBB ,W1 ,0. ,2. )
COVAL(OBB ,KE ,0. ,0.012093 )
COVAL(OBB ,EP ,0. ,9.469614E-03 )
COVAL(OBB ,H1 ,0. ,3.85E+05 )
PATCH(OBC ,HIGH , 1, 1, 1, 20, 30, 30, 1, 1)
COVAL(OBC ,P1 ,1. ,0. )
COVAL(OBC ,V1 ,0. ,0. )
COVAL(OBC ,W1 ,0. ,0. )
COVAL(OBC ,KE ,0. , SAME )
COVAL(OBC ,EP ,0. , SAME )
COVAL(OBC ,H1 ,0. ,0. )
PATCH(OBD ,WEST , 1, 1, 1, 20, 15, 25, 1, 1)
COVAL(OBD ,P1 , FIXFLU ,4.398777 )
COVAL(OBD ,V1 ,0. ,0. )
COVAL(OBD ,W1 ,0. ,0. )
COVAL(OBD ,KE ,0. ,0.048373 )
COVAL(OBD ,EP ,0. ,0.049269 )
COVAL(OBD ,H1 ,0. ,3.85E+05 )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 90 ;ISWC1 = 1
LITHYD = 10 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(V1)=10
LITER(W1)=10 ;LITER(KE)=20
LITER(EP)=20 ;LITER(H1)=20
LITER(SMOK)=20 ;LITER(FCL)=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(SMOK)=1.0E-03 ;ENDIT(FCL)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,0.175)
RELAX(V1,FALSDT,0.1)
RELAX(W1,FALSDT,0.1)
RELAX(KE,FALSDT,0.1)
RELAX(EP,FALSDT,0.1)
RELAX(H1,FALSDT,1.0E+09)
RELAX(PRPS,LINRLX,1.)
RELAX(EPKE,LINRLX,1.)
RELAX(YSUM,LINRLX,1.)
RELAX(GOPA,LINRLX,1.)
RELAX(GOFU,LINRLX,1.)
RELAX(GH,LINRLX,1.)
RELAX(GC,LINRLX,1.)
RELAX(GO,LINRLX,1.)
RELAX(FRAC,LINRLX,1.)
RELAX(FLIM,LINRLX,1.)
RELAX(DEN1,LINRLX,0.1)
RELAX(YH2O,LINRLX,1.)
RELAX(YCO2,LINRLX,1.)
RELAX(YCO,LINRLX,1.)
RELAX(YO2,LINRLX,1.)
RELAX(YH2,LINRLX,1.)
RELAX(YN2,LINRLX,1.)
RELAX(TEMP,LINRLX,1.)
RELAX(HSUB,LINRLX,1.)
RELAX(RMIX,LINRLX,1.)
RELAX(SMOK,FALSDT,1.0E-02)
RELAX(FCL,FALSDT,0.1)
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(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
VARMAX(YSUM)=1.0E+10 ;VARMIN(YSUM)=-1.0E+10
VARMAX(HPOR)=1.0E+10 ;VARMIN(HPOR)=-1.0E+10
VARMAX(NPOR)=1.0E+10 ;VARMIN(NPOR)=-1.0E+10
VARMAX(VPOR)=1.0E+10 ;VARMIN(VPOR)=-1.0E+10
VARMAX(GOPA)=1.0E+10 ;VARMIN(GOPA)=-1.0E+10
VARMAX(GOFU)=1.0E+10 ;VARMIN(GOFU)=-1.0E+10
VARMAX(GH)=1.0E+10 ;VARMIN(GH)=-1.0E+10
VARMAX(GC)=1.0E+10 ;VARMIN(GC)=-1.0E+10
VARMAX(GO)=1.0E+10 ;VARMIN(GO)=-1.0E+10
VARMAX(FRAC)=1.0E+10 ;VARMIN(FRAC)=-1.0E+10
VARMAX(FLIM)=1.0E+10 ;VARMIN(FLIM)=-1.0E+10
VARMAX(DEN1)=3. ;VARMIN(DEN1)=1.0E-03
VARMAX(YH2O)=1.0E+10 ;VARMIN(YH2O)=-1.0E+10
VARMAX(YCO2)=1.0E+10 ;VARMIN(YCO2)=-1.0E+10
VARMAX(YCO)=1.0E+10 ;VARMIN(YCO)=-1.0E+10
VARMAX(YO2)=1.0E+10 ;VARMIN(YO2)=-1.0E+10
VARMAX(YH2)=1.0E+10 ;VARMIN(YH2)=-1.0E+10
VARMAX(YN2)=1.0E+10 ;VARMIN(YN2)=-1.0E+10
VARMAX(TEMP)=5000. ;VARMIN(TEMP)=350.
VARMAX(HSUB)=1.0E+10 ;VARMIN(HSUB)=-1.0E+10
VARMAX(RMIX)=1.0E+10 ;VARMIN(RMIX)=-1.0E+10
VARMAX(SMOK)=1.0E+10 ;VARMIN(SMOK)=-1.0E+10
VARMAX(FCL)=1. ;VARMIN(FCL)=0.
************************************************************
Group 19. Data transmitted to GROUND
NAMSAT =MOSG
GENK = T
PARSOL = F
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,OBJNAM,!OB1,C,Z1)
SPEDAT(SET,OBJTYP,!OB1,C,NULL)
SPEDAT(SET,OBJNAM,!OB2,C,POROS)
SPEDAT(SET,OBJTYP,!OB2,C,BLOCKAGE)
SPEDAT(SET,POROS,MATERIAL,R,-1.)
SPEDAT(SET,OBJNAM,!OV2,C,POROS)
SPEDAT(SET,OBJTYP,!OV2,C,BLOCKAGE)
SPEDAT(SET,OBJNAM,^OB3,C,SCHIMNEY)
SPEDAT(SET,OBJTYP,^OB3,C,BLOCKAGE)
SPEDAT(SET,SCHIMNEY,MATERIAL,R,198.)
SPEDAT(SET,OBJNAM,^OB4,C,NCHIMNEY)
SPEDAT(SET,OBJTYP,^OB4,C,BLOCKAGE)
SPEDAT(SET,NCHIMNEY,MATERIAL,R,198.)
SPEDAT(SET,OBJNAM,!OB5,C,POR1)
SPEDAT(SET,OBJTYP,!OB5,C,PLATE)
SPEDAT(SET,OBJNAM,!OB5-L,C,POR1)
SPEDAT(SET,OBJTYP,!OB5-L,C,PLATE)
SPEDAT(SET,OBJNAM,!OB5-H,C,POR1)
SPEDAT(SET,OBJTYP,!OB5-H,C,PLATE)
SPEDAT(SET,OBJNAM,!OB6,C,POR2)
SPEDAT(SET,OBJTYP,!OB6,C,PLATE)
SPEDAT(SET,OBJNAM,!OB6-L,C,POR2)
SPEDAT(SET,OBJTYP,!OB6-L,C,PLATE)
SPEDAT(SET,OBJNAM,!OB6-H,C,POR2)
SPEDAT(SET,OBJTYP,!OB6-H,C,PLATE)
SPEDAT(SET,OBJNAM,!OB7,C,POR3)
SPEDAT(SET,OBJTYP,!OB7,C,PLATE)
SPEDAT(SET,OBJNAM,!OB7-L,C,POR3)
SPEDAT(SET,OBJTYP,!OB7-L,C,PLATE)
SPEDAT(SET,OBJNAM,!OB7-H,C,POR3)
SPEDAT(SET,OBJTYP,!OB7-H,C,PLATE)
SPEDAT(SET,OBJNAM,!OB8,C,POR4)
SPEDAT(SET,OBJTYP,!OB8,C,PLATE)
SPEDAT(SET,OBJNAM,!OB8-L,C,POR4)
SPEDAT(SET,OBJTYP,!OB8-L,C,PLATE)
SPEDAT(SET,OBJNAM,!OB8-H,C,POR4)
SPEDAT(SET,OBJTYP,!OB8-H,C,PLATE)
SPEDAT(SET,OBJNAM,!OB9,C,POR5)
SPEDAT(SET,OBJTYP,!OB9,C,PLATE)
SPEDAT(SET,OBJNAM,!OB9-L,C,POR5)
SPEDAT(SET,OBJTYP,!OB9-L,C,PLATE)
SPEDAT(SET,OBJNAM,!OB9-H,C,POR5)
SPEDAT(SET,OBJTYP,!OB9-H,C,PLATE)
SPEDAT(SET,OBJNAM,!OBA,C,POR6)
SPEDAT(SET,OBJTYP,!OBA,C,PLATE)
SPEDAT(SET,OBJNAM,!OBA-L,C,POR6)
SPEDAT(SET,OBJTYP,!OBA-L,C,PLATE)
SPEDAT(SET,OBJNAM,!OBA-H,C,POR6)
SPEDAT(SET,OBJTYP,!OBA-H,C,PLATE)
SPEDAT(SET,OBJNAM,!OBB,C,BOTAIR)
SPEDAT(SET,OBJTYP,!OBB,C,INLET)
SPEDAT(SET,OBJNAM,!OBC,C,OUTLET)
SPEDAT(SET,OBJTYP,!OBC,C,OUTLET)
SPEDAT(SET,ARATIO,!OBC,R,1.)
SPEDAT(SET,OBJNAM,!OBD,C,FLAMAIR)
SPEDAT(SET,OBJTYP,!OBD,C,INLET)
SPEDAT(SET,FACETDAT,NUMOBJ,I,83)
SPEDAT(SET,MATERIAL,198,L,T)
RG( 1) =32.
RG( 2) =28.
RG( 3) =2.
RG( 4) =44.
RG( 5) =18.
RG( 6) =28.
RG( 7) =12.
RG( 8) =0.232
RG( 9) =8314.299805
RG( 10) =1500.
RG( 11) =0.05
RG( 12) =0.95
RG( 14) =9.
RG( 15) =0.154804
RG( 16) =1.010657E+07
RG( 17) =1.209E+06
RG( 18) =2.991049E+07
RG( 19) =0.768
************************************************************
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(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(EPKE,Y,N,Y,N,N,N)
OUTPUT(YSUM,Y,N,Y,N,N,N)
OUTPUT(HPOR,Y,N,Y,N,N,N)
OUTPUT(NPOR,Y,N,Y,N,N,N)
OUTPUT(VPOR,Y,N,Y,N,N,N)
OUTPUT(GOPA,Y,N,Y,N,N,N)
OUTPUT(GOFU,Y,N,Y,N,N,N)
OUTPUT(GH,Y,N,Y,N,N,N)
OUTPUT(GC,Y,N,Y,N,N,N)
OUTPUT(GO,Y,N,Y,N,N,N)
OUTPUT(FRAC,Y,N,Y,N,N,N)
OUTPUT(FLIM,Y,N,Y,N,N,N)
OUTPUT(DEN1,Y,N,Y,N,N,N)
OUTPUT(YH2O,Y,N,Y,N,N,N)
OUTPUT(YCO2,Y,N,Y,N,N,N)
OUTPUT(YCO,Y,N,Y,N,N,N)
OUTPUT(YO2,Y,N,Y,N,N,N)
OUTPUT(YH2,Y,N,Y,N,N,N)
OUTPUT(YN2,Y,N,Y,N,N,N)
OUTPUT(TEMP,Y,N,Y,N,N,N)
OUTPUT(HSUB,Y,N,Y,N,N,N)
OUTPUT(RMIX,Y,N,Y,N,N,N)
OUTPUT(SMOK,Y,N,Y,Y,Y,Y)
OUTPUT(FCL,Y,N,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 4 ;IZMON = 30
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 30
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
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP