TALK=T;RUN(1,1) DISPLAY 1D Laminar Premixed Burner-Stabilised H2-Air Flame ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ The case considered is a steady, 1D, laminar, burner-stabilised, pre-mixed, stoichiometric H2-air flame at atmospheric pressure. The reactants enter the domain at 298 K with a velocity of 120 cm/s. The flame is stabilised on a porous-plug burner in the sense that the heat loss necessary to stabilise the flame is achieved through conductive heat transfer to the inlet boundary. The present case involves 8 species and 11 reactions, and Ficks law is used to represent species diffusion. The computation may be performed using the PHOENICS point-by-point (PBP) solver, or alternatively the more efficient CHEMKIN TWOPNT solver which solves the species and TEM1 equations simultaneously in a PBP manner. ENDDIS AUTOPLOT USE file PHI 5 clear d 1 MH2;d 1 MO2;d 1 MH2O;d 1 TEM1; mult y 0.0001 4 plot 1 3;col3 4;colf 4;blb3 3;scale x 0. 0.2 msg H2, O2, H2O mole-fraction & Temp(/1000) profiles msg Pressto continue pause clear d 1 MH;d 1 MO;d 1 MOH;d 1 MHO2 mult y 10. 2;mult y 10. 3;mult y 1000. 4 col3 1 4;col6 2;col9 3;colf 4;scale x 0. 0.2 msg H, O(*10), OH (*10) & HO2(*1000) mole fraction profiles msg Press to continue pause clear msg Heat-release rate profile in erg/s/cm**3 d 1 QDOT;scale x 0. 0.2;plot 1;blb1 1 enduse ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT( CHEMKIN - 1DY Premixed H2-Air Flame ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 204 ************************************************************ 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 = 40 YVLAST =0.6 YFRAC(1)=2.733405E-03 ;YFRAC(2)=8.286134E-03 YFRAC(3)=0.015853 ;YFRAC(4)=0.025119 YFRAC(5)=0.035897 ;YFRAC(6)=0.048056 YFRAC(7)=0.061498 ;YFRAC(8)=0.076146 YFRAC(9)=0.091937 ;YFRAC(10)=0.108819 YFRAC(11)=0.126745 ;YFRAC(12)=0.145678 YFRAC(13)=0.165582 ;YFRAC(14)=0.186427 YFRAC(15)=0.208185 ;YFRAC(16)=0.230832 YFRAC(17)=0.254344 ;YFRAC(18)=0.278702 YFRAC(19)=0.303885 ;YFRAC(20)=0.329877 YFRAC(21)=0.35666 ;YFRAC(22)=0.38422 YFRAC(23)=0.412542 ;YFRAC(24)=0.441613 YFRAC(25)=0.47142 ;YFRAC(26)=0.501951 YFRAC(27)=0.533195 ;YFRAC(28)=0.565141 YFRAC(29)=0.597779 ;YFRAC(30)=0.6311 YFRAC(31)=0.665093 ;YFRAC(32)=0.699752 YFRAC(33)=0.735066 ;YFRAC(34)=0.771028 YFRAC(35)=0.807631 ;YFRAC(36)=0.844866 YFRAC(37)=0.882728 ;YFRAC(38)=0.921208 YFRAC(39)=0.960301 ;YFRAC(40)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 1 ZWLAST =1. ZFRAC(1)=1. ************************************************************ Group 6. Body-Fitted Coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(1)=P1 ;NAME(5)=V1 NAME(16)=H2 ;NAME(17)=H NAME(18)=O2 ;NAME(19)=O NAME(20)=OH ;NAME(21)=HO2 NAME(22)=H2O ;NAME(23)=N2 NAME(133)=MHO2 ;NAME(134)=MO NAME(135)=MH ;NAME(136)=MOH NAME(137)=MH2O ;NAME(138)=MO2 NAME(139)=MH2 ;NAME(140)=ELN NAME(141)=ELO ;NAME(142)=PRPS NAME(143)=ELH ;NAME(144)=SPH1 NAME(145)=KOND ;NAME(146)=ENTH NAME(147)=QDOT ;NAME(148)=DEN1 NAME(149)=VISL ;NAME(150)=TEM1 * 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,N,N,N,Y) SOLUTN(V1,Y,Y,N,N,N,Y) SOLUTN(H2,Y,Y,Y,N,N,N) SOLUTN(H,Y,Y,Y,N,N,N) SOLUTN(O2,Y,Y,Y,N,N,N) SOLUTN(O,Y,Y,Y,N,N,N) SOLUTN(OH,Y,Y,Y,N,N,N) SOLUTN(HO2,Y,Y,Y,N,N,N) SOLUTN(H2O,Y,Y,Y,N,N,N) SOLUTN(N2,Y,N,N,N,N,Y) SOLUTN(MHO2,Y,N,N,N,N,Y) SOLUTN(MO,Y,N,N,N,N,Y) SOLUTN(MH,Y,N,N,N,N,Y) SOLUTN(MOH,Y,N,N,N,N,Y) SOLUTN(MH2O,Y,N,N,N,N,Y) SOLUTN(MO2,Y,N,N,N,N,Y) SOLUTN(MH2,Y,N,N,N,N,Y) SOLUTN(ELN,Y,N,N,N,N,Y) SOLUTN(ELO,Y,N,N,N,N,Y) SOLUTN(PRPS,Y,N,N,N,N,Y) SOLUTN(ELH,Y,N,N,N,N,Y) SOLUTN(SPH1,Y,N,N,N,N,Y) SOLUTN(KOND,Y,N,N,N,N,Y) SOLUTN(ENTH,Y,N,N,N,N,Y) SOLUTN(QDOT,Y,N,N,N,N,Y) SOLUTN(DEN1,Y,N,N,N,N,Y) SOLUTN(VISL,Y,N,N,N,N,Y) SOLUTN(TEM1,Y,Y,N,N,N,N) DEN1 = 148 VISL = 149 PRPS = 142 ************************************************************ 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(H2,N,Y,Y,Y,Y,Y) TERMS(H,N,Y,Y,Y,N,Y) TERMS(O2,N,Y,Y,Y,Y,Y) TERMS(O,N,Y,Y,Y,N,Y) TERMS(OH,N,Y,Y,Y,Y,Y) TERMS(HO2,N,Y,Y,Y,N,Y) TERMS(H2O,N,Y,Y,Y,Y,Y) TERMS(TEM1,N,Y,Y,Y,Y,N) DIFCUT =0. ;ZDIFAC =1. GALA = F ;ADDDIF = F NEWRH1 = T NEWENL = 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 = GRND9 ;TMP1 =0. ;EL1 =0. TSURR =0. ;TEMP0 =0. ;PRESS0 =0. DVO1DT =0. ;DRH1DP =0. RHO1A =0. ;RHO1B =0. ;RHO1C =0. EMISS =0. ;SCATT =0. RADIA =0. ;RADIB =0. ENUL = GRND9 ;ENUT =0. ENULA =0. ;ENULB =0. ;ENULC =0. ;ENULD =0. ENULE =0. ;ENULF =0. ;ENULG =0. IENULA = 0 ;IENULB = 0 PRNDTL(V1)=1. ;PRNDTL(H2)=1. PRNDTL(H)=1. ;PRNDTL(O2)=1. PRNDTL(O)=1. ;PRNDTL(OH)=1. PRNDTL(HO2)=1. ;PRNDTL(H2O)=1. PRNDTL(TEM1)=-GRND9 PRLH1A =0. ;PRLH1B =0. ;PRLH1C =0. PRT(V1)=1. ;PRT(H2)=1. PRT(H)=1. ;PRT(O2)=1. PRT(O)=1. ;PRT(OH)=1. PRT(HO2)=1. ;PRT(H2O)=1. PRT(TEM1)=1. CP1 = GRND9 ;CP2 =1. CP1A =0. ;CP1B =0. ;CP1C =0. ;CP1D =0. * List of user-defined materials to be read by EARTH MATFLG=T;IMAT=1 * Name *Ind. Dens. Viscos. Spec.heat Conduct. Expans. Compr. * 71 GRND9 GRND9 GRND9 GRND9 0. * constants for GRND option no 1 0.0 * constants for GRND option no 2 0.0 * constants for GRND option no 3 0.0 * constants for GRND option no 4 0.0 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(P1)=1.0E-10 ;FIINIT(V1)=120. FIINIT(H2)=1.0E-04 ;FIINIT(H)=1.0E-03 FIINIT(O2)=1.0E-04 ;FIINIT(O)=0.02 FIINIT(OH)=0.015 ;FIINIT(HO2)=1.0E-04 FIINIT(H2O)=0.129368 ;FIINIT(N2)=0.834232 FIINIT(MHO2)=1.0E-10 ;FIINIT(MO)=1.0E-10 FIINIT(MH)=1.0E-10 ;FIINIT(MOH)=1.0E-10 FIINIT(MH2O)=1.0E-10 ;FIINIT(MO2)=1.0E-10 FIINIT(MH2)=1.0E-10 ;FIINIT(ELN)=1.0E-10 FIINIT(ELO)=1.0E-10 ;FIINIT(PRPS)=71. FIINIT(ELH)=1.0E-10 ;FIINIT(SPH1)=1.0E-10 FIINIT(KOND)=1.0E-10 ;FIINIT(ENTH)=1.0E-10 FIINIT(QDOT)=1.0E-10 ;FIINIT(DEN1)=1.0E-10 FIINIT(VISL)=1.0E-10 ;FIINIT(TEM1)=2250. PATCH(ICHEMK1 ,INIVAL, 1, 1, 1, 40, 1, 1, 1, 1) INIT(ICHEMK1 ,DEN1,0. , GRND1 ) 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(NOCPCK1 ,SOUTH , 1, 1, 1, 1, 1, 1, 1, 1) COVAL(NOCPCK1 ,P1 , FIXFLU , GRND9 ) COVAL(NOCPCK1 ,V1 ,0. ,120. ) COVAL(NOCPCK1 ,H2 ,0. ,0.028522 ) COVAL(NOCPCK1 ,H ,0. ,0. ) COVAL(NOCPCK1 ,O2 ,0. ,0.226364 ) COVAL(NOCPCK1 ,O ,0. ,0. ) COVAL(NOCPCK1 ,OH ,0. ,0. ) COVAL(NOCPCK1 ,HO2 ,0. ,0. ) COVAL(NOCPCK1 ,H2O ,0. ,0. ) COVAL(NOCPCK1 ,N2 ,0. ,0. ) COVAL(NOCPCK1 ,TEM1,0. , GRND9 ) PATCH(FLUXIN ,SWALL , 1, 1, 1, 1, 1, 1, 1, 1) COVAL(FLUXIN ,TEM1,1. ,298. ) PATCH(OUT1 ,NORTH , 1, 1, 40, 40, 1, 1, 1, 1) COVAL(OUT1 ,P1 ,1000. ,0. ) COVAL(OUT1 ,V1 ,0. ,0. ) COVAL(OUT1 ,H2 ,0. , SAME ) COVAL(OUT1 ,H ,0. , SAME ) COVAL(OUT1 ,O2 ,0. , SAME ) COVAL(OUT1 ,O ,0. , SAME ) COVAL(OUT1 ,OH ,0. , SAME ) COVAL(OUT1 ,HO2 ,0. , SAME ) COVAL(OUT1 ,H2O ,0. , SAME ) COVAL(OUT1 ,TEM1,0. , SAME ) 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-05 ************************************************************ Group 16. Terminate Iterations LITER(P1)=20 ;LITER(V1)=10 LITER(H2)=20 ;LITER(H)=20 LITER(O2)=20 ;LITER(O)=20 LITER(OH)=20 ;LITER(HO2)=20 LITER(H2O)=20 ;LITER(TEM1)=20 ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03 ENDIT(H2)=1.0E-06 ;ENDIT(H)=1.0E-06 ENDIT(O2)=1.0E-06 ;ENDIT(O)=1.0E-06 ENDIT(OH)=1.0E-06 ;ENDIT(HO2)=1.0E-06 ENDIT(H2O)=1.0E-06 ;ENDIT(TEM1)=1.0E-06 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,1.) RELAX(V1,FALSDT,6.0E-05) RELAX(H2,FALSDT,0.06) RELAX(H,FALSDT,0.06) RELAX(O2,FALSDT,0.06) RELAX(O,FALSDT,0.06) RELAX(OH,FALSDT,0.06) RELAX(HO2,FALSDT,0.06) RELAX(H2O,FALSDT,0.06) RELAX(N2,LINRLX,1.) RELAX(MHO2,LINRLX,1.) RELAX(MO,LINRLX,1.) RELAX(MH,LINRLX,1.) RELAX(MOH,LINRLX,1.) RELAX(MH2O,LINRLX,1.) RELAX(MO2,LINRLX,1.) RELAX(MH2,LINRLX,1.) RELAX(ELN,LINRLX,1.) RELAX(ELO,LINRLX,1.) RELAX(PRPS,LINRLX,1.) RELAX(ELH,LINRLX,1.) RELAX(SPH1,LINRLX,1.) RELAX(KOND,LINRLX,1.) RELAX(ENTH,LINRLX,1.) RELAX(QDOT,LINRLX,1.) RELAX(DEN1,LINRLX,1.) RELAX(VISL,LINRLX,1.) RELAX(TEM1,FALSDT,1.0E-05) 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(H2)=1. ;VARMIN(H2)=0. VARMAX(H)=1. ;VARMIN(H)=0. VARMAX(O2)=1. ;VARMIN(O2)=0. VARMAX(O)=1. ;VARMIN(O)=0. VARMAX(OH)=1. ;VARMIN(OH)=0. VARMAX(HO2)=1. ;VARMIN(HO2)=0. VARMAX(H2O)=1. ;VARMIN(H2O)=0. VARMAX(N2)=1. ;VARMIN(N2)=0. VARMAX(MHO2)=1.0E+10 ;VARMIN(MHO2)=-1.0E+10 VARMAX(MO)=1.0E+10 ;VARMIN(MO)=-1.0E+10 VARMAX(MH)=1.0E+10 ;VARMIN(MH)=-1.0E+10 VARMAX(MOH)=1.0E+10 ;VARMIN(MOH)=-1.0E+10 VARMAX(MH2O)=1.0E+10 ;VARMIN(MH2O)=-1.0E+10 VARMAX(MO2)=1.0E+10 ;VARMIN(MO2)=-1.0E+10 VARMAX(MH2)=1.0E+10 ;VARMIN(MH2)=-1.0E+10 VARMAX(ELN)=1.0E+10 ;VARMIN(ELN)=-1.0E+10 VARMAX(ELO)=1.0E+10 ;VARMIN(ELO)=-1.0E+10 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10 VARMAX(ELH)=1.0E+10 ;VARMIN(ELH)=-1.0E+10 VARMAX(SPH1)=1.0E+10 ;VARMIN(SPH1)=-1.0E+10 VARMAX(KOND)=1.0E+10 ;VARMIN(KOND)=-1.0E+10 VARMAX(ENTH)=1.0E+10 ;VARMIN(ENTH)=-1.0E+10 VARMAX(QDOT)=1.0E+10 ;VARMIN(QDOT)=-1.0E+10 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10 VARMAX(VISL)=1.0E+10 ;VARMIN(VISL)=-1.0E+10 VARMAX(TEM1)=1.0E+10 ;VARMIN(TEM1)=273. ************************************************************ Group 19. Data transmitted to GROUND PARSOL = F CHMKIN = T ISG62 = 1 CHSOA = GRND9 ;CHSOB =16. CHSOC =1. CSG4 ='ho11' CSG10 ='q1' SPEDAT(SET,NOCPCK1,TINLET,R,298.) SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,MATERIAL,71,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=13.3 ;EX(V1)=692.599976 EX(H2)=4.629E-03 ;EX(H)=4.477E-04 EX(O2)=0.03609 ;EX(O)=1.525E-03 EX(OH)=5.686E-03 ;EX(HO2)=3.197E-06 EX(H2O)=0.2065 ;EX(N2)=0.7451 EX(MHO2)=2.134E-06 ;EX(MO)=2.219E-03 EX(MH)=0.0102 ;EX(MOH)=7.992E-03 EX(MH2O)=0.2747 ;EX(MO2)=0.02486 EX(MH2)=0.05084 ;EX(ELN)=0.7451 EX(ELO)=0.2264 ;EX(PRPS)=0. EX(ELH)=0.02852 ;EX(SPH1)=1.232E+07 EX(KOND)=1.691E+04 ;EX(ENTH)=1.91E+05 EX(QDOT)=1.003E+10 ;EX(DEN1)=1.822E-04 EX(VISL)=4.671 ;EX(TEM1)=1975. ************************************************************ 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(H2,Y,N,Y,Y,Y,Y) OUTPUT(H,Y,N,Y,Y,Y,Y) OUTPUT(O2,Y,N,Y,Y,Y,Y) OUTPUT(O,Y,N,Y,Y,Y,Y) OUTPUT(OH,Y,N,Y,Y,Y,Y) OUTPUT(HO2,Y,N,Y,Y,Y,Y) OUTPUT(H2O,Y,N,Y,Y,Y,Y) OUTPUT(N2,Y,N,Y,N,Y,N) OUTPUT(MHO2,Y,N,Y,N,N,N) OUTPUT(MO,Y,N,Y,N,N,N) OUTPUT(MH,Y,N,Y,N,N,N) OUTPUT(MOH,Y,N,Y,N,N,N) OUTPUT(MH2O,Y,N,Y,N,N,N) OUTPUT(MO2,Y,N,Y,N,N,N) OUTPUT(MH2,Y,N,Y,N,N,N) OUTPUT(ELN,Y,N,Y,N,N,N) OUTPUT(ELO,Y,N,Y,N,N,N) OUTPUT(PRPS,N,N,Y,N,N,N) OUTPUT(ELH,Y,N,Y,N,N,N) OUTPUT(SPH1,Y,N,Y,N,N,N) OUTPUT(KOND,Y,N,Y,N,N,N) OUTPUT(ENTH,Y,N,Y,N,N,N) OUTPUT(QDOT,Y,N,Y,N,N,N) OUTPUT(DEN1,Y,Y,Y,N,Y,N) OUTPUT(VISL,Y,N,Y,N,N,N) OUTPUT(TEM1,Y,N,Y,Y,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 39 ;IZMON = 1 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1 UWATCH = T ;USTEER = T HIGHLO = F ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 250 ;NUMCLS = 5 NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 40 IPLTF = 1 ;IPLTL = -1 ;NPLT = 25 ISWPRF = 1 ;ISWPRL = 100000 ITABL = 3 ;IPROF = 1 ABSIZ =0.5 ;ORSIZ =0.4 NTZPRF = 1 ;NCOLPF = 50 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20 PATCH(PROFSY ,PROFIL, 1, 1, 1, 40, 1, 1, 1, 1) PLOT(PROFSY ,H2 ,0. ,0. ) PLOT(PROFSY ,O2 ,0. ,0. ) PATCH(PROFSYB ,PROFIL, 1, 1, 1, 40, 1, 1, 1, 1) PLOT(PROFSYB ,N2 ,0. ,0. ) PLOT(PROFSYB ,TEM1,0. ,0. ) PATCH(PROFSYA ,PROFIL, 1, 1, 1, 40, 1, 1, 1, 1) PLOT(PROFSYA ,H ,0. ,0. ) PLOT(PROFSYA ,O ,0. ,0. ) PLOT(PROFSYA ,OH ,0. ,0. ) PLOT(PROFSYA ,HO2 ,0. ,0. ) PLOT(PROFSYA ,H2O ,0. ,0. ) ************************************************************ Group 24. Dumps For Restarts SAVE = T ;NOWIPE = F NSAVE =CHAM STOP