TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2018, Date 08/01/18 CPVNAM=VDI; SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY Library case 642: Tutorial: Combustion model - SCRS ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 642 ************************************************************ Group 1. Run Title TEXT(Combustion SCRS ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,15,15,35) * Cylindrical-polar grid CARTES=F ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(141)=PRPS ;NAME(142)=FUEL NAME(143)=TMP1 ;NAME(144)=PROD NAME(145)=OXID ;NAME(146)=MIXF NAME(147)=EPKE ;NAME(148)=DEN1 NAME(149)=EL1 ;NAME(150)=ENUT * Solved variables list SOLVE(P1,U1,V1,W1,H1,FUEL,MIXF) * Stored variables list STORE(ENUT,EL1,DEN1,EPKE,OXID,PROD,TMP1,PRPS) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(H1,Y,Y,Y,N,N,Y) TURMOD(KECHEN) ************************************************************ Group 8. Terms & Devices TERMS (H1 ,N,Y,Y,Y,Y,Y) NEWRH1 = T ************************************************************ Group 9. Properties PRESS0 =8.0E+05 RHO1 = GRND6 RHO1A =16. ;RHO1B =29. RHO1C =28. ENUL =1.544E-05 CP1 = GRND10 CP1A =1500. ;CP1B =1500. CP1C =1500. ;CP1D =0. TMP1 = GRND8 TMP2A =0.054825 ;TMP2B =4.9E+07 TMP2C =0. DVO1DT =3.41E-03 PRNDTL(H1)=-0.0258 PRT(KE)=0.75 ;PRT(EP)=1.15 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(W1)=10. ;FIINIT(KE)=11.25 FIINIT(EP)=353. ;FIINIT(H1)=1.1595E+04 FIINIT(PRPS)=-1. ;FIINIT(FUEL)=0.0548 FIINIT(MIXF)=0.0548 ;FIINIT(EPKE)=0.09 FIINIT(DEN1)=3.6 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 PATCH(CHSO, PHASEM, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(CHSO, FUEL, GRND9, GRND9) EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 500 RESFAC =1.0E-04 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ;LITER(KE)=10 LITER(EP)=10 ;LITER(H1)=10 LITER(FUEL)=10 ;LITER(MIXF)=10 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,0.5 ) RELAX(U1 ,LINRLX,0.5 ) RELAX(V1 ,LINRLX,0.5 ) RELAX(W1 ,LINRLX,0.5 ) RELAX(KE ,LINRLX,0.5 ) RELAX(EP ,LINRLX,0.5 ) RELAX(H1 ,LINRLX,0.5 ) RELAX(FUEL,LINRLX,0.5 ) RELAX(MIXF,LINRLX,0.5 ) RELAX(DEN1,LINRLX,0.5 ) KELIN = 3 ************************************************************ Group 18. Limits VARMAX(P1)=1.0E+05 ;VARMIN(P1)=-7.9999E+05 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=1. VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=273. VARMAX(DEN1)=4. ;VARMIN(DEN1)=-1.0E+10 ************************************************************ Group 19. EARTH Calls To GROUND Station GENK = T PARSOL = F CONWIZ = T IENUTA = 2 ISG62 = 1 CHSOA =0.054825 ;CHSOB =1. SPEDAT(SET,FACETDAT,NUMOBJ,I,8) SPEDAT(SET,RLXFAC,REFVEL,R,100.) SPEDAT(SET,RLXFAC,REFLEN,R,0.2) SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,MAXINC,H1,R,1.0E+05) SPEDAT(SET,MAXINC,MIXF,R,0.1) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=705.799988 ;EX(U1)=1.808 EX(V1)=8.685 ;EX(W1)=14.44 EX(KE)=19.76 ;EX(EP)=2.526E+04 EX(H1)=1.6E+06 ;EX(PRPS)=0.8381 EX(FUEL)=1.653E-03 ;EX(TMP1)=1013. EX(PROD)=0.2036 ;EX(OXID)=0.6329 EX(MIXF)=0.01281 ;EX(EPKE)=659.900024 EX(DEN1)=2.064 ;EX(EL1)=1.008E-03 EX(ENUT)=2.557E-03 ************************************************************ Group 21. Print-out of Variables ************************************************************ Group 22. Monitor Print-Out IXMON = 9 ;IYMON = 6 ;IZMON = 15 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,0.849838,0.03053,0.526159) GVIEW(UP,0.010462,0.997147,-0.074757) GVIEW(VDIS,0.088888) GVIEW(CENTRE,9.516163E-04,0.030113,0.151897) > DOM, SIZE, 6.280000E-01, 6.500000E-02, 2.450000E-01 > DOM, MONIT, 3.390000E-01, 2.383333E-02, 6.945834E-02 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > GRID, MINCELL, 0.025 ,1.0E-02 ,1.0E-02 > GRID, RSET_X_1, 3, 1.000000E+00 > GRID, RSET_X_2, 5, 1.000000E+00 > GRID, RSET_X_3, 2, 1.000000E+00 > GRID, RSET_X_4, 5, 1.000000E+00 > GRID, RSET_Y_1, 3, 1.000000E+00 > GRID, RSET_Y_2, 12, 1.000000E+00 > GRID, RSET_Z_1, 3, 1.000000E+00 > GRID, RSET_Z_2, 3, 1.000000E+00 > GRID, RSET_Z_3, 3, 1.000000E+00 > GRID, RSET_Z_4, 6, 1.000000E+00 > GRID, RSET_Z_5, 3, 1.000000E+00 > GRID, RSET_Z_6, 6, 1.000000E+00 > GRID, RSET_Z_7, 3, 1.000000E+00 > GRID, RSET_Z_8, 8, 1.000000E+00 > OBJ, NAME, CMP1 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, 1.300000E-02, 2.500000E-02 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, polcu8 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 198,Solid with smooth-wall friction > OBJ, NAME, CMP2 > OBJ, POSITION, 0.000000E+00, 1.300000E-02, 0.000000E+00 > OBJ, SIZE, TO_END, TO_END, 2.500000E-02 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, polwe2 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 198,Solid with smooth-wall friction > OBJ, NAME, INLET > OBJ, POSITION, 0.000000E+00, 1.300000E-02, 9.000000E-03 > OBJ, SIZE, TO_END, 0.000000E+00, 4.000000E-03 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, polcu5t > OBJ, TYPE, INLET > OBJ, DENSITY, 3.606 > OBJ, VELOCITY, 0. ,100. ,0. > OBJ, TEMPERATURE, 773. > OBJ, INLET_FUEL, 0.0548 > OBJ, INLET_MIXF, 0.0548 > OBJ, TURB-INTENS, 5. > OBJ, OBJECT-SIDE, HIGH > OBJ, NAME, SIN1 > OBJ, POSITION, 0.000000E+00, AT_END, 7.350000E-02 > OBJ, SIZE, 1.000000E-01, 0.000000E+00, 2.450000E-02 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, polcu5t > OBJ, TYPE, INLET > OBJ, PRESSURE, 0. > OBJ, VELOCITY, 0. ,-40. ,0. > OBJ, TEMPERATURE, 773. > OBJ, TURB-INTENS, 5. > OBJ, NAME, SIN2 > OBJ, POSITION, 3.140000E-01, AT_END, 7.350000E-02 > OBJ, SIZE, 1.000000E-01, 0.000000E+00, 2.450000E-02 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, polcu5t > OBJ, TYPE, INLET > OBJ, PRESSURE, 0. > OBJ, VELOCITY, 0. ,-40. ,0. > OBJ, TEMPERATURE, 773. > OBJ, TURB-INTENS, 5. > OBJ, NAME, DIN > OBJ, POSITION, 0.000000E+00, AT_END, 1.500000E-01 > OBJ, SIZE, 1.000000E-01, 0.000000E+00, 2.450000E-02 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, polcu5t > OBJ, TYPE, INLET > OBJ, PRESSURE, 0. > OBJ, VELOCITY, 0. ,-40. ,0. > OBJ, TEMPERATURE, 773. > OBJ, TURB-INTENS, 5. > OBJ, NAME, OUTLET > 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, 0. > OBJ, TEMPERATURE, SAME > OBJ, COEFFICIENT, 1. > OBJ, TURBULENCE, SAME , SAME > OBJ, NAME, WALL > OBJ, POSITION, 0.000000E+00, AT_END, 2.500000E-02 > OBJ, SIZE, TO_END, 0.000000E+00, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, polcu10 > OBJ, VISIBLE, NO > OBJ, TYPE, PLATE STOP