TALK=T;RUN(1,1) DISPLAY The problem considered is turbulent reacting flow in a cylindrical combustion chamber with central and annular coaxial feed injectors supplying premixed methane and air. The geometry is axisymmetric with an outer radius of 0.2m and an axial length of 2m. The reactants are supplied at a temperature of 300K, and the composition of each inlet stream is identical. The default fuel/air equivalence ratio is 1.15, and the operating pressure is 1 bar. ENDDIS ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(Premixed 2Step CH4 Combustion ) ************************************************************ ************************************************************ 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.2073 YFRAC(1)=0.046068 ;YFRAC(2)=0.092137 YFRAC(3)=0.161119 ;YFRAC(4)=0.230101 YFRAC(5)=0.362277 ;YFRAC(6)=0.494452 YFRAC(7)=0.747226 ;YFRAC(8)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 10 ZWLAST =2.743 ZFRAC(1)=0.034725 ;ZFRAC(2)=0.06945 ZFRAC(3)=0.104174 ;ZFRAC(4)=0.138899 ZFRAC(5)=0.22274 ;ZFRAC(6)=0.34534 ZFRAC(7)=0.488615 ;ZFRAC(8)=0.647217 ZFRAC(9)=0.818288 ;ZFRAC(10)=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(136)=ELH NAME(137)=ELO ;NAME(138)=ELN NAME(139)=ELC ;NAME(140)=MMWT NAME(141)=S2RS ;NAME(142)=S1RS NAME(143)=P1RS ;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,Y,N,N,N,N) SOLUTN(O2,Y,N,N,N,N,N) SOLUTN(H2,Y,Y,N,N,N,N) SOLUTN(CO,Y,Y,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(ELH,Y,N,N,N,N,Y) SOLUTN(ELO,Y,N,N,N,N,Y) SOLUTN(ELN,Y,N,N,N,N,Y) SOLUTN(ELC,Y,N,N,N,N,Y) SOLUTN(MMWT,Y,N,N,N,N,Y) SOLUTN(S2RS,Y,N,N,N,N,Y) SOLUTN(S1RS,Y,N,N,N,N,Y) SOLUTN(P1RS,Y,N,N,N,N,Y) 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(CH4,N,Y,Y,Y,Y,Y) TERMS(H2,N,Y,Y,Y,Y,Y) TERMS(CO,N,Y,Y,Y,N,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 =1.45 ;RADIB =0. EL1A =0. ;EL1B =0. EL1C =0. ENUL =4.2E-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(CH4)=1. PRNDTL(H2)=1. ;PRNDTL(CO)=1. PRNDTL(F)=1. ;PRNDTL(SRAD)=-GRND1 PRT(V1)=1. ;PRT(W1)=1. PRT(KE)=1. ;PRT(EP)=1.314 PRT(H1)=1. ;PRT(CH4)=1. PRT(H2)=1. ;PRT(CO)=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)=20. ;FIINIT(KE)=4.999999E-03 FIINIT(EP)=0.022981 ;FIINIT(H1)=-1.39E+05 FIINIT(CH4)=0.07 ;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(ELH)=1.0E-10 FIINIT(ELO)=1.0E-10 ;FIINIT(ELN)=1.0E-10 FIINIT(ELC)=1.0E-10 ;FIINIT(MMWT)=1.0E-10 FIINIT(S2RS)=1.0E-10 ;FIINIT(S1RS)=1.0E-10 FIINIT(P1RS)=1.0E-10 ;FIINIT(F)=1.0E-10 FIINIT(SRAD)=248.049393 ;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 No PATCHes yet used for this Group 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(RADISO ,FREEVL, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(RADISO ,H1 , GRND1 , GRND1 ) COVAL(RADISO ,SRAD, GRND1 , GRND1 ) PATCH(SCRSPRSO,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(SCRSPRSO,CH4 , GRND2 , GRND2 ) COVAL(SCRSPRSO,H2 , FIXFLU , GRND2 ) COVAL(SCRSPRSO,CO , FIXFLU , GRND2 ) PATCH(SCRSSRSO,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(SCRSSRSO,H2 , GRND5 , GRND5 ) COVAL(SCRSSRSO,CO , GRND5 , GRND5 ) PATCH(SCRSF ,LOW , 1, 1, 3, 4, 1, 1, 1, 1) COVAL(SCRSF ,P1 , FIXFLU , GRND1 ) COVAL(SCRSF ,V1 ,0. ,0. ) COVAL(SCRSF ,W1 ,0. ,115. ) COVAL(SCRSF ,KE ,0. ,14.77505 ) COVAL(SCRSF ,EP ,0. ,9822.168945 ) COVAL(SCRSF ,H1 ,0. , GRND3 ) COVAL(SCRSF ,CH4 ,0. ,0.07 ) COVAL(SCRSF ,H2 ,0. ,0. ) COVAL(SCRSF ,CO ,0. ,0. ) COVAL(SCRSF ,F ,0. ,1. ) COVAL(SCRSF ,SRAD,0. ,0. ) PATCH(SCRSO ,LOW , 1, 1, 7, 8, 1, 1, 1, 1) COVAL(SCRSO ,P1 , FIXFLU , GRND1 ) COVAL(SCRSO ,V1 ,0. ,0. ) COVAL(SCRSO ,W1 ,0. ,10. ) COVAL(SCRSO ,KE ,0. ,0.763848 ) COVAL(SCRSO ,EP ,0. ,21.937037 ) COVAL(SCRSO ,H1 ,0. , GRND3 ) COVAL(SCRSO ,CH4 ,0. ,0.0275 ) COVAL(SCRSO ,H2 ,0. ,0. ) COVAL(SCRSO ,CO ,0. ,0. ) COVAL(SCRSO ,F ,0. ,0. ) COVAL(SCRSO ,SRAD,0. ,0. ) PATCH(OUT ,HIGH , 1, 1, 1, 8, 10, 10, 1, 1) COVAL(OUT ,P1 ,10. ,0. ) COVAL(OUT ,H1 ,0. , SAME ) COVAL(OUT ,F ,0. , SAME ) PATCH(NWALL3 ,NWALL , 1, 1, 8, 8, 1, 10, 1, 1) COVAL(NWALL3 ,W1 , GRND2 ,0. ) COVAL(NWALL3 ,KE , GRND2 , GRND2 ) COVAL(NWALL3 ,EP , GRND2 , GRND2 ) PATCH(NWALL3R ,NORTH , 1, 1, 8, 8, 1, 10, 1, 1) COVAL(NWALL3R ,SRAD,1. ,1097.476929 ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 250 ;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(CH4)=20 ;LITER(H2)=20 LITER(CO)=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(CH4)=1.0E-03 ;ENDIT(H2)=1.0E-03 ENDIT(CO)=1.0E-03 ;ENDIT(F)=1.0E-03 ENDIT(SRAD)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,1.) RELAX(V1,FALSDT,5.0E-03) RELAX(W1,FALSDT,5.0E-03) RELAX(KE,FALSDT,5.0E-03) RELAX(EP,FALSDT,5.0E-03) RELAX(H1,FALSDT,1.) RELAX(CH4,FALSDT,1.0E-03) RELAX(O2,LINRLX,1.) RELAX(H2,FALSDT,1.0E-03) RELAX(CO,FALSDT,1.0E-03) RELAX(H2O,LINRLX,1.) RELAX(CO2,LINRLX,1.) RELAX(N2,LINRLX,1.) RELAX(ELH,LINRLX,1.) RELAX(ELO,LINRLX,1.) RELAX(ELN,LINRLX,1.) RELAX(ELC,LINRLX,1.) RELAX(MMWT,LINRLX,1.) RELAX(S2RS,LINRLX,1.) RELAX(S1RS,LINRLX,1.) RELAX(P1RS,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 = 0 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(ELH)=1.0E+10 ;VARMIN(ELH)=-1.0E+10 VARMAX(ELO)=1.0E+10 ;VARMIN(ELO)=-1.0E+10 VARMAX(ELN)=1.0E+10 ;VARMIN(ELN)=-1.0E+10 VARMAX(ELC)=1.0E+10 ;VARMIN(ELC)=-1.0E+10 VARMAX(MMWT)=1.0E+10 ;VARMIN(MMWT)=-1.0E+10 VARMAX(S2RS)=1.0E+10 ;VARMIN(S2RS)=-1.0E+10 VARMAX(S1RS)=1.0E+10 ;VARMIN(S1RS)=-1.0E+10 VARMAX(P1RS)=1.0E+10 ;VARMIN(P1RS)=-1.0E+10 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,2) SPEDAT(SET,SCRS,NRSYS,I,-6) SPEDAT(SET,SCRS,NREAC,I,3) SPEDAT(SET,SCRS,NELEM,I,4) SPEDAT(SET,SCRS,IPDF,I,0) SPEDAT(SET,SCRS,CEBP,R,4.) SPEDAT(SET,SCRS,CEBS1,R,4.) SPEDAT(SET,SCRS,CEBS2,R,4.) 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,CFUEL,C,CH4) SPEDAT(SET,SCRS,COXID,C,O2) SPEDAT(SET,SCRS,CFP1,C,CO) SPEDAT(SET,SCRS,CFP2,C,H2) SPEDAT(SET,SCRS,YFU1,R,0.07) SPEDAT(SET,SCRS,YFU2,R,0.1764) SPEDAT(SET,SCRS,YFU3,R,0.) SPEDAT(SET,SCRS,YFU4,R,0.) SPEDAT(SET,SCRS,YFU5,R,0.02118) SPEDAT(SET,SCRS,YFU6,R,0.0263) SPEDAT(SET,SCRS,YFU7,R,0.70612) SPEDAT(SET,SCRS,TEMPFU,R,470.) SPEDAT(SET,SCRS,YOX1,R,0.0275) SPEDAT(SET,SCRS,YOX2,R,0.1914) SPEDAT(SET,SCRS,YOX3,R,0.) SPEDAT(SET,SCRS,YOX4,R,0.) SPEDAT(SET,SCRS,YOX5,R,0.01862) SPEDAT(SET,SCRS,YOX6,R,0.02353) SPEDAT(SET,SCRS,YOX7,R,0.73895) SPEDAT(SET,SCRS,TEMPOX,R,503.) SPEDAT(SET,CHEM,SIUNIT,L,T) SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ 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(CH4,Y,N,Y,Y,Y,Y) OUTPUT(O2,Y,N,Y,Y,N,N) OUTPUT(H2,Y,N,Y,Y,Y,Y) OUTPUT(CO,Y,N,Y,Y,Y,Y) 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(ELH,Y,N,Y,N,N,N) OUTPUT(ELO,Y,N,Y,N,N,N) OUTPUT(ELN,Y,N,Y,N,N,N) OUTPUT(ELC,Y,N,Y,N,N,N) OUTPUT(MMWT,Y,N,Y,N,N,N) OUTPUT(S2RS,Y,N,Y,N,N,N) OUTPUT(S1RS,Y,N,Y,N,N,N) OUTPUT(P1RS,Y,N,Y,N,N,N) OUTPUT(F,Y,N,Y,Y,Y,Y) OUTPUT(SRAD,Y,N,Y,Y,Y,Y) OUTPUT(YSUM,Y,N,Y,N,N,N) 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 = 4 ;IZMON = 6 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 =CHAM STOP