TALK=T;RUN(1,1) R401

 TALK=T;RUN( 1, 1)
  DISPLAY

  A RADIANT HEATER PANEL

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

  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 well insulated wall. The panel
  has dimensions 30x150 cm and is located 30 cm above a workpiece
  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; the workpiece floor and side walls are
  assumed black at 300K.

  The task is to calculate the operating temperatures of the heater
  elements.

  ENDDIS
  PHOTON USE
  p
 
 
 
  con tem1 z 1 fi;0.001 ; pause
  con t3   z 1 fi;0.001 ; pause
  con qrx  z 1 fi;0.001 ; pause
  con qry  z 1 fi;0.001 ; pause
  
  con #an1 z 1 fi;0.001 ; pause
  con #ae1 z 1 fi;0.001 ; pause
  con #ap1 z 1 fi;0.001 ; pause
  con #rs1 z 1 fi;0.001 ; pause
  con #cr1 z 1 fi;0.001 ; pause
  
  
  con #an3 z 1 fi;0.001  ; pause
  con #ae3 z 1 fi;0.001  ; pause
  con #ap3 z 1 fi;0.001  ; pause
  con #rs3 z 1 fi;0.001  ; pause
  con #cr3 z 1 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(A Radiant Heater Panel                  )
 ************************************************************
  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,1)
    * 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 )
store(rhoin,gtin,wvelin)
    * Additional solver options
 SOLUTN(TEM1,Y,Y,Y,N,N,Y)

 ************************************************************
  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
   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  =    1000
 RESFAC  = 1.000000E-03
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(TEM1,linrlx,0.1)
 RELAX(T3,linrlx,0.1)
 ************************************************************
  Group 18. Limits
VARMIN(T3)=300.   ;VARMAX(T3)=3000.
VARMIN(TEM1)=300. ;VARMAX(TEM1)=3000.
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 ASAP    =    T
 RADIA   = 1.000000E-09 ;RADIB  = 1.000000E-09
 RSG41   = 1.000000E+00
 ************************************************************
  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   =       16 ;IYMON  =      16 ;IZMON  =       1
 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,0.000000E+00,0.000000E+00,-1.000000E+00)
 GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)

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

> OBJ1,   NAME,        B1
> OBJ1,   POSITION,    3.100000E-01, 0.000000E+00, 0.000000E+00
> OBJ1,   SIZE,        1.000000E-02, 3.200000E-01, 1.500000E+00
> OBJ1,   CLIPART,     cube4
> OBJ1,   ROTATION24,        1
> OBJ1,   GRID,        2.000000E+00
> OBJ1,   TYPE,        BLOCKAGE
> OBJ1,   MATERIAL,      100
> OBJ1,   FIXED_TMP,     0.000000E+00, 3.000000E+02
> OBJ1,   EMISSIVITY,    1.000000E+00

> OBJ2,   NAME,        B2
> OBJ2,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ2,   SIZE,        1.000000E-02, 3.200000E-01, 1.500000E+00
> OBJ2,   CLIPART,     cube4
> OBJ2,   ROTATION24,        1
> OBJ2,   GRID,        2.000000E+00
> OBJ2,   TYPE,        BLOCKAGE
> OBJ2,   MATERIAL,      100
> OBJ2,   FIXED_TMP,     0.000000E+00, 3.000000E+02
> OBJ2,   EMISSIVITY,    1.000000E+00

> OBJ3,   NAME,        B3
> OBJ3,   POSITION,    0.000000E+00, 3.100000E-01, 0.000000E+00
> OBJ3,   SIZE,        3.200000E-01, 1.000000E-02, 1.500000E+00
> OBJ3,   CLIPART,     cube14
> OBJ3,   ROTATION24,        1
> OBJ3,   GRID,        2.000000E+00
> OBJ3,   TYPE,        BLOCKAGE
> OBJ3,   MATERIAL,      198

> OBJ4,   NAME,        B4
> OBJ4,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ4,   SIZE,        3.200000E-01, 1.000000E-02, 1.500000E+00
> OBJ4,   CLIPART,     cube4
> OBJ4,   ROTATION24,        1
> OBJ4,   GRID,        2.000000E+00
> OBJ4,   TYPE,        BLOCKAGE
> OBJ4,   MATERIAL,      100
> OBJ4,   FIXED_TMP,     0.000000E+00, 3.000000E+02
> OBJ4,   EMISSIVITY,    1.000000E+00

> OBJ5,   NAME,        B5
> OBJ5,   POSITION,    0.000000E+00, 3.100000E-01, 0.000000E+00
> OBJ5,   SIZE,        3.200000E-01, 0.000000E+00, 1.500000E+00
> OBJ5,   CLIPART,     cube11
> OBJ5,   ROTATION24,        1
> OBJ5,   GRID,        2.000000E+00
> OBJ5,   TYPE,        PLATE
> OBJ5,   POROSITY,      0.000000E+00
> OBJ5,   SIDE,        BOTH
> OBJ5,   ADIABATIC_L,   0.000000E+00, 0.000000E+00
> OBJ5,   EMI_LOW,       8.000000E-01
> OBJ5,   ADIABATIC_H,   0.000000E+00, 0.000000E+00
> OBJ5,   EMI_HIGH,      8.000000E-01

> OBJ6,   NAME,        B15
> OBJ6,   POSITION,    5.000000E-02, 2.800000E-01, 0.000000E+00
> OBJ6,   SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ6,   CLIPART,     cube
> 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,        B16
> OBJ7,   POSITION,    8.000002E-02, 2.800000E-01, 0.000000E+00
> OBJ7,   SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ7,   CLIPART,     cube
> 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,        B17
> OBJ8,   POSITION,    1.100000E-01, 2.800000E-01, 0.000000E+00
> OBJ8,   SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ8,   CLIPART,     cube
> 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,        B18
> OBJ9,   POSITION,    1.400000E-01, 2.800000E-01, 0.000000E+00
> OBJ9,   SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ9,   CLIPART,     cube
> 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,        B19
> OBJ10,  POSITION,    1.700000E-01, 2.800000E-01, 0.000000E+00
> OBJ10,  SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ10,  CLIPART,     cube
> 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,        B20
> OBJ11,  POSITION,    2.000000E-01, 2.800000E-01, 0.000000E+00
> OBJ11,  SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ11,  CLIPART,     cube
> 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,        B21
> OBJ12,  POSITION,    2.300000E-01, 2.800000E-01, 0.000000E+00
> OBJ12,  SIZE,        1.000000E-02, 1.000000E-02, 1.500000E+00
> OBJ12,  CLIPART,     cube
> OBJ12,  ROTATION24,        1
> OBJ12,  GRID,        2.000000E+00
> OBJ12,  TYPE,        BLOCKAGE
> OBJ12,  MATERIAL,      100
> OBJ12,  HEAT_FLUX,     0.000000E+00, 5.000000E+03
> OBJ12,  EMISSIVITY,    9.000000E-01

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

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

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

(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))