TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.6, Date 08/09/04
CPVNAM=VDI;SPPNAM=FLAIR
************************************************************
Echo DISPLAY / USE settings
PHOTON USE
p
xyz
0.20443E+04 0.15633E+04 CR
red
MSG Geometry of the problem
use patgeo
msg
msg Surface of smoke concentration of 0.3
surf c1 x 0.3 col 12
surf c1 y 0.3 col 12
surf c1 z 0.3 col 12
msg
msg Press return to plot velocities
pause
vec z 1 sh; vec z m sh
vec x 5 sh
msg
msg Type E to end
ENDUSE
DISPLAY
Library case I103:
FIRE IN A BUILDING
This case represents three-dimensional section of a
building, described with the use of Body Fitted
Coordinates (BFC). There are two air-condition vents and
several partitions on both floors. Fire is burning in one
corner of the atrium-like space and it can spread to both
floors. Turbulence is taken into account with the help of
k-e model.
(Use PHOTON to display results. Enter USE Q1 in PHOTON
command mode, or create a file called U, which contains
the one line: USE Q1 before typing RUNPHO.)
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 103
************************************************************
Group 1. Run Title
TEXT(Fire in a building: I103 )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,10,10,20)
* Set overall domain extent:
* xulast yvlast zwlast
name
* Set overall domain extent:
* xulast yvlast zwlast
name
XSI= 6.000000E+00; YSI= 5.000000E+00; ZSI= 1.500000E+01
RSET(D,CHAM )
* Set objects: x0 y0 z0
* dx dy dz
name
XPO= 0.000000E+00; YPO= 2.500000E+00; ZPO= 0.000000E+00
XSI= 6.000000E+00; YSI= 0.000000E+00; ZSI= 9.000000E+00
RSET(B,FLOOR )
XPO= 0.000000E+00; YPO= 2.500000E+00; ZPO= 4.500000E+00
XSI= 4.200000E+00; YSI= 2.500000E+00; ZSI= 0.000000E+00
RSET(B,PART1 )
XPO= 1.800000E+00; YPO= 0.000000E+00; ZPO= 4.500000E+00
XSI= 4.200000E+00; YSI= 2.000000E+00; ZSI= 0.000000E+00
RSET(B,PART2 )
XPO= 1.800000E+00; YPO= 0.000000E+00; ZPO= 9.000000E+00
XSI= 4.200000E+00; YSI= 2.000000E+00; ZSI= 0.000000E+00
RSET(B,PART3 )
XPO= 3.000000E+00; YPO= 0.000000E+00; ZPO= 1.350000E+01
XSI= 1.200000E+00; YSI= 1.500000E+00; ZSI= 1.499999E+00
RSET(B,FIRE )
XPO= 3.000000E+00; YPO= 2.000000E+00; ZPO= 0.000000E+00
XSI= 6.000001E-01; YSI= 5.000000E-01; ZSI= 0.000000E+00
RSET(B,FN1 )
XPO= 2.400000E+00; YPO= 5.000000E+00; ZPO= 1.350000E+01
XSI= 5.999999E-01; YSI= 0.000000E+00; ZSI= 7.500000E-01
RSET(B,FN2 )
XPO= 1.800000E+00; YPO= 4.000000E+00; ZPO= 0.000000E+00
XSI= 2.400000E+00; YSI= 1.000000E+00; ZSI= 0.000000E+00
RSET(B,WN1 )
XPO= 1.800000E+00; YPO= 2.500000E+00; ZPO= 1.500000E+01
XSI= 2.400000E+00; YSI= 2.500000E+00; ZSI= 0.000000E+00
RSET(B,WN2 )
************************************************************
Group 6. Body-Fitted coordinates
BFC=T
* Set points
XPO= 0.000000E+00;YPO= 0.000000E+00;ZPO= 1.500000E+01;GSET(P,A )
XPO= 0.000000E+00;YPO= 2.500000E+00;ZPO= 1.500000E+01;GSET(P,B )
XPO= 6.000000E+00;YPO= 2.500000E+00;ZPO= 1.500000E+01;GSET(P,C )
XPO= 6.000000E+00;YPO= 0.000000E+00;ZPO= 1.500000E+01;GSET(P,D )
* Set lines/arcs
GSET(L,AB,A,B,10,1.0)
GSET(L,BC,B,C,10,1.0)
GSET(L,CD,C,D,10,1.0)
GSET(L,DA,D,A,10,1.0)
* Set frames
GSET(F,ABCD,A,-,B,-,C,-,D,-)
* Match a grid mesh
GSET(M,ABCD,+J+I,1,1,21,TRANS)
* Copy/Transfer/Block grid planes
GSET(C,I4,F,I1,1,10,1,20,+,2,0,0,INC,1)
GSET(C,I5,F,I4,1,10,1,20,+,.5,0,0,INC,1)
GSET(C,I6,F,I5,1,10,1,20,+,.5,0,0,INC,1)
GSET(C,I7,F,I6,1,10,1,20,+,.5,0,0,INC,1)
GSET(C,I8,F,I7,1,10,1,20,+,.5,0,0,INC,1)
GSET(C,I11,F,I8,1,10,1,20,+,2,0,0,INC,1)
GSET(C,J4,F,J1,1,10,1,20,+,0,1.5,0,INC,1)
GSET(C,J5,F,J4,1,10,1,20,+,0,.5,0,INC,1)
GSET(C,J6,F,J5,1,10,1,20,+,0,.5,0,INC,1)
GSET(C,J9,F,J6,1,10,1,20,+,0,1.5,0,INC,1)
GSET(C,J11,F,J9,1,10,1,20,+,0,1,0,INC,1)
GSET(C,K7,F,K1,1,10,1,10,+,0,0,5,INC,1)
GSET(C,K13,F,K7,1,10,1,10,+,0,0,5,INC,1)
GSET(C,K19,F,K13,1,10,1,10,+,0,0,4,INC,1)
GSET(C,K20,F,K19,1,10,1,10,+,0,0,.5,INC,1)
GSET(C,K21,F,K20,1,10,1,10,+,0,0,.5,INC,1)
GSET(T,K21,F,K13,1,10,1,10,1.000000E+00)
**********
NONORT = T
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME( 38) =HPOR ; NAME( 39) =NPOR
NAME( 42) =ABSV ; NAME( 43) =BLOK
NAME( 44) =PRPS ; NAME( 46) =DEN1
NAME( 47) =TEM1 ; NAME( 48) =WCRT
NAME( 49) =VCRT ; NAME( 50) =UCRT
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,W1 ,C1 ,TEM1)
* Stored variables list
STORE(UCRT,VCRT,WCRT,DEN1,PRPS,BLOK,ABSV,NPOR)
STORE(HPOR)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,N)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
IVARBK = -1 ;ISOLBK = 1
TURMOD(KEMODL)
************************************************************
Group 8. Terms & Devices
NEWRH1 = T
NEWENL = T
ISOLX = 0 ;ISOLY = 0 ;ISOLZ = 0
DENPCO = T
************************************************************
Group 9. Properties
PRESS0 = 1.000000E+05 ;TEMP0 = 2.730000E+02
* Domain material is: 2 Air using Ideal Gas Law, STP
SETPRPS(1, 2)
DRH1DP = GRND5
DVO1DT = 3.330000E-03
PRNDTL(TEM1) = -2.630000E-02
* Non-default SOLUTN setting for KE
SOLUTN(KE ,Y,Y,N,N,N,Y)
* Non-default SOLUTN setting for EP
SOLUTN(EP ,Y,Y,N,N,N,Y)
PRT (EP ) = 1.314000E+00
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(KE ) = 7.225000E-03 ;FIINIT(EP ) = 1.007000E-02
FIINIT(C1 ) = 1.001000E-10 ;FIINIT(HPOR) = 1.000000E+00
FIINIT(NPOR) = 1.000000E+00 ;FIINIT(ABSV) = 1.001000E-10
FIINIT(BLOK) = 1.000000E+00 ;FIINIT(PRPS) = -1.000000E+00
FIINIT(DEN1) = 1.161000E+00 ;FIINIT(TEM1) = 2.000000E+01
FIINIT(WCRT) = 1.001000E-10 ;FIINIT(VCRT) = 1.001000E-10
FIINIT(UCRT) = 1.001000E-10
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 )
BUOYA = 0.000000E+00 ; BUOYB =-9.810000E+00
BUOYC = 0.000000E+00
BUOYD = 1.189000E+00
BFCA = 1.189000E+00
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 100
RESFAC = 1.000000E-02
************************************************************
Group 16. Terminate Iterations
LITER (TEM1) = 15
ENDIT (TEM1) = GRND1
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000000E+00)
RELAX(U1 ,FALSDT, 7.500000E-02)
RELAX(V1 ,FALSDT, 7.500000E-02)
RELAX(W1 ,FALSDT, 7.500000E-02)
RELAX(KE ,FALSDT, 7.500000E-02)
RELAX(EP ,FALSDT, 7.500000E-02)
RELAX(DEN1,LINRLX, 7.000000E-01)
RELAX(TEM1,FALSDT, 1.000000E+02)
KELIN = 1
************************************************************
Group 18. Limits
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
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
NAMGRD =FLAR
GENK = T
PARSOL = T
LSG57 = T
SPEDAT(SET,FLAIR,CLOTHING,R,9.30000E-02)
SPEDAT(SET,FLAIR,RELHUMID,R,5.00000E+01)
SPEDAT(SET,FLAIR,TRAD,R,2.50000E+01)
SPEDAT(SET,FACETDAT,NUMOBJ,I,10)
SPEDAT(SET,IMMERSOL,EMCON,R,5.00000E+00)
************************************************************
Group 20. Preliminary Printout
ECHO = T
NULLPR = F
NDST = 0
DSTTOL = 1.000000E-02
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 7 ;IZMON = 20
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NXPRIN = 1
NYPRIN = 1
NZPRIN = 4
NPLT = 1
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,5.773503E-01,5.773503E-01,5.773503E-01)
GVIEW(UP,-4.082483E-01,8.164966E-01,-4.082483E-01)
> DOM, SIZE, 1.000000E+01, 1.000000E+01, 2.000000E+01
> DOM, MONIT, 1.000000E+00, 7.000000E+00, 2.000000E+01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> OBJ, NAME, FLOOR
> OBJ, POSITION, 0.000000E+00, 5.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+01, 0.000000E+00, 1.200000E+01
> OBJ, GEOMETRY, FLOOR
> OBJ, TYPE, PLATE
> OBJ, POROSITY, 0.000000E+00
> OBJ, SIDE, BOTH
> OBJ, NAME, PART1
> OBJ, POSITION, 0.000000E+00, 5.000000E+00, 6.000000E+00
> OBJ, SIZE, 7.000000E+00, 5.000000E+00, 0.000000E+00
> OBJ, GEOMETRY, PART1
> OBJ, TYPE, PLATE
> OBJ, POROSITY, 0.000000E+00
> OBJ, SIDE, BOTH
> OBJ, NAME, PART2
> OBJ, POSITION, 3.000000E+00, 0.000000E+00, 6.000000E+00
> OBJ, SIZE, 7.000000E+00, 4.000000E+00, 0.000000E+00
> OBJ, GEOMETRY, PART2
> OBJ, TYPE, PLATE
> OBJ, POROSITY, 0.000000E+00
> OBJ, SIDE, BOTH
> OBJ, NAME, PART3
> OBJ, POSITION, 3.000000E+00, 0.000000E+00, 1.200000E+01
> OBJ, SIZE, 7.000000E+00, 4.000000E+00, 0.000000E+00
> OBJ, GEOMETRY, PART3
> OBJ, TYPE, PLATE
> OBJ, POROSITY, 0.000000E+00
> OBJ, SIDE, BOTH
> OBJ, NAME, FIRE
> OBJ, POSITION, 5.000000E+00, 0.000000E+00, 1.800000E+01
> OBJ, SIZE, 2.000000E+00, 3.000000E+00, 2.000000E+00
> OBJ, GEOMETRY, FIRE
> OBJ, TYPE, BLOCKAGE
> OBJ, COLOR-MODE, DEFAULT2
> OBJ, COLOR-VAL, 0
> OBJ, OPAQUE, 50
> OBJ, MATERIAL, DOMAIN
> OBJ, HEAT_FLVX, 0.000000E+00, 7.500000E+05
> OBJ, SCAL_FIXF, 2.000000E+00
> OBJ, FIX_VAL_C1, FIXVAL , 1.000000E+00
> OBJ, INI_C1, 0.000000E+00
> OBJ, NAME, FN1
> OBJ, POSITION, 5.000000E+00, 4.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ, GEOMETRY, FN1
> OBJ, TYPE, INLET
> OBJ, COLOR-MODE, DEFAULT2
> OBJ, COLOR-VAL, 0
> OBJ, OPAQUE, 50
> OBJ, DENSITY, 1.189000E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 2.000000E+00
> OBJ, TEMPERATURE, 2.000000E+01
> OBJ, KE_IN, 1.000000E-02
> OBJ, EP_IN, 1.640000E-02
> OBJ, NAME, FN2
> OBJ, POSITION, 4.000000E+00, 1.000000E+01, 1.800000E+01
> OBJ, SIZE, 1.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ, GEOMETRY, FN2
> OBJ, TYPE, INLET
> OBJ, COLOR-MODE, DEFAULT2
> OBJ, COLOR-VAL, 0
> OBJ, OPAQUE, 50
> OBJ, DENSITY, 1.189000E+00
> OBJ, VELOCITY, 0.000000E+00,-1.700000E+00, 0.000000E+00
> OBJ, TEMPERATURE, 2.000000E+01
> OBJ, KE_IN, 7.225000E-03
> OBJ, EP_IN, 1.007000E-02
> OBJ, NAME, WN1
> OBJ, POSITION, 3.000000E+00, 8.000000E+00, 0.000000E+00
> OBJ, SIZE, 4.000000E+00, 2.000000E+00, 0.000000E+00
> OBJ, GEOMETRY, WN1
> OBJ, TYPE, OPENING
> OBJ, COLOR-MODE, DEFAULT2
> OBJ, COLOR-VAL, 0
> OBJ, OPAQUE, 50
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, 2.000000E+01
> OBJ, COEFFICIENT, 1.000000E+00
> OBJ, TURBULENCE, SAME , SAME
> OBJ, NAME, WN2
> OBJ, POSITION, 3.000000E+00, 5.000000E+00, 2.000000E+01
> OBJ, SIZE, 4.000000E+00, 5.000000E+00, 0.000000E+00
> OBJ, GEOMETRY, WN2
> OBJ, TYPE, OPENING
> OBJ, COLOR-MODE, DEFAULT2
> OBJ, COLOR-VAL, 0
> OBJ, OPAQUE, 50
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, 2.000000E+01
> OBJ, COEFFICIENT, 1.000000E+00
> OBJ, TURBULENCE, SAME , SAME
STOP