TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2020, Date 07/04/20
CPVNAM=VDI; SPPNAM=FLAIR
************************************************************
Echo DISPLAY / USE settings
PHOTON USE
p
f1
xyz
0.20443E+04 0.15633E+04 CR
red
msg Geometry of the problem
use patgeo
set prop on
msg
msg Smoke concentration contours at the end
msg of the first time-step (100 sec)
cont c1 z 1 fil;.001
msg
msg Press return to plot velocity field
pause
vec z 1 sh
msg
msg Press return to load results after second time-step
pause
p
f2
xyz
0.20443E+04 0.15633E+04 CR
red
msg Geometry of the problem
use patgeo
set prop on
msg
msg Smoke concentration contours at the end
msg of the second time-step (200 sec)
msg Note smoke development
cont c1 z 1 fil;.001
msg
msg Press return to plot velocity field
pause
vec z 1 sh
msg
msg Press return to load results after third time-step
pause
p
f3
xyz
0.20443E+04 0.15633E+04 CR
red
msg Geometry of the problem
use patgeo
set prop on
msg
msg Smoke concentration contours at the end
msg of the third time-step (300 sec)
msg Note further smoke development
cont c1 z 1 fil;.001
msg
msg Press return to plot velocity field
pause
vec z 1 sh
msg
msg Press return to load results after last time-step
pause
p
f4
xyz
0.20443E+04 0.15633E+04 CR
red
msg Geometry of the problem
use patgeo
set prop on
msg
msg Smoke concentration contours at the end
msg of simulated period (400 sec)
msg Note that smoke fills the chimney
cont c1 z 1 fil;.001
msg
msg Press return to plot velocity field
pause
vec z 1 sh
msg
msg Type E to end
ENDUSE
DISPLAY
Library case I102:
SMOKE IN A CHIMNEY
This case models a simplifed two-dimensional chimney,
described with the aid of Body Fitted Coordinates (BFC).
There is a small opening at the lower end of the chimney
and the top is fully opened. Fire starts near the bottom
end; it can be clearly seen how the smoke develops with
time and finally fills the chimney.
The flow is treated as laminar.
(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 = 102
************************************************************
Group 1. Run Title
TEXT(Smoke in a chimney: I102 )
************************************************************
Group 2. Transience
STEADY=F
* Set overall time and no. of steps
RSET(U,0.,400.,4)
* Cut regions
RSET(T,L,0,100.)
RSET(T,L,0,200.)
RSET(T,L,0,300.)
* Modify regions
RSET(T,1,1,1.)
RSET(T,2,1,1.)
RSET(T,3,1,1.)
RSET(T,4,1,1.)
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,30,8,1)
* Set overall domain extent:
* xulast yvlast zwlast
name
* Set overall domain extent:
* xulast yvlast zwlast
name
XSI= 3.000000E+01; YSI= 2.000000E+00; ZSI= 1.000000E+00
RSET(D,CHAM )
* Set objects: x0 y0 z0
* dx dy dz
name
XPO= 1.500000E+01; YPO= 7.500000E-01; ZPO= 0.000000E+00
XSI= 5.000000E+00; YSI= 5.000000E-01; ZSI= 1.000000E+00
RSET(B,BFL )
XPO= 4.000000E+00; YPO= 0.000000E+00; ZPO= 0.000000E+00
XSI= 1.000000E+00; YSI= 7.500000E-01; ZSI= 1.000000E+00
RSET(B,TM1FIRE )
XPO= 0.000000E+00; YPO= 7.500000E-01; ZPO= 0.000000E+00
XSI= 0.000000E+00; YSI= 5.000000E-01; ZSI= 1.000000E+00
RSET(B,FANIN )
XPO= 3.000000E+01; YPO= 0.000000E+00; ZPO= 0.000000E+00
XSI= 0.000000E+00; YSI= 2.000000E+00; ZSI= 1.000000E+00
RSET(B,EXIT )
************************************************************
Group 6. Body-Fitted coordinates
BFC=T
* Copy/Transfer/Block grid planes
GSET(C,J4,F,J1,1,30,1,1,+,0,.75,0,INC,1)
GSET(C,J6,F,J4,1,30,1,1,+,0,.5,0,INC,1)
GSET(C,I5,F,I1,1,8,1,1,+,4,0,0,INC,1)
GSET(C,I6,F,I5,1,8,1,1,+,1,0,0,INC,1)
GSET(C,I11,F,I6,1,8,1,1,RZ,1.5708,5,6,INC,1)
GSET(C,I31,F,I11,1,8,1,1,+,0,14,0,INC,1)
GSET(C,J1,F,J1,16,20,1,1,+,2,0,0)
GSET(C,J9,F,J1,16,20,1,1,+,-6,0,0,INC,1)
GSET(T,I16,F,I11,1,8,1,1,1.000000E+00)
GSET(T,I26,F,I21,1,8,1,1,1.000000E+00)
**********
NONORT = T
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(42)=ABSV
NAME(44)=PRPS ;NAME(45)=DEN1
NAME(46)=WCRT ;NAME(47)=TEM1
NAME(48)=SPH1 ;NAME(49)=VCRT
NAME( 50) =UCRT
* Solved variables list
SOLVE(P1,U1,V1,C1,TEM1)
* Stored variables list
STORE(UCRT,VCRT,SPH1,WCRT,DEN1,PRPS,ABSV)
* Additional solver options
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
NEWRH1 = T
NEWENL = T
ISOLX = 0 ;ISOLY = 0 ;ISOLZ = 0
DENPCO = T
************************************************************
Group 9. Properties
PRESS0 =1.0E+05
* Domain material index is 2 signifying:
* Air using Ideal Gas Law, STP
SETPRPS(1, 2)
ENUT =0.
DRH1DP = GRND5
DVO1DT =3.33E-03
PRNDTL(TEM1)=-0.0263
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(C1)=1.001E-10 ;FIINIT(ABSV)=1.001E-10
FIINIT(PRPS)=-1.
FIINIT(DEN1)=1.161 ;FIINIT(WCRT)=1.001E-10
FIINIT(TEM1)=292. ;FIINIT(SPH1)=1.001E-10
FIINIT(VCRT)=1.001E-10 ;FIINIT(UCRT)=1.001E-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, 4)
COVAL(BUOYANCY, U1, FIXFLU, GRND2)
COVAL(BUOYANCY, V1, FIXFLU, GRND2)
BUOYA =0. ; BUOYB =-9.81 ; BUOYC =0.
BUOYD =1.161
BFCA =1.161
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 450
RESFAC =1.0E-02
************************************************************
Group 16. Terminate Iterations
LITER(TEM1)=15
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,0.3 )
RELAX(U1 ,FALSDT,0.15 )
RELAX(V1 ,FALSDT,0.15 )
RELAX(DEN1,LINRLX,0.3 )
RELAX(TEM1,FALSDT,1. )
************************************************************
Group 18. Limits
VARMAX(TEM1)=1.0E+10 ;VARMIN(TEM1)=-30.
************************************************************
Group 19. EARTH Calls To GROUND Station
NAMGRD =FLAR
GENK = T
PARSOL = F
ISG62 = 1
SPEDAT(SET,FACETDAT,NUMOBJ,I,5)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
RG( 41) =3.0E+05
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 25 ;IYMON = 1 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 10000
NXPRIN = 6
NYPRIN = 0
NPLT = 1
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0
CSG1 ='F'
GVIEW(P,0.,0.,1.)
GVIEW(UP,0.,1.,0.)
> DOM, SIZE, 3.000000E+01, 8.000000E+00, 1.000000E+00
> DOM, MONIT, 2.500000E+01, 1.000000E+00, 1.000000E+00
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID, RSET_X_1, 4, 1.000000E+00
> GRID, RSET_X_2, 1, 1.000000E+00
> GRID, RSET_X_3, 5, 1.000000E+00
> GRID, RSET_X_4, 5, 1.000000E+00
> GRID, RSET_X_5, 5, 1.000000E+00
> GRID, RSET_X_6, 5, 1.000000E+00
> GRID, RSET_X_7, 5, 1.000000E+00
> GRID, RSET_Y_1, 3, 1.000000E+00
> GRID, RSET_Y_2, 2, 1.000000E+00
> GRID, RSET_Y_3, 3, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> OBJ, NAME, BFL
> OBJ, POSITION, 1.500000E+01, 3.000000E+00, 0.000000E+00
> OBJ, SIZE, 5.000000E+00, 2.000000E+00, 1.000000E+00
> OBJ, GEOMETRY, BFL
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 117,Cast concrete (dense)
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, NAME, TM1FIRE
> OBJ, POSITION, 4.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 3.000000E+00, 1.000000E+00
> OBJ, GEOMETRY, TM1FIRE
> OBJ, TYPE, BLOCKAGE
> OBJ, COLOR-MODE, FROM_FIL
> OBJ, COLOR-VAL, 0
> OBJ, OPAQUE, 50
> OBJ, MATERIAL, DOMAIN
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, LIN_T_HVA, 0. ,3.0E+05
> OBJ, SCAL_FIXF, Per-unit-volume
> OBJ, FIX_VAL_C1, FIXVAL ,1.
> OBJ, INI_C1, 0.
> OBJ, NAME, FANIN
> OBJ, POSITION, 0.000000E+00, 3.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 2.000000E+00, 1.000000E+00
> OBJ, GEOMETRY, FANIN
> OBJ, TYPE, INLET
> OBJ, COLOR-MODE, FROM_FIL
> OBJ, COLOR-VAL, 0
> OBJ, OPAQUE, 50
> OBJ, DENSITY, 1.161
> OBJ, VELOCITY, 1.2 ,0. ,0.
> OBJ, TEMPERATURE, 293.
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, NAME, EXIT
> OBJ, POSITION, 3.000000E+01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 8.000000E+00, 1.000000E+00
> OBJ, GEOMETRY, EXIT
> OBJ, TYPE, OPENING
> OBJ, COLOR-MODE, FROM_FIL
> OBJ, COLOR-VAL, 0
> OBJ, OPAQUE, 50
> OBJ, PRESSURE, 0.
> OBJ, TEMPERATURE, 293.
> OBJ, COEFFICIENT, 1.
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
DISTIL=T
RESFAC=1E-6
SPEDAT(SET,SELREF,RESFC,R,1E-6)
STOP