TALK=T;RUN(1,1) DISPLAY The case considered is the turbulent combustion of CO and H2 with O2 in a pipe of 32cm diameter. Concentric 'fuel' and 'oxidiser' streams enter the pipe axially, and then having mixed and burnt, exit axially 2m downstream of the inlet. The fuel stream, which comprises H2, CO, H2O, CO2 and N2, enters through an inner injector of 8cm diameter, and the oxidiser stream, which comprises O2, N2, H2O and CO2 enters through an outer annulus. The tempera- tures of the fuel and oxidiser streams are 1754 and 1518 Deg C, and the corresponding mass flow rates are 0.0726 and 0.144 kg/s. The fast-chemistry assumption is invoked for main combustion calculation in terms of two competing single-step infinite-rate reactions: 2H2 + O2 > 2H2O and CO + O2 > CO2. The NOX calculation is peformed as a post-processing run with a frozen hydrodynamic and temperature fields. The NO distributions are calculated by CHEMKIN via a reaction mechanism involving the Zeldovich reactions and fuel kinetics. ENDDIS ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(RADIAN: 2D Turbulent Diffusion Flame ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 14 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = F NX = 1 XULAST =0.1 XFRAC(1)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 8 YVLAST =0.16 YFRAC(1)=0.0625 ;YFRAC(2)=0.125 YFRAC(3)=0.1875 ;YFRAC(4)=0.25 YFRAC(5)=0.4375 ;YFRAC(6)=0.625 YFRAC(7)=0.8125 ;YFRAC(8)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 8 ZWLAST =2. ZFRAC(1)=0.044194 ;ZFRAC(2)=0.125 ZFRAC(3)=0.22964 ;ZFRAC(4)=0.353553 ZFRAC(5)=0.494106 ;ZFRAC(6)=0.649519 ZFRAC(7)=0.818488 ;ZFRAC(8)=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(16)=CH4 ;NAME(17)=O2 NAME(18)=H2 ;NAME(19)=CO NAME(20)=H2O ;NAME(21)=CO2 NAME(22)=N2 ;NAME(144)=F NAME(145)=SRAD ;NAME(146)=YSUM NAME(147)=SPH1 ;NAME(148)=TMP1 NAME(149)=DEN1 ;NAME(150)=VIST * 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,N) SOLUTN(W1,Y,Y,N,N,N,N) 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(CH4,Y,N,N,N,N,N) SOLUTN(O2,Y,N,N,N,N,N) SOLUTN(H2,Y,N,N,N,N,N) SOLUTN(CO,Y,N,N,N,N,N) SOLUTN(H2O,Y,N,N,N,N,N) SOLUTN(CO2,Y,N,N,N,N,N) SOLUTN(N2,Y,N,N,N,N,N) SOLUTN(F,Y,Y,Y,N,N,Y) SOLUTN(SRAD,Y,Y,Y,N,N,Y) SOLUTN(YSUM,Y,N,N,N,N,Y) SOLUTN(SPH1,Y,N,N,N,N,Y) SOLUTN(TMP1,Y,N,N,N,N,Y) SOLUTN(DEN1,Y,N,N,N,N,Y) SOLUTN(VIST,Y,N,N,N,N,Y) DEN1 = 149 VIST = 150 TEMP1 = 148 ************************************************************ 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,Y,Y,Y) TERMS(F,N,Y,Y,Y,Y,Y) TERMS(SRAD,N,N,Y,N,Y,N) 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 = GRND9 ;TMP1 = GRND9 EL1 = GRND4 TSURR =0. ;TEMP0 =0. PRESS0 =1.0E+05 DVO1DT =0. ;DRH1DP =0. RHO1A = GRND9 ;RHO1B =0. RHO1C =0. TMP1A =0. ;TMP1B =0. TMP1C =0. TMP2A =0. ;TMP2B =0. TMP2C =0. EMISS =0. ;SCATT =0. RADIA =0.35 ;RADIB =0. EL1A =0. ;EL1B =0. EL1C =0. ENUL =3.0E-05 ;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(F)=1. PRNDTL(SRAD)=-GRND1 PRT(V1)=1. ;PRT(W1)=1. PRT(KE)=1. ;PRT(EP)=1.314 PRT(H1)=1. ;PRT(F)=1. PRT(SRAD)=1.0E+10 CP1 = GRND7 ;CP2 =1. CP1A =0. ;CP1B =0. CP1C =0. ;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)=23.040001 FIINIT(EP)=4542.565918 ;FIINIT(H1)=1.267E+05 FIINIT(CH4)=1.0E-10 ;FIINIT(O2)=1.0E-10 FIINIT(H2)=1.0E-10 ;FIINIT(CO)=1.0E-10 FIINIT(H2O)=1.0E-10 ;FIINIT(CO2)=1.0E-10 FIINIT(N2)=1.0E-10 ;FIINIT(F)=1.0E-10 FIINIT(SRAD)=4.083574E+04 ;FIINIT(YSUM)=1.0E-10 FIINIT(SPH1)=1.0E-10 ;FIINIT(TMP1)=1.0E-10 FIINIT(DEN1)=1.0E-10 ;FIINIT(VIST)=1.0E-10 PATCH(INIT ,INIVAL, 1, 1, 1, 4, 1, 8, 1, 1) INIT(INIT ,W1 ,0. ,96. ) 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(RADISO ,FREEVL, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(RADISO ,H1 , GRND1 , GRND1 ) COVAL(RADISO ,SRAD, GRND1 , GRND1 ) PATCH(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(KESOURCE,KE , GRND4 , GRND4 ) COVAL(KESOURCE,EP , GRND4 , GRND4 ) PATCH(SCRSF ,LOW , 1, 1, 1, 4, 1, 1, 1, 1) COVAL(SCRSF ,P1 , FIXFLU , GRND1 ) COVAL(SCRSF ,V1 ,0. ,0. ) COVAL(SCRSF ,W1 ,0. ,96. ) COVAL(SCRSF ,KE ,0. ,23.040001 ) COVAL(SCRSF ,EP ,0. ,4542.565918 ) COVAL(SCRSF ,H1 ,0. , GRND1 ) COVAL(SCRSF ,F ,0. ,1. ) COVAL(SCRSF ,SRAD,0. ,0. ) PATCH(SCRSO ,LOW , 1, 1, 5, 8, 1, 1, 1, 1) COVAL(SCRSO ,P1 , FIXFLU , GRND1 ) COVAL(SCRSO ,V1 ,0. ,0. ) COVAL(SCRSO ,W1 ,0. ,9.9 ) COVAL(SCRSO ,KE ,0. ,0.245025 ) COVAL(SCRSO ,EP ,0. ,1.660626 ) COVAL(SCRSO ,H1 ,0. , GRND1 ) COVAL(SCRSO ,F ,0. ,0. ) COVAL(SCRSO ,SRAD,0. ,0. ) PATCH(NWALL1 ,NWALL , 1, 1, 8, 8, 1, 8, 1, 1) COVAL(NWALL1 ,W1 , GRND2 ,0. ) COVAL(NWALL1 ,KE , GRND2 , GRND2 ) COVAL(NWALL1 ,EP , GRND2 , GRND2 ) PATCH(RWALLN ,NORTH , 1, 1, 8, 8, 1, 8, 1, 1) COVAL(RWALLN ,SRAD,1. ,1097.476929 ) PATCH(RADINF ,LOW , 1, 1, 1, 4, 1, 1, 1, 1) COVAL(RADINF ,SRAD,0.5 ,3.349987E+05 ) PATCH(RADINO ,LOW , 1, 1, 5, 8, 1, 1, 1, 1) COVAL(RADINO ,SRAD,0.5 ,2.041787E+05 ) PATCH(OUT ,HIGH , 1, 1, 1, 8, 8, 8, 1, 1) COVAL(OUT ,P1 ,10. ,0. ) COVAL(OUT ,V1 ,0. ,0. ) COVAL(OUT ,W1 ,0. ,0. ) COVAL(OUT ,KE ,0. , SAME ) COVAL(OUT ,EP ,0. , SAME ) COVAL(OUT ,H1 ,0. , SAME ) COVAL(OUT ,F ,0. ,0. ) COVAL(OUT ,SRAD,0. , SAME ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 200 ;ISWC1 = 1 LITHYD = 1 ;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(F)=20 ;LITER(SRAD)=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(F)=1.0E-03 ;ENDIT(SRAD)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,1.) RELAX(V1,FALSDT,5.208333E-03) RELAX(W1,FALSDT,5.208333E-03) RELAX(KE,LINRLX,0.4) RELAX(EP,LINRLX,0.4) RELAX(H1,FALSDT,1.) RELAX(CH4,LINRLX,1.) RELAX(O2,LINRLX,1.) RELAX(H2,LINRLX,1.) RELAX(CO,LINRLX,1.) RELAX(H2O,LINRLX,1.) RELAX(CO2,LINRLX,1.) RELAX(N2,LINRLX,1.) RELAX(F,LINRLX,0.8) RELAX(SRAD,FALSDT,1.) RELAX(YSUM,LINRLX,1.) RELAX(SPH1,LINRLX,1.) RELAX(TMP1,LINRLX,1.) RELAX(DEN1,LINRLX,0.5) RELAX(VIST,LINRLX,1.) KELIN = 1 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(CH4)=1. ;VARMIN(CH4)=0. VARMAX(O2)=1. ;VARMIN(O2)=0. VARMAX(H2)=1. ;VARMIN(H2)=0. VARMAX(CO)=1. ;VARMIN(CO)=0. VARMAX(H2O)=1. ;VARMIN(H2O)=0. VARMAX(CO2)=1. ;VARMIN(CO2)=0. VARMAX(N2)=1. ;VARMIN(N2)=0. VARMAX(F)=1.0E+10 ;VARMIN(F)=-1.0E+10 VARMAX(SRAD)=1.0E+10 ;VARMIN(SRAD)=-1.0E+10 VARMAX(YSUM)=1.0E+10 ;VARMIN(YSUM)=-1.0E+10 VARMAX(SPH1)=1.0E+10 ;VARMIN(SPH1)=-1.0E+10 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=1.0E-10 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=1.0E-10 VARMAX(VIST)=1.0E+10 ;VARMIN(VIST)=-1.0E+10 ************************************************************ Group 19. Data transmitted to GROUND GENK = T PARSOL = F ISG62 = 1 CHSOB =16. CSG4 ='scrs' SPEDAT(SET,SCRS,IREACT,I,1) SPEDAT(SET,SCRS,NRSYS,I,-3) SPEDAT(SET,SCRS,NREAC,I,2) SPEDAT(SET,SCRS,NELEM,I,4) SPEDAT(SET,SCRS,IPDF,I,0) SPEDAT(SET,SCRS,CEBP,R,0.) SPEDAT(SET,SCRS,CEBS1,R,0.) SPEDAT(SET,SCRS,CEBS2,R,0.) SPEDAT(SET,SCRS,IFULIN,I,0) SPEDAT(SET,SCRS,IPRMOD,I,0) SPEDAT(SET,SCRS,ISRMOD,I,0) SPEDAT(SET,SCRS,IPCHEM,I,0) SPEDAT(SET,SCRS,ISCHEM,I,0) SPEDAT(SET,SCRS,PREXP,R,1.585E+13) SPEDAT(SET,SCRS,FUEXP,R,0.7) SPEDAT(SET,SCRS,OXEXP,R,0.8) SPEDAT(SET,SCRS,EACTP,R,4.8E+04) SPEDAT(SET,SCRS,PREXS1,R,3.98E+14) SPEDAT(SET,SCRS,FUEXS1,R,1.) SPEDAT(SET,SCRS,OXEXS1,R,0.25) SPEDAT(SET,SCRS,EACTS1,R,4.0E+04) SPEDAT(SET,SCRS,PREXS2,R,0.) SPEDAT(SET,SCRS,FUEXS2,R,1.) SPEDAT(SET,SCRS,OXEXS2,R,1.) SPEDAT(SET,SCRS,EACTS2,R,0.) SPEDAT(SET,SCRS,NSPEC,I,7) SPEDAT(SET,SCRS,COXID,C,O2) SPEDAT(SET,SCRS,CFP1,C,CO) SPEDAT(SET,SCRS,CFP2,C,H2) SPEDAT(SET,SCRS,YFU1,R,0.) SPEDAT(SET,SCRS,YFU2,R,0.) SPEDAT(SET,SCRS,YFU3,R,-0.089) SPEDAT(SET,SCRS,YFU4,R,0.091) SPEDAT(SET,SCRS,YFU5,R,0.164) SPEDAT(SET,SCRS,YFU6,R,0.036) SPEDAT(SET,SCRS,YFU7,R,0.62) SPEDAT(SET,SCRS,TEMPFU,R,2027.) SPEDAT(SET,SCRS,YOX1,R,0.) SPEDAT(SET,SCRS,YOX2,R,0.086) SPEDAT(SET,SCRS,YOX3,R,0.) SPEDAT(SET,SCRS,YOX4,R,0.) SPEDAT(SET,SCRS,YOX5,R,0.111) SPEDAT(SET,SCRS,YOX6,R,0.055) SPEDAT(SET,SCRS,YOX7,R,0.748) SPEDAT(SET,SCRS,TEMPOX,R,1791.) SPEDAT(SET,CHEM,SIUNIT,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=28. ;EX(V1)=0.6527 EX(W1)=33.299999 ;EX(KE)=96.269997 EX(EP)=1.521E+04 ;EX(H1)=3.221E+05 EX(CH4)=0. ;EX(O2)=0.02489 EX(H2)=1.093E-03 ;EX(CO)=0.01553 EX(H2O)=0.1093 ;EX(CO2)=0.121 EX(N2)=0.7281 ;EX(F)=0.4214 EX(SRAD)=8.462E+04 ;EX(YSUM)=1. EX(SPH1)=1482. ;EX(TMP1)=2024. EX(DEN1)=0.1637 ;EX(VIST)=0.07197 ************************************************************ 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(CH4,N,N,Y,N,N,N) OUTPUT(O2,Y,N,Y,Y,N,N) OUTPUT(H2,Y,N,Y,Y,N,N) OUTPUT(CO,Y,N,Y,Y,N,N) OUTPUT(H2O,Y,N,Y,Y,N,N) OUTPUT(CO2,Y,N,Y,Y,N,N) OUTPUT(N2,Y,N,Y,Y,N,N) OUTPUT(F,Y,N,Y,Y,Y,Y) OUTPUT(SRAD,Y,N,Y,Y,Y,Y) OUTPUT(YSUM,Y,N,Y,N,Y,Y) OUTPUT(SPH1,Y,N,Y,N,N,N) OUTPUT(TMP1,Y,N,Y,N,Y,Y) OUTPUT(DEN1,Y,N,Y,N,Y,Y) OUTPUT(VIST,Y,N,Y,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 6 ;IZMON = 7 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 = F IPLTF = 1 ;IPLTL = -1 ;NPLT = 10 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 =coh2 STOP