************************************************************ Q1 created by VDI menu, Version 3.3, Date 27/01/00 CPVNAM=VDI;SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings PHOTON USE P up z use patgeo msg Boiler geometry msg Hit Enter to continue pause cl;red vec x 3 sh use patgeo msg Gas velocities msg Hit Enter to continue pause vec cl;red set vec ph 2 vec x 3 sh use patgeo msg Coal fines velocities msg Hit Enter to continue pause vec cl;red con r2 x 3 fil;.001 use patgeo msg Fines volume fractions msg Hit Enter to continue pause con cl;red con tmp1 x 3 fil;.001 use patgeo msg Gas temperature msg Hit Enter to continue pause con cl;red con tmp2 x 3 fil;.001 use patgeo msg Particle temperatures msg Hit Enter to continue pause con cl;red con yO2 x 3 fil;.001 use patgeo msg Oxygen mass fractions msg Hit Enter to continue pause con cl;red con yN2 x 3 fil;.001 use patgeo msg Nitrogen mass fractions msg Hit Enter to continue pause con cl;red con yH2 x 3 fil;.001 use patgeo msg Hydrogen mass fractions msg Hit Enter to continue pause con cl;red con yH2O x 3 fil;.001 use patgeo msg Water vapour mass fractions msg Hit Enter to continue pause con cl;red con yCO x 3 fil;.001 use patgeo msg Carbon monoxide mass fractions msg Hit Enter to continue pause con cl;red con yCO2 x 3 fil;.001 use patgeo msg Carbon dioxide mass fractions msg Type e and hit Enter to quit ENDUSE DISPLAY A model for the simulation of the flow, heat and mass transfer in the combustion and convection chambers is presented aimed at the simulation of the relevant physical phenomena taking place in a coal-fired utility boilers. ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 302 ************************************************************ Group 1. Run Title TEXT( Coal-Fired Boiler ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,5,10,20) * Set overall domain extent: * xulast yvlast zwlast name XSI= 1.000000E+00; YSI= 2.000000E+00; ZSI= 1.000000E+01 RSET(D,CHAM ) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = F * Non-default variable names NAME( 11) =SHAD ; NAME( 16) =MIXF NAME(132) =PRPS ; NAME(133) =VPOR NAME(134) =HPOR ; NAME(135) =EPOR NAME(136) =NPOR ; NAME(137) =CMDO NAME(138) =TMP2 ; NAME(139) =TMP1 NAME(140) =RHO2 ; NAME(141) =RHO1 NAME(142) =YH2O ; NAME(143) =YH2 NAME(144) =YN2 ; NAME(145) =YCO2 NAME(146) =YO2 ; NAME(147) =YCO NAME(148) =EL1 ; NAME(149) =ENUT NAME(150) =EPKE * Solved variables list SOLVE(P1 ,U1 ,U2 ,V1 ,V2 ,W1 ,W2 ,R1 ) SOLVE(R2 ,SHAD,H1 ,H2 ,MIXF) * Stored variables list STORE(EPKE,ENUT,EL1 ,YCO ,YO2 ,YCO2,YN2 ,YH2 ) STORE(YH2O,RHO1,RHO2,TMP1,TMP2,CMDO,NPOR,EPOR) STORE(HPOR,VPOR,PRPS) * Additional solver options SOLUTN(P1 ,Y,Y,Y,N,N,Y) SOLUTN(R1 ,Y,Y,N,Y,N,Y) SOLUTN(R2 ,Y,Y,N,Y,N,Y) TURMOD(KEMODL) ************************************************************ Group 8. Terms & Devices TERMS (SHAD,Y,Y,Y,Y,N,N) TERMS (H1 ,N,Y,Y,Y,Y,Y) TERMS (H2 ,N,Y,Y,Y,N,Y) NEWRH1 = T ************************************************************ Group 9. Properties RHO1 = GRND7 RHO1A = 9.500000E-01 ;RHO1B = 5.000000E-02 RHO1C = 0.000000E+00 PRESS0 = 1.000000E+05 CP1 = 1.100000E+03 RHO2 = 1.000000E+03 RHO2A = 1.000000E+00 ;RHO2B = 0.000000E+00 RHO2C = 0.000000E+00 CP2 = 1.100000E+03 ENUL = 1.000000E-05 PHINT (H1 ) = GRND7 ;PHINT (H2 ) = GRND7 PHINT (MIXF) = 1.000000E+00 PHNH1A = 0.000000E+00 ;PHNH1B = 0.000000E+00 PHNH1C = 0.000000E+00 PHNH2A = 0.000000E+00 ;PHNH2B = 0.000000E+00 PHNH2C = 0.000000E+00 PRT (EP ) = 1.314000E+00 ************************************************************ Group 10.Inter-Phase Transfer Processes CFIPS = GRND7 RLOLIM = 1.000000E-09 ;CMDOT = GRND3 CFIPA = 1.000000E-09 ;CFIPB = 1.000000E-04 CFIPC = 1.000000E+02 ;CFIPD = 0.000000E+00 CMDTA = 1.400000E+05 ;CMDTB = 0.000000E+00 CMDTC = 1.482112E-01 ;CMDTD = 0.000000E+00 CINT (MIXF) = 0.000000E+00 ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(R1 ) = 9.999000E-01 ;FIINIT(R2 ) = 1.000000E-04 FIINIT(SHAD) = 1.000000E-04 ;FIINIT(VPOR) = 1.000000E+00 FIINIT(HPOR) = 1.000000E+00 ;FIINIT(EPOR) = 1.000000E+00 FIINIT(NPOR) = 1.000000E+00 ;FIINIT(EPKE) = 1.000000E+00 PATCH (STEAMPOR,INIVAL,10,0,0,0,0,0,1,1) INIT(STEAMPOR,VPOR, 0.000000E+00, 5.000000E-01) INIT(STEAMPOR,HPOR, 0.000000E+00, 5.000000E-01) INIT(STEAMPOR,EPOR, 0.000000E+00, 5.000000E-01) INIT(STEAMPOR,NPOR, 0.000000E+00, 5.000000E-01) PATCH (HEATPOR ,INIVAL,12,0,0,0,0,0,1,1) INIT(HEATPOR ,VPOR, 0.000000E+00, 5.000000E-01) INIT(HEATPOR ,HPOR, 0.000000E+00, 5.000000E-01) INIT(HEATPOR ,EPOR, 0.000000E+00, 5.000000E-01) INIT(HEATPOR ,NPOR, 0.000000E+00, 5.000000E-01) INIADD = F ************************************************************ Group 12. Convection and diffusion adjustments No PATCHes used for this Group ************************************************************ Group 13. Boundary & Special Sources PATCH (AGRAV ,PHASEM,0,0,0,0,0,0,1,1) COVAL (AGRAV ,W1 , FIXFLU ,-9.810000E+00) COVAL (AGRAV ,W2 , FIXFLU ,-9.810000E+00) PATCH (STEAM ,VOLUME,10,0,0,0,0,0,1,1) COVAL (STEAM ,H1 , 5.000000E-01, 1.000000E+05) PATCH (AIRHEATE,VOLUME,12,0,0,0,0,0,1,1) COVAL (AIRHEATE,H1 , 4.000000E-01, 1.000000E+05) ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 400 RESFAC = 1.000000E-03 ************************************************************ Group 16. Terminate Iterations LITER (R1 ) = 1 ;LITER (R2 ) = 1 LITER (SHAD) = 1 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX, 5.000000E-01) RELAX(U1 ,FALSDT, 5.000000E-02) RELAX(U2 ,FALSDT, 5.000000E-02) RELAX(V1 ,FALSDT, 5.000000E-02) RELAX(V2 ,FALSDT, 5.000000E-02) RELAX(W1 ,FALSDT, 5.000000E-02) RELAX(W2 ,FALSDT, 5.000000E-02) RELAX(R1 ,LINRLX, 5.000000E-01) RELAX(R2 ,LINRLX, 5.000000E-01) RELAX(SHAD,FALSDT, 5.000000E-02) RELAX(KE ,FALSDT, 5.000000E-02) RELAX(EP ,FALSDT, 5.000000E-02) RELAX(H1 ,FALSDT, 5.000000E-02) RELAX(H2 ,FALSDT, 5.000000E-02) RELAX(MIXF,FALSDT, 5.000000E-02) RELAX(RHO1,LINRLX, 3.000000E-01) KELIN = 3 ************************************************************ Group 18. Limits VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06 VARMAX(U2 ) = 1.000000E+06 ;VARMIN(U2 ) =-1.000000E+06 VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06 VARMAX(V2 ) = 1.000000E+06 ;VARMIN(V2 ) =-1.000000E+06 VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06 VARMAX(W2 ) = 1.000000E+06 ;VARMIN(W2 ) =-1.000000E+06 VARMAX(R1 ) = 1.000000E+00 ;VARMIN(R1 ) = 1.000000E-06 VARMAX(R2 ) = 1.000000E+00 ;VARMIN(R2 ) = 1.000000E-09 VARMAX(H1 ) = 1.000000E+10 ;VARMIN(H1 ) = 0.000000E+00 VARMAX(H2 ) = 1.000000E+10 ;VARMIN(H2 ) = 0.000000E+00 VARMAX(MIXF) = 1.500000E-01 ;VARMIN(MIXF) = 1.000000E-06 VARMAX(TMP2) = 1.000000E+10 ;VARMIN(TMP2) = 3.000000E+02 VARMAX(TMP1) = 1.000000E+10 ;VARMIN(TMP1) = 3.000000E+02 VARMAX(RHO2) = 1.000000E+03 ;VARMIN(RHO2) = 1.000000E-03 VARMAX(RHO1) = 2.000000E+00 ;VARMIN(RHO1) = 1.000000E-03 VARMAX(YH2O) = 1.000000E+00 ;VARMIN(YH2O) = 0.000000E+00 VARMAX(YH2 ) = 1.000000E+00 ;VARMIN(YH2 ) = 0.000000E+00 VARMAX(YN2 ) = 1.000000E+00 ;VARMIN(YN2 ) = 0.000000E+00 VARMAX(YCO2) = 1.000000E+00 ;VARMIN(YCO2) = 0.000000E+00 VARMAX(YO2 ) = 1.000000E+00 ;VARMIN(YO2 ) = 0.000000E+00 VARMAX(YCO ) = 1.000000E+00 ;VARMIN(YCO ) = 0.000000E+00 VARMAX(EL1 ) = 1.000000E+10 ;VARMIN(EL1 ) = 0.000000E+00 VARMAX(ENUT) = 1.000000E+10 ;VARMIN(ENUT) = 0.000000E+00 VARMAX(EPKE) = 1.000000E+10 ;VARMIN(EPKE) = 1.000000E-06 ************************************************************ Group 19. EARTH Calls To GROUND Station USEGRD = T ;USEGRX = T GENK = T ASAP = T ************************************************************ Group 20. Preliminary Printout ECHO = T ************************************************************ Group 21. Print-out of Variables OUTPUT(P1 ,Y,Y,Y,Y,Y,Y) OUTPUT(U1 ,Y,Y,Y,N,Y,Y) OUTPUT(U2 ,Y,Y,Y,N,Y,Y) OUTPUT(V1 ,Y,Y,Y,N,Y,Y) OUTPUT(V2 ,Y,Y,Y,N,Y,Y) OUTPUT(W1 ,Y,Y,Y,Y,Y,Y) OUTPUT(W2 ,Y,Y,Y,Y,Y,Y) OUTPUT(R1 ,Y,Y,Y,Y,Y,Y) OUTPUT(R2 ,Y,Y,Y,Y,Y,Y) OUTPUT(SHAD,Y,N,Y,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 1 ;IZMON = 4 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 NZPRIN = 1 ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts NOWIPE = T GVIEW(P,-9.493308E-01,2.914123E-01,1.176854E-01) GVIEW(UP,1.166439E-01,-2.100885E-02,9.929516E-01) > DOM, SIZE, 1.000000E+00, 2.000000E+00, 1.000000E+01 > DOM, MONIT, 1.000000E-01, 1.000000E-01, 1.750000E+00 > DOM, SCALE, 3.000000E+00, 3.000000E+00, 1.000000E+00 > DOM, SNAPSIZE, 1.000000E-02 > GRID, RSET_Y_1, 10, 1.000000E+00 > DOM, RELAX, 3.000000E-01 > OBJ1, NAME, OUTLET > OBJ1, POSITION, 0.000000E+00, 2.000001E+00, 0.000000E+00 > OBJ1, SIZE, 1.000000E+00, 0.000000E+00, 1.000000E+00 > OBJ1, CLIPART, cube12t > OBJ1, ROTATION24, 1 > OBJ1, GRID, 2.000000E+00 > OBJ1, TYPE, OUTLET > OBJ1, PRESSURE, 0.000000E+00 > OBJ1, TEMPERATURE, -1.026000E+04 > OBJ1, COEFFICIENT, 1.000000E+03 > OBJ1, TEMPERATURE-2,-1.026000E+04 > OBJ1, COEFFICIENT-2, 1.000000E+06 > OBJ1, TURBULENCE, -1.026000E+04,-1.026000E+04 > OBJ2, NAME, INWALL > OBJ2, POSITION, 0.000000E+00, 1.000000E+00, 0.000000E+00 > OBJ2, SIZE, 1.000000E+00, 0.000000E+00, 8.000000E+00 > OBJ2, CLIPART, cube11 > OBJ2, ROTATION24, 1 > OBJ2, GRID, 2.000000E+00 > OBJ2, TYPE, PLATE > OBJ2, POROSITY, 0.000000E+00 > OBJ2, SIDE, BOTH > OBJ3, NAME, SOUWALL > OBJ3, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ3, SIZE, 1.000000E+00, 0.000000E+00, 1.000000E+01 > OBJ3, CLIPART, cube11 > OBJ3, ROTATION24, 1 > OBJ3, TYPE, PLATE > OBJ3, ADIABATI1, 0.000000E+00, 0.000000E+00 > OBJ3, ADIABATI2, 0.000000E+00, 0.000000E+00 > OBJ4, NAME, BURNER1 > OBJ4, POSITION, 6.000000E-01, 0.000000E+00, 5.000000E-01 > OBJ4, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01 > OBJ4, GEOMETRY, cylinder > OBJ4, ROTATION24, 5 > OBJ4, GRID, 2.000000E+00 > OBJ4, TYPE, INLET > OBJ4, DENSITY, 9.999000E-01 > OBJ4, DENSITY-2, 1.000000E-01 > OBJ4, VELOCITY, 0.000000E+00, 2.000000E+01, 0.000000E+00 > OBJ4, VELOCITY-2, 0.000000E+00, 2.000000E+01, 0.000000E+00 > OBJ4, TEMPERATURE, 3.300000E+05 > OBJ4, TEMPERATURE-2, 3.300000E+05 > OBJ4, TURB-INTENS, 5.000000E+00 > OBJ5, NAME, BURNER2 > OBJ5, POSITION, 2.000000E-01, 0.000000E+00, 5.000000E-01 > OBJ5, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01 > OBJ5, GEOMETRY, cylinder > OBJ5, ROTATION24, 5 > OBJ5, GRID, 2.000000E+00 > OBJ5, TYPE, INLET > OBJ5, DENSITY, 9.999000E-01 > OBJ5, DENSITY-2, 1.000000E-01 > OBJ5, VELOCITY, 0.000000E+00, 2.000000E+01, 0.000000E+00 > OBJ5, VELOCITY-2, 0.000000E+00, 2.000000E+01, 0.000000E+00 > OBJ5, TEMPERATURE, 3.300000E+05 > OBJ5, TEMPERATURE-2, 3.300000E+05 > OBJ5, TURB-INTENS, 5.000000E+00 > OBJ6, NAME, TOPWALL > OBJ6, POSITION, 0.000000E+00, 0.000000E+00, 1.000000E+01 > OBJ6, SIZE, 1.000000E+00, 2.000000E+00, 0.000000E+00 > OBJ6, CLIPART, cube11 > OBJ6, ROTATION24, 1 > OBJ6, TYPE, PLATE > OBJ6, ADIABATI1, 0.000000E+00, 0.000000E+00 > OBJ6, ADIABATI2, 0.000000E+00, 0.000000E+00 > OBJ7, NAME, BOTWALL > OBJ7, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ7, SIZE, 1.000000E+00, 2.000000E+00, 0.000000E+00 > OBJ7, CLIPART, cube11 > OBJ7, ROTATION24, 1 > OBJ7, TYPE, PLATE > OBJ7, ADIABATI1, 0.000000E+00, 0.000000E+00 > OBJ7, ADIABATI2, 0.000000E+00, 0.000000E+00 > OBJ8, NAME, NORWALL > OBJ8, POSITION, 0.000000E+00, 2.000000E+00, 1.000000E+00 > OBJ8, SIZE, 1.000000E+00, 0.000000E+00, 9.000000E+00 > OBJ8, CLIPART, cube11 > OBJ8, ROTATION24, 1 > OBJ8, TYPE, PLATE > OBJ8, ADIABATI1, 0.000000E+00, 0.000000E+00 > OBJ8, ADIABATI2, 0.000000E+00, 0.000000E+00 > OBJ9, NAME, EASWALL > OBJ9, POSITION, 1.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ9, SIZE, 0.000000E+00, 2.000000E+00, 1.000000E+01 > OBJ9, CLIPART, cube11 > OBJ9, ROTATION24, 1 > OBJ9, TYPE, PLATE > OBJ9, ADIABATI1, 0.000000E+00, 0.000000E+00 > OBJ9, ADIABATI2, 0.000000E+00, 0.000000E+00 > OBJ10, NAME, B10 > OBJ10, POSITION, 0.000000E+00, 6.000000E-01, 8.600000E+00 > OBJ10, SIZE, 1.000000E+00, 8.666667E-01, 1.400000E+00 > OBJ10, CLIPART, cube12t > OBJ10, ROTATION24, 1 > OBJ10, VISIBLE, -1.000000E+00 > OBJ10, GRID, 2.000000E+00 > OBJ10, TYPE, USER_DEFINED > OBJ11, NAME, B11 > OBJ11, POSITION, 0.000000E+00, 6.000000E-01, 8.600000E+00 > OBJ11, SIZE, 1.000000E+00, 8.670000E-01, 1.400000E+00 > OBJ11, CLIPART, 14 > OBJ11, ROTATION24, 9 > OBJ11, TYPE, CELLTYPE > OBJ12, NAME, B12 > OBJ12, POSITION, 0.000000E+00, 1.200000E+00, 2.400000E+00 > OBJ12, SIZE, 1.000000E+00, 8.000002E-01, 3.600000E+00 > OBJ12, CLIPART, cube12t > OBJ12, ROTATION24, 1 > OBJ12, VISIBLE, -1.000000E+00 > OBJ12, GRID, 2.000000E+00 > OBJ12, TYPE, USER_DEFINED > OBJ13, NAME, B13 > OBJ13, POSITION, 0.000000E+00, 1.200000E+00, 2.400000E+00 > OBJ13, SIZE, 1.000000E+00, 8.000000E-01, 3.600000E+00 > OBJ13, CLIPART, 14 > OBJ13, ROTATION24, 12 > OBJ13, GRID, 2.000000E+00 > OBJ13, TYPE, USER_DEFINED STOP distil=t EX(P1 )= 4.785E+01 EX(U1 )= 4.764E-01 EX(U2 )= 3.964E-01 EX(V1 )= 2.364E+00 EX(V2 )= 2.383E+00 EX(W1 )= 1.278E+01 EX(W2 )= 1.259E+01 EX(R1 )= 1.000E+00 EX(R2 )= 3.466E-06 EX(SHAD)= 2.763E-05 EX(KE )= 2.438E+00 EX(EP )= 1.433E+01 EX(H1 )= 1.717E+06 EX(H2 )= 3.786E+07 EX(MIXF)= 7.561E-02 EX(PRPS)= 1.000E+00 EX(VPOR)= 9.000E-01 EX(HPOR)= 9.000E-01 EX(EPOR)= 9.000E-01 EX(NPOR)= 8.200E-01 EX(CMDO)= 2.983E-04 EX(TMP2)= 1.101E+03 EX(TMP1)= 1.321E+03 EX(RHO2)= 1.000E+03 EX(RHO1)= 2.821E-01 EX(YH2O)= 2.986E-02 EX(YH2 )= 4.626E-04 EX(YN2 )= 7.099E-01 EX(YCO2)= 2.312E-01 EX(YO2 )= 8.085E-03 EX(YCO )= 2.051E-02 EX(EL1 )= 5.266E-02 EX(ENUT)= 4.102E-02 EX(EPKE)= 4.299E+00