TALK=T;RUN(1,1)
DISPLAY
This case is based on a standard fire model test case used for
the validation of CFD field models. Its primary purpose is to
test the fire models predictive capability in predicting
temperature and flow distributions in a small compartment
subject to a steady non-spreading fire. This case was studied
experimentally by Steckler et al (NBSIR 82-2520, NBS, 1982).
A steady non-spreading fire of 100kW is located centrally in a
compartment measuring 2.8m square by 2.18m high. The walls and
ceiling were 0.1m thick and covered with a ceramic fibre
insulation board to establish near steady-state conditions within
30 minutes. The door measures 0.74m wide and 1.83m high and is
located centrally in one of the walls. The CFD model assumes all
the walls are adiabatic and perfect radiative reflectors. The
fire is simulated by means of a simple volumetric heat source and
radiative transfer is modelled by means of the P-1 radiation
model.
ix=11 locates mid-plane of doorway
The total sweep number is 800 and the case is run for 600 sweeps,
and therafter GENTRA is activated so as to represent the spray of
a steady sprinkler system located on the ceiling above the fire.
ENDDIS
VRV USE
* Start of frame
* Setting display switches
AXIS ON
POSITION TITLE 2.975420E-01 8.998145E-01
POSITION PROBE 7.671410E-01 9.276438E-02
* Settings for current slice
PROBE 5.040420E-01 1.389167E+00 9.259090E-01; PROBE ON
SLICE Y
SLICE OUTLINE ON
* View and up directions
VIEW -3.537522E-01 -9.002048E-01 2.539546E-01
UP 6.479594E-02 2.472723E-01 9.667771E-01
* View perspective
VIEW DEPTH 3.000000E+00;VIEW TILT 0.8
* Setting object visibility and painting status
OBJECT SHOW TYPE BLOCKAGE
OBJECT HIDE LIST
LIST FIRE
OBJECT SHOW TYPE OPENING
OBJECT HIDE LIST
LIST OPN2
OBJECT PAINT TYPE OPENING OFF
OBJECT SHOW TYPE SPRAY-HEAD
OBJECT PAINT TYPE SPRAY-HEAD OFF
OBJECT HIDE TYPE PLATE
OBJECT PAINT TYPE PLATE OFF
VARIABLE Pressure;VECTOR ON
msg Velocity vectors in Y-plane
PAUSE
* Start of frame
VARIABLE Temperature; VARIABLE RANGE 3.027381E+02 2.564551E+03
VECTOR OFF
CONTOUR ON
CONTOUR OPAQUENESS 100
msg Temperature contours in Y-plane
PAUSE
* Start of frame
VARIABLE Pressure; VECTOR OFF;CONTOUR ON
msg Pressure contours in Y-plane
PAUSE
* Start of frame
VARIABLE Temperature; VARIABLE RANGE 3.027381E+02 2.564551E+03
VECTOR ON
CONTOUR ON
msg Temprature contours and velocity vectors
PAUSE
* Start of frame
VARIABLE Temperature; VARIABLE RANGE 2.938006E+02 2.564551E+03
VECTOR OFF
CONTOUR OFF
* Generating streamlines
STREAM CLEAR
STREAM ANIMATE OFF
HISTORY READ
GHIS
msg Spray tracks
PAUSE
* Start of frame
PROBE 1.314041E+00 2.800000E+00 0.000000E+00; PROBE ON
SLICE X
SLICE OUTLINE ON
* View and up directions
VIEW -2.541715E-01 -9.079202E-01 3.332854E-01
UP 3.135502E-02 3.366860E-01 9.410949E-01
* View centre
VIEW CENTRE 1.905373E+00 1.588076E+00 1.068990E+00
* View size
VIEW SIZE 2.946434E+00
* View perspective
VIEW DEPTH 3.000000E+00;VIEW TILT 0.8
VARIABLE Temperature; VARIABLE RANGE 2.938006E+02 2.564551E+03
VECTOR OFF
CONTOUR ON
HISTORY READ
GHIS
msg Temperature contours and spray tracks in X-plane
PAUSE
* Start of frame
POSITION CONTOURKEY 2.587322E-02 9.276438E-02
POSITION TITLE 2.975420E-01 8.998145E-01
POSITION PROBE 7.671410E-01 9.276438E-02
* Settings for current slice
PROBE 1.140410E-01 1.320001E+00 0.000000E+00; PROBE ON
SLICE Y
SLICE OUTLINE ON
* View and up directions
VIEW -8.648933E-01 -2.411327E-01 4.402422E-01
UP 3.938385E-01 2.177698E-01 8.930104E-01
* View centre
VIEW CENTRE 1.905373E+00 1.588076E+00 1.068990E+00
* View size
VIEW SIZE 2.946434E+00
* View perspective
VIEW DEPTH 3.000000E+00;VIEW TILT 0.8
VARIABLE Temperature; VARIABLE RANGE 3.051951E+02 2.650010E+03
VECTOR OFF
CONTOUR ON
* Generating streamlines
HISTORY READ
GHIS
msg Temperature contours and spray tracks in Y-plane
PAUSE
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Fire-Spray in a compartment )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 1
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 31
XULAST =5.
XFRAC(1)=0.029692 ;XFRAC(2)=0.058729
XFRAC(3)=0.087049 ;XFRAC(4)=0.114568
XFRAC(5)=0.141181 ;XFRAC(6)=0.166745
XFRAC(7)=0.19105 ;XFRAC(8)=0.213761
XFRAC(9)=0.234226 ;XFRAC(10)=0.25
XFRAC(11)=0.268 ;XFRAC(12)=0.281333
XFRAC(13)=0.294667 ;XFRAC(14)=0.308
XFRAC(15)=0.31 ;XFRAC(16)=0.325774
XFRAC(17)=0.346239 ;XFRAC(18)=0.36895
XFRAC(19)=0.393255 ;XFRAC(20)=0.418819
XFRAC(21)=0.445432 ;XFRAC(22)=0.472951
XFRAC(23)=0.50127 ;XFRAC(24)=0.530308
XFRAC(25)=0.56 ;XFRAC(26)=0.58
XFRAC(27)=0.617566 ;XFRAC(28)=0.686253
XFRAC(29)=0.775198 ;XFRAC(30)=0.880527
XFRAC(31)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 20
YVLAST =2.8
YFRAC(1)=0.06131 ;YFRAC(2)=0.122619
YFRAC(3)=0.183929 ;YFRAC(4)=0.245238
YFRAC(5)=0.306548 ;YFRAC(6)=0.367857
YFRAC(7)=0.407143 ;YFRAC(8)=0.446429
YFRAC(9)=0.478571 ;YFRAC(10)=0.514286
YFRAC(11)=0.55 ;YFRAC(12)=0.553571
YFRAC(13)=0.589286 ;YFRAC(14)=0.625
YFRAC(15)=0.6875 ;YFRAC(16)=0.75
YFRAC(17)=0.8125 ;YFRAC(18)=0.875
YFRAC(19)=0.9375 ;YFRAC(20)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 17
ZWLAST =2.2
ZFRAC(1)=0.045455 ;ZFRAC(2)=0.090909
ZFRAC(3)=0.136364 ;ZFRAC(4)=0.199587
ZFRAC(5)=0.26281 ;ZFRAC(6)=0.326033
ZFRAC(7)=0.389256 ;ZFRAC(8)=0.452479
ZFRAC(9)=0.515702 ;ZFRAC(10)=0.578926
ZFRAC(11)=0.642149 ;ZFRAC(12)=0.705372
ZFRAC(13)=0.768595 ;ZFRAC(14)=0.831818
ZFRAC(15)=0.893182 ;ZFRAC(16)=0.954545
ZFRAC(17)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(7)=W1
NAME(12)=KE ;NAME(13)=EP
NAME(127)=MH2O ;NAME(128)=PRPS
NAME(129)=QRZ ;NAME(130)=ENUL
NAME(131)=QRY ;NAME(132)=TDEC
NAME(133)=SPH1 ;NAME(134)=QRX
NAME(135)=LTLS ;NAME(136)=WDIS
NAME(137)=SCAT ;NAME(138)=EMIS
NAME(139)=WGAP ;NAME(140)=T3
NAME(141)=MASS ;NAME(142)=HEAT
NAME(143)=MOMZ ;NAME(144)=MOMY
NAME(145)=MOMX ;NAME(146)=TEM1
NAME(147)=EPKE ;NAME(148)=DEN1
NAME(149)=EL1 ;NAME(150)=ENUT
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(P1,Y,Y,Y,N,N,Y)
SOLUTN(U1,Y,Y,N,N,N,Y)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(W1,Y,Y,N,N,N,Y)
SOLUTN(KE,Y,Y,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
SOLUTN(MH2O,Y,Y,Y,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,N)
SOLUTN(QRZ,Y,N,N,N,N,Y)
SOLUTN(ENUL,Y,N,N,N,N,Y)
SOLUTN(QRY,Y,N,N,N,N,Y)
SOLUTN(TDEC,Y,N,N,N,N,Y)
SOLUTN(SPH1,Y,N,N,N,N,Y)
SOLUTN(QRX,Y,N,N,N,N,Y)
SOLUTN(LTLS,Y,Y,Y,N,N,Y)
SOLUTN(WDIS,Y,N,N,N,N,Y)
SOLUTN(SCAT,Y,N,N,N,N,Y)
SOLUTN(EMIS,Y,N,N,N,N,Y)
SOLUTN(WGAP,Y,N,N,N,N,Y)
SOLUTN(T3,Y,Y,N,N,N,Y)
SOLUTN(MASS,Y,N,N,N,N,Y)
SOLUTN(HEAT,Y,N,N,N,N,Y)
SOLUTN(MOMZ,Y,N,N,N,N,Y)
SOLUTN(MOMY,Y,N,N,N,N,Y)
SOLUTN(MOMX,Y,N,N,N,N,Y)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
SOLUTN(EPKE,Y,N,N,N,N,Y)
SOLUTN(DEN1,Y,N,N,N,N,Y)
SOLUTN(EL1,Y,N,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
DEN1 = 148
VISL = 130
VIST = 150
LEN1 = 149
PRPS = 128
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
TERMS(P1,Y,Y,Y,N,Y,Y)
TERMS(U1,Y,Y,Y,Y,Y,Y)
TERMS(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(KE,N,Y,Y,Y,Y,N)
TERMS(EP,N,Y,Y,Y,Y,N)
TERMS(MH2O,N,Y,Y,Y,N,Y)
TERMS(LTLS,N,N,Y,N,Y,Y)
TERMS(T3,N,N,Y,N,Y,N)
TERMS(TEM1,N,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWRH1 = T
NEWENT = T
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 = GRND5 ;TMP1 =0. ;EL1 = GRND4
TSURR =0. ;TEMP0 =0. ;PRESS0 =1.0E+05
DVO1DT =3.33E-03 ;DRH1DP = GRND5
RHO1A =287.591705 ;RHO1B =461.744507 ;RHO1C =0.71429
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0. ;EL1B =0. ;EL1C =0.
ENUL =1.0E-05 ;ENUT = GRND3
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(W1)=1. ;PRNDTL(KE)=1.
PRNDTL(EP)=1. ;PRNDTL(MH2O)=1.
PRNDTL(LTLS)=1. ;PRNDTL(T3)=-0.0263
PRNDTL(TEM1)=-0.0263
PRT(U1)=1. ;PRT(V1)=1.
PRT(W1)=1. ;PRT(KE)=1.
PRT(EP)=1.314 ;PRT(MH2O)=1.
PRT(LTLS)=1. ;PRT(T3)=1.
PRT(TEM1)=1.
CP1 =1004. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1. ;FIINIT(U1)=1.0E-02
FIINIT(V1)=1.0E-02 ;FIINIT(W1)=1.0E-02
FIINIT(KE)=1.0E-10 ;FIINIT(EP)=1.0E-10
FIINIT(MH2O)=0. ;FIINIT(PRPS)=-1.
FIINIT(QRZ)=1.0E-10 ;FIINIT(ENUL)=1.0E-10
FIINIT(QRY)=1.0E-10 ;FIINIT(TDEC)=1.0E-10
FIINIT(SPH1)=1.0E-10 ;FIINIT(QRX)=1.0E-10
FIINIT(LTLS)=1.0E-10 ;FIINIT(WDIS)=0.22
FIINIT(SCAT)=0. ;FIINIT(EMIS)=0.
FIINIT(WGAP)=1.0E-10 ;FIINIT(T3)=298.
FIINIT(MASS)=1.0E-10 ;FIINIT(HEAT)=1.0E-10
FIINIT(MOMZ)=1.0E-10 ;FIINIT(MOMY)=1.0E-10
FIINIT(MOMX)=1.0E-10 ;FIINIT(TEM1)=293.
FIINIT(EPKE)=1.0E-10 ;FIINIT(DEN1)=1.0E-10
FIINIT(EL1)=1.0E-10 ;FIINIT(ENUT)=1.0E-10
PATCH(OB1 ,INIVAL, 26, 26, 1, 6, 1, 17, 1, 1)
INIT(OB1 ,PRPS,0. ,199. )
PATCH(OB2 ,INIVAL, 26, 26, 15, 20, 1, 17, 1, 1)
INIT(OB2 ,PRPS,0. ,199. )
PATCH(OB3 ,INIVAL, 26, 26, 7, 14, 15, 17, 1, 1)
INIT(OB3 ,PRPS,0. ,199. )
PATCH(OBC ,INIVAL, 11, 15, 9, 12, 1, 3, 1, 1)
INIT(OBC ,P1 ,0. ,0. )
INIT(OBC ,MH2O,0. ,0. )
INIT(OBC ,PRPS,0. ,-1. )
INIADD = F
FSWEEP = 1
NAMFI =stdy
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources
PATCH(BUOYANCY,PHASEM, 1, 31, 1, 20, 1, 17, 1, 1)
COVAL(BUOYANCY,U1 , FIXFLU , GRND2 )
COVAL(BUOYANCY,V1 , FIXFLU , GRND2 )
COVAL(BUOYANCY,W1 , FIXFLU , GRND2 )
PATCH(GENPAT ,CELL , 1, 31, 1, 20, 1, 17, 1, 1)
COVAL(GENPAT ,U1 , FIXFLU , GRND )
COVAL(GENPAT ,V1 , FIXFLU , GRND )
COVAL(GENPAT ,W1 , FIXFLU , GRND )
COVAL(GENPAT ,MH2O, FIXFLU , GRND )
COVAL(GENPAT ,TEM1, FIXFLU , GRND )
PATCH(GENMAS ,CELL , 1, 31, 1, 20, 1, 17, 1, 1)
COVAL(GENMAS ,P1 , FIXFLU , GRND )
PATCH(KESOURCE,PHASEM, 1, 31, 1, 20, 1, 17, 1, 1)
COVAL(KESOURCE,KE , GRND4 , GRND4 )
COVAL(KESOURCE,EP , GRND4 , GRND4 )
PATCH(GX4 ,NORTH , 27, 31, 20, 20, 1, 17, 1, 1)
COVAL(GX4 ,P1 ,1. ,0. )
COVAL(GX4 ,U1 ,0. ,0. )
COVAL(GX4 ,V1 ,0. ,0. )
COVAL(GX4 ,W1 ,0. ,0. )
COVAL(GX4 ,KE ,0. ,0. )
COVAL(GX4 ,EP ,0. ,0. )
COVAL(GX4 ,MH2O,0. ,0. )
COVAL(GX4 ,TEM1,0. ,303. )
PATCH(*RAD4 ,NORTH , 27, 31, 20, 20, 1, 17, 1, 1)
COVAL(*RAD4 ,T3 ,5.6697E-08 ,300. )
PATCH(GX5 ,SOUTH , 27, 31, 1, 1, 1, 17, 1, 1)
COVAL(GX5 ,P1 ,1. ,0. )
COVAL(GX5 ,U1 ,0. ,0. )
COVAL(GX5 ,V1 ,0. ,0. )
COVAL(GX5 ,W1 ,0. ,0. )
COVAL(GX5 ,KE ,0. ,0. )
COVAL(GX5 ,EP ,0. ,0. )
COVAL(GX5 ,MH2O,0. ,0. )
COVAL(GX5 ,TEM1,0. ,303. )
PATCH(*RAD5 ,SOUTH , 27, 31, 1, 1, 1, 17, 1, 1)
COVAL(*RAD5 ,T3 ,5.6697E-08 ,300. )
PATCH(GX6 ,EAST , 31, 31, 1, 20, 1, 17, 1, 1)
COVAL(GX6 ,P1 ,1. ,0. )
COVAL(GX6 ,U1 ,0. ,0. )
COVAL(GX6 ,V1 ,0. ,0. )
COVAL(GX6 ,W1 ,0. ,0. )
COVAL(GX6 ,KE ,0. ,0. )
COVAL(GX6 ,EP ,0. ,0. )
COVAL(GX6 ,MH2O,0. ,0. )
COVAL(GX6 ,TEM1,0. ,303. )
PATCH(*RAD6 ,EAST , 31, 31, 1, 20, 1, 17, 1, 1)
COVAL(*RAD6 ,T3 ,5.6697E-08 ,300. )
PATCH(GX7 ,HIGH , 27, 31, 1, 20, 17, 17, 1, 1)
COVAL(GX7 ,P1 ,1. ,0. )
COVAL(GX7 ,U1 ,0. ,0. )
COVAL(GX7 ,V1 ,0. ,0. )
COVAL(GX7 ,W1 ,0. ,0. )
COVAL(GX7 ,KE ,0. ,0. )
COVAL(GX7 ,EP ,0. ,0. )
COVAL(GX7 ,MH2O,0. ,0. )
COVAL(GX7 ,TEM1,0. ,303. )
PATCH(*RAD7 ,HIGH , 27, 31, 1, 20, 17, 17, 1, 1)
COVAL(*RAD7 ,T3 ,5.6697E-08 ,300. )
PATCH(IM8 ,NWALL , 1, 25, 20, 20, 1, 17, 1, 1)
COVAL(IM8 ,U1 , GRND2 ,0. )
COVAL(IM8 ,W1 , GRND2 ,0. )
COVAL(IM8 ,KE , GRND2 , GRND2 )
COVAL(IM8 ,EP , GRND2 , GRND2 )
COVAL(IM8 ,T3 , GRND4 ,300. )
COVAL(IM8 ,TEM1, GRND2 ,300. )
PATCH(IM9 ,SWALL , 1, 25, 1, 1, 1, 17, 1, 1)
COVAL(IM9 ,U1 , GRND2 ,0. )
COVAL(IM9 ,W1 , GRND2 ,0. )
COVAL(IM9 ,KE , GRND2 , GRND2 )
COVAL(IM9 ,EP , GRND2 , GRND2 )
COVAL(IM9 ,T3 , GRND4 ,300. )
COVAL(IM9 ,TEM1, GRND2 ,300. )
PATCH(IMA ,WWALL , 1, 1, 1, 20, 1, 17, 1, 1)
COVAL(IMA ,V1 , GRND2 ,0. )
COVAL(IMA ,W1 , GRND2 ,0. )
COVAL(IMA ,KE , GRND2 , GRND2 )
COVAL(IMA ,EP , GRND2 , GRND2 )
COVAL(IMA ,T3 , GRND4 ,300. )
COVAL(IMA ,TEM1, GRND2 ,300. )
PATCH(IMB ,LWALL , 1, 31, 1, 20, 1, 1, 1, 1)
COVAL(IMB ,U1 , GRND2 ,0. )
COVAL(IMB ,V1 , GRND2 ,0. )
COVAL(IMB ,KE , GRND2 , GRND2 )
COVAL(IMB ,EP , GRND2 , GRND2 )
COVAL(IMB ,T3 , GRND4 ,300. )
COVAL(IMB ,TEM1, GRND2 ,300. )
PATCH(OCC ,VOLUME, 11, 15, 9, 12, 1, 3, 1, 1)
COVAL(OCC ,TEM1,In-Form:source - see Grp 19)
XCYCLE = F
EGWF = T
WALLCO = GRND2
BUOYA =0. ; BUOYB =0. ; BUOYC =-9.81
BUOYD =1.2763
BUOYE =27.
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 800 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
LITER(P1)=50 ;LITER(U1)=10
LITER(V1)=10 ;LITER(W1)=10
LITER(KE)=10 ;LITER(EP)=10
LITER(MH2O)=20 ;LITER(LTLS)=20
LITER(T3)=20 ;LITER(TEM1)=50
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03
ENDIT(MH2O)=1.0E-03 ;ENDIT(LTLS)=1.0E-03
ENDIT(T3)=1.0E-03 ;ENDIT(TEM1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(U1,FALSDT,1.)
RELAX(V1,FALSDT,1.)
RELAX(W1,FALSDT,1.)
RELAX(KE,LINRLX,0.5)
RELAX(EP,LINRLX,0.5)
RELAX(MH2O,LINRLX,0.3)
RELAX(PRPS,LINRLX,1.)
RELAX(QRZ,LINRLX,1.)
RELAX(ENUL,LINRLX,1.)
RELAX(QRY,LINRLX,1.)
RELAX(TDEC,LINRLX,0.3)
RELAX(SPH1,LINRLX,1.)
RELAX(QRX,LINRLX,1.)
RELAX(LTLS,FALSDT,1.0E+09)
RELAX(WDIS,LINRLX,1.)
RELAX(SCAT,LINRLX,1.)
RELAX(EMIS,LINRLX,1.)
RELAX(WGAP,LINRLX,1.)
RELAX(T3,FALSDT,1.0E+09)
RELAX(MASS,LINRLX,0.3)
RELAX(HEAT,LINRLX,0.3)
RELAX(MOMZ,LINRLX,0.3)
RELAX(MOMY,LINRLX,0.3)
RELAX(MOMX,LINRLX,0.3)
RELAX(TEM1,FALSDT,1.0E+09)
RELAX(EPKE,LINRLX,1.)
RELAX(DEN1,LINRLX,0.3)
RELAX(EL1,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
KELIN = 3
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-7.5E+04
VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
VARMAX(MH2O)=1. ;VARMIN(MH2O)=0.
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
VARMAX(QRZ)=1.0E+10 ;VARMIN(QRZ)=-1.0E+10
VARMAX(ENUL)=1.0E+10 ;VARMIN(ENUL)=-1.0E+10
VARMAX(QRY)=1.0E+10 ;VARMIN(QRY)=-1.0E+10
VARMAX(TDEC)=1.0E+10 ;VARMIN(TDEC)=-1.0E+10
VARMAX(SPH1)=1.0E+10 ;VARMIN(SPH1)=-1.0E+10
VARMAX(QRX)=1.0E+10 ;VARMIN(QRX)=-1.0E+10
VARMAX(LTLS)=1.0E+10 ;VARMIN(LTLS)=-1.0E+10
VARMAX(WDIS)=1.0E+10 ;VARMIN(WDIS)=-1.0E+10
VARMAX(SCAT)=1.0E+10 ;VARMIN(SCAT)=-1.0E+10
VARMAX(EMIS)=1.0E+10 ;VARMIN(EMIS)=-1.0E+10
VARMAX(WGAP)=1.0E+10 ;VARMIN(WGAP)=-1.0E+10
VARMAX(T3)=1.0E+10 ;VARMIN(T3)=-30.
VARMAX(MASS)=1.0E+10 ;VARMIN(MASS)=-1.0E+10
VARMAX(HEAT)=1.0E+10 ;VARMIN(HEAT)=-1.0E+10
VARMAX(MOMZ)=1.0E+10 ;VARMIN(MOMZ)=-1.0E+10
VARMAX(MOMY)=1.0E+10 ;VARMIN(MOMY)=-1.0E+10
VARMAX(MOMX)=1.0E+10 ;VARMIN(MOMX)=-1.0E+10
VARMAX(TEM1)=3000. ;VARMIN(TEM1)=-204.75
VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
VARMAX(EL1)=1.0E+10 ;VARMIN(EL1)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
NAMGRD =FLAR
GENK = T
PARSOL = F
CONWIZ = T
ISG62 = 1
SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,2)
SPEDAT(SET,FACETDAT,NUMOBJ,I,13)
SPEDAT(SET,GENTRA,RESTART,C,NONE)
SPEDAT(SET,MAXINC,T3,R,100.)
SPEDAT(SET,MAXINC,TEM1,R,100.)
SPEDAT(SET,STORED,TDEC,C,=TEM1-273.!ZSLFIN)
SPEDAT(SET,PROPERTY,RHO1,C,=(P1+1.0E+05)/((175.2*MH2O+286.7)*TEM1$)
SPEDAT(SET,PROPERTY,RHO1,C,)!IMAT<100)
SPEDAT(SET,SOURCE,TEM1!FIRE,C,=2.E6*EXP(-10.*MH2O*TIM)!VOLU)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,OBJNAM,!OB1,C,BLK1)
SPEDAT(SET,OBJTYP,!OB1,C,BLOCKAGE)
SPEDAT(SET,BLK1,MATERIAL,R,199.)
SPEDAT(SET,OBJNAM,!OB2,C,BLK2)
SPEDAT(SET,OBJTYP,!OB2,C,BLOCKAGE)
SPEDAT(SET,BLK2,MATERIAL,R,199.)
SPEDAT(SET,OBJNAM,!OB3,C,BLK3)
SPEDAT(SET,OBJTYP,!OB3,C,BLOCKAGE)
SPEDAT(SET,BLK3,MATERIAL,R,199.)
SPEDAT(SET,OBJNAM,!GX4,C,OPN1)
SPEDAT(SET,OBJTYP,!GX4,C,OPENING)
SPEDAT(SET,OBJNAM,!*RAD4,C,OPN1)
SPEDAT(SET,OBJTYP,!*RAD4,C,OPENING)
SPEDAT(SET,ARATIO,!GX4,R,1.)
SPEDAT(SET,OBJNAM,!GX5,C,OPN2)
SPEDAT(SET,OBJTYP,!GX5,C,OPENING)
SPEDAT(SET,OBJNAM,!*RAD5,C,OPN2)
SPEDAT(SET,OBJTYP,!*RAD5,C,OPENING)
SPEDAT(SET,ARATIO,!GX5,R,1.)
SPEDAT(SET,OBJNAM,!GX6,C,OPN3)
SPEDAT(SET,OBJTYP,!GX6,C,OPENING)
SPEDAT(SET,OBJNAM,!*RAD6,C,OPN3)
SPEDAT(SET,OBJTYP,!*RAD6,C,OPENING)
SPEDAT(SET,ARATIO,!GX6,R,1.)
SPEDAT(SET,OBJNAM,!GX7,C,OPN4)
SPEDAT(SET,OBJTYP,!GX7,C,OPENING)
SPEDAT(SET,OBJNAM,!*RAD7,C,OPN4)
SPEDAT(SET,OBJTYP,!*RAD7,C,OPENING)
SPEDAT(SET,ARATIO,!GX7,R,1.)
SPEDAT(SET,OBJNAM,!OB8,C,WALL1)
SPEDAT(SET,OBJTYP,!OB8,C,PLATE)
SPEDAT(SET,EMISSIVITY,OF!IM8,R,1.)
SPEDAT(SET,OBJNAM,!IM8,C,WALL1)
SPEDAT(SET,OBJTYP,!IM8,C,PLATE)
SPEDAT(SET,OBJNAM,!OB9,C,WALL2)
SPEDAT(SET,OBJTYP,!OB9,C,PLATE)
SPEDAT(SET,EMISSIVITY,OF!IM9,R,1.)
SPEDAT(SET,OBJNAM,!IM9,C,WALL2)
SPEDAT(SET,OBJTYP,!IM9,C,PLATE)
SPEDAT(SET,OBJNAM,!OBA,C,WALL3)
SPEDAT(SET,OBJTYP,!OBA,C,PLATE)
SPEDAT(SET,EMISSIVITY,OF!IMA,R,1.)
SPEDAT(SET,OBJNAM,!IMA,C,WALL3)
SPEDAT(SET,OBJTYP,!IMA,C,PLATE)
SPEDAT(SET,OBJNAM,!OBB,C,FLOOR)
SPEDAT(SET,OBJTYP,!OBB,C,PLATE)
SPEDAT(SET,EMISSIVITY,OF!IMB,R,1.)
SPEDAT(SET,OBJNAM,!IMB,C,FLOOR)
SPEDAT(SET,OBJTYP,!IMB,C,PLATE)
SPEDAT(SET,OBJNAM,!OBC,C,FIRE)
SPEDAT(SET,OBJTYP,!OBC,C,BLOCKAGE)
SPEDAT(SET,FIRE,MATERIAL,R,-1.)
SPEDAT(SET,OBJNAM,!OCC,C,FIRE)
SPEDAT(SET,OBJTYP,!OCC,C,BLOCKAGE)
SPEDAT(SET,OBJNAM,!OBD,C,SPRAY)
SPEDAT(SET,OBJTYP,!OBD,C,USER_DEFINED)
SPEDAT(SET,MATERIAL,199,L,T)
LG( 20) = T
IG( 9) = 1
IG( 10) = 1
IG( 14) = 1
IG( 15) = 100
IG( 16) = 5
IG( 17) = 1
IG( 18) = 600
IG( 19) = 1
IG( 20) = 60
RG( 71) = GRND1
RG( 72) =4131.799805
RG( 73) =8.617739E-05
RG( 74) = GRND1
RG( 75) = GRND1
RG( 76) = GRND1
RG( 77) =1.
RG( 80) =1.
RG( 82) =4131.799805
RG( 83) =100.
RG( 84) =1.
RG( 86) =18.
RG( 87) =28.9
RG( 88) = GRND1
RG( 89) =1007.
RG( 90) = GRND1
RG( 91) =273.149994
RG( 92) =273.149994
RG( 93) =1.
RG( 94) =3.335E+05
RG( 95) = GRND1
RG( 96) =0.0263
RG( 97) = GRND1
RG( 98) =-9.81
CG( 6) ='GHIS'
CG( 8) ='spin'
CG( 9) ='NONE'
CG( 10) ='NONE'
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
OUTPUT(P1,Y,N,Y,Y,Y,Y)
OUTPUT(U1,Y,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(W1,Y,N,Y,Y,Y,Y)
OUTPUT(KE,Y,N,Y,Y,Y,Y)
OUTPUT(EP,Y,N,Y,Y,Y,Y)
OUTPUT(MH2O,Y,N,Y,Y,Y,Y)
OUTPUT(PRPS,Y,N,Y,N,N,N)
OUTPUT(QRZ,Y,N,Y,N,N,N)
OUTPUT(ENUL,Y,N,N,N,N,N)
OUTPUT(QRY,Y,N,Y,N,N,N)
OUTPUT(TDEC,Y,N,Y,N,N,N)
OUTPUT(SPH1,Y,N,Y,N,N,N)
OUTPUT(QRX,Y,N,Y,N,N,N)
OUTPUT(LTLS,Y,N,Y,Y,Y,Y)
OUTPUT(WDIS,Y,N,N,N,N,N)
OUTPUT(SCAT,Y,N,Y,N,N,N)
OUTPUT(EMIS,Y,N,Y,N,N,N)
OUTPUT(WGAP,Y,N,Y,N,N,N)
OUTPUT(T3,Y,N,Y,Y,Y,Y)
OUTPUT(MASS,Y,N,Y,N,N,N)
OUTPUT(HEAT,Y,N,Y,N,N,N)
OUTPUT(MOMZ,Y,N,Y,N,N,N)
OUTPUT(MOMY,Y,N,Y,N,N,N)
OUTPUT(MOMX,Y,N,Y,N,N,N)
OUTPUT(TEM1,Y,N,Y,Y,Y,Y)
OUTPUT(EPKE,Y,N,Y,N,N,N)
OUTPUT(DEN1,Y,N,Y,N,N,N)
OUTPUT(EL1,Y,N,Y,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 4 ;IYMON = 3 ;IZMON = 8
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP