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:
- Ventilation air is heated by the fire source and rises due
to buoyancy and subsequently leaves the room through the top of
the door.
- Ambient air enters through the lower portion of the doorway;
brings in an extra oxidant and is heated up by combustion,
rises due to buoyancy and leaves the room through the higher
part of the door.
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