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;;;;
msg Geometry of the problem
use patgeo
msg
msg Pressure distribution.
con p1 z 1 fi;.1
con p1 y 1 fi;.1
con p1 x 8 fi;.1
msg Press return to continue.
pause
con cl
red
msg Pollutant distribution.
con c1 y 1 fi;.1
con c1 x 8 fi;.1
msg
msg Press return to plot velocity field
pause
con cl
red
vec x 8 sh;vec y 1 sh;vec z 1 sh
msg
msg Type E to end.
ENDUSE
DISPLAY
Library case I107:
EXHAUST DISPERSION NEAR A TOWER BLOCK
This Q1 is an example of modelling urban pollution. The
wind flows around a 31-storey tower block with a low rise
block on the other side of the street.
The source of pollution comes from the vehicular exhaust.
The road passes between two residential blocks.
Although FLAIR does not have capability of automatically
setting wind profile, this example shows that the same
effect can be achieved by defining several 'fans' with
height-dependent inflow. The only non-standard setting is
switching-off buoyancy Patch which has to be done from a
command-mode level by typing: BUOYANCY=SKIP.
(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 = 107
************************************************************
Group 1. Run Title
TEXT(Exhaust dispersion near tower block:I107)
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,13,15,21)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME( 45) =HPOR ; NAME( 46) =NPOR
NAME( 47) =EPOR ; NAME( 48) =VPOR
NAME( 49) =DEN1 ; NAME( 50) =PRPS
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,W1 ,C1 )
* Stored variables list
STORE(PRPS,DEN1,VPOR,EPOR,NPOR,HPOR)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,N)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
* Domain material is: 0 Air at 20 deg C, 1 atm
SETPRPS(1, 0)
ENUT = 0.000000E+00
DVO1DT = 3.410000E-03
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(C1 ) = 1.001000E-10 ;FIINIT(HPOR) = 1.000000E+00
FIINIT(NPOR) = 1.000000E+00 ;FIINIT(EPOR) = 1.000000E+00
FIINIT(VPOR) = 1.000000E+00 ;FIINIT(DEN1) = 1.001000E-10
FIINIT(PRPS) = -1.000000E+00
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
No PATCHes used for this Group
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 100
RESFAC = 1.000000E-02
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 5.000000E-01)
RELAX(U1 ,FALSDT, 4.139000E-01)
RELAX(V1 ,FALSDT, 4.139000E-01)
RELAX(W1 ,FALSDT, 4.139000E-01)
RELAX(C1 ,FALSDT, 6.667000E+00)
************************************************************
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
ASAP = 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,23)
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 = 6 ;IYMON = 3 ;IZMON = 20
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NXPRIN = 2
NYPRIN = 3
NZPRIN = 4
NPLT = 1
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,8.564056E-01,4.528628E-01,-2.479614E-01)
GVIEW(UP,-4.338343E-01,8.915725E-01,1.299470E-01)
> DOM, SIZE, 1.300000E+02, 1.500000E+02, 2.100000E+02
> DOM, MONIT, 5.500000E+01, 2.500000E+01, 1.950000E+02
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> GRID, RSET_X_1, 5, 1.000000E+00
> GRID, RSET_X_2, 3, 1.000000E+00
> GRID, RSET_X_3, 5, 1.000000E+00
> GRID, RSET_Y_1, 1, 1.000000E+00
> GRID, RSET_Y_2, 1, 1.000000E+00
> GRID, RSET_Y_3, 1, 1.000000E+00
> GRID, RSET_Y_4, 1, 1.000000E+00
> GRID, RSET_Y_5, 1, 1.000000E+00
> GRID, RSET_Y_6, 1, 1.000000E+00
> GRID, RSET_Y_7, 1, 1.000000E+00
> GRID, RSET_Y_8, 1, 1.000000E+00
> GRID, RSET_Y_9, 1, 1.000000E+00
> GRID, RSET_Y_10, 1, 1.000000E+00
> GRID, RSET_Y_11, 1, 1.000000E+00
> GRID, RSET_Y_12, 1, 1.000000E+00
> GRID, RSET_Y_13, 1, 1.000000E+00
> GRID, RSET_Y_14, 1, 1.000000E+00
> GRID, RSET_Y_15, 1, 1.000000E+00
> GRID, RSET_Z_1, 5, 1.000000E+00
> GRID, RSET_Z_2, 3, 1.000000E+00
> GRID, RSET_Z_3, 1, 1.000000E+00
> GRID, RSET_Z_4, 4, 1.000000E+00
> GRID, RSET_Z_5, 1, 1.000000E+00
> GRID, RSET_Z_6, 2, 1.000000E+00
> GRID, RSET_Z_7, 5, 1.000000E+00
> OBJ, NAME, TOWGR
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 6.000000E+01
> OBJ, SIZE, 1.300000E+02, 1.000000E+01, 3.000000E+01
> OBJ, GEOMETRY, cube14
> OBJ, ROTATION24, 1
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, TOWER
> OBJ, POSITION, 5.000000E+01, 1.000000E+01, 6.000001E+01
> OBJ, SIZE, 3.000000E+01, 9.000000E+01, 2.999999E+01
> OBJ, GEOMETRY, cube14
> OBJ, ROTATION24, 1
> OBJ, VISIBLE, NO
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, REST
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.400000E+02
> OBJ, SIZE, 1.300000E+02, 3.000000E+01, 2.000000E+01
> OBJ, GEOMETRY, cube14
> OBJ, ROTATION24, 1
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, INL1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 1.000000E+01, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.188929E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 2.186000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, INL2
> OBJ, POSITION, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 1.000000E+01, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.189086E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 2.877000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, INL3
> OBJ, POSITION, 0.000000E+00, 2.000000E+01, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 9.999998E+00, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.189049E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 3.269000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, INL4
> OBJ, POSITION, 0.000000E+00, 3.000000E+01, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 1.000000E+01, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.188976E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 3.556000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, INL5
> OBJ, POSITION, 0.000000E+00, 4.000000E+01, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 9.999996E+00, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.189068E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 3.787000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, INL6
> OBJ, POSITION, 0.000000E+00, 5.000000E+01, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 1.000000E+01, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.189101E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 3.982000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, INL7
> OBJ, POSITION, 0.000000E+00, 6.000000E+01, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 1.000000E+01, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.189065E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 4.152000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, INL8
> OBJ, POSITION, 0.000000E+00, 7.000000E+01, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 1.000001E+01, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.188938E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 4.303000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, INL9
> OBJ, POSITION, 0.000000E+00, 8.000002E+01, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 9.999985E+00, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.189006E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 4.439000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, INL10
> OBJ, POSITION, 0.000000E+00, 9.000000E+01, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 1.000000E+01, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.189047E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 4.565000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, INL11
> OBJ, POSITION, 0.000000E+00, 1.000000E+02, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 1.000000E+01, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.189103E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 4.680000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, INL12
> OBJ, POSITION, 0.000000E+00, 1.100000E+02, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 1.000000E+01, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.189014E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 4.788000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, INL13
> OBJ, POSITION, 0.000000E+00, 1.200000E+02, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 1.000000E+01, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.188996E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 4.889000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, INL14
> OBJ, POSITION, 0.000000E+00, 1.300000E+02, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 1.000000E+01, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.189005E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 4.984000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, INL15
> OBJ, POSITION, 0.000000E+00, 1.400000E+02, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 1.000000E+01, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.189003E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 5.074000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.547000E-08
> OBJ, NAME, OUT
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 2.100000E+02
> OBJ, SIZE, 1.300000E+02, 1.500000E+02, 0.000000E+00
> OBJ, GEOMETRY, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OPENING
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, SAME
> OBJ, COEFFICIENT, 1.000000E+00
> OBJ, NAME, TOP
> OBJ, POSITION, 0.000000E+00, 1.500000E+02, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 0.000000E+00, 2.100000E+02
> OBJ, GEOMETRY, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OPENING
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, SAME
> OBJ, COEFFICIENT, 1.000000E+00
> OBJ, NAME, GROUN
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.300000E+02, 0.000000E+00, 2.100000E+02
> OBJ, GEOMETRY, cube11
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, NAME, SRCE
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 9.400000E+01
> OBJ, SIZE, 1.300000E+02, 0.000000E+00, 4.200000E+01
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.189000E+00
> OBJ, VELOCITY, 0.000000E+00, 5.000000E-01, 0.000000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, INLET_C1, 1.800000E-08
STOP