DISPLAY Library case 643: Tutorial: Combustion model - ESCRS ENDDIS TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 3.6, Date 02/11/04 CPVNAM=VDI;SPPNAM=Core ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(ESCRS - Fast Chemistry ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,1,30,40) * Cylindrical-polar grid CARTES=F ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T * Non-default variable names NAME( 16) =CH4 ; NAME( 17) =O2 NAME( 18) =H2 ; NAME( 19) =CO NAME( 20) =H2O ; NAME( 21) =CO2 NAME( 22) =N2 ; NAME(140) =TMP1 NAME(141) =ENUL ; NAME(142) =MCO2 NAME(143) =MCO ; NAME(144) =MMWT NAME(145) =SPH1 ; NAME(146) =F NAME(147) =EPKE ; NAME(148) =DEN1 NAME(149) =EL1 ; NAME(150) =ENUT * Solved variables list SOLVE(P1 ,V1 ,W1 ,H1 ,F ) * Stored variables list STORE(ENUT,EL1 ,DEN1,EPKE,SPH1,MMWT,MCO ,MCO2) STORE(ENUL,TMP1,N2 ,CO2 ,H2O ,CO ,H2 ,O2 ) STORE(CH4 ) * Additional solver options SOLUTN(P1 ,Y,Y,Y,N,N,Y) SOLUTN(H1 ,Y,Y,Y,N,N,Y) SOLUTN(F ,Y,Y,Y,N,N,Y) TURMOD(KEMODL) ************************************************************ Group 8. Terms & Devices TERMS (H1 ,N,Y,Y,Y,Y,Y) NEWRH1 = T DENPCO = T ************************************************************ Group 9. Properties PRESS0 = 3.800000E+05 RHO1 = GRND9 RHO1A = GRND9 ;RHO1B = 0.000000E+00 RHO1C = 0.000000E+00 ENUL = 4.250000E-05 CP1 = GRND7 TMP1 = GRND9 PRT (EP ) = 1.314000E+00 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1 ) = -2.900000E+01 ;FIINIT(W1 ) = 5.500000E+00 FIINIT(KE ) = 9.458000E+00 ;FIINIT(EP ) = 4.725000E+03 FIINIT(H1 ) = -4.865000E+05 ;FIINIT(F ) = 1.860000E-01 No PATCHes used for this Group INIADD = F ************************************************************ Group 12. Convection and diffusion adjustments No PATCHes used for this Group ************************************************************ Group 13. Boundary & Special Sources INLET (SCRSO ,LOW ,1,0,0,0,0,0,1,1) VALUE (SCRSO ,P1 , GRND1 ) VALUE (SCRSO ,V1 , 0.000000E+00) VALUE (SCRSO ,W1 , 1.960000E+01) VALUE (SCRSO ,KE , 1.540000E+01) VALUE (SCRSO ,EP , 6.190000E+03) VALUE (SCRSO ,H1 , 0.000000E+00) VALUE (SCRSO ,F , 0.000000E+00) INLET (SCRSF ,LOW ,2,0,0,0,0,0,1,1) VALUE (SCRSF ,P1 , GRND1 ) VALUE (SCRSF ,V1 , 0.000000E+00) VALUE (SCRSF ,W1 , 9.100000E-01) VALUE (SCRSF ,KE , 3.300000E-02) VALUE (SCRSF ,EP , 3.100000E-01) VALUE (SCRSF ,H1 , 0.000000E+00) VALUE (SCRSF ,CH4 , 1.000000E+00) VALUE (SCRSF ,F , 1.000000E+00) EGWF = F ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 500 RESFAC = 1.000000E-03 ************************************************************ Group 16. Terminate Iterations LITER (P1 ) = 200 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX, 5.000000E-01) RELAX(V1 ,FALSDT, 1.000000E-03) RELAX(W1 ,FALSDT, 1.000000E-03) RELAX(KE ,LINRLX, 5.000000E-01) RELAX(EP ,LINRLX, 5.000000E-01) RELAX(H1 ,LINRLX, 5.000000E-03) RELAX(F ,LINRLX, 5.000000E-03) RELAX(DEN1,LINRLX, 1.000000E-02) KELIN = 3 ************************************************************ Group 18. Limits VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06 VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06 VARMAX(CH4 ) = 1.000000E+00 ;VARMIN(CH4 ) = 0.000000E+00 VARMAX(O2 ) = 1.000000E+00 ;VARMIN(O2 ) = 0.000000E+00 VARMAX(H2 ) = 1.000000E+00 ;VARMIN(H2 ) = 0.000000E+00 VARMAX(CO ) = 1.000000E+00 ;VARMIN(CO ) = 0.000000E+00 VARMAX(H2O ) = 1.000000E+00 ;VARMIN(H2O ) = 0.000000E+00 VARMAX(CO2 ) = 1.000000E+00 ;VARMIN(CO2 ) = 0.000000E+00 VARMAX(N2 ) = 1.000000E+00 ;VARMIN(N2 ) = 0.000000E+00 ************************************************************ Group 19. EARTH Calls To GROUND Station USEGRD = T ;USEGRX = T GENK = T ASAP = T PARSOL = T CHSOB = 1.600000E+01 CSG4 ='sch4' SPEDAT(SET,SCRS,IREACT,I,1) SPEDAT(SET,SCRS,NRSYS,I,-2) SPEDAT(SET,SCRS,NREAC,I,1) SPEDAT(SET,SCRS,NELEM,I,4) SPEDAT(SET,SCRS,IPDF,I,0) SPEDAT(SET,SCRS,CEBP,R,0.00000E+00) SPEDAT(SET,SCRS,CEBS1,R,0.00000E+00) SPEDAT(SET,SCRS,CEBS2,R,0.00000E+00) 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.58500E+13) SPEDAT(SET,SCRS,FUEXP,R,7.00000E-01) SPEDAT(SET,SCRS,OXEXP,R,8.00000E-01) SPEDAT(SET,SCRS,EACTP,R,4.80000E+04) SPEDAT(SET,SCRS,PREXS1,R,3.98000E+14) SPEDAT(SET,SCRS,FUEXS1,R,1.00000E+00) SPEDAT(SET,SCRS,OXEXS1,R,2.50000E-01) SPEDAT(SET,SCRS,EACTS1,R,4.00000E+04) SPEDAT(SET,SCRS,PREXS2,R,0.00000E+00) SPEDAT(SET,SCRS,FUEXS2,R,1.00000E+00) SPEDAT(SET,SCRS,OXEXS2,R,1.00000E+00) SPEDAT(SET,SCRS,EACTS2,R,0.00000E+00) SPEDAT(SET,SCRS,NSPEC,I,7) SPEDAT(SET,SCRS,CFUEL,C,CH4) SPEDAT(SET,SCRS,COXID,C,O2) SPEDAT(SET,SCRS,CFP1,C,CO2) SPEDAT(SET,SCRS,CFP2,C,H2O) SPEDAT(SET,CHEM,SIUNIT,L,T) SPEDAT(SET,SCRS,YFU1,R,1.00000E+00) SPEDAT(SET,SCRS,YFU2,R,0.00000E+00) SPEDAT(SET,SCRS,YFU3,R,0.00000E+00) SPEDAT(SET,SCRS,YFU4,R,0.00000E+00) SPEDAT(SET,SCRS,YFU5,R,0.00000E+00) SPEDAT(SET,SCRS,YFU6,R,0.00000E+00) SPEDAT(SET,SCRS,YFU7,R,0.00000E+00) SPEDAT(SET,SCRS,TEMPFU,R,2.88000E+02) SPEDAT(SET,SCRS,YOX1,R,0.00000E+00) SPEDAT(SET,SCRS,YOX2,R,2.32000E-01) SPEDAT(SET,SCRS,YOX3,R,0.00000E+00) SPEDAT(SET,SCRS,YOX4,R,0.00000E+00) SPEDAT(SET,SCRS,YOX5,R,0.00000E+00) SPEDAT(SET,SCRS,YOX6,R,0.00000E+00) SPEDAT(SET,SCRS,YOX7,R,7.68000E-01) SPEDAT(SET,SCRS,TEMPOX,R,7.50000E+02) SPEDAT(SET,FACETDAT,NUMOBJ,I,5) ************************************************************ Group 20. Preliminary Printout ECHO = T NULLPR = F NDST = 0 DSTTOL = 1.000000E-02 ************************************************************ Group 21. Print-out of Variables OUTPUT(CH4 ,Y,N,Y,Y,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) ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 13 ;IZMON = 20 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts GVIEW(P,-9.999875E-01,0.000000E+00,4.999979E-03) GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00) > DOM, SIZE, 1.000000E-01, 6.100000E-02, 5.000000E-01 > DOM, MONIT, 5.000000E-02, 2.625000E-02, 1.829000E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, SNAPSIZE, 1.000000E-02 > GRID, RSET_Z_1, 40, 1.400000E+00 > OBJ, NAME, SCRSO > OBJ, POSITION, 0.000000E+00, 3.150000E-02, 0.000000E+00 > OBJ, SIZE, 1.000000E-01, 1.600000E-02, 0.000000E+00 > OBJ, GEOMETRY, poldef > OBJ, ROTATION24, 1 > OBJ, TYPE, USER_DEFINED > OBJ, NAME, SCRSF > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E-01, 3.150000E-02, 0.000000E+00 > OBJ, GEOMETRY, poldef > OBJ, ROTATION24, 1 > OBJ, TYPE, USER_DEFINED > OBJ, NAME, OUTLET > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 5.000000E-01 > OBJ, SIZE, 1.000000E-01, 6.100000E-02, 0.000000E+00 > OBJ, GEOMETRY, polcubet > OBJ, ROTATION24, 1 > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0.000000E+00 > OBJ, TEMPERATURE, SAME > OBJ, COEFFICIENT, 1.000000E+03 > OBJ, TURBULENCE, SAME , SAME > OBJ, NAME, NWALL > OBJ, POSITION, 0.000000E+00, 6.100000E-02, 0.000000E+00 > OBJ, SIZE, 1.000000E-01, 0.000000E+00, 5.000000E-01 > OBJ, GEOMETRY, polcu10 > OBJ, ROTATION24, 1 > OBJ, TYPE, PLATE STOP