TALK=T;RUN(1,1) PHOTON USE p up z con temp x 1 z 1 16 fil;.001 gr ou x 1 z 1 16 msg SAFIR sub-model: Combustion of descending coke msg --------------------------------------------- msg msg Gas temperature pause con cl;red msg SAFIR sub-model: Combustion of descending coke msg --------------------------------------------- msg msg Fines volume fractions con rf x 1 z 1 16 fil;.001 *enduse DISPLAY Fines-and-coke-combustion-driven solid flow: One-phase, 2D, two spaces, scalar treatment of coal fines combustion, lump coke combustion, provision for raceway. ENDDIS ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(Fines-and-coke-combustion-driven solid f) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 14 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 1 XULAST =0.1 XFRAC(1)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 15 YVLAST =8. YFRAC(1)=0.066667 ;YFRAC(2)=0.133333 YFRAC(3)=0.2 ;YFRAC(4)=0.266667 YFRAC(5)=0.333333 ;YFRAC(6)=0.4 YFRAC(7)=0.466667 ;YFRAC(8)=0.533333 YFRAC(9)=0.6 ;YFRAC(10)=0.666667 YFRAC(11)=0.733333 ;YFRAC(12)=0.8 YFRAC(13)=0.866667 ;YFRAC(14)=0.933333 YFRAC(15)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 33 ZWLAST =57. ZFRAC(1)=0.030303 ;ZFRAC(2)=0.060606 ZFRAC(3)=0.090909 ;ZFRAC(4)=0.121212 ZFRAC(5)=0.151515 ;ZFRAC(6)=0.181818 ZFRAC(7)=0.212121 ;ZFRAC(8)=0.242424 ZFRAC(9)=0.272727 ;ZFRAC(10)=0.30303 ZFRAC(11)=0.333333 ;ZFRAC(12)=0.363636 ZFRAC(13)=0.393939 ;ZFRAC(14)=0.424242 ZFRAC(15)=0.454545 ;ZFRAC(16)=0.484848 ZFRAC(17)=0.515152 ;ZFRAC(18)=0.545455 ZFRAC(19)=0.575758 ;ZFRAC(20)=0.606061 ZFRAC(21)=0.636364 ;ZFRAC(22)=0.666667 ZFRAC(23)=0.69697 ;ZFRAC(24)=0.727273 ZFRAC(25)=0.757576 ;ZFRAC(26)=0.787879 ZFRAC(27)=0.818182 ;ZFRAC(28)=0.848485 ZFRAC(29)=0.878788 ;ZFRAC(30)=0.909091 ZFRAC(31)=0.939394 ;ZFRAC(32)=0.969697 ZFRAC(33)=1. ************************************************************ Group 6. Body-Fitted Coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(1)=P1 ;NAME(5)=V1 NAME(7)=W1 ;NAME(14)=H1 NAME(120)=PRPS ;NAME(121)=FURN NAME(122)=VVPO ;NAME(123)=BURN NAME(124)=YSUM ;NAME(125)=GOPA NAME(126)=GOFU ;NAME(127)=GH NAME(128)=GC ;NAME(129)=GO NAME(130)=FRAC ;NAME(131)=FLIM NAME(132)=YH2O ;NAME(133)=YCO2 NAME(134)=YCO ;NAME(135)=YO2 NAME(136)=YH2 ;NAME(137)=YN2 NAME(138)=TEMP ;NAME(139)=HSUB NAME(140)=RMIX ;NAME(141)=RF NAME(142)=POR ;NAME(143)=FCL NAME(144)=MARK ;NAME(145)=SPOR NAME(146)=VPOR ;NAME(147)=NPOR NAME(148)=HPOR ;NAME(149)=VISL NAME(150)=RHO1 * 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(H1,Y,Y,N,N,N,Y) SOLUTN(PRPS,Y,N,N,N,N,N) SOLUTN(FURN,Y,N,N,N,N,Y) SOLUTN(VVPO,Y,N,N,N,N,Y) SOLUTN(BURN,Y,N,N,N,N,Y) SOLUTN(YSUM,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(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(RF,Y,Y,N,N,N,Y) SOLUTN(POR,Y,Y,N,N,N,Y) SOLUTN(FCL,Y,Y,N,N,N,Y) SOLUTN(MARK,Y,N,N,N,N,Y) SOLUTN(SPOR,Y,N,N,N,N,Y) SOLUTN(VPOR,Y,N,N,N,N,Y) SOLUTN(NPOR,Y,N,N,N,N,Y) SOLUTN(HPOR,Y,N,N,N,N,Y) SOLUTN(VISL,Y,N,N,N,N,Y) SOLUTN(RHO1,Y,N,N,N,N,Y) DEN1 = 150 VISL = 149 EPOR = 0 ;HPOR = 148 ;NPOR = 147 ;VPOR = 146 PRPS = 120 ************************************************************ 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(H1,N,Y,N,Y,Y,Y) TERMS(RF,N,Y,N,Y,N,Y) TERMS(POR,N,Y,N,Y,Y,Y) TERMS(FCL,N,Y,N,Y,N,Y) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = F NEWRH1 = T NEWENL = 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 = GRND ;TMP1 =0. EL1 =0. 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. ENUL = GRND ;ENUT =0. ENULA =0. ;ENULB =0. ENULC =0. PRNDTL(V1)=1. ;PRNDTL(W1)=1. PRT(V1)=1. ;PRT(W1)=1. CP1 =1. ;CP2 =1. ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(P1)=1.0E-10 ;FIINIT(V1)=1.0E-10 FIINIT(W1)=1.0E-10 ;FIINIT(H1)=1.0E-10 FIINIT(PRPS)=-1. ;FIINIT(FURN)=1.0E-10 FIINIT(VVPO)=1.0E-10 ;FIINIT(BURN)=1.0E-10 FIINIT(YSUM)=1.0E-10 ;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(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)=350. ;FIINIT(HSUB)=1.0E-10 FIINIT(RMIX)=1.0E-10 ;FIINIT(RF)=1.0E-10 FIINIT(POR)=1.0E-04 ;FIINIT(FCL)=0.154804 FIINIT(MARK)=1.0E-10 ;FIINIT(SPOR)=1.0E-10 FIINIT(VPOR)=1. ;FIINIT(NPOR)=1. FIINIT(HPOR)=1. ;FIINIT(VISL)=1.0E-10 FIINIT(RHO1)=1.0E-10 PATCH(ISOLAT ,INIVAL, 1, 1, 1, 15, 17, 17, 1, 1) INIT(ISOLAT ,PRPS,0. ,199. ) INIT(ISOLAT ,NPOR,0. ,0. ) INIT(ISOLAT ,HPOR,0. ,0. ) PATCH(ISOL ,INIVAL, 1, 1, 1, 15, 16, 16, 1, 1) INIT(ISOL ,HPOR,0. ,0. ) PATCH(LOWMARK ,INIVAL, 1, 1, 1, 15, 1, 16, 1, 1) INIT(LOWMARK ,MARK,0. ,1. ) PATCH(UPPMARK ,INIVAL, 1, 1, 1, 15, 18, 33, 1, 1) INIT(UPPMARK ,MARK,0. ,2. ) PATCH(INILOW ,INIVAL, 1, 1, 1, 15, 1, 16, 1, 1) INIT(INILOW ,SPOR,0. ,0.25 ) INIT(INILOW ,VPOR,0. ,0.25 ) INIT(INILOW ,NPOR,0. ,0.25 ) INIT(INILOW ,HPOR,0. ,0.25 ) INIT(INILOW ,RHO1,0. ,1. ) PATCH(INIUP ,INIVAL, 1, 1, 1, 15, 18, 33, 1, 1) INIT(INIUP ,SPOR,0. ,0.75 ) INIT(INIUP ,VPOR,0. ,0.75 ) INIT(INIUP ,NPOR,0. ,0.75 ) INIT(INIUP ,HPOR,0. ,0.75 ) INIT(INIUP ,RHO1,0. ,1000. ) 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(GASINL ,NORTH , 1, 1, 15, 15, 2, 2, 1, 1) COVAL(GASINL ,P1 , FIXFLU ,20. ) COVAL(GASINL ,V1 ,0. ,-20. ) COVAL(GASINL ,H1 ,0. ,3.85E+05 ) COVAL(GASINL ,RF ,0. ,1.0E-04 ) COVAL(GASINL ,POR ,0. ,0. ) COVAL(GASINL ,FCL ,0. ,0. ) PATCH(TOPLOW ,HIGH , 1, 1, 1, 15, 16, 16, 1, 1) COVAL(TOPLOW ,P1 ,1. ,0. ) COVAL(TOPLOW ,RF ,0. , SAME ) COVAL(TOPLOW ,POR ,0. , SAME ) PATCH(SOLINL ,HIGH , 1, 1, 1, 15, 33, 33, 1, 1) COVAL(SOLINL ,P1 ,5000. ,0. ) COVAL(SOLINL ,RF ,0. , SAME ) COVAL(SOLINL ,POR ,0. , SAME ) COVAL(SOLINL ,FCL ,0. , SAME ) PATCH(FRIC ,VOLUME, 1, 1, 1, 15, 1, 16, 1, 1) COVAL(FRIC ,V1 , GRND ,0. ) COVAL(FRIC ,W1 , GRND ,0. ) PATCH(GASFCOKE,VOLUME, 1, 1, 1, 15, 1, 16, 1, 1) COVAL(GASFCOKE,P1 , FIXFLU , GRND ) COVAL(GASFCOKE,H1 ,0. , GRND ) COVAL(GASFCOKE,POR ,0. , SAME ) COVAL(GASFCOKE,FCL ,0. ,1. ) PATCH(COKE2GAS,VOLUME, 1, 1, 1, 15, 18, 33, 1, 1) COVAL(COKE2GAS,P1 ,1.0E-20 , GRND ) COVAL(COKE2GAS,RF ,0. , SAME ) COVAL(COKE2GAS,POR ,0. , SAME ) PATCH(GASFFINE,VOLUME, 1, 1, 1, 15, 1, 16, 1, 1) COVAL(GASFFINE,P1 , FIXFLU , GRND ) COVAL(GASFFINE,RF ,0. ,0. ) COVAL(GASFFINE,POR ,0. , SAME ) COVAL(GASFFINE,FCL ,0. ,1. ) PATCH(RACEWAY ,VOLUME, 1, 1, 1, 15, 1, 16, 1, 1) COVAL(RACEWAY ,POR , GRND , GRND ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 700 ;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(H1)=20 LITER(RF)=20 ;LITER(POR)=20 LITER(FCL)=20 ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03 ENDIT(W1)=1.0E-03 ;ENDIT(H1)=1.0E-03 ENDIT(RF)=1.0E-03 ;ENDIT(POR)=1.0E-03 ENDIT(FCL)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,0.15) RELAX(V1,FALSDT,1.0E-02) RELAX(W1,FALSDT,1.0E-02) RELAX(H1,FALSDT,1.0E-02) RELAX(PRPS,LINRLX,1.) RELAX(FURN,LINRLX,1.) RELAX(VVPO,LINRLX,1.) RELAX(BURN,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(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(RF,FALSDT,1.0E-02) RELAX(POR,FALSDT,1.0E-02) RELAX(FCL,FALSDT,1.0E-02) RELAX(MARK,LINRLX,1.) RELAX(SPOR,LINRLX,1.) RELAX(VISL,LINRLX,1.) RELAX(RHO1,LINRLX,0.15) 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(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10 VARMAX(FURN)=1.0E+10 ;VARMIN(FURN)=-1.0E+10 VARMAX(VVPO)=1.0E+10 ;VARMIN(VVPO)=-1.0E+10 VARMAX(BURN)=1.0E+10 ;VARMIN(BURN)=-1.0E+10 VARMAX(YSUM)=1.0E+10 ;VARMIN(YSUM)=-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(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)=3000. ;VARMIN(TEMP)=350. VARMAX(HSUB)=1.0E+10 ;VARMIN(HSUB)=-1.0E+10 VARMAX(RMIX)=1.0E+10 ;VARMIN(RMIX)=-1.0E+10 VARMAX(RF)=1.0E-04 ;VARMIN(RF)=0. VARMAX(POR)=1. ;VARMIN(POR)=0. VARMAX(FCL)=0.154804 ;VARMIN(FCL)=0. VARMAX(MARK)=1.0E+10 ;VARMIN(MARK)=-1.0E+10 VARMAX(SPOR)=1.0E+10 ;VARMIN(SPOR)=-1.0E+10 VARMAX(VPOR)=1.0E+10 ;VARMIN(VPOR)=-1.0E+10 VARMAX(NPOR)=1.0E+10 ;VARMIN(NPOR)=-1.0E+10 VARMAX(HPOR)=1.0E+10 ;VARMIN(HPOR)=-1.0E+10 VARMAX(VISL)=1.0E+10 ;VARMIN(VISL)=-1.0E+10 VARMAX(RHO1)=1000. ;VARMIN(RHO1)=1.0E-03 ************************************************************ Group 19. Data transmitted to GROUND NAMSAT =MOSG PARSOL = F ISG62 = 1 SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,MATERIAL,199,L,T) IG( 1) = 17 RG( 1) =1. RG( 2) =1000. ************************************************************ 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(H1,Y,N,Y,Y,Y,Y) OUTPUT(PRPS,Y,N,Y,N,N,N) OUTPUT(FURN,Y,N,Y,N,N,N) OUTPUT(VVPO,Y,N,Y,N,N,N) OUTPUT(BURN,Y,N,Y,N,N,N) OUTPUT(YSUM,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(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(RF,Y,N,Y,Y,Y,Y) OUTPUT(POR,Y,N,Y,Y,Y,Y) OUTPUT(FCL,Y,N,Y,Y,Y,Y) OUTPUT(MARK,Y,N,Y,N,N,N) OUTPUT(SPOR,Y,N,Y,N,N,N) OUTPUT(VPOR,Y,N,Y,N,N,N) OUTPUT(NPOR,Y,N,Y,N,N,N) OUTPUT(HPOR,Y,N,Y,N,N,N) OUTPUT(VISL,Y,N,Y,N,N,N) OUTPUT(RHO1,Y,N,Y,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 7 ;IZMON = 16 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 = 10000 NZPRIN = 1 ;IZPRF = 1 ;IZPRL = 10000 XZPR = F ;YZPR = T 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