TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2019, Date 02/08/19
CPVNAM=VDI; SPPNAM=FLAIR
************************************************************
Echo DISPLAY / USE settings
DISPLAY
Library Case I407: Aerosol dispersion and deposition in a small
room.
The case considered is steady, isothermal airflow with aerosol
transport and deposition of three particle sizes in a small room.
The Eulerian drift flux model is used to simulate aerosol
dispersion with surface deposition via the 3-layer deposition
model of Chen & Lai (2004), which accounts for the deposition
mechanisms of gravity, Brownian & turbulent diffusion.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(I407: Aerosol deposition in a small room)
************************************************************
Echo save-block settings for Group 1
save1begin
Aerosol particles of 2, 6 and 10 microns and each with a
density of 1000 kg/m^3 enter a 2m cube room vertically
with a central supply of air at floor level. The air supply
has a Reynolds number of 1000, an extraction fan is located
on the ceiling with a rating of 0.006 m^3/s, and the room
is naturally ventilated by an opening located near the
floor of the north wall. The flow is isothermal with a
temperature of 20degC. In the simulations, each particle
concentration is normalised by its inlet value.
References
----------
Chen, F.Z. & Lai, A.C.K, "An Eulerian model for particle
deposition under electrostatic and turbulent conditions",
J.Aerosol Science, Vol.35, p47-62, (2004).
save1end
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,25,39,58)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(126)=EPKE ;NAME(127)=V3DH
NAME(128)=VD3 ;NAME(129)=VD2
NAME(130)=VS3 ;NAME(131)=VS1
NAME(132)=VDPH ;NAME(133)=DI1
NAME(134)=VP1 ;NAME(135)=VSTR
NAME(136)=TP1 ;NAME(137)=V2DH
NAME(138)=PRPS ;NAME(139)=VABS
NAME(140)=YPLS ;NAME(141)=VD1
NAME(142)=STRS ;NAME(143)=V1DH
NAME(144)=ENUL ;NAME(145)=SPH1
NAME(146)=VS2 ;NAME(147)=TEM1
NAME(148)=DEN1 ;NAME(149)=EL1
NAME(150) =ENUT
* Solved variables list
SOLVE(P1,U1,V1,W1,C6,C7,C8,TEM1)
* Stored variables list
STORE(ENUT,EL1,DEN1,VS2,SPH1,ENUL,V1DH,STRS)
STORE(VD1,YPLS,VABS,PRPS,V2DH,TP1,VSTR,VP1)
STORE(DI1,VDPH,VS1,VS3,VD2,VD3,V3DH,EPKE)
* Additional solver options
SOLUTN(P1,Y,Y,Y,N,N,Y)
SOLUTN(C6,Y,Y,Y,N,N,Y)
SOLUTN(C7,Y,Y,Y,N,N,Y)
SOLUTN(C8,Y,Y,Y,N,N,Y)
TURMOD(KERNG)
************************************************************
Group 8. Terms & Devices
NEWRH1 = T
************************************************************
Group 9. Properties
PRESS0 =1.01325E+05 ;TEMP0 =273.
* Domain material index is 2 signifying:
* Air using Ideal Gas Law, STP
SETPRPS(1, 2)
ENUL =1.0E-05
ENULA =-4.94679E-06 ;ENULB =4.58394E-08 ;ENULC =8.0974E-11
DRH1DP = GRND5
DVO1DT =3.33E-03
PRNDTL(C6)=0.7 ;PRNDTL(TEM1)=-0.0263
PRT(KE)=0.7194 ;PRT(EP)=0.7194
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1)=0. ;FIINIT(W1)=0.225
FIINIT(C6)=0. ;FIINIT(C7)=0.
FIINIT(C8)=0. ;FIINIT(PRPS)=-1.
FIINIT(TEM1)=20.
No PATCHes used for this Group
INIADD = F
NAMFI ='phida'
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
PATCH(DFLUX, CELL, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(DFLUX, C6, GRND4, GRND4)
COVAL(DFLUX, C7, GRND4, GRND4)
COVAL(DFLUX, C8, GRND4, GRND4)
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. ; BUOYB =0. ; BUOYC =-9.81
BUOYD =1.204938
EGWF = T
************************************************************
Echo save-block settings for Group 13
save13begin
** compute average floor deposition velocity in m/hr & floor-averaged C
(make1 SUMAF is 0.0)
(store1 SUMAF at PW4 is SUM(AHIGH))
(print Floor_area is SUMAF)
INTEGER(KK,NPART);KK=6;NPART=3
DO JJ=1,NPART
** deposition velocity in m/hr
(STORED of V:JJ:DH is VD:JJ:*3600.0)
** compute particle deposition fraction
(MAKE1 P:JJ:DFR is 0.0)
(make1 SUMVD:JJ: is 0.0)
(make1 SUMC:KK: is 0.0)
(STORE1 P:JJ:DFR is ABS(NETS(C:KK:,PW4))/NETS(C:KK:,INL))
(store1 SUMVD:JJ: at PW4 is SUM(V:JJ:DH*AHIGH))
(store1 SUMC:KK: at PW4 is SUM(C:KK:*AHIGH))
(PRINT Par:JJ:_dfrac is P:JJ:DFR)
(print Floor_avg_Vdep:JJ: is SUMVD:JJ:/SUMAF)
(print Floor_avg_C:KK: is SUMC:KK:/SUMAF)
KK=KK+1
ENDDO
save13end
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 2000
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
LITER(KE)=10 ;LITER(EP)=10
LITER(C6)=50 ;LITER(C7)=50
LITER(C8)=50 ;LITER(TEM1)=50
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
RELAX(U1 ,FALSDT,10. )
RELAX(V1 ,FALSDT,10. )
RELAX(W1 ,FALSDT,10. )
RELAX(KE ,LINRLX,0.5 )
RELAX(EP ,LINRLX,0.5 )
RELAX(C6 ,FALSDT,10. )
RELAX(C7 ,FALSDT,10. )
RELAX(C8 ,FALSDT,10. )
RELAX(TEM1,LINRLX,0.5 )
KELIN = 3
************************************************************
Group 18. Limits
VARMAX(C6)=1.1 ;VARMIN(C6)=0.
VARMAX(C7)=1.1 ;VARMIN(C7)=0.
VARMAX(C8)=1.1 ;VARMIN(C8)=0.
VARMAX(TEM1)=3000. ;VARMIN(TEM1)=-30.
************************************************************
Group 19. EARTH Calls To GROUND Station
NAMGRD =FLAR
GENK = T
PARSOL = F
IENUTA = 1
ISG62 = 1
SPEDAT(SET,DFLUX,CONSLP,L,F)
SPEDAT(SET,DFLUX,JWDEP,I,2)
SPEDAT(SET,OUTPUT,NOFIELD,L,T)
SPEDAT(SET,DFLUX,DFMODL,L,T)
SPEDAT(SET,DFLUX,DEPOMOD,I,4)
SPEDAT(SET,DFLUX,DENP1,R,1000.)
SPEDAT(SET,DFLUX,DIAP1,R,1.0E-05)
SPEDAT(SET,DFLUX,DENP2,R,1000.)
SPEDAT(SET,DFLUX,DIAP2,R,6.0E-06)
SPEDAT(SET,DFLUX,DENP3,R,1000.)
SPEDAT(SET,DFLUX,DIAP3,R,2.0E-06)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=0.0208 ;EX(U1)=1.545E-03
EX(V1)=2.758E-03 ;EX(W1)=4.951E-03
EX(KE)=1.535E-05 ;EX(EP)=3.061E-06
EX(C6)=0.09899 ;EX(C7)=0.1749
EX(C8)=0.3205 ;EX(EPKE)=0.
EX(V3DH)=0.01089 ;EX(VD3)=3.025E-06
EX(VD2)=2.583E-05 ;EX(VS3)=1.957E-04
EX(VS1)=4.593E-03 ;EX(VDPH)=1.0E-10
EX(DI1)=2.433E+05 ;EX(VP1)=0.08505
EX(VSTR)=1.885E-04 ;EX(TP1)=1.639E-05
EX(V2DH)=0.09299 ;EX(PRPS)=1.
EX(VABS)=7.045E-03 ;EX(YPLS)=0.4537
EX(VD1)=7.098E-05 ;EX(STRS)=3.501E-07
EX(V1DH)=0.2555 ;EX(ENUL)=1.0E-05
EX(SPH1)=0. ;EX(VS2)=1.671E-03
EX(TEM1)=20. ;EX(DEN1)=1.205
EX(EL1)=0.02831 ;EX(ENUT)=4.428E-05
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 16 ;IYMON = 15 ;IZMON = 47
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
IDISPA = 2000 ;IDISPB = 0 ;IDISPC = 0
CSG1 ='SW'
GVIEW(P,0.815082,0.368271,0.447234)
GVIEW(UP,-0.407564,-0.184146,0.894417)
GVIEW(WINDOW,976,1056)
GVIEW(DEPTH,3.0E+04)
GVIEW(VDIS,2.026717)
GVIEW(CENTRE,0.848345,0.988749,1.383435)
> DOM, SIZE, 2.000000E+00, 2.000000E+00, 2.800000E+00
> DOM, MONIT, 1.067742E+00, 8.026798E-01, 2.410216E+00
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID, AUTO, T T T
> GRID, MINCELL, 0.034674 ,0.018621 ,4.570882E-03
> GRID, RSET_X_1, -9, 1.100000E+00,G
> GRID, RSET_X_2, -3, 1.100000E+00,G
> GRID, RSET_X_3, 1, 1.000000E+00
> GRID, RSET_X_4, -3, 1.100000E+00,G
> GRID, RSET_X_5, 1, 1.000000E+00
> GRID, RSET_X_6, -2, 1.100000E+00,G
> GRID, RSET_X_7, -6, 1.100000E+00,G
> GRID, RSET_Y_1, -17, 1.100000E+00,G
> GRID, RSET_Y_2, -1, 1.100000E+00,G
> GRID, RSET_Y_3, -3, 1.100000E+00,G
> GRID, RSET_Y_4, -1, 1.100000E+00,G
> GRID, RSET_Y_5, -17, 1.100000E+00,G
> GRID, RSET_Z_1, -16, 1.100000E+00,G
> GRID, RSET_Z_2, -42, 1.100000E+00,G
> DOM, INI_AMB, YES
> DOM, INI_BUOY, YES
> OBJ, NAME, INL
> OBJ, POSITION, 8.500000E-01, 9.500000E-01, 0.000000E+00
> OBJ, SIZE, 3.000000E-01, 1.000000E-01, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, TYPE, INLET
> OBJ, DENSITY, 1.176822
> OBJ, VELOCITY, 0. ,0. ,0.1
> OBJ, TEMPERATURE, T_AMBIENT
> OBJ, INLET_C6, 1.
> OBJ, INLET_C7, 1.
> OBJ, INLET_C8, 1.
> OBJ, KE_IN, 1.000000E-03
> OBJ, EP_IN, 9.000000E-03
> OBJ, NAME, EXHAUST
> OBJ, POSITION, 9.000000E-01, 9.000000E-01, AT_END
> OBJ, SIZE, 2.000000E-01, 2.000000E-01, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, cube3t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, VOLUFLOW, -6.0E-03
> OBJ, TEMPERATURE, T_AMBIENT
> OBJ, TURB-INTENS, 5.
> OBJ, NAME, NW
> OBJ, POSITION, 0.000000E+00, AT_END, 0.000000E+00
> OBJ, SIZE, TO_END, 0.000000E+00, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, cube11
> OBJ, VISIBLE, NO
> OBJ, TYPE, PLATE
> OBJ, NAME, LW1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, TO_END, TO_END, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, cube11
> OBJ, VISIBLE, NO
> OBJ, TYPE, PLATE
> OBJ, NAME, HW1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, AT_END
> OBJ, SIZE, TO_END, TO_END, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, cube11
> OBJ, VISIBLE, NO
> OBJ, TYPE, PLATE
> OBJ, NAME, EW
> OBJ, POSITION, AT_END, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, TO_END, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, cube11
> OBJ, VISIBLE, NO
> OBJ, TYPE, PLATE
> OBJ, NAME, WW
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, TO_END, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, cube11
> OBJ, VISIBLE, NO
> OBJ, TYPE, PLATE
> OBJ, ROUGH, 0.
> OBJ, WALL_CO, LOG-LAW
> OBJ, NAME, NOPEN
> OBJ, POSITION, 7.000000E-01, 2.000000E+00, 0.000000E+00
> OBJ, SIZE, 6.000000E-01, 0.000000E+00, 3.000000E-01
> OBJ, GEOMETRY, cube12t
> OBJ, TYPE, OPENING
> OBJ, PRESSURE, P_AMBIENT
> OBJ, TEMPERATURE, T_AMBIENT
> OBJ, COEFFICIENT, 1.
> OBJ, TURBULENCE, SAME , SAME
> OBJ, NAME, SW
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, TO_END, 0.000000E+00, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, cube11
> OBJ, VISIBLE, NO
> OBJ, TYPE, PLATE
STOP