TALK=T;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2012, Date 19/10/12
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
  Echo DISPLAY / USE settings
  DISPLAY
 
   A direct-fired furnace for heat treating solid materials is a
   rectangular 3D box in shape. The side walls and roof are
   naturally cooled by ambient air with an estimated heat transfer
   coefficient. The have an emmitance of 0.9. The furnace floor is
   refractory and covered by three beams of load made from brick
   material with an emmitance of 0.9.
 
   The preheated methane gas premixed with air enters the furnace
   with an angle to the floor through closely spaced ports at the
   side wall. It ignites and steady kinetically-controled combustion
   takes place producing hot non-scattering furnace gases with an
   absorption coefficient equal to 0.8. The exhaust opening is at
   the opposite side wall.
 
    The task is to calculate the operating temperatures of the
    load and furnace gases along with all related field
    distributions.
 
  ENDDIS
  VRV USE
  SLICE Z
  * View and up directions
  VIEW  0.000000E+00  0.000000E+00 -1.000000E+00
  UP  0.000000E+00  1.000000E+00  0.000000E+00
  * View centre
  VIEW CENTRE  5.000000E-01  1.600000E-01  5.000000E-01
  * View size
  VIEW SIZE  6.524845E-01
  * View perspective
  VIEW DEPTH  3.000000E+00;VIEW TILT 0.8
  VARIABLE TEMPERATURE
  VECTOR OFF
  CONTOUR ON
  SURFACE OFF
  PAUSE
  VARIABLE T3
  PAUSE
  VARIABLE QRY
  PAUSE
  VARIABLE QRX
 
  ENDUSE
  PHOTON USE
  p; ; ; ; ;
 
  msg Temperature contours (Tmp1)
  con tmp1 z 3 fi;0.001 ; pause;cl
  msg Radiation temperature contours (T3)
  con t3   z 3 fi;0.001 ; pause;cl
  msg x-direction heat flux (qrx) contours
  con qrx  z 3 fi;0.001 ; pause;cl
  msg y-direction heat flux (qry) contours
  con qry  z 3 fi;0.001 ; pause;cl
 
  MSG Solution coefficients for H1
  con #an1 z 3 fi;0.001 ; pause
  con #ae1 z 3 fi;0.001 ; pause
  con #ap1 z 3 fi;0.001 ; pause
  con #rs1 z 3 fi;0.001 ; pause
  con #cr1 z 3 fi;0.001 ; pause;cl
 
  MSG Solution coefficients for T3
  con #an3 z 3 fi;0.001  ; pause
  con #ae3 z 3 fi;0.001  ; pause
  con #ap3 z 3 fi;0.001  ; pause
  con #rs3 z 3 fi;0.001  ; pause
  con #cr3 z 3 fi;0.001  ; pause
 
  ENDUSE
 ************************************************************
 IRUNN = 1 ;LIBREF = 0
 ************************************************************
  Group 1. Run Title
 TEXT(Heat treatment in a direct-fired furnace)
 ************************************************************
  Group 2. Transience
 STEADY = T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,20,10,5)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
    * Non-default variable names
 NAME(132)=CP1 ;NAME(133)=SCAT
 NAME(134)=EMIS ;NAME(135)=EPKE
 NAME(136)=ENUT ;NAME(137)=GEN1
 NAME(138)=WDIS ;NAME(139)=WGAP
 NAME(140)=PROD ;NAME(141)=OXID
 NAME(142)=PRPS ;NAME(143)=QRZ
 NAME(144)=QRY ;NAME(145)=QRX
 NAME(146)=DEN1 ;NAME(147)=TMP1
 NAME(148)=FUEL ;NAME(149)=T3
 NAME(150) =MIXF
    * Solved variables list
 SOLVE(P1,U1,V1,W1,H1,FUEL,T3,MIXF)
    * Stored variables list
 STORE(TMP1,DEN1,QRX,QRY,QRZ,PRPS,OXID,PROD)
 STORE(WGAP,WDIS,GEN1,ENUT,EPKE,EMIS,SCAT,CP1)
 TURMOD(KEMODL)
 
 ************************************************************
  Group 8. Terms & Devices
 TERMS (H1  ,N,Y,Y,N,Y,N)
 TERMS (T3  ,N,N,Y,N,Y,N)
 NEWRH1 = T
 ************************************************************
  Group 9. Properties
 PRESS0 =1.0E+05
    * Domain material index is  30 signifying:
    * Ideal gas Law for mixture of 3 gases (ENUL=3.6E-5/GRHO)
 SETPRPS(1, 30)
 TMP1 = GRND8
 TMP1A =1500. ;TMP1B =1500.
 TMP1C =1500.
 TMP2A =0.054825 ;TMP2B =4.9E+07
 TMP2C =0.
 DISWAL
 DRH1DP =3.41E-03
 DVO1DT =3.33E-03
 PRNDTL(H1)=-0.0263 ;PRNDTL(T3)=-0.0263
 PRT(EP)=1.314
 EMISS =0.8 ;SCATT =0.
 RADIA =0.8 ;RADIB =0.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(U1)=2. ;FIINIT(KE)=0.45
 FIINIT(EP)=0.57404 ;FIINIT(H1)=1.1595E+06
 FIINIT(EMIS)=0.8 ;FIINIT(EPKE)=1.
 FIINIT(ENUT)=1.001E-10 ;FIINIT(GEN1)=1.001E-10
 FIINIT(WDIS)=0.032 ;FIINIT(PROD)=1.001E-10
 FIINIT(OXID)=1.001E-10 ;FIINIT(PRPS)=30.
 FIINIT(DEN1)=0.3 ;FIINIT(TMP1)=300.
 FIINIT(FUEL)=0.1 ;FIINIT(MIXF)=0.1
   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 = 1000
 RESFAC =1.0E-06
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1 ,LINRLX,0.9 )
 RELAX(KE ,LINRLX,0.5 )
 RELAX(EP ,LINRLX,0.5 )
 RELAX(H1 ,LINRLX,0.5 )
 RELAX(T3 ,LINRLX,0.5 )
 RELAX(DEN1,LINRLX,0.9 )
 RELAX(MIXF,LINRLX,0.9 )
 KELIN = 3
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-5.0E-06
 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=300.
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 GENK = T
 CONWIZ = T
 RELAX(DEN1,LINRLX,0.1)
 CHSOA =0.054825 ;CHSOB =1.
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(W1  ,Y,Y,Y,N,Y,Y)
 OUTPUT(WDIS,Y,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 10 ;IYMON = 1 ;IZMON = 3
 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.,0.,-1.)
 GVIEW(UP,0.,1.,0.)
 
> DOM,    SIZE,        1.000000E+00, 3.200000E-01, 1.000000E+00
> DOM,    MONIT,       4.750000E-01, 1.600000E-02, 5.000000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID,   RSET_X_1,     20, 1.000000E+00
> GRID,   RSET_Y_1,     10, 1.000000E+00
  > GRID,   RSET_Z_1,      5, 1.000000E+00
 
> OBJ,    NAME,        IMSWL1
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    cube11
> OBJ,    TYPE,        PLATE
> OBJ,    EMISSIVITY,  1.
 
> OBJ,    NAME,        IMSWL2
> OBJ,    POSITION,    0.000000E+00, 3.200000E-01, 0.000000E+00
> OBJ,    SIZE,        1.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    cube13
> OBJ,    TYPE,        PLATE
> OBJ,    SURF_HEAT,   1. ,0.
> OBJ,    PARA_HEAT,   1. ,1. ,0.
> OBJ,    T_EXT,       300.
> OBJ,    EMISSIVITY,  0.9
 
> OBJ,    NAME,        B5
> OBJ,    POSITION,    2.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.000000E-01, 9.600000E-02, 6.000000E-01
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    cube
> OBJ,    GRID,        N,N,N
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    102,BRICK at 20 deg c
> OBJ,    EMISSIVITY,  1.
 
> OBJ,    NAME,        B6
> OBJ,    POSITION,    1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 3.200000E-01, 1.000000E+00
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    cube13
> OBJ,    TYPE,        PLATE
> OBJ,    SURF_HEAT,   0. ,0.
> OBJ,    PARA_HEAT,   1. ,1. ,0.
> OBJ,    T_EXT,       300.
> OBJ,    EMISSIVITY,  0.9
 
> OBJ,    NAME,        B7
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 3.200000E-01, 1.000000E+00
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    cube13
> OBJ,    TYPE,        PLATE
> OBJ,    SURF_HEAT,   0. ,0.
> OBJ,    PARA_HEAT,   1. ,1. ,0.
> OBJ,    T_EXT,       300.
> OBJ,    EMISSIVITY,  0.9
 
> OBJ,    NAME,        OUTLET
> OBJ,    POSITION,    1.000000E+00, 2.400000E-01, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 8.000000E-02, 1.000000E+00
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    cube12t
> OBJ,    GRID,        N,N,N
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,    0.
> OBJ,    TEMPERATURE, 0.
> OBJ,    COEFFICIENT, 1000.
> OBJ,    TURBULENCE,  SAME , SAME
> OBJ,    T_EXT,         3.000000E+02
 
> OBJ,    NAME,        INLET
> OBJ,    POSITION,    0.000000E+00, 1.200000E-01, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 8.000000E-02, 1.000000E+00
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    cube3t
> OBJ,    GRID,        N,N,N
> OBJ,    TYPE,        INLET
> OBJ,    PRESSURE,    3.606
> OBJ,    VELOCITY,    10. ,-10. ,0.
> OBJ,    TEMPERATURE, 700.
> OBJ,    INLET_FUEL,  0.1
> OBJ,    INLET_MIXF,  0.1
> OBJ,    KE_IN,         1.000000E-02
> OBJ,    EP_IN,         9.400417E-02
> OBJ,    T_EXT,         3.000000E+02
 
> OBJ,    NAME,        B8
> OBJ,    POSITION,    5.000001E-01, 0.000000E+00, 4.000000E-01
> OBJ,    SIZE,        1.000000E-01, 9.600000E-02, 6.000000E-01
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        N,N,N
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    102,BRICK at 20 deg c
> OBJ,    EMISSIVITY,  1.
 
> OBJ,    NAME,        B9
> OBJ,    POSITION,    7.999998E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.000000E-01, 9.600000E-02, 6.000000E-01
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        N,N,N
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    102,BRICK at 20 deg c
> OBJ,    EMISSIVITY,  1.
 
> OBJ,    NAME,        B19
> OBJ,    POSITION,    0.000000E+00, 1.200000E-01, 2.800000E-01
> OBJ,    SIZE,        1.999998E-02, 7.200000E-02, 1.000000E-01
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    port_vca
> OBJ,    ROTATION24,  9
> OBJ,    TYPE,        CELLTYPE
 
> OBJ,    NAME,        B11
> OBJ,    POSITION,    0.000000E+00, 1.200000E-01, 0.000000E+00
> OBJ,    SIZE,        1.999998E-02, 7.200000E-02, 1.000000E-01
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    port_vca
> OBJ,    ROTATION24,  9
> OBJ,    TYPE,        CELLTYPE
 
> OBJ,    NAME,        B12
> OBJ,    POSITION,    0.000000E+00, 1.200000E-01, 6.000000E-01
> OBJ,    SIZE,        1.999998E-02, 7.200000E-02, 1.000000E-01
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    port_vca
> OBJ,    ROTATION24,  9
> OBJ,    TYPE,        CELLTYPE
 
> OBJ,    NAME,        B13
> OBJ,    POSITION,    0.000000E+00, 1.200000E-01, 1.400000E-01
> OBJ,    SIZE,        1.999998E-02, 7.200000E-02, 1.000000E-01
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    port_vca
> OBJ,    ROTATION24,  9
> OBJ,    TYPE,        CELLTYPE
 
> OBJ,    NAME,        B14
> OBJ,    POSITION,    0.000000E+00, 1.200000E-01, 7.799999E-01
> OBJ,    SIZE,        1.999998E-02, 7.200000E-02, 1.000000E-01
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    port_vca
> OBJ,    ROTATION24,  9
> OBJ,    TYPE,        CELLTYPE
 
> OBJ,    NAME,        B15
> OBJ,    POSITION,    0.000000E+00, 1.200000E-01, 4.400001E-01
> OBJ,    SIZE,        1.999998E-02, 7.200000E-02, 1.000000E-01
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    port_vca
> OBJ,    ROTATION24,  9
> OBJ,    TYPE,        CELLTYPE
 
> OBJ,    NAME,        B16
> OBJ,    POSITION,    9.880000E-01, 2.400000E-01, 2.000000E-01
> OBJ,    SIZE,        1.200000E-02, 8.000000E-02, 6.000000E-01
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    port_vca
> OBJ,    ROTATION24,  21
> OBJ,    TYPE,        CELLTYPE
 
> OBJ,    NAME,        B17
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ,    SIZE,        1.000000E+00, 3.200000E-01, 0.000000E+00
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    cube13
> OBJ,    TYPE,        PLATE
> OBJ,    SURF_HEAT,   0. ,0.
> OBJ,    PARA_HEAT,   1. ,1. ,0.
> OBJ,    T_EXT,       300.
> OBJ,    EMISSIVITY,  0.9
 
> OBJ,    NAME,        B18
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.000000E+00, 3.200000E-01, 0.000000E+00
> OBJ,    ROT-MODE,    OLD
> OBJ,    GEOMETRY,    cube13
> OBJ,    VISIBLE,     NO
> OBJ,    TYPE,        PLATE
> OBJ,    SURF_HEAT,   0. ,0.
> OBJ,    PARA_HEAT,   1. ,1. ,0.
> OBJ,    T_EXT,       300.
> OBJ,    EMISSIVITY,  0.9
DISTIL=T
STOP
  inform7begin
(stored var #an1 is anco(H1))
(stored var #ae1 is aeco(H1))
(stored var #ap1 is apco(H1))
(stored var #cr1 is corr(H1))
(stored var #rs1 is resi(H1))
 
 
(stored var #an3 is anco(t3))
(stored var #ae3 is aeco(t3))
(stored var #ap3 is apco(t3))
(stored var #cr3 is corr(t3))
(stored var #rs3 is resi(t3))
  inform7end