TALK=T;RUN( 1, 1)

  DISPLAY

  
    RADIATIVE OVEN

  An IMMERSOL application is presented aimed at the demonstration
  of the method for the   3D situation when the thermal radiation
  is dominant.

  A 3D oven for heat-treating samples is 0.32mx0.32m in cross
  section and 1.5 m long.

  A radiant heater panel consists of a row of cylindrical electrical
  heating elements 1 cm in diameter, 150 cm long, spaced at a 3 cm
  pitch, and backed by a wall convectively cooled from outside.
  The panel has dimensions 30x150 cm and is located 30 cm above a
  well insulated floor, which is also 30x150 cm. The heater elements
  are rated at 5 kW each. The emittances of the elements and back wall
  are 0.9 and 0.8, respectively.

  Two pieces of metal and brick loads are placed on the asbestos sheets.

  The heat transfer coefficient of the ceiling and back/front walls to
  the ambient air of   300K are 20 and 10 W/m^2K, respectively. The heat
  transfer from black side walls is   rated at 300 W/m^2K.

  The task is to calculate the distribution of temperatures and radiant
  heat fluxes.

  ENDDIS
  PHOTON USE
  p; ; ; ; ;
 
  msg Temperature contours (TEM1)
  do kk = 1, 10,3
  con tem1 k kk  fi;0.001
  enddo
  pause;cl
  msg Radiation temperature contours (T3)
  do kk = 1, 10,3
  con t3 k kk  fi;0.001
  enddo
  pause;cl
  msg x-direction heat flux (qrx) contours
  con qrx  z 5 fi;0.001 ; pause;cl
  msg y-direction heat flux (qry) contours
  con qry  z 5 fi;0.001 ; pause;cl
 
  MSG Solution coefficients for Tem1
  con #an1 z 5 fi;0.001 ; pause
  con #ae1 z 5 fi;0.001 ; pause
  con #ap1 z 5 fi;0.001 ; pause
  con #rs1 z 5 fi;0.001 ; pause
  con #cr1 z 5 fi;0.001 ; pause;cl
 
  MSG Solution coefficients for T3
  con #an3 z 5 fi;0.001  ; pause
  con #ae3 z 5 fi;0.001  ; pause
  con #ap3 z 5 fi;0.001  ; pause
  con #rs3 z 5 fi;0.001  ; pause
  con #cr3 z 5 fi;0.001  ; pause

  ENDUSE

 ************************************************************
   Q1 created by VDI menu, Version 3.3, Date 24/07/00
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
 IRUNN   =       1 ;LIBREF =     197
 ************************************************************
  Group 1. Run Title
 TEXT(Radiative oven )
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,32,32,15)
    * Set overall domain extent:
    *        xulast  yvlast  zwlast
            name
 XSI= 3.200000E-01; YSI= 3.200000E-01; ZSI= 1.500000E+00
 RSET(D,CHAM    )
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(144) =QRY  ; NAME(145) =QRX
 NAME(147) =WDIS ; NAME(148) =WGAP
 NAME(149) =T3   ; NAME(150) =TEM1
    * Solved variables list
 SOLVE(T3  ,TEM1)
    * Stored variables list
 STORE(WGAP,WDIS,QRX ,QRY,QRZ )

 ************************************************************
  Group 8. Terms & Devices
 TERMS (T3  ,N,N,Y,N,Y,N)
 ************************************************************
  Group 9. Properties
 SETPRPS(1,  0)
 RHO1    = 1.189000E+00
 PRESS0  = 1.000000E+05
 CP1     = 1.005000E+03
 DISWAL
 ENUL    = 1.544000E-05 ;ENUT   = 0.000000E+00
 DVO1DT  = 3.410000E-03
 PRNDTL(T3  ) = -2.580000E-02 ;PRNDTL(TEM1) = -2.580000E-02
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(WDIS) =  3.200000E-02 ;FIINIT(T3  ) =  5.000000E+02
 FIINIT(TEM1) =  5.000000E+02
   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

 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  = 2000
 RESFAC  = 1.000000E-03
 ************************************************************
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(T3  ,LINRLX, 5.000000E-01)
 RELAX(TEM1,LINRLX, 5.000000E-01)
 ************************************************************
  Group 18. Limits
 VARMAX(T3  ) = 3.000000E+03 ;VARMIN(T3  ) = 3.000000E+02
 VARMAX(TEM1) = 3.000000E+03 ;VARMIN(TEM1) = 3.000000E+02
  ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 ASAP    =    T
 RADIA   = 1.000000E-09 ;RADIB  = 1.000000E-09
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(WDIS,Y,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =      17 ;IYMON  =       1 ;IZMON  =       8
 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.862918E-02,8.898216E-02,9.920822E-01)
 GVIEW(UP,7.917824E-03,9.960333E-01,-8.862919E-02)

> DOM,    SIZE,        3.200000E-01, 3.200000E-01, 1.500000E+00
> DOM,    MONIT,       1.650000E-01, 5.000000E-03, 7.500000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   RSET_X_1,     32, 1.000000E+00
> GRID,   RSET_Y_1,     32, 1.000000E+00
> DOM,    RELAX,       5.000000E-01

> OBJ1,   NAME,        B5
> OBJ1,   POSITION,    0.000000E+00, 3.200000E-01, 0.000000E+00
> OBJ1,   SIZE,        3.200000E-01, 0.000000E+00, 1.500000E+00
> OBJ1,   CLIPART,     cube13
> OBJ1,   ROTATION24,        1
> OBJ1,   VISIBLE,    -1.000000E+00
> OBJ1,   GRID,        2.000000E+00
> OBJ1,   TYPE,        PLATE
> OBJ1,   SURF_HEAT,     0.000000E+00, 0.000000E+00
> OBJ1,   PARA_HEAT,     2.000000E+01, 1.000000E+00, 0.000000E+00
> OBJ1,   T_EXT,         3.000000E+02
> OBJ1,   EMISSIVITY,    8.000000E-01

> OBJ2,   NAME,        B15
> OBJ2,   POSITION,    5.000000E-02, 2.800000E-01, 0.000000E+00
> OBJ2,   SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ2,   CLIPART,     cylinder
> OBJ2,   ROTATION24,        1
> OBJ2,   GRID,        2.000000E+00
> OBJ2,   TYPE,        BLOCKAGE
> OBJ2,   MATERIAL,      100
> OBJ2,   HEAT_FLUX,     0.000000E+00, 5.000000E+03
> OBJ2,   EMISSIVITY,    9.000000E-01

> OBJ3,   NAME,        B16
> OBJ3,   POSITION,    8.000002E-02, 2.800000E-01, 0.000000E+00
> OBJ3,   SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ3,   CLIPART,     cylinder
> OBJ3,   ROTATION24,        1
> OBJ3,   GRID,        2.000000E+00
> OBJ3,   TYPE,        BLOCKAGE
> OBJ3,   MATERIAL,      100
> OBJ3,   HEAT_FLUX,     0.000000E+00, 5.000000E+03
> OBJ3,   EMISSIVITY,    9.000000E-01

> OBJ4,   NAME,        B17
> OBJ4,   POSITION,    1.100000E-01, 2.800000E-01, 0.000000E+00
> OBJ4,   SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ4,   CLIPART,     cylinder
> OBJ4,   ROTATION24,        1
> OBJ4,   GRID,        2.000000E+00
> OBJ4,   TYPE,        BLOCKAGE
> OBJ4,   MATERIAL,      100
> OBJ4,   HEAT_FLUX,     0.000000E+00, 5.000000E+03
> OBJ4,   EMISSIVITY,    9.000000E-01

> OBJ5,   NAME,        B18
> OBJ5,   POSITION,    1.400000E-01, 2.800000E-01, 0.000000E+00
> OBJ5,   SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ5,   CLIPART,     cylinder
> OBJ5,   ROTATION24,        1
> OBJ5,   GRID,        2.000000E+00
> OBJ5,   TYPE,        BLOCKAGE
> OBJ5,   MATERIAL,      100
> OBJ5,   HEAT_FLUX,     0.000000E+00, 5.000000E+03
> OBJ5,   EMISSIVITY,    9.000000E-01

> OBJ6,   NAME,        B19
> OBJ6,   POSITION,    1.700000E-01, 2.800000E-01, 0.000000E+00
> OBJ6,   SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ6,   CLIPART,     cylinder
> OBJ6,   ROTATION24,        1
> OBJ6,   GRID,        2.000000E+00
> OBJ6,   TYPE,        BLOCKAGE
> OBJ6,   MATERIAL,      100
> OBJ6,   HEAT_FLUX,     0.000000E+00, 5.000000E+03
> OBJ6,   EMISSIVITY,    9.000000E-01

> OBJ7,   NAME,        B20
> OBJ7,   POSITION,    2.000000E-01, 2.800000E-01, 0.000000E+00
> OBJ7,   SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ7,   CLIPART,     cylinder
> OBJ7,   ROTATION24,        1
> OBJ7,   GRID,        2.000000E+00
> OBJ7,   TYPE,        BLOCKAGE
> OBJ7,   MATERIAL,      100
> OBJ7,   HEAT_FLUX,     0.000000E+00, 5.000000E+03
> OBJ7,   EMISSIVITY,    9.000000E-01

> OBJ8,   NAME,        B21
> OBJ8,   POSITION,    2.300000E-01, 2.800000E-01, 0.000000E+00
> OBJ8,   SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ8,   CLIPART,     cylinder
> OBJ8,   ROTATION24,        1
> OBJ8,   GRID,        2.000000E+00
> OBJ8,   TYPE,        BLOCKAGE
> OBJ8,   MATERIAL,      100
> OBJ8,   HEAT_FLUX,     0.000000E+00, 5.000000E+03
> OBJ8,   EMISSIVITY,    9.000000E-01

> OBJ9,   NAME,        B22
> OBJ9,   POSITION,    2.600000E-01, 2.800000E-01, 0.000000E+00
> OBJ9,   SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ9,   CLIPART,     cylinder
> OBJ9,   ROTATION24,        1
> OBJ9,   GRID,        2.000000E+00
> OBJ9,   TYPE,        BLOCKAGE
> OBJ9,   MATERIAL,      100
> OBJ9,   HEAT_FLUX,     0.000000E+00, 5.000000E+03
> OBJ9,   EMISSIVITY,    9.000000E-01

> OBJ10,  NAME,        B14
> OBJ10,  POSITION,    2.000000E-02, 2.800000E-01, 0.000000E+00
> OBJ10,  SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ10,  CLIPART,     cylinder
> OBJ10,  ROTATION24,        1
> OBJ10,  GRID,        2.000000E+00
> OBJ10,  TYPE,        BLOCKAGE
> OBJ10,  MATERIAL,      100
> OBJ10,  HEAT_FLUX,     0.000000E+00, 5.000000E+03
> OBJ10,  EMISSIVITY,    9.000000E-01

> OBJ11,  NAME,        B23
> OBJ11,  POSITION,    2.900002E-01, 2.800000E-01, 0.000000E+00
> OBJ11,  SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ11,  CLIPART,     cylinder
> OBJ11,  ROTATION24,        1
> OBJ11,  GRID,        2.000000E+00
> OBJ11,  TYPE,        BLOCKAGE
> OBJ11,  MATERIAL,      100
> OBJ11,  HEAT_FLUX,     0.000000E+00, 5.000000E+03
> OBJ11,  EMISSIVITY,    9.000000E-01

> OBJ12,  NAME,        B24
> OBJ12,  POSITION,    3.200000E-01, 0.000000E+00, 0.000000E+00
> OBJ12,  SIZE,        0.000000E+00, 1.400000E-01, 1.500000E+00
> OBJ12,  CLIPART,     cube13
> OBJ12,  ROTATION24,        1
> OBJ12,  GRID,        2.000000E+00
> OBJ12,  TYPE,        PLATE
> OBJ12,  SURF_HEAT,     0.000000E+00, 0.000000E+00
> OBJ12,  PARA_HEAT,     3.000000E+02, 1.000000E+00, 0.000000E+00
> OBJ12,  T_EXT,         3.000000E+02
> OBJ12,  EMISSIVITY,    1.000000E+00

> OBJ13,  NAME,        B25
> OBJ13,  POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ13,  SIZE,        0.000000E+00, 3.200000E-01, 1.500000E+00
> OBJ13,  CLIPART,     cube13
> OBJ13,  ROTATION24,        1
> OBJ13,  TYPE,        PLATE
> OBJ13,  SURF_HEAT,     0.000000E+00, 0.000000E+00
> OBJ13,  PARA_HEAT,     3.000000E+02, 1.000000E+00, 0.000000E+00
> OBJ13,  T_EXT,         3.000000E+02
> OBJ13,  EMISSIVITY,    1.000000E+00

> OBJ14,  NAME,        B26
> OBJ14,  POSITION,    4.000000E-02, 1.000000E-01, 2.000000E-01
> OBJ14,  SIZE,        1.000000E-01, 3.000000E-02, 5.000000E-01
> OBJ14,  CLIPART,     cylinder
> OBJ14,  ROTATION24,        5
> OBJ14,  GRID,        2.000000E+00
> OBJ14,  TYPE,        BLOCKAGE
> OBJ14,  MATERIAL,      111
> OBJ14,  ADIABATIC,     0.000000E+00, 0.000000E+00
> OBJ14,  EMISSIVITY,    5.000000E-01

> OBJ15,  NAME,        B27
> OBJ15,  POSITION,    1.800000E-01, 1.500000E-01, 8.000000E-01
> OBJ15,  SIZE,        1.000000E-01, 5.000000E-02, 5.000000E-01
> OBJ15,  CLIPART,     cube
> OBJ15,  ROTATION24,        1
> OBJ15,  GRID,        2.000000E+00
> OBJ15,  TYPE,        BLOCKAGE
> OBJ15,  MATERIAL,      102
> OBJ15,  ADIABATIC,     0.000000E+00, 0.000000E+00
> OBJ15,  EMISSIVITY,    1.000000E+00

> OBJ16,  NAME,        B28
> OBJ16,  POSITION,    1.600000E-01, 1.400000E-01, 0.000000E+00
> OBJ16,  SIZE,        1.600000E-01, 1.000000E-02, 1.500000E+00
> OBJ16,  CLIPART,     cube14
> OBJ16,  ROTATION24,        1
> OBJ16,  GRID,        2.000000E+00
> OBJ16,  TYPE,        BLOCKAGE
> OBJ16,  MATERIAL,      113
> OBJ16,  ADIABATIC,     0.000000E+00, 0.000000E+00
> OBJ16,  EMISSIVITY,    5.000000E-01

> OBJ17,  NAME,        B29
> OBJ17,  POSITION,    2.000000E-02, 9.000000E-02, 0.000000E+00
> OBJ17,  SIZE,        1.400000E-01, 1.000000E-02, 1.500000E+00
> OBJ17,  CLIPART,     cube14
> OBJ17,  ROTATION24,        1
> OBJ17,  GRID,        2.000000E+00
> OBJ17,  TYPE,        BLOCKAGE
> OBJ17,  MATERIAL,      113
> OBJ17,  ADIABATIC,     0.000000E+00, 0.000000E+00
> OBJ17,  EMISSIVITY,    5.000000E-01

> OBJ18,  NAME,        B30
> OBJ18,  POSITION,    3.200000E-01, 1.500000E-01, 0.000000E+00
> OBJ18,  SIZE,        0.000000E+00, 1.700000E-01, 1.500000E+00
> OBJ18,  CLIPART,     cube13
> OBJ18,  ROTATION24,        1
> OBJ18,  GRID,        2.000000E+00
> OBJ18,  TYPE,        PLATE
> OBJ18,  SURF_HEAT,     0.000000E+00, 0.000000E+00
> OBJ18,  PARA_HEAT,     3.000000E+02, 1.000000E+00, 0.000000E+00
> OBJ18,  T_EXT,         3.000000E+02
> OBJ18,  EMISSIVITY,    1.000000E+00

> OBJ19,  NAME,        B31
> OBJ19,  POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ19,  SIZE,        3.200000E-01, 0.000000E+00, 1.500000E+00
> OBJ19,  CLIPART,     cube11
> OBJ19,  ROTATION24,        1
> OBJ19,  GRID,        2.000000E+00
> OBJ19,  TYPE,        PLATE
> OBJ19,  ADIABATIC,     0.000000E+00, 0.000000E+00
> OBJ19,  EMISSIVITY,    8.000000E-01

> OBJ20,  NAME,        B32
> OBJ20,  POSITION,    0.000000E+00, 0.000000E+00, 1.500000E+00
> OBJ20,  SIZE,        3.200000E-01, 3.200000E-01, 0.000000E+00
> OBJ20,  CLIPART,     cube13
> OBJ20,  ROTATION24,        1
> OBJ20,  VISIBLE,    -1.000000E+00
> OBJ20,  GRID,        2.000000E+00
> OBJ20,  TYPE,        PLATE
> OBJ20,  SURF_HEAT,     0.000000E+00, 0.000000E+00
> OBJ20,  PARA_HEAT,     1.000000E+01, 1.000000E+00, 0.000000E+00
> OBJ20,  T_EXT,         3.000000E+02
> OBJ20,  EMISSIVITY,    5.000000E-01

> OBJ21,  NAME,        B33
> OBJ21,  POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ21,  SIZE,        3.200000E-01, 3.200000E-01, 0.000000E+00
> OBJ21,  CLIPART,     cube13
> OBJ21,  ROTATION24,        1
> OBJ21,  GRID,        2.000000E+00
> OBJ21,  TYPE,        PLATE
> OBJ21,  SURF_HEAT,     0.000000E+00, 0.000000E+00
> OBJ21,  PARA_HEAT,     1.000000E+01, 1.000000E+00, 0.000000E+00
> OBJ21,  T_EXT,         3.000000E+02
> OBJ21,  EMISSIVITY,    5.000000E-01
STOP
  inform7begin
(stored var #an1 is anco(tem1))
(stored var #ae1 is aeco(tem1))
(stored var #ap1 is apco(tem1))
(stored var #cr1 is corr(tem1))
(stored var #rs1 is resi(tem1))
 
 
(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