************************************************************
   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