DISPLAY Library case 643: Tutorial: Combustion model - ESCRS ENDDIS TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2018, Date 08/01/18 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 * 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(141)=TMP1 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,F) * Stored variables list STORE(ENUT,EL1,DEN1,EPKE,SPH1,MMWT,MCO,MCO2) STORE(TMP1,N2,CO2,H2O,CO,H2,O2,CH4) STORE(H1) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(F,Y,Y,Y,N,N,Y) TURMOD(KECHEN) ************************************************************ Group 8. Terms & Devices NEWRH1 = T DENPCO = T ************************************************************ Group 9. Properties PRESS0 =3.8E+05 RHO1 = GRND9 RHO1A = GRND9 ;RHO1B =0. RHO1C =0. ENUL =4.2E-05 CP1 = GRND7 TMP1 = GRND9 PRT(KE)=0.75 ;PRT(EP)=1.15 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1)=0. ;FIINIT(W1)=5.5 FIINIT(KE)=9.458 ;FIINIT(EP)=4725. FIINIT(F)=0.186 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 This patch is attached to object SCRSO INLET(SCRSO, LOW, -1, 0, 0, 0, 0, 0, 1, 1) VALUE(SCRSO, P1, GRND1) VALUE(SCRSO, V1, 0.) VALUE(SCRSO, W1, 19.6) VALUE(SCRSO, KE, 15.4) VALUE(SCRSO, EP, 6190.) VALUE(SCRSO, CH4, 0.) VALUE(SCRSO, F, 0.) This patch is attached to object SCRSF INLET(SCRSF, LOW, -1, 0, 0, 0, 0, 0, 1, 1) VALUE(SCRSF, P1, GRND1) VALUE(SCRSF, V1, 0.) VALUE(SCRSF, W1, 0.91) VALUE(SCRSF, KE, 0.033) VALUE(SCRSF, EP, 0.31) VALUE(SCRSF, CH4, 1.) VALUE(SCRSF, F, 1.) EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 400 SARAH =1. RESFAC =1.0E-04 ************************************************************ Group 16. Terminate Iterations LITER(KE)=10 ;LITER(EP)=10 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,1. ) RELAX(KE ,LINRLX,0.5 ) RELAX(EP ,LINRLX,0.5 ) RELAX(F ,LINRLX,0.5 ) RELAX(DEN1,LINRLX,1.0E-02 ) KELIN = 3 ************************************************************ Group 18. Limits 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. ************************************************************ Group 19. EARTH Calls To GROUND Station GENK = T PARSOL = F IENUTA = 2 ISG62 = 1 CHSOB =16. 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.) 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,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.) SPEDAT(SET,SCRS,YFU2,R,0.) SPEDAT(SET,SCRS,YFU3,R,0.) SPEDAT(SET,SCRS,YFU4,R,0.) SPEDAT(SET,SCRS,YFU5,R,0.) SPEDAT(SET,SCRS,YFU6,R,0.) SPEDAT(SET,SCRS,YFU7,R,0.) SPEDAT(SET,SCRS,TEMPFU,R,288.) SPEDAT(SET,SCRS,YOX1,R,0.) SPEDAT(SET,SCRS,YOX2,R,0.232) SPEDAT(SET,SCRS,YOX3,R,0.) SPEDAT(SET,SCRS,YOX4,R,0.) SPEDAT(SET,SCRS,YOX5,R,0.) SPEDAT(SET,SCRS,YOX6,R,0.) SPEDAT(SET,SCRS,YOX7,R,0.768) SPEDAT(SET,SCRS,TEMPOX,R,750.) SPEDAT(SET,OUTPUT,NOFIELD,L,T) SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=31.299999 ;EX(V1)=0.5365 EX(W1)=12.21 ;EX(KE)=6.236 EX(EP)=1453. ;EX(H1)=7.251E+05 EX(CH4)=0.1186 ;EX(O2)=0.09841 EX(H2)=0. ;EX(CO)=0. EX(H2O)=0.05644 ;EX(CO2)=0.06894 EX(N2)=0.6576 ;EX(TMP1)=1410. EX(MCO2)=0.04011 ;EX(MCO)=0. EX(MMWT)=26.35 ;EX(SPH1)=1672. EX(F)=0.1437 ;EX(EPKE)=177.5 EX(DEN1)=1.023 ;EX(EL1)=3.651E-03 EX(ENUT)=4.547E-03 ************************************************************ Group 21. Print-out of Variables OUTPUT(CH4 ,Y,N,Y,Y,Y,Y) OUTPUT(O2 ,Y,N,Y,Y,N,N) OUTPUT(H2 ,Y,N,Y,Y,N,N) 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(DEN1,Y,N,Y,N,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 10 ;IZMON = 25 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 YZPR = T ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts GVIEW(P,-1.,0.,0.) GVIEW(UP,0.,1.,0.) GVIEW(VDIS,0.236081) GVIEW(CENTRE,3.044919E-03,0.0305,0.244) > DOM, SIZE, 1.000000E-01, 6.100000E-02, 5.000000E-01 > DOM, MONIT, 5.000000E-02, 2.992500E-02, 2.517498E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02 > GRID, AUTO, T F F > GRID, RSET_X_1, 1, 1.000000E+00 > GRID, RSET_Y_1, 10, 1.000000E+00 > GRID, RSET_Y_2, 10, 1.000000E+00 > GRID, RSET_Y_3, 10, 1.000000E+00 > GRID, RSET_Z_1, 40, 1.400000E+00 > DOM, INI_AMB, YES > DOM, INI_BUOY, YES > OBJ, NAME, SCRSO > OBJ, POSITION, 0.000000E+00, 3.150000E-02, 0.000000E+00 > OBJ, SIZE, TO_END, 1.600000E-02, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, poldef > OBJ, TYPE, USER_DEFINED > OBJ, PATCHES, SCRSO > OBJ, NAME, SCRSF > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, 3.150000E-02, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, poldef > OBJ, TYPE, USER_DEFINED > OBJ, PATCHES, SCRSF > OBJ, NAME, OUTL3 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, AT_END > OBJ, SIZE, TO_END, TO_END, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, polcubet > OBJ, TYPE, OUTLET > OBJ, PRESSURE, P_AMBIENT > OBJ, TEMPERATURE, T_AMBIENT > OBJ, COEFFICIENT, 1000. > OBJ, TURBULENCE, SAME , SAME > OBJ, NAME, NWAL > OBJ, POSITION, 0.000000E+00, AT_END, 0.000000E+00 > OBJ, SIZE, TO_END, 0.000000E+00, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, polcu10 > OBJ, TYPE, PLATE STOP