TALK=T;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 3.6, Date 23/09/04
 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 =         0
 ************************************************************
  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.000000E+00,4.000000E+02,4)
    * Cut regions
 RSET(T,L,0,1.000000E+02)
 RSET(T,L,0,2.000000E+02)
 RSET(T,L,0,3.000000E+02)
    * Modify regions
 RSET(T,1,1,1.000000E+00)
 RSET(T,2,1,1.000000E+00)
 RSET(T,3,1,1.000000E+00)
 RSET(T,4,1,1.000000E+00)
 ************************************************************
  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
 ONEPHS  =    T
    * Non-default variable names
 NAME( 42) =ABSV ; NAME( 43) =BLOK
 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,BLOK,ABSV)
    * Additional solver options
 SOLUTN(TEM1,Y,Y,Y,N,N,Y)
 IVARBK  =        -1 ;ISOLBK =         1
 
 ************************************************************
  Group 8. Terms & Devices
 NEWRH1  =    T
 NEWENL  =    T
 ISOLX   =         0 ;ISOLY  =         0 ;ISOLZ  =         0
 DENPCO  =    T
 ************************************************************
  Group 9. Properties
 PRESS0  = 1.000000E+05
    * Domain material is: 2   Air using Ideal Gas Law, STP
 SETPRPS(1,  2)
 ENUT    = 0.000000E+00
 DRH1DP  = GRND5
 DVO1DT  = 3.330000E-03
 PRNDTL(TEM1) = -2.630000E-02
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(C1  ) =  1.001000E-10 ;FIINIT(ABSV) =  1.001000E-10
 FIINIT(BLOK) =  1.000000E+00 ;FIINIT(PRPS) = -1.000000E+00
 FIINIT(DEN1) =  1.161000E+00 ;FIINIT(WCRT) =  1.001000E-10
 FIINIT(TEM1) =  2.920000E+02 ;FIINIT(SPH1) =  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,4)
 COVAL (BUOYANCY,U1  , FIXFLU      , GRND2       )
 COVAL (BUOYANCY,V1  , FIXFLU      , GRND2       )
 
 BUOYA   = 0.000000E+00 ; BUOYB =-9.810000E+00
 BUOYC   = 0.000000E+00
 BUOYD   = 1.161000E+00
 BFCA    = 1.161000E+00
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =       450
 RESFAC  = 1.000000E-02
 ************************************************************
  Group 16. Terminate Iterations
 LITER (TEM1) =       15
 ENDIT (TEM1) =  GRND1
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 3.000000E-01)
 RELAX(U1  ,FALSDT, 1.500000E-01)
 RELAX(V1  ,FALSDT, 1.500000E-01)
 RELAX(DEN1,LINRLX, 3.000000E-01)
 RELAX(TEM1,FALSDT, 1.000000E+00)
 ************************************************************
  Group 18. Limits
 VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 NAMGRD  =FLAR
 GENK    =    T
 PARSOL  =    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,5)
 SPEDAT(SET,IMMERSOL,EMCON,R,5.00000E+00)
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 RG( 41) = 3.000000E+05
 ************************************************************
  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   =        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.000000E+00,0.000000E+00,1.000000E+00)
 GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
 
> 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,    SNAPSIZE,    1.000000E-02
 
> 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,    INI_BLOK,            2.000000E+00
 
> 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,  DEFAULT2
> OBJ,    COLOR-VAL,          0
> OBJ,    OPAQUE,            50
> OBJ,    MATERIAL,    DOMAIN
> OBJ,    TIME_LIMITS,   ALWAYS_ACTIVE
> OBJ,    LIN_T_HVA,     0.000000E+00, 3.000000E+05
> OBJ,    SCAL_FIXF,           2.000000E+00
> OBJ,    FIX_VAL_C1,    FIXVAL      , 1.000000E+00
> OBJ,    INI_C1,              0.000000E+00
 
> 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,  DEFAULT2
> OBJ,    COLOR-VAL,          0
> OBJ,    OPAQUE,            50
> OBJ,    DENSITY,             1.161000E+00
> OBJ,    VELOCITY,      1.200000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    TEMPERATURE,         2.930000E+02
> 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,  DEFAULT2
> OBJ,    COLOR-VAL,          0
> OBJ,    OPAQUE,            50
> OBJ,    PRESSURE,            0.000000E+00
> OBJ,    TEMPERATURE,         2.930000E+02
> OBJ,    COEFFICIENT,         1.000000E+00
> OBJ,    TIME_LIMITS,   ALWAYS_ACTIVE
STOP