TALK=T;RUN( 1, 1)
  DISPLAY

  A SECONDARY COMBUSTION CHAMBER

  This case models three-dimensional secondary combustion
  chamber of incinerator unit.

  SCRS model of carbon monoxide combustion is extended here
  to handle five gases mixture flow with heat and mass
  transfer resulting from three feeding streams of different
  composition entering the chamber.

  ENDDIS
 ************************************************************
   Q1 created by VDI menu, Version 3.3, Date 27/01/00
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
     PHOTON USE
     p




     use patgeo
     msg Chamber geometry
     msg Press RETURN to continue
     pause
     cl
     con tmp1 x m fil;.001
     con tmp1 y m fil;.001
     con tmp1 z m fil;.001
     vec x m
     vec y m
     vec z m
     msg Outer wall velocity vectors and temperature distribution
     use patgeo
     msg Press RETURN to continue
     pause
     cl
     con MO2 x m fil;.001
     con MO2 y m fil;.001
     con MO2 z m fil;.001
     msg Outer wall oxygen, O2, mass fraction distribution
     use patgeo
     msg Press RETURN to continue
     pause
     cl
     con MN2 x m fil;.001
     con MN2 y m fil;.001
     con MN2 z m fil;.001
     msg Outer wall nytrogen, N2, mass fraction distribution
     use patgeo
     msg Press RETURN to continue
     pause
     cl
     con MCO x m fil;.001
     con MCO y m fil;.001
     con MCO z m fil;.001
     msg Outer wall carbon monoxide, CO, mass fraction distribution
     use patgeo
     msg Press RETURN to continue
     pause
     cl
     con MCO2 x m fil;.001
     con MCO2 y m fil;.001
     con MCO2 z m fil;.001
     msg Outer wall carbon dioxide, CO2, mass fraction distribution
     use patgeo
     msg Press RETURN to continue
     pause
     cl
     con MH2O x m fil;.001
     con MH2O y m fil;.001
     con MH2O z m fil;.001
     msg Outer wall water vapour, H2O, mass fraction distribution
     use patgeo
     msg Press E to end
     enduse
  
 ************************************************************
 IRUNN   =       1 ;LIBREF =     302
 ************************************************************
  Group 1. Run Title
 TEXT( Secondary Combustion Chamber
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,14,26,22)
    * Set overall domain extent:
    *        xulast  yvlast  zwlast
            name
 XSI= 9.200000E-01; YSI= 2.100000E+00; ZSI= 2.400000E+00
 RSET(D,CHAM    )
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(134) =MO2  ; NAME(135) =MN2
 NAME(136) =MH2O ; NAME(137) =MCO2
 NAME(138) =MCO  ; NAME(139) =RHO1
 NAME(140) =MOXI ; NAME(141) =MCOO
 NAME(142) =MFUE ; NAME(143) =TMP1
 NAME(144) =MPRO ; NAME(145) =MDIL
 NAME(146) =EL1  ; NAME(147) =ENUT
 NAME(148) =EPKE ; NAME(149) =PRPS
 NAME(150) =FMIX
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  ,W1  ,H1  ,MN2 ,MCOO,MDIL)
 SOLVE(FMIX)
    * Stored variables list
 STORE(PRPS,EPKE,ENUT,EL1 ,MPRO,TMP1,MFUE,MOXI)
 STORE(RHO1,MCO ,MCO2,MH2O,MO2 )
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 TURMOD(KEMODL)

 ************************************************************
  Group 8. Terms & Devices
 TERMS (H1  ,N,Y,Y,Y,Y,Y)
 NEWRH1  =    T
 ************************************************************
  Group 9. Properties
 RHO1    = GRND
 PRESS0  = 1.000000E+05
 CP1     = 1.004000E+03
 ENUL    = 1.000000E-05
 ENULA   =-4.946790E-06 ;ENULB  = 4.583940E-08
 ENULC   = 8.097400E-11
 DRH1DP  = GRND5
 DVO1DT  = 3.330000E-03
 PRNDTL(H1  ) = -2.630000E-02
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Echo PLANT settings
    PLANTBEGIN
    DEN1=MCO/28.+MO2/32.+MN2/28.+MCO2/44.+MH2O/18.
    DEN1= (PRESS0+P1)/(8314.*DEN1*TMP1+TINY)

   FMIX=AMAX1(0.,AMIN1(1.,FMIX))
   MOXI=AMAX1(0.,(0.803-FMIX)*0.232/0.803)
   MFUE=AMAX1(0.,(FMIX-0.803)*0.232/(0.5174*0.803))
   MPRO=1.-MOXI-MFUE-MDIL
   TMP1=(H1-1.e7*MFUE)/1100.
   MCO2=MCOO+MPRO
   MO2=MOXI
   MCO=MFUE
   MH2O=MDIL-MCOO-MN2
    PLANTEND
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(RHO1) =  1.000000E+00 ;FIINIT(MFUE) =  1.000000E+00
 FIINIT(TMP1) =  1.000000E+03 ;FIINIT(MPRO) =  0.000000E+00
 FIINIT(MDIL) =  0.000000E+00 ;FIINIT(ENUT) =  1.000000E-05
 FIINIT(EPKE) =  1.000000E+00 ;FIINIT(FMIX) =  1.000000E+00
   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
   No PATCHes used for this Group

 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =     500
 RESFAC  = 1.000000E-03
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 3.000000E-01)
 RELAX(U1  ,FALSDT, 2.000000E-02)
 RELAX(V1  ,FALSDT, 2.000000E-02)
 RELAX(W1  ,FALSDT, 2.000000E-02)
 RELAX(KE  ,FALSDT, 2.000000E-02)
 RELAX(EP  ,FALSDT, 2.000000E-02)
 RELAX(H1  ,FALSDT, 2.000000E-02)
 RELAX(MN2 ,FALSDT, 2.000000E-02)
 RELAX(RHO1,LINRLX, 2.000000E-01)
 RELAX(MCOO,FALSDT, 2.000000E-02)
 RELAX(MDIL,FALSDT, 2.000000E-02)
 RELAX(FMIX,FALSDT, 2.000000E-02)
 KELIN   =       3
 ************************************************************
  Group 18. Limits
 VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
 VARMAX(W1  ) = 1.000000E+06 ;VARMIN(W1  ) =-1.000000E+06
 VARMAX(EP  ) = 1.000000E+10 ;VARMIN(EP  ) = 0.000000E+00
 VARMAX(H1  ) = 1.000000E+10 ;VARMIN(H1  ) = 1.000000E+05
 VARMAX(MO2 ) = 1.000000E+00 ;VARMIN(MO2 ) = 0.000000E+00
 VARMAX(MN2 ) = 1.000000E+00 ;VARMIN(MN2 ) = 0.000000E+00
 VARMAX(MH2O) = 1.000000E+00 ;VARMIN(MH2O) = 0.000000E+00
 VARMAX(MCO2) = 1.000000E+00 ;VARMIN(MCO2) = 0.000000E+00
 VARMAX(MCO ) = 1.000000E+00 ;VARMIN(MCO ) = 0.000000E+00
 VARMAX(RHO1) = 5.000000E+00 ;VARMIN(RHO1) = 3.000000E-02
 VARMAX(MOXI) = 1.000000E+00 ;VARMIN(MOXI) = 0.000000E+00
 VARMAX(MCOO) = 1.000000E+00 ;VARMIN(MCOO) = 0.000000E+00
 VARMAX(MFUE) = 1.000000E+00 ;VARMIN(MFUE) = 0.000000E+00
 VARMAX(TMP1) = 1.000000E+04 ;VARMIN(TMP1) = 2.890000E+02
 VARMAX(MPRO) = 1.000000E+00 ;VARMIN(MPRO) = 0.000000E+00
 VARMAX(MDIL) = 1.000000E+00 ;VARMIN(MDIL) = 0.000000E+00
 VARMAX(ENUT) = 1.000000E+03 ;VARMIN(ENUT) = 1.000000E-05
 VARMAX(FMIX) = 1.000000E+00 ;VARMIN(FMIX) = 0.000000E+00
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 NAMSAT  =MOSG
 GENK    =    T
 ASAP    =    T
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =       8 ;IYMON  =      11 ;IZMON  =      17
 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
 NOWIPE  =    T

 GVIEW(P,8.829637E-01,4.309091E-01,1.862594E-01)
 GVIEW(UP,-4.216301E-01,9.023954E-01,-8.894201E-02)

> DOM,    SIZE,        9.200000E-01, 2.100000E+00, 2.400000E+00
> DOM,    MONIT,       4.925000E-01, 7.879170E-01, 1.602500E+00
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> DOM,    RELAX,       1.400000E-01

> OBJ1,   NAME,        PRIMINL
> OBJ1,   POSITION,    8.500000E-02, 0.000000E+00, 4.000000E-02
> OBJ1,   SIZE,        7.500000E-01, 0.000000E+00, 7.500000E-01
> OBJ1,   CLIPART,     cylinder
> OBJ1,   ROTATION24,        5
> OBJ1,   TYPE,        INLET
> OBJ1,   DENSITY,       2.620000E-01
> OBJ1,   VOLUFLOW,      1.420000E+00
> OBJ1,   TEMPERATURE,   2.330000E+06
> OBJ1,   INLET_MN2,     6.000000E-01
> OBJ1,   INLET_MCOO,    2.000000E-02
> OBJ1,   INLET_MDIL,    8.900000E-01
> OBJ1,   INLET_FMIX,    1.000000E+00
> OBJ1,   TURB-INTENS,   5.000000E+00

> OBJ2,   NAME,        SECAIR1
> OBJ2,   POSITION,    1.250000E-01, 3.150000E-01, 0.000000E+00
> OBJ2,   SIZE,        1.300000E-01, 1.300000E-01, 0.000000E+00
> OBJ2,   CLIPART,     cylinder
> OBJ2,   ROTATION24,        1
> OBJ2,   TYPE,        INLET
> OBJ2,   DENSITY,       1.195000E+00
> OBJ2,   VOLUFLOW,      9.330000E-02
> OBJ2,   TEMPERATURE,   3.180000E+05
> OBJ2,   INLET_MN2,     7.680000E-01
> OBJ2,   INLET_MDIL,    7.680000E-01
> OBJ2,   TURB-INTENS,   5.000000E+00

> OBJ3,   NAME,        SECAIR2
> OBJ3,   POSITION,    3.950000E-01, 3.150000E-01, 0.000000E+00
> OBJ3,   SIZE,        1.300000E-01, 1.300000E-01, 0.000000E+00
> OBJ3,   CLIPART,     cylinder
> OBJ3,   ROTATION24,        1
> OBJ3,   TYPE,        INLET
> OBJ3,   DENSITY,       1.195000E+00
> OBJ3,   VOLUFLOW,      9.300000E-02
> OBJ3,   TEMPERATURE,   3.180000E+05
> OBJ3,   INLET_MN2,     7.680000E-01
> OBJ3,   INLET_MDIL,    7.680000E-01
> OBJ3,   TURB-INTENS,   5.000000E+00

> OBJ4,   NAME,        SECAIR3
> OBJ4,   POSITION,    6.650000E-01, 3.150000E-01, 0.000000E+00
> OBJ4,   SIZE,        1.300000E-01, 1.300000E-01, 0.000000E+00
> OBJ4,   CLIPART,     cylinder
> OBJ4,   ROTATION24,        1
> OBJ4,   TYPE,        INLET
> OBJ4,   DENSITY,       1.195000E+00
> OBJ4,   VOLUFLOW,      9.300000E-02
> OBJ4,   TEMPERATURE,   3.180000E+05
> OBJ4,   INLET_MN2,     7.680000E-01
> OBJ4,   INLET_MDIL,    7.680000E-01
> OBJ4,   TURB-INTENS,   5.000000E+00

> OBJ5,   NAME,        BURNPROD
> OBJ5,   POSITION,    0.000000E+00, 2.550000E-01, 2.810000E-01
> OBJ5,   SIZE,        0.000000E+00, 2.500000E-01, 3.381250E-01
> OBJ5,   CLIPART,     cylinder
> OBJ5,   ROTATION24,        9
> OBJ5,   TYPE,        INLET
> OBJ5,   DENSITY,       1.990000E-01
> OBJ5,   VOLUFLOW,      5.659000E-01
> OBJ5,   TEMPERATURE,   1.910000E+06
> OBJ5,   INLET_MN2,     7.500000E-01
> OBJ5,   INLET_MCOO,    1.300000E-01
> OBJ5,   INLET_MDIL,    9.195000E-01
> OBJ5,   INLET_FMIX,    5.175000E-01
> OBJ5,   TURB-INTENS,   5.000000E+00

> OBJ6,   NAME,        OUTLET
> OBJ6,   POSITION,    8.500000E-02, 2.100000E+00, 4.000000E-02
> OBJ6,   SIZE,        7.500000E-01, 0.000000E+00, 7.500000E-01
> OBJ6,   CLIPART,     cylinder
> OBJ6,   ROTATION24,       20
> OBJ6,   TYPE,        OUTLET
> OBJ6,   PRESSURE,      0.000000E+00
> OBJ6,   TEMPERATURE,   2.890000E+02
> OBJ6,   COEFFICIENT,   1.000000E+00
> OBJ6,   TURBULENCE,   -1.026000E+04,-1.026000E+04

> OBJ7,   NAME,        PLATE010
> OBJ7,   POSITION,    0.000000E+00, 1.115000E+00, 0.000000E+00
> OBJ7,   SIZE,        9.200000E-01, 4.000000E-02, 1.530000E+00
> OBJ7,   CLIPART,     cube14
> OBJ7,   ROTATION24,        1
> OBJ7,   TYPE,        BLOCKAGE
> OBJ7,   MATERIAL,      198

> OBJ8,   NAME,        B9
> OBJ8,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ8,   SIZE,        9.200000E-01, 2.100000E+00, 0.000000E+00
> OBJ8,   CLIPART,     default
> OBJ8,   ROTATION24,        1
> OBJ8,   TYPE,        CELLTYPE

> OBJ9,   NAME,        B10
> OBJ9,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ9,   SIZE,        0.000000E+00, 2.100000E+00, 2.400000E+00
> OBJ9,   CLIPART,     default
> OBJ9,   ROTATION24,        1
> OBJ9,   TYPE,        CELLTYPE

> OBJ10,  NAME,        B11
> OBJ10,  POSITION,    9.200000E-01, 0.000000E+00, 0.000000E+00
> OBJ10,  SIZE,        0.000000E+00, 2.100000E+00, 2.400000E+00
> OBJ10,  CLIPART,     default
> OBJ10,  ROTATION24,        1
> OBJ10,  VISIBLE,    -1.000000E+00
> OBJ10,  TYPE,        CELLTYPE

> OBJ11,  NAME,        B12
> OBJ11,  POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ11,  SIZE,        9.200000E-01, 0.000000E+00, 2.400000E+00
> OBJ11,  CLIPART,     default
> OBJ11,  ROTATION24,        1
> OBJ11,  TYPE,        CELLTYPE

> OBJ12,  NAME,        B13
> OBJ12,  POSITION,    0.000000E+00, 2.100000E+00, 0.000000E+00
> OBJ12,  SIZE,        9.200000E-01, 0.000000E+00, 2.400000E+00
> OBJ12,  CLIPART,     default
> OBJ12,  ROTATION24,        1
> OBJ12,  TYPE,        CELLTYPE

> OBJ13,  NAME,        B14
> OBJ13,  POSITION,    0.000000E+00, 0.000000E+00, 2.400000E+00
> OBJ13,  SIZE,        9.200000E-01, 2.100000E+00, 0.000000E+00
> OBJ13,  CLIPART,     default
> OBJ13,  ROTATION24,        1
> OBJ13,  TYPE,        CELLTYPE

> OBJ14,  NAME,        B15
> OBJ14,  POSITION,    0.000000E+00, 9.900000E-01, 1.376000E+00
> OBJ14,  SIZE,        9.200000E-01, 1.250000E-01, 1.540000E-01
> OBJ14,  CLIPART,     cube14
> OBJ14,  ROTATION24,        1
> OBJ14,  TYPE,        BLOCKAGE
> OBJ14,  MATERIAL,      198

> OBJ15,  NAME,        B16
> OBJ15,  POSITION,    0.000000E+00, 9.900000E-01, 0.000000E+00
> OBJ15,  SIZE,        9.200000E-01, 1.250000E-01, 1.280000E-01
> OBJ15,  CLIPART,     cube14
> OBJ15,  ROTATION24,        1
> OBJ15,  TYPE,        BLOCKAGE
> OBJ15,  MATERIAL,      198

> OBJ16,  NAME,        B17
> OBJ16,  POSITION,    0.000000E+00, 9.900000E-01, 7.000000E-01
> OBJ16,  SIZE,        9.200000E-01, 1.250000E-01, 9.000000E-02
> OBJ16,  CLIPART,     cube14
> OBJ16,  ROTATION24,        1
> OBJ16,  TYPE,        BLOCKAGE
> OBJ16,  MATERIAL,      198
M
STOP