TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2018, Date 25/01/18
CPVNAM=VDI; SPPNAM=FLAIR
************************************************************
Echo DISPLAY / USE settings
DISPLAY
Library Case I406: Deposition, downpipe, Re=5.04E4
--------------------------------------------------
The case considered is 2d steady, isothermal, turbulent flow
of air and particles down a pipe with particle deposition on
the wall of the pipe. The Eulerian drift flux model is used to
simulate particle transport with deposition by means of the
3-layer model of Bin & Zhou(2006), which accounts for deposition
via gravity, Brownian and turbulent diffusion, and turbophoresis.
This case has been studied experimentally by Liu & Agrawal(1974),
who for fully-developed flow, measured deposition rates of particles
to the walls of the pipe with nominal sizes ranging from 1.4 to 21
microns in a 0.0127m pipe at air speeds of 11.85m/s and 60 m/s,
corresponding to Reynolds numbers of 9,894 & 50,000, respectively.
The present Q1 file considers the air speed of 60m/s at a temperature
of 22degC so as to replicate the experimental deposition of 13
different particles sizes onto the pipe wall in the fully-developed
region of the flow.
ENDDIS
PHOTON USE
p
up -z
view y
con C6 x 1 fi;.1;pa
vec x 1;pa
ENDUSE
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(I406: Deposition, downpipe, Re=5E4 )
************************************************************
Echo save-block settings for Group 1
save1begin
Nineteen different aerosol particle sizes with a density of
920 kg/m^3 enter a vertical down pipe duct in an air stream of
stream of 60 m/s at 22degC. The pipe is 80 diameters long, at
which location the flow is fully developed. In the simulations
the particle concentration is normalized by the inlet concentration.
The case has been studied numerically by Zhao and Wu (2006),
and experimentally by Liu & Agrawal(1974). The particle
characteristics are listed below, and cover the 13 experimental
sizes (runs 1 to 15), plus 6 additional smaller sizes( runs a
to f) so as to cover the diffusion-dominated deposition regime.
Run diameter Tr+ Vd+ Tr+ Vd+
(microns) measured predicted
a 0.00375 - - 1.37E-03 4.20E-04
b 0.025 - - 9.67E-03 1.21E-04
c 0.15 - - 0.08 4.18E-05
d 0.25 - - 0.17 3.51E-05
e 0.5 - - 0.55 4.18E-05
f 0.7 - - 1.0 7.02E-05
1 1.4 3.7 9.70E-3 3.6 2.67E-3
2 2.0 7.8 6.70E-2 7.2 0.121
3 2.7 13.6 0.11 12.8 0.209
4 3.6 22.7 0.13 22.4 0.251
5 5.4 53.3 0.14 49.7 0.321
7 5.6 58.4 0.12 53.4 0.326
8 6.0 63.2 0.12 61.2 0.336
9 7.2 110.0 0.12 10.3 0.365
10 8.4 124.0 0.12 119.0 0.372
11 10.5 193.0 0.13 185.2 0.388
12 14.0 344.0 0.11 328.0 0.399
14 16.8 526.0 0.097 471.4 0.405
15 21.0 774.0 0.089 735.1 0.409
where Tr+ (=Tr*Ust^2/enul_g) is the dimensional particle relaxation
time, Vd+ (=Vd/Ust)is the dimensionless deposition velocity, enul_g
is the kinematic laminar viscosity, Ust (=3 m/s) is the friction
velocity, Tr (=rho_p*diam_p^2*C/(18.*rho_g*enul_g) is the particle
relaxation time and C is the Cunningham slip coefficient.
References
----------.
Zhao, B. & Wu, J., "Modelling particle deposition from fully
developed flow in ventilation duct", J. Atmospheric
Environment, Vol.40, p457-466, (2006).
Liu, B.Y.H.,& Agarwal,J.K.(1974). "Experimental observation
of aerosol deposition in turbulent flow", Journal of Aerosol
Science, 5, Vol.5, p145-155, (1974).
Zhao, B. & Wu, J., "Particle deposition in indoor
environments: Analysis of influencing factors", Journal of
Hazardous Materials 147, p439-448, (2007).
Lai, C.K., Nazaroff, W.W., "Modeling indoor particle
deposition from turbulent flow onto smooth surfaces", Journal
of Aerosol Science 31, p463-476, (2000)
save1end
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,1,20,80)
* Cylindrical-polar grid
CARTES=F
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(56)=VP19 ;NAME(57)=TP19
NAME(58)=VP18 ;NAME(59)=TP18
NAME(60)=VP17 ;NAME(61)=TP17
NAME(62)=VP16 ;NAME(63)=TP16
NAME(64)=VP15 ;NAME(65)=TP15
NAME(66)=VP14 ;NAME(67)=TP14
NAME(69)=VP13 ;NAME(70)=TP13
NAME(74)=VP12 ;NAME(75)=TP12
NAME(79)=VP11 ;NAME(80)=TP11
NAME(83)=VS10 ;NAME(84)=VP10
NAME(85)=TP10 ;NAME(89)=VP9
NAME(90)=TP9 ;NAME(94)=VP8
NAME(95)=TP8 ;NAME(99)=VP7
NAME(100)=TP7 ;NAME(101)=VD7
NAME(102)=TR7 ;NAME(103)=VS6
NAME(104)=VP6 ;NAME(105)=TP6
NAME(108)=VP5 ;NAME(109)=TP5
NAME(110)=VP4 ;NAME(111)=TP4
NAME(112)=VP3 ;NAME(113)=TP3
NAME(114)=VP2 ;NAME(115)=TP2
NAME(116)=VP1 ;NAME(117)=TP1
NAME(121)=DEP5 ;NAME(125)=DEP4
NAME(129)=DEP3 ;NAME(133)=DEP2
NAME(137)=VSTR ;NAME(140)=DEP1
NAME(144)=ENUL ;NAME(145)=YPLS
NAME(146)=STRS ;NAME(147)=EPKE
NAME(148)=DEN1 ;NAME(149)=EL1
NAME(150) =ENUT
* Solved variables list
SOLVE(P1,V1,W1,C6,C7,C8,C9,C10)
SOLVE(C11,C12,C13,C14,C15,C16,C17,C18)
SOLVE(C19,C20,C21,C22,C23,C24)
* Stored variables list
STORE(ENUT,EL1,DEN1,EPKE,STRS,YPLS,ENUL,DEP1)
STORE(VSTR,DEP2,DEP3,DEP4,DEP5,TP1,VP1,TP2)
STORE(VP2,TP3,VP3,TP4,VP4,TP5,VP5,TP6)
STORE(VP6,VS6,TR7,VD7,TP7,VP7,TP8,VP8)
STORE(TP9,VP9,TP10,VP10,VS10,TP11,VP11,TP12)
STORE(VP12,TP13,VP13,TP14,VP14,TP15,VP15,TP16)
STORE(VP16,TP17,VP17,TP18,VP18,TP19,VP19)
* 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)
SOLUTN(C9,Y,Y,Y,N,N,Y)
SOLUTN(C10,Y,Y,Y,N,N,Y)
SOLUTN(C11,Y,Y,Y,N,N,Y)
SOLUTN(C12,Y,Y,Y,N,N,Y)
SOLUTN(C13,Y,Y,Y,N,N,Y)
SOLUTN(C14,Y,Y,Y,N,N,Y)
SOLUTN(C15,Y,Y,Y,N,N,Y)
SOLUTN(C16,Y,Y,Y,N,N,Y)
SOLUTN(C17,Y,Y,Y,N,N,Y)
SOLUTN(C18,Y,Y,Y,N,N,Y)
SOLUTN(C19,Y,Y,Y,N,N,Y)
SOLUTN(C20,Y,Y,Y,N,N,Y)
SOLUTN(C21,Y,Y,Y,N,N,Y)
SOLUTN(C22,Y,Y,Y,N,N,Y)
SOLUTN(C23,Y,Y,Y,N,N,Y)
SOLUTN(C24,Y,Y,Y,N,N,Y)
TURMOD(KEMODL)
************************************************************
Echo save-block settings for Group 7
save7begin
save7end
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
PRESS0 =1.01325E+05 ;TEMP0 =273.
* Domain material index is 0 signifying:
* Air at 20 deg C, 1 atm, treated as incompressible
SETPRPS(1, 0)
RHO1 =1.198
ENUL =1.524E-05
TMP1 = GRND1
TMP1A =22. ;TMP1B =0.
TMP1C =0.
DVO1DT =3.41E-03
PRT(EP)=1.314
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1)=0. ;FIINIT(W1)=60.
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(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)
COVAL(DFLUX, C9, GRND4, GRND4)
COVAL(DFLUX, C10, GRND4, GRND4)
COVAL(DFLUX, C11, GRND4, GRND4)
COVAL(DFLUX, C12, GRND4, GRND4)
COVAL(DFLUX, C13, GRND4, GRND4)
COVAL(DFLUX, C14, GRND4, GRND4)
COVAL(DFLUX, C15, GRND4, GRND4)
COVAL(DFLUX, C16, GRND4, GRND4)
COVAL(DFLUX, C17, GRND4, GRND4)
COVAL(DFLUX, C18, GRND4, GRND4)
COVAL(DFLUX, C19, GRND4, GRND4)
COVAL(DFLUX, C20, GRND4, GRND4)
COVAL(DFLUX, C21, GRND4, GRND4)
COVAL(DFLUX, C22, GRND4, GRND4)
COVAL(DFLUX, C23, GRND4, GRND4)
COVAL(DFLUX, C24, GRND4, GRND4)
BUOYA =0. ; BUOYB =0.
BUOYC =9.81
EGWF = T
************************************************************
Echo save-block settings for Group 13
save13begin
** estimate for turbulent inlet intensity
for fully developed duct flow
real(din,win,tint,kein,epin,mixl,fric,temk)
real(reyn,ust,tvis,tvisn,enugas,rhogas)
temk=22.0+temp0
din=0.0127
win=60.0
rhogas=press0/(287.0*temk)
rhogas
enugas=1.46e-6*temk**1.5/(110.+temk)/rhogas
enugas
reyn=win*din/enugas
reyn
fric=1./(1.82*log10(reyn)-1.64)**2
fric
ust=win*(fric/8.)**0.5
ust
kein=fric*win*win/4.
mixl=0.09*0.5*din
mixl
epin=0.1643*kein**1.5/mixl
kein
epin
tint=kein**0.5/win
tint
tvis=0.09*kein*kein/epin
tvis
tvisn=tvis/(ust*din)
tvisn
real(vslip,rho_p,dpart,grava,vdep,rey,tpart,tpl)
** no of phases set by menu limit
integer(npart);npart=5
** add more phases
integer(npext,nptot)
npext=14 ! no of extra phases above menu default of 5
nptot=npart+npext
nptot
** Particle diameter & density
array(pdia,real,nptot)
array(pden,real,nptot)
rho_p=920.
** additional 6 phases for full coverage
of diffusion-dominated regime
pdia(1)=0.00375;pdia(2)=0.025;pdia(3)=0.15
pdia(4)=0.25;pdia(5)=0.5;pdia(6)=0.7
pdia(7)=1.4;pdia(8)=2.0;pdia(9)=2.7
pdia(10)=3.6;pdia(11)=5.4
do ii=1,nptot
+ pden(ii)=rho_p
enddo
pdia(12)=5.6;pdia(13)=6.0;pdia(14)=7.8;pdia(15)=8.4
pdia(16)=10.5;pdia(17)=14.0;pdia(18)=16.8
pdia(19)=21.0
integer(kk)
kk=10
set densities & diameters of extra phases
DO ii=npart+1,nptot
+ kk=kk+1
+ SOLVE(C:kk:)
+ SOLUTN(C:kk:,y,y,y,p,p,p)
+ RELAX(C:kk:,FALSDT,10.0)
+ COVAL(DFLUX ,C:kk: , GRND4 , GRND4 )
+ spedat(set,dflux,diap:ii:,r,pdia(ii)*1.e-6)
+ spedat(set,dflux,denp:ii:,r,pden(ii))
+ store(tr:ii:)
+ store(vd:ii:)
+ store(tp:ii:)
+ store(vp:ii:)
+ store(vs:ii:)
ENDDO
do ii=1,nptot
dpart=pdia(ii)*1.e-6
rho_p=pden(ii)
** Stokes-flow slip velocity
grava=9.81
vslip=rho_p*dpart*dpart*grava/(18.*rhogas*enugas)
rey=vslip*dpart/enugas
ii
vdep=vslip
vslip
rey
tpart=rho_p*dpart**2/(18.*rhogas*enugas)
tpart
tpl=tpart*ust*ust/enugas
tpl
enddo
** print to the "inforout" file the wall deposition
results in the fully-developed flow region.
integer(izp)
izp=75
(make1 vs_ny is 0)
(store1 vs_ny is VSTR[1,NY,:izp:])
(print V* is vs_ny)
(make1 yp_ny is 0)
(store1 yp_ny is YPLS[1,NY,:izp:])
(print y+ is yp_ny)
do jj=1,nptot
(make1 tr_:jj: is 0)
(store1 tr_:jj: is TR:jj:[1,NY,:izp:,1])
(print Tr_:jj: is tr_:jj:)
(make1 vd_:jj: is 0)
(store1 vd_:jj: is VD:jj:[1,NY,:izp:])
(print Vd_:jj: is vd_:jj:)
(make1 tp_:jj: is 0)
(store1 tp_:jj: is TP:jj:[1,NY,:izp:])
(print Tr+_:jj: is tp_:jj:)
(make1 vp_:jj: is 0)
(store1 vp_:jj: is VP:jj:[1,NY,:izp:])
(print Vd+_:jj: is vp_:jj:)
enddo
save13end
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 140
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
LITER(P1)=50
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
RELAX(KE ,LINRLX,0.5 )
RELAX(EP ,LINRLX,0.5 )
RELAX(C6 ,FALSDT,10. )
RELAX(C11 ,FALSDT,10. )
RELAX(C12 ,FALSDT,10. )
RELAX(C13 ,FALSDT,10. )
RELAX(C14 ,FALSDT,10. )
RELAX(C15 ,FALSDT,10. )
RELAX(C16 ,FALSDT,10. )
RELAX(C17 ,FALSDT,10. )
RELAX(C18 ,FALSDT,10. )
RELAX(C19 ,FALSDT,10. )
RELAX(C20 ,FALSDT,10. )
RELAX(C21 ,FALSDT,10. )
RELAX(C22 ,FALSDT,10. )
RELAX(C23 ,FALSDT,10. )
RELAX(C24 ,FALSDT,10. )
KELIN = 3
************************************************************
Group 18. Limits
VARMAX(C6)=1. ;VARMIN(C6)=0.
************************************************************
Group 19. EARTH Calls To GROUND Station
NAMGRD =FLAR
GENK = T
PARSOL = F
ISG62 = 0
SPEDAT(SET,DFLUX,DIAP6,R,5.6)
SPEDAT(SET,DFLUX,DENP6,R,920.)
SPEDAT(SET,DFLUX,DIAP7,R,6.0E-06)
SPEDAT(SET,DFLUX,DENP7,R,920.)
SPEDAT(SET,DFLUX,DIAP8,R,7.8E-06)
SPEDAT(SET,DFLUX,DENP8,R,920.)
SPEDAT(SET,DFLUX,DIAP9,R,8.399999E-06)
SPEDAT(SET,DFLUX,DENP9,R,920.)
SPEDAT(SET,DFLUX,DIAP10,R,1.05E-05)
SPEDAT(SET,DFLUX,DENP10,R,920.)
SPEDAT(SET,DFLUX,DIAP11,R,1.4E-05)
SPEDAT(SET,DFLUX,DENP11,R,920.)
SPEDAT(SET,DFLUX,DIAP12,R,1.68E-05)
SPEDAT(SET,DFLUX,DENP12,R,920.)
SPEDAT(SET,DFLUX,DIAP13,R,2.1E-05)
SPEDAT(SET,DFLUX,DENP13,R,920.)
SPEDAT(SET,DFLUX,DIAP14,R,3.75E-09)
SPEDAT(SET,DFLUX,DENP14,R,920.)
SPEDAT(SET,DFLUX,DIAP15,R,2.5E-08)
SPEDAT(SET,DFLUX,DENP15,R,920.)
SPEDAT(SET,DFLUX,DIAP16,R,1.5E-07)
SPEDAT(SET,DFLUX,DENP16,R,920.)
SPEDAT(SET,DFLUX,DIAP17,R,2.5E-07)
SPEDAT(SET,DFLUX,DENP17,R,920.)
SPEDAT(SET,DFLUX,DIAP18,R,5.0E-07)
SPEDAT(SET,DFLUX,DENP18,R,920.)
SPEDAT(SET,DFLUX,DIAP19,R,7.0E-07)
SPEDAT(SET,DFLUX,DENP19,R,920.)
SPEDAT(SET,DFLUX,DFMODL,L,T)
SPEDAT(SET,DFLUX,DEPOMOD,I,3)
SPEDAT(SET,DFLUX,DENP1,R,920.)
SPEDAT(SET,DFLUX,DIAP1,R,3.75E-09)
SPEDAT(SET,DFLUX,DENP2,R,920.)
SPEDAT(SET,DFLUX,DIAP2,R,2.5E-08)
SPEDAT(SET,DFLUX,DENP3,R,920.)
SPEDAT(SET,DFLUX,DIAP3,R,1.5E-07)
SPEDAT(SET,DFLUX,DENP4,R,920.)
SPEDAT(SET,DFLUX,DIAP4,R,2.5E-07)
SPEDAT(SET,DFLUX,DENP5,R,920.)
SPEDAT(SET,DFLUX,DIAP5,R,5.0E-07)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=1713. ;EX(V1)=0.01308
EX(W1)=63.790001 ;EX(KE)=15.97
EX(EP)=4.871E+04 ;EX(C6)=0.997
EX(C7)=0.9991 ;EX(C8)=0.9997
EX(C9)=0.9997 ;EX(C10)=0.9997
EX(C11)=0.9994 ;EX(C12)=0.9691
EX(C13)=0.6524 ;EX(C14)=0.5959
EX(C15)=0.5791 ;EX(C16)=0.5599
EX(C17)=0.5588 ;EX(C18)=0.5568
EX(C19)=0.5513 ;EX(C20)=0.5502
EX(C21)=0.5478 ;EX(C22)=0.5458
EX(C23)=0.5451 ;EX(C24)=0.5444
EX(VP19)=0.02045 ;EX(TP19)=37.290001
EX(VP18)=0.02023 ;EX(TP18)=23.91
EX(VP17)=0.01998 ;EX(TP17)=16.639999
EX(VP16)=0.01935 ;EX(TP16)=9.396
EX(VP15)=0.0186 ;EX(TP15)=6.037
EX(VP14)=0.01828 ;EX(TP14)=5.213
EX(VP13)=0.01683 ;EX(TP13)=3.104
EX(VP12)=0.01636 ;EX(TP12)=2.709
EX(VP11)=0.01609 ;EX(TP11)=2.522
EX(VS10)=3.718E-04 ;EX(VP10)=0.01258
EX(TP10)=1.137 ;EX(VP9)=0.01047
EX(TP9)=0.6492 ;EX(VP8)=6.156E-03
EX(TP8)=0.3634 ;EX(VP7)=1.727E-04
EX(TP7)=0.1839 ;EX(VD7)=5.48E-04
EX(TR7)=6.129E-06 ;EX(VS6)=1.663E-05
EX(VP6)=3.626E-06 ;EX(TP6)=0.05086
EX(VP5)=2.123E-06 ;EX(TP5)=0.02797
EX(VP4)=1.762E-06 ;EX(TP4)=8.862E-03
EX(VP3)=2.094E-06 ;EX(TP3)=4.166E-03
EX(VP2)=6.029E-06 ;EX(TP2)=4.904E-04
EX(VP1)=2.099E-05 ;EX(TP1)=6.948E-05
EX(DEP5)=7.706E-06 ;EX(DEP4)=6.381E-06
EX(DEP3)=7.581E-06 ;EX(DEP2)=2.179E-05
EX(VSTR)=0.1511 ;EX(DEP1)=7.552E-05
EX(ENUL)=1.524E-05 ;EX(YPLS)=1.574
EX(STRS)=0.4573 ;EX(EPKE)=1980.
EX(DEN1)=1.198 ;EX(EL1)=6.605E-04
EX(ENUT)=1.269E-03
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 20 ;IZMON = 75
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
GVIEW(P,0.767594,0.440068,0.465982)
GVIEW(UP,-0.444126,0.88939,-0.108339)
GVIEW(NEARPLANE,1.524E-03)
GVIEW(VDIS,0.490557)
GVIEW(CENTRE,3.169711E-04,3.175E-03,0.51)
> DOM, SIZE, 1.000000E-01, 6.350000E-03, 1.020000E+00
> DOM, MONIT, 5.000000E-02, 6.191250E-03, 9.498749E-01
> DOM, SCALE, 1.000000E+01, 1.000000E+01, 1.000000E+00
> DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID, AUTO, T F F
> GRID, RSET_X_1, 1, 1.000000E+00
> GRID, RSET_Y_1, 20, 1.000000E+00
> GRID, RSET_Z_1, 80, 1.000000E+00
> DOM, INI_AMB, YES
> DOM, INI_BUOY, YES
> OBJ, NAME, WALLN
> OBJ, POSITION, 0.000000E+00, AT_END, 0.000000E+00
> OBJ, SIZE, TO_END, 0.000000E+00, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu10
> OBJ, VISIBLE, NO
> OBJ, TYPE, PLATE
> OBJ, NAME, INL
> 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, polcu5t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, VELOCITY, 0. ,0. ,60.
> OBJ, INLET_C6, 1.
> OBJ, INLET_C7, 1.
> OBJ, INLET_C8, 1.
> OBJ, INLET_C9, 1.
> OBJ, INLET_C10, 1.
> OBJ, INLET_C11, 1.
> OBJ, INLET_C12, 1.
> OBJ, INLET_C13, 1.
> OBJ, INLET_C14, 1.
> OBJ, INLET_C15, 1.
> OBJ, INLET_C16, 1.
> OBJ, INLET_C17, 1.
> OBJ, INLET_C18, 1.
> OBJ, INLET_C19, 1.
> OBJ, INLET_C20, 1.
> OBJ, INLET_C21, 1.
> OBJ, INLET_C22, 1.
> OBJ, INLET_C23, 1.
> OBJ, INLET_C24, 1.
> OBJ, TURB-INTENS, 5.
> OBJ, NAME, OUT
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, AT_END
> OBJ, SIZE, TO_END, TO_END, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcubet
> OBJ, TYPE, OPENING
> OBJ, PRESSURE, P_AMBIENT
> OBJ, COEFFICIENT, 1000.
> OBJ, TURBULENCE, SAME , SAME
STOP
kk=10
** set deposition boundary conditions
for additional aerosol phases
DO ii=npart+1,nptot
+ kk=kk+1
+ COVAL(PW1,C:kk:, GRND3 ,0. )
ENDDO