TALK=T;RUN( 1, 1)
  DISPLAY

  THERMAL RADIATION IN A COMPARTMENT FIRE


    An IMMERSOL application is presented aimed at the demonstration
    of the method for the simulation of fire in a 3D enclosure when
    the convection, chemical reactions, thermal radiation and
    conjugate heat transfer in the internal structures have to be
    considered simultaneously.

    The demonstration case considered consists of rectangular room
    with an internal wall of finite thickness and door opening at
    the front end of the room. The dimensions of the room are those
    of the Test Compartment used by Swedish National Testing
    Institute for calibration experiments.

    It is supposed that a fuel (methane) leaks from the
    communication system and enters the room vertically upwards
    through an aperture located at the floor in the bottom of the
    rear wall. In the vicinity of fuel leak, there is a vent
    supplying the air for heating and ventilation purposes.

    The fuel is ignited on entry and steady combustion is in
    progress producing the high temperature combustion products.
    Their movement is greatly influenced by buoyancy with the
    salient features as follows:
  
  The task is to calculate the temperatures of the internal wall
  structure and combustion gases along with all related field
  distributions.


  ENDDIS
  VRV USE
  SLICE X
  SLICE OUTLINE ON
  * View and up directions
  VIEW  1.000000E+00  0.000000E+00  0.000000E+00
  UP  0.000000E+00  1.000000E+00  0.000000E+00
  * View centre
  VIEW CENTRE  4.980295E-01  1.129838E+00  1.500320E+00
  * View size
  VIEW SIZE  2.066437E+00
  * View perspective
  VIEW DEPTH  3.000000E+00;VIEW TILT 0.8
  VARIABLE TEMPERATURE
  VECTOR OFF
  CONTOUR ON
  SURFACE OFF
  PAUSE
  VARIABLE H1
  PAUSE
  VARIABLE PROD
  PAUSE
  VARIABLE OXID
  PAUSE
  VARIABLE FUEL
  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 Enthalpy (H1) contours
  con h1  z 3 fi;0.001 ; pause;cl
  msg Mass fraction of fuel  (FUEL) contours
  con fuel  z 3 fi;0.001 ; pause;cl
  msg Mass fraction of product  (PROD) contours
  con prod  z 3 fi;0.001 ; pause;cl
  msg Mass fraction of oxidant  (PROD) contours
  con oxid  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
 ************************************************************
 CPVNAM=VDI;SPPNAM=Core
 IRUNN   =       1 ;LIBREF =      84
 ************************************************************
  Group 1. Run Title
 TEXT(Thermal radiation in a compartment fire )
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,6,12,15)
 ************************************************************
  Group 6. Body-Fitted coordinates

 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(137) =CP1  ; NAME(138) =SCAT
 NAME(139) =EMIS ; NAME(140) =EPKE
 NAME(141) =FUEL ; NAME(142) =TMP1
 NAME(143) =PROD ; NAME(144) =OXID
 NAME(146) =WDIS ; NAME(147) =WGAP
 NAME(148) =DEN1 ; NAME(149) =T3
 NAME(150) =MIXF
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  ,W1  ,H1  ,T3  ,MIXF)
    * Stored variables list
 STORE(DEN1,WGAP,WDIS,OXID,PROD,TMP1,FUEL,EPKE)
 STORE(EMIS,SCAT,CP1,TWAL,ENUT )
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 TURMOD(LVEL)
 
 ************************************************************
  Group 8. Terms & Devices
 TERMS (H1  ,N,Y,Y,Y,Y,Y)
 TERMS (T3  ,N,N,Y,N,Y,N)
 NEWRH1  =    T
 ************************************************************
  Group 9. Properties
 PRESS0  = 1.000000E+05
 RHO1    = GRND6
 RHO1A   = 1.600000E+01 ;RHO1B  = 2.900000E+01
 RHO1C   = 2.800000E+01
 ENUL    = 1.000000E-05
 ENULA   =-4.946790E-06 ;ENULB  = 4.583940E-08
 ENULC   = 8.097400E-11
 CP1     = GRND10
 CP1A    = 1.500000E+03 ;CP1B   = 1.500000E+03
 CP1C    = 1.500000E+03 ;CP1D   = 0.000000E+00
 TMP1    = GRND7
 TMP2A   = 5.482456E-02 ;TMP2B  = 4.900000E+07
 TMP2C   = 0.000000E+00
 DISWAL
 DVO1DT  = 3.330000E-03
 PRNDTL(H1  ) =  7.000000E-01 ;PRNDTL(T3  ) =  7.000000E-01
 PRNDTL(MIXF) =  7.000000E-01
 PRT   (EP  ) =  1.314000E+00
 EMISS   = 8.000000E-01 ;SCATT  = 2.000000E-01
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(TMP1) =  3.000000E+02 ;FIINIT(WDIS) =  1.200000E-01
 FIINIT(T3  ) =  3.000000E+02 ;fiinit(h1)=300*cp1b
 FIINIT(EMIS) =  8.0E-01  ; FIINIT(SCAT) =  2.0E-01
   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 (BUOYANCY,PHASEM,0,0,0,0,0,0,1,1)
 COVAL (BUOYANCY,U1  , FIXFLU      , GRND2       )
 COVAL (BUOYANCY,V1  , FIXFLU      , GRND2       )
 COVAL (BUOYANCY,W1  , FIXFLU      , GRND2       )
 
 PATCH (CHSO    ,VOLUME,0,0,0,0,0,0,1,1)
 
 BUOYA   = 0.000000E+00 ; BUOYB =-9.810000E+00
 BUOYC   = 0.000000E+00
 BUOYD   = 1.200000E+00
 BUOYE   = 2.000000E+01
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =   10000
 RESFAC  = 1.000000E-03
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 9.000000E-01)
 RELAX(U1  ,LINRLX, 1.000000E-01)
 RELAX(V1  ,LINRLX, 1.000000E-01)
 RELAX(W1  ,LINRLX, 1.000000E-01)
 RELAX(KE  ,LINRLX, 0.3)
 RELAX(EP  ,LINRLX, 0.3)
 RELAX(H1  ,LINRLX, 1.000000E-01)
 RELAX(DEN1,LINRLX, 1.000000E-02)
 RELAX(T3  ,LINRLX, 1.000000E-01)
 RELAX(MIXF,LINRLX, 1.000000E-01)
 KELIN   =       3
 ************************************************************
  Group 18. Limits
 VARMIN(P1  ) = -PRESS0+100
 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(KE  ) = 1.000000E+10 ;VARMIN(KE  ) = 0.000000E+00
 VARMAX(EP  ) = 1.000000E+10 ;VARMIN(EP  ) = 0.000000E+00
 VARMAX(H1  ) = 1.000000E+10 ;VARMIN(H1  ) = 1.000000E+05
 VARMAX(FUEL) = 1.000000E+00 ;VARMIN(FUEL) = 0.000000E+00
 VARMAX(TMP1) = 3.000000E+03 ;VARMIN(TMP1) = 3.000000E+02
 VARMAX(PROD) = 1.000000E+00 ;VARMIN(PROD) = 0.000000E+00
 VARMAX(OXID) = 1.000000E+00 ;VARMIN(OXID) = 0.000000E+00
 VARMAX(DEN1) = 2.000000E+00 ;VARMIN(DEN1) = 1.000000E-02
 VARMAX(T3  ) = 3.000000E+03 ;VARMIN(T3  ) = 3.000000E+02
 VARMAX(MIXF) = 1.000000E+00 ;VARMIN(MIXF) = 0.000000E+00
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 GENK    =    T
 SPEDAT(SET,MAXINC,T3,R,20)
 SPEDAT(SET,MACINC,H1,R,20*1500)
 conwiz=t
  ************************************************************
  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   =       4 ;IYMON  =       3 ;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
 NOWIPE  =    T
 
 GVIEW(P,8.301357E-01,3.522742E-01,4.321778E-01)
 GVIEW(UP,-3.124654E-01,9.358968E-01,-1.626729E-01)
 
> DOM,    SIZE,        1.200000E+00, 2.400000E+00, 3.600000E+00
> DOM,    MONIT,       7.000000E-01, 5.000000E-01, 5.200000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   RSET_Z_1,      2, 1.000000E+00
> GRID,   RSET_Z_2,      5, 1.000000E+00
> GRID,   RSET_Z_3,      2, 1.000000E+00
> GRID,   RSET_Z_4,      6, 1.000000E+00
 
> OBJ,    NAME,        DOOR
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 3.600000E+00
> OBJ,    SIZE,        4.000000E-01, 2.000000E+00, 0.000000E+00
> OBJ,    CLIPART,     box
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    TEMPERATURE,   3.000000E+02
> OBJ,    COEFFICIENT,   1.000000E+03
> OBJ,    TURBULENCE,   0.001, 0.009
> OBJ,    T_EXT,         3.000000E+02
 
> OBJ,    NAME,        VENT
> OBJ,    POSITION,    1.000000E+00, 4.000000E-01, 0.000000E+00
> OBJ,    SIZE,        2.000000E-01, 4.000000E-01, 0.000000E+00
> OBJ,    CLIPART,     box
> OBJ,    ROTATION24,       24
> OBJ,    TYPE,        INLET
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    VELOCITY,      0.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ,    TEMPERATURE,   3.000000E+02
> OBJ,    TURB-INTENS,   5.000000E+00
> OBJ,    T_EXT,         3.000000E+02
 
> OBJ,    NAME,        FUELLEAK
> OBJ,    POSITION,    1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.000000E-01, 0.000000E+00, 4.000000E-01
> OBJ,    CLIPART,     box
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        INLET
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    VELOCITY,      0.000000E+00, 3.000000E+00, 0.000000E+00
> OBJ,    TEMPERATURE,   3.000000E+02
> OBJ,    INLET_MIXF,    1.000000E+00
> OBJ,    TURB-INTENS,   5.000000E+00
> OBJ,    T_EXT,         3.000000E+02
 
> OBJ,    NAME,        WALL
> OBJ,    POSITION,    8.000001E-01, 0.000000E+00, 1.600000E+00
> OBJ,    SIZE,        4.000000E-01, 2.400000E+00, 2.000000E-01
> OBJ,    CLIPART,     box
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      102
> OBJ,    EMISSIVITY,    8.000000E-01
 
> OBJ,    NAME,        WALL2
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 1.600000E+00
> OBJ,    SIZE,        4.000000E-01, 2.400000E+00, 2.000000E-01
> OBJ,    CLIPART,     box
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      102
> OBJ,    EMISSIVITY,    8.000000E-01
 
> OBJ,    NAME,        WALL3
> OBJ,    POSITION,    4.000000E-01, 2.000000E+00, 1.600000E+00
> OBJ,    SIZE,        4.000000E-01, 4.000000E-01, 2.000000E-01
> OBJ,    CLIPART,     box
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      102
> OBJ,    EMISSIVITY,    8.000000E-01
 
> OBJ,    NAME,        B7
> OBJ,    POSITION,    1.200000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 2.400000E+00, 3.600000E+00
> OBJ,    CLIPART,     box
> OBJ,    ROTATION24,        1
> OBJ,    VISIBLE,          -1
> OBJ,    TYPE,        PLATE
> OBJ,    EMISSIVITY,    9.000000E-01
 
> OBJ,    NAME,        B8
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 2.400000E+00, 3.600000E+00
> OBJ,    CLIPART,     box
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
> OBJ,    EMISSIVITY,    9.000000E-01
STOP
  inform7begin
(stored var #an1 is anco(h1))
(stored var #ae1 is aeco(h1))
(stored var #ah1 is ahco(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