TALK=T;RUN(1,1)
DISPLAY
Library Case I404: Deposition, horizontal duct @ 2.2m/s
The case considered is 3d steady, isothermal, turbulent air flow
in a horizontal square ventilation duct with the transport
and deposition of particles. The Eulerian drift flux model is
used to simulate particles with deposition by means of the
3-layer deposition model of Chen & Lai (2004), which accounts
for deposition by gravity, Brownian diffusion and turbophoresis.
This case has been studied experimentally by Sippola and Nazaroff
(2004), who for fully-developed flow, measured deposition rates
of particles to the walls of the duct with nominal sizes ranging
from 1 to 16 microns in a 0.152m square duct at each of three
nominal air speeds: 2.2m/s, 5.3m/s, 9m/s, corresponding to
Reynolds numbers of 21,957, 53,894 & 87,733. This particular Q1
file employs the lowest air speed of 2.2m/s at a temperature of
22degC so as to replicate Runs 1 to 5 of the experimental series,
which consider deposition of 5 different particle sizes onto the
four walls of the duct in the fully-developed region of the flow.
ENDDIS
PHOTON USE
p
up x
view y
con C6 y 1 fi;.1;pa
vec y 1;pa
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(I404: Deposition, horizontal [email protected]/s)
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 60
XULAST =3.04
XFRAC(1)=0.016667 ;XFRAC(2)=0.033333
XFRAC(3)=0.05 ;XFRAC(4)=0.066667
XFRAC(5)=0.083333 ;XFRAC(6)=0.1
XFRAC(7)=0.116667 ;XFRAC(8)=0.133333
XFRAC(9)=0.15 ;XFRAC(10)=0.166667
XFRAC(11)=0.183333 ;XFRAC(12)=0.2
XFRAC(13)=0.216667 ;XFRAC(14)=0.233333
XFRAC(15)=0.25 ;XFRAC(16)=0.266667
XFRAC(17)=0.283333 ;XFRAC(18)=0.3
XFRAC(19)=0.316667 ;XFRAC(20)=0.333333
XFRAC(21)=0.35 ;XFRAC(22)=0.366667
XFRAC(23)=0.383333 ;XFRAC(24)=0.4
XFRAC(25)=0.416667 ;XFRAC(26)=0.433333
XFRAC(27)=0.45 ;XFRAC(28)=0.466667
XFRAC(29)=0.483333 ;XFRAC(30)=0.5
XFRAC(31)=0.516667 ;XFRAC(32)=0.533333
XFRAC(33)=0.55 ;XFRAC(34)=0.566667
XFRAC(35)=0.583333 ;XFRAC(36)=0.6
XFRAC(37)=0.616667 ;XFRAC(38)=0.633333
XFRAC(39)=0.65 ;XFRAC(40)=0.666667
XFRAC(41)=0.683333 ;XFRAC(42)=0.7
XFRAC(43)=0.716667 ;XFRAC(44)=0.733333
XFRAC(45)=0.75 ;XFRAC(46)=0.766667
XFRAC(47)=0.783333 ;XFRAC(48)=0.8
XFRAC(49)=0.816667 ;XFRAC(50)=0.833333
XFRAC(51)=0.85 ;XFRAC(52)=0.866667
XFRAC(53)=0.883333 ;XFRAC(54)=0.9
XFRAC(55)=0.916667 ;XFRAC(56)=0.933333
XFRAC(57)=0.95 ;XFRAC(58)=0.966667
XFRAC(59)=0.983333 ;XFRAC(60)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 30
YVLAST =0.152
YFRAC(1)=0.033333 ;YFRAC(2)=0.066667
YFRAC(3)=0.1 ;YFRAC(4)=0.133333
YFRAC(5)=0.166667 ;YFRAC(6)=0.2
YFRAC(7)=0.233333 ;YFRAC(8)=0.266667
YFRAC(9)=0.3 ;YFRAC(10)=0.333333
YFRAC(11)=0.366667 ;YFRAC(12)=0.4
YFRAC(13)=0.433333 ;YFRAC(14)=0.466667
YFRAC(15)=0.5 ;YFRAC(16)=0.533333
YFRAC(17)=0.566667 ;YFRAC(18)=0.6
YFRAC(19)=0.633333 ;YFRAC(20)=0.666667
YFRAC(21)=0.7 ;YFRAC(22)=0.733333
YFRAC(23)=0.766667 ;YFRAC(24)=0.8
YFRAC(25)=0.833333 ;YFRAC(26)=0.866667
YFRAC(27)=0.9 ;YFRAC(28)=0.933333
YFRAC(29)=0.966667 ;YFRAC(30)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 30
ZWLAST =0.152
ZFRAC(1)=0.033333 ;ZFRAC(2)=0.066667
ZFRAC(3)=0.1 ;ZFRAC(4)=0.133333
ZFRAC(5)=0.166667 ;ZFRAC(6)=0.2
ZFRAC(7)=0.233333 ;ZFRAC(8)=0.266667
ZFRAC(9)=0.3 ;ZFRAC(10)=0.333333
ZFRAC(11)=0.366667 ;ZFRAC(12)=0.4
ZFRAC(13)=0.433333 ;ZFRAC(14)=0.466667
ZFRAC(15)=0.5 ;ZFRAC(16)=0.533333
ZFRAC(17)=0.566667 ;ZFRAC(18)=0.6
ZFRAC(19)=0.633333 ;ZFRAC(20)=0.666667
ZFRAC(21)=0.7 ;ZFRAC(22)=0.733334
ZFRAC(23)=0.766667 ;ZFRAC(24)=0.8
ZFRAC(25)=0.833334 ;ZFRAC(26)=0.866667
ZFRAC(27)=0.9 ;ZFRAC(28)=0.933334
ZFRAC(29)=0.966667 ;ZFRAC(30)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(7)=W1
NAME(12)=KE ;NAME(13)=EP
NAME(21)=C6 ;NAME(22)=C7
NAME(23)=C8 ;NAME(24)=C9
NAME(25)=C10 ;NAME(106)=ENUN
NAME(107)=DI5 ;NAME(108)=DI4
NAME(109)=DI3 ;NAME(110)=DI2
NAME(111)=VP5 ;NAME(112)=TP5
NAME(113)=VP4 ;NAME(114)=TP4
NAME(115)=VP3 ;NAME(116)=TP3
NAME(117)=VP2 ;NAME(118)=TP2
NAME(119)=VS5 ;NAME(120)=TR5
NAME(121)=VD5 ;NAME(122)=DEP5
NAME(123)=VS4 ;NAME(124)=TR4
NAME(125)=VD4 ;NAME(126)=DEP4
NAME(127)=VS3 ;NAME(128)=TR3
NAME(129)=VP1 ;NAME(130)=TP1
NAME(131)=DI1 ;NAME(132)=VS1
NAME(133)=TR1 ;NAME(134)=VD1
NAME(135)=VD3 ;NAME(136)=VSTR
NAME(137)=DEP1 ;NAME(138)=DEP3
NAME(139)=VS2 ;NAME(140)=ENUL
NAME(141)=TR2 ;NAME(142)=VD2
NAME(143)=DEP2 ;NAME(145)=YPLS
NAME(146)=STRS ;NAME(147)=EPKE
NAME(148)=DEN1 ;NAME(149)=EL1
NAME(150)=ENUT
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(P1,Y,Y,Y,N,N,Y)
SOLUTN(U1,Y,Y,N,N,N,Y)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(W1,Y,Y,N,N,N,Y)
SOLUTN(KE,Y,Y,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
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(ENUN,Y,N,N,N,N,Y)
SOLUTN(DI5,Y,N,N,N,N,Y)
SOLUTN(DI4,Y,N,N,N,N,Y)
SOLUTN(DI3,Y,N,N,N,N,Y)
SOLUTN(DI2,Y,N,N,N,N,Y)
SOLUTN(VP5,Y,N,N,N,N,Y)
SOLUTN(TP5,Y,N,N,N,N,Y)
SOLUTN(VP4,Y,N,N,N,N,Y)
SOLUTN(TP4,Y,N,N,N,N,Y)
SOLUTN(VP3,Y,N,N,N,N,Y)
SOLUTN(TP3,Y,N,N,N,N,Y)
SOLUTN(VP2,Y,N,N,N,N,Y)
SOLUTN(TP2,Y,N,N,N,N,Y)
SOLUTN(VS5,Y,N,N,N,N,Y)
SOLUTN(TR5,Y,N,N,N,N,Y)
SOLUTN(VD5,Y,N,N,N,N,Y)
SOLUTN(DEP5,Y,N,N,N,N,Y)
SOLUTN(VS4,Y,N,N,N,N,Y)
SOLUTN(TR4,Y,N,N,N,N,Y)
SOLUTN(VD4,Y,N,N,N,N,Y)
SOLUTN(DEP4,Y,N,N,N,N,Y)
SOLUTN(VS3,Y,N,N,N,N,Y)
SOLUTN(TR3,Y,N,N,N,N,Y)
SOLUTN(VP1,Y,N,N,N,N,Y)
SOLUTN(TP1,Y,N,N,N,N,Y)
SOLUTN(DI1,Y,N,N,N,N,Y)
SOLUTN(VS1,Y,N,N,N,N,Y)
SOLUTN(TR1,Y,N,N,N,N,Y)
SOLUTN(VD1,Y,N,N,N,N,Y)
SOLUTN(VD3,Y,N,N,N,N,Y)
SOLUTN(VSTR,Y,N,N,N,N,Y)
SOLUTN(DEP1,Y,N,N,N,N,Y)
SOLUTN(DEP3,Y,N,N,N,N,Y)
SOLUTN(VS2,Y,N,N,N,N,Y)
SOLUTN(ENUL,Y,N,N,N,N,Y)
SOLUTN(TR2,Y,N,N,N,N,Y)
SOLUTN(VD2,Y,N,N,N,N,Y)
SOLUTN(DEP2,Y,N,N,N,N,Y)
SOLUTN(YPLS,Y,N,N,N,N,Y)
SOLUTN(STRS,Y,N,N,N,N,Y)
SOLUTN(EPKE,Y,N,N,N,N,Y)
SOLUTN(DEN1,Y,N,N,N,N,Y)
SOLUTN(EL1,Y,N,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
DEN1 = 148
VISL = 140
VIST = 150
LEN1 = 149
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
TERMS(P1,Y,Y,Y,N,Y,Y)
TERMS(U1,Y,Y,Y,Y,Y,Y)
TERMS(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(KE,N,Y,Y,Y,Y,N)
TERMS(EP,N,Y,Y,Y,Y,N)
TERMS(C6,N,Y,Y,Y,N,Y)
TERMS(C7,N,Y,Y,Y,Y,Y)
TERMS(C8,N,Y,Y,Y,N,Y)
TERMS(C9,N,Y,Y,Y,Y,Y)
TERMS(C10,N,Y,Y,Y,N,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWRH1 = T
NEWENL = T
NEWENT = T
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 = GRND5 ;TMP1 = GRND1 ;EL1 = GRND4
TSURR =0. ;TEMP0 =273. ;PRESS0 =1.01325E+05
DVO1DT =3.33E-03 ;DRH1DP = GRND5
RHO1A =0. ;RHO1B =3.4843E-03 ;RHO1C =0.71429
TMP1A =22. ;TMP1B =0. ;TMP1C =0.
TMP2A =0. ;TMP2B =0. ;TMP2C =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0. ;EL1B =0. ;EL1C =0.
ENUL = GRND6 ;ENUT = GRND3
ENULA =1.46E-06 ;ENULB =110. ;ENULC =8.0974E-11 ;ENULD =0.
ENULE =0. ;ENULF =0. ;ENULG =0.
IENULA = 0 ;IENULB = 0
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(W1)=1. ;PRNDTL(KE)=1.
PRNDTL(EP)=1. ;PRNDTL(C6)=1.
PRNDTL(C7)=1. ;PRNDTL(C8)=1.
PRNDTL(C9)=1. ;PRNDTL(C10)=1.
PRT(U1)=1. ;PRT(V1)=1.
PRT(W1)=1. ;PRT(KE)=1.
PRT(EP)=1.314 ;PRT(C6)=1.
PRT(C7)=1. ;PRT(C8)=1.
PRT(C9)=1. ;PRT(C10)=1.
CP1 =1004. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=0. ;FIINIT(U1)=2.2
FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.0E-10
FIINIT(KE)=0.02238 ;FIINIT(EP)=0.07239
FIINIT(C6)=0. ;FIINIT(C7)=0.
FIINIT(C8)=0. ;FIINIT(C9)=0.
FIINIT(C10)=0. ;FIINIT(ENUN)=1.0E-10
FIINIT(DI5)=1.0E-10 ;FIINIT(DI4)=1.0E-10
FIINIT(DI3)=1.0E-10 ;FIINIT(DI2)=1.0E-10
FIINIT(VP5)=1.0E-10 ;FIINIT(TP5)=1.0E-10
FIINIT(VP4)=1.0E-10 ;FIINIT(TP4)=1.0E-10
FIINIT(VP3)=1.0E-10 ;FIINIT(TP3)=1.0E-10
FIINIT(VP2)=1.0E-10 ;FIINIT(TP2)=1.0E-10
FIINIT(VS5)=1.0E-10 ;FIINIT(TR5)=1.0E-10
FIINIT(VD5)=1.0E-10 ;FIINIT(DEP5)=1.0E-10
FIINIT(VS4)=1.0E-10 ;FIINIT(TR4)=1.0E-10
FIINIT(VD4)=1.0E-10 ;FIINIT(DEP4)=1.0E-10
FIINIT(VS3)=1.0E-10 ;FIINIT(TR3)=1.0E-10
FIINIT(VP1)=1.0E-10 ;FIINIT(TP1)=1.0E-10
FIINIT(DI1)=1.0E-10 ;FIINIT(VS1)=1.0E-10
FIINIT(TR1)=1.0E-10 ;FIINIT(VD1)=1.0E-10
FIINIT(VD3)=1.0E-10 ;FIINIT(VSTR)=1.0E-10
FIINIT(DEP1)=1.0E-10 ;FIINIT(DEP3)=1.0E-10
FIINIT(VS2)=1.0E-10 ;FIINIT(ENUL)=1.0E-10
FIINIT(TR2)=1.0E-10 ;FIINIT(VD2)=1.0E-10
FIINIT(DEP2)=1.0E-10 ;FIINIT(YPLS)=1.0E-10
FIINIT(STRS)=1.0E-10 ;FIINIT(EPKE)=1.0E-10
FIINIT(DEN1)=1.0E-10 ;FIINIT(EL1)=1.0E-10
FIINIT(ENUT)=1.0E-10
No PATCHes yet used for this Group
INIADD = F
FSWEEP = 1
NAMFI =CHAM
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources
PATCH(DFLUX ,CELL , 1, 60, 1, 30, 1, 30, 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 )
PATCH(KESOURCE,PHASEM, 1, 60, 1, 30, 1, 30, 1, 1)
COVAL(KESOURCE,KE , GRND4 , GRND4 )
COVAL(KESOURCE,EP , GRND4 , GRND4 )
PATCH(OB1 ,WEST , 1, 1, 1, 30, 1, 30, 1, 1)
COVAL(OB1 ,P1 , FIXFLU ,2.632891 )
COVAL(OB1 ,U1 ,0. ,2.2 )
COVAL(OB1 ,V1 ,0. ,0. )
COVAL(OB1 ,W1 ,0. ,0. )
COVAL(OB1 ,KE ,0. ,0.02238 )
COVAL(OB1 ,EP ,0. ,0.07239 )
COVAL(OB1 ,C6 ,0. ,1. )
COVAL(OB1 ,C7 ,0. ,1. )
COVAL(OB1 ,C8 ,0. ,1. )
COVAL(OB1 ,C9 ,0. ,1. )
COVAL(OB1 ,C10 ,0. ,1. )
PATCH(OB2 ,EAST , 60, 60, 1, 30, 1, 30, 1, 1)
COVAL(OB2 ,P1 ,1000. ,0. )
COVAL(OB2 ,U1 ,0. ,0. )
COVAL(OB2 ,V1 ,0. ,0. )
COVAL(OB2 ,W1 ,0. ,0. )
COVAL(OB2 ,KE ,0. , SAME )
COVAL(OB2 ,EP ,0. , SAME )
PATCH(OB3 ,HWALL , 1, 60, 1, 30, 30, 30, 1, 1)
COVAL(OB3 ,U1 , GRND2 ,0. )
COVAL(OB3 ,V1 , GRND2 ,0. )
COVAL(OB3 ,KE , GRND2 , GRND2 )
COVAL(OB3 ,EP , GRND2 , GRND2 )
PATCH(PW3 ,HIGH , 1, 60, 1, 30, 30, 30, 1, 1)
COVAL(PW3 ,C6 , GRND3 ,0. )
COVAL(PW3 ,C7 , GRND3 ,0. )
COVAL(PW3 ,C8 , GRND3 ,0. )
COVAL(PW3 ,C9 , GRND3 ,0. )
COVAL(PW3 ,C10 , GRND3 ,0. )
PATCH(OB4 ,SWALL , 1, 60, 1, 1, 1, 30, 1, 1)
COVAL(OB4 ,U1 , GRND2 ,0. )
COVAL(OB4 ,W1 , GRND2 ,0. )
COVAL(OB4 ,KE , GRND2 , GRND2 )
COVAL(OB4 ,EP , GRND2 , GRND2 )
PATCH(PW4 ,SOUTH , 1, 60, 1, 1, 1, 30, 1, 1)
COVAL(PW4 ,C6 , GRND3 ,0. )
COVAL(PW4 ,C7 , GRND3 ,0. )
COVAL(PW4 ,C8 , GRND3 ,0. )
COVAL(PW4 ,C9 , GRND3 ,0. )
COVAL(PW4 ,C10 , GRND3 ,0. )
PATCH(OB5 ,NWALL , 1, 60, 30, 30, 1, 30, 1, 1)
COVAL(OB5 ,U1 , GRND2 ,0. )
COVAL(OB5 ,W1 , GRND2 ,0. )
COVAL(OB5 ,KE , GRND2 , GRND2 )
COVAL(OB5 ,EP , GRND2 , GRND2 )
PATCH(PW5 ,NORTH , 1, 60, 30, 30, 1, 30, 1, 1)
COVAL(PW5 ,C6 , GRND3 ,0. )
COVAL(PW5 ,C7 , GRND3 ,0. )
COVAL(PW5 ,C8 , GRND3 ,0. )
COVAL(PW5 ,C9 , GRND3 ,0. )
COVAL(PW5 ,C10 , GRND3 ,0. )
PATCH(OB6 ,LWALL , 1, 60, 1, 30, 1, 1, 1, 1)
COVAL(OB6 ,U1 , GRND2 ,0. )
COVAL(OB6 ,V1 , GRND2 ,0. )
COVAL(OB6 ,KE , GRND2 , GRND2 )
COVAL(OB6 ,EP , GRND2 , GRND2 )
PATCH(PW6 ,LOW , 1, 60, 1, 30, 1, 1, 1, 1)
COVAL(PW6 ,C6 , GRND3 ,0. )
COVAL(PW6 ,C7 , GRND3 ,0. )
COVAL(PW6 ,C8 , GRND3 ,0. )
COVAL(PW6 ,C9 , GRND3 ,0. )
COVAL(PW6 ,C10 , GRND3 ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
BUOYA =0. ; BUOYB =0. ; BUOYC =-9.81
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 500 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(P1)=200 ;LITER(U1)=10
LITER(V1)=10 ;LITER(W1)=10
LITER(KE)=20 ;LITER(EP)=20
LITER(C6)=20 ;LITER(C7)=20
LITER(C8)=20 ;LITER(C9)=20
LITER(C10)=20
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03
ENDIT(C6)=1.0E-03 ;ENDIT(C7)=1.0E-03
ENDIT(C8)=1.0E-03 ;ENDIT(C9)=1.0E-03
ENDIT(C10)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(U1,FALSDT,1.)
RELAX(V1,FALSDT,1.)
RELAX(W1,FALSDT,1.)
RELAX(KE,LINRLX,0.5)
RELAX(EP,LINRLX,0.5)
RELAX(C6,LINRLX,0.5)
RELAX(C7,LINRLX,0.5)
RELAX(C8,LINRLX,0.5)
RELAX(C9,LINRLX,0.5)
RELAX(C10,LINRLX,0.5)
RELAX(ENUN,LINRLX,1.)
RELAX(DI5,LINRLX,1.)
RELAX(DI4,LINRLX,1.)
RELAX(DI3,LINRLX,1.)
RELAX(DI2,LINRLX,1.)
RELAX(VP5,LINRLX,1.)
RELAX(TP5,LINRLX,1.)
RELAX(VP4,LINRLX,1.)
RELAX(TP4,LINRLX,1.)
RELAX(VP3,LINRLX,1.)
RELAX(TP3,LINRLX,1.)
RELAX(VP2,LINRLX,1.)
RELAX(TP2,LINRLX,1.)
RELAX(VS5,LINRLX,1.)
RELAX(TR5,LINRLX,1.)
RELAX(VD5,LINRLX,1.)
RELAX(DEP5,LINRLX,1.)
RELAX(VS4,LINRLX,1.)
RELAX(TR4,LINRLX,1.)
RELAX(VD4,LINRLX,1.)
RELAX(DEP4,LINRLX,1.)
RELAX(VS3,LINRLX,1.)
RELAX(TR3,LINRLX,1.)
RELAX(VP1,LINRLX,1.)
RELAX(TP1,LINRLX,1.)
RELAX(DI1,LINRLX,1.)
RELAX(VS1,LINRLX,1.)
RELAX(TR1,LINRLX,1.)
RELAX(VD1,LINRLX,1.)
RELAX(VD3,LINRLX,1.)
RELAX(VSTR,LINRLX,1.)
RELAX(DEP1,LINRLX,1.)
RELAX(DEP3,LINRLX,1.)
RELAX(VS2,LINRLX,1.)
RELAX(ENUL,LINRLX,1.)
RELAX(TR2,LINRLX,1.)
RELAX(VD2,LINRLX,1.)
RELAX(DEP2,LINRLX,1.)
RELAX(YPLS,LINRLX,1.)
RELAX(STRS,LINRLX,1.)
RELAX(EPKE,LINRLX,1.)
RELAX(DEN1,LINRLX,1.)
RELAX(EL1,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
KELIN = 3
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
VARMAX(C6)=1. ;VARMIN(C6)=0.
VARMAX(C7)=1. ;VARMIN(C7)=0.
VARMAX(C8)=1. ;VARMIN(C8)=0.
VARMAX(C9)=1. ;VARMIN(C9)=0.
VARMAX(C10)=1. ;VARMIN(C10)=0.
VARMAX(ENUN)=1.0E+10 ;VARMIN(ENUN)=-1.0E+10
VARMAX(DI5)=1.0E+10 ;VARMIN(DI5)=-1.0E+10
VARMAX(DI4)=1.0E+10 ;VARMIN(DI4)=-1.0E+10
VARMAX(DI3)=1.0E+10 ;VARMIN(DI3)=-1.0E+10
VARMAX(DI2)=1.0E+10 ;VARMIN(DI2)=-1.0E+10
VARMAX(VP5)=1.0E+10 ;VARMIN(VP5)=-1.0E+10
VARMAX(TP5)=1.0E+10 ;VARMIN(TP5)=-1.0E+10
VARMAX(VP4)=1.0E+10 ;VARMIN(VP4)=-1.0E+10
VARMAX(TP4)=1.0E+10 ;VARMIN(TP4)=-1.0E+10
VARMAX(VP3)=1.0E+10 ;VARMIN(VP3)=-1.0E+10
VARMAX(TP3)=1.0E+10 ;VARMIN(TP3)=-1.0E+10
VARMAX(VP2)=1.0E+10 ;VARMIN(VP2)=-1.0E+10
VARMAX(TP2)=1.0E+10 ;VARMIN(TP2)=-1.0E+10
VARMAX(VS5)=1.0E+10 ;VARMIN(VS5)=-1.0E+10
VARMAX(TR5)=1.0E+10 ;VARMIN(TR5)=-1.0E+10
VARMAX(VD5)=1.0E+10 ;VARMIN(VD5)=-1.0E+10
VARMAX(DEP5)=1.0E+10 ;VARMIN(DEP5)=-1.0E+10
VARMAX(VS4)=1.0E+10 ;VARMIN(VS4)=-1.0E+10
VARMAX(TR4)=1.0E+10 ;VARMIN(TR4)=-1.0E+10
VARMAX(VD4)=1.0E+10 ;VARMIN(VD4)=-1.0E+10
VARMAX(DEP4)=1.0E+10 ;VARMIN(DEP4)=-1.0E+10
VARMAX(VS3)=1.0E+10 ;VARMIN(VS3)=-1.0E+10
VARMAX(TR3)=1.0E+10 ;VARMIN(TR3)=-1.0E+10
VARMAX(VP1)=1.0E+10 ;VARMIN(VP1)=-1.0E+10
VARMAX(TP1)=1.0E+10 ;VARMIN(TP1)=-1.0E+10
VARMAX(DI1)=1.0E+10 ;VARMIN(DI1)=-1.0E+10
VARMAX(VS1)=1.0E+10 ;VARMIN(VS1)=-1.0E+10
VARMAX(TR1)=1.0E+10 ;VARMIN(TR1)=-1.0E+10
VARMAX(VD1)=1.0E+10 ;VARMIN(VD1)=-1.0E+10
VARMAX(VD3)=1.0E+10 ;VARMIN(VD3)=-1.0E+10
VARMAX(VSTR)=1.0E+10 ;VARMIN(VSTR)=-1.0E+10
VARMAX(DEP1)=1.0E+10 ;VARMIN(DEP1)=-1.0E+10
VARMAX(DEP3)=1.0E+10 ;VARMIN(DEP3)=-1.0E+10
VARMAX(VS2)=1.0E+10 ;VARMIN(VS2)=-1.0E+10
VARMAX(ENUL)=1.0E+10 ;VARMIN(ENUL)=-1.0E+10
VARMAX(TR2)=1.0E+10 ;VARMIN(TR2)=-1.0E+10
VARMAX(VD2)=1.0E+10 ;VARMIN(VD2)=-1.0E+10
VARMAX(DEP2)=1.0E+10 ;VARMIN(DEP2)=-1.0E+10
VARMAX(YPLS)=1.0E+10 ;VARMIN(YPLS)=-1.0E+10
VARMAX(STRS)=1.0E+10 ;VARMIN(STRS)=-1.0E+10
VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
VARMAX(EL1)=1.0E+10 ;VARMIN(EL1)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
NAMGRD =FLAR
GENK = T
PARSOL = F
ISG62 = 0
SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,2)
SPEDAT(SET,DFLUX,DFMODL,L,T)
SPEDAT(SET,DFLUX,DEPOMOD,I,3)
SPEDAT(SET,DFLUX,DENP1,R,1400.)
SPEDAT(SET,DFLUX,DIAP1,R,1.0E-06)
SPEDAT(SET,DFLUX,DENP2,R,1200.)
SPEDAT(SET,DFLUX,DIAP2,R,2.8E-06)
SPEDAT(SET,DFLUX,DENP3,R,1200.)
SPEDAT(SET,DFLUX,DIAP3,R,5.2E-06)
SPEDAT(SET,DFLUX,DENP4,R,1000.)
SPEDAT(SET,DFLUX,DIAP4,R,9.1E-06)
SPEDAT(SET,DFLUX,DENP5,R,1000.)
SPEDAT(SET,DFLUX,DIAP5,R,1.6E-05)
SPEDAT(SET,STORED,ENUN,C,=ENUT/(0.152*0.43)!ZSLFIN)
SPEDAT(SET,MAKE,VS_NX,C,=0)
SPEDAT(SET,STORED,VS_NX,C,=VSTR[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,V*,C,=VS_NX)
SPEDAT(SET,MAKE,YP_NX,C,=0)
SPEDAT(SET,STORED,YP_NX,C,=YPLS[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,Y+,C,=YP_NX)
SPEDAT(SET,MAKE,TR_1,C,=0)
SPEDAT(SET,STORED,TR_1,C,=TR1[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,TR_1,C,=TR_1)
SPEDAT(SET,MAKE,VD_1,C,=0)
SPEDAT(SET,STORED,VD_1,C,=VD1[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,VD_1,C,=VD_1)
SPEDAT(SET,MAKE,TP_1,C,=0)
SPEDAT(SET,STORED,TP_1,C,=TP1[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_1,C,=TP_1)
SPEDAT(SET,MAKE,VP_1,C,=0)
SPEDAT(SET,STORED,VP_1,C,=VP1[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_1,C,=VP_1)
SPEDAT(SET,MAKE,TR_2,C,=0)
SPEDAT(SET,STORED,TR_2,C,=TR2[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,TR_2,C,=TR_2)
SPEDAT(SET,MAKE,VD_2,C,=0)
SPEDAT(SET,STORED,VD_2,C,=VD2[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,VD_2,C,=VD_2)
SPEDAT(SET,MAKE,TP_2,C,=0)
SPEDAT(SET,STORED,TP_2,C,=TP2[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_2,C,=TP_2)
SPEDAT(SET,MAKE,VP_2,C,=0)
SPEDAT(SET,STORED,VP_2,C,=VP2[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_2,C,=VP_2)
SPEDAT(SET,MAKE,TR_3,C,=0)
SPEDAT(SET,STORED,TR_3,C,=TR3[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,TR_3,C,=TR_3)
SPEDAT(SET,MAKE,VD_3,C,=0)
SPEDAT(SET,STORED,VD_3,C,=VD3[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,VD_3,C,=VD_3)
SPEDAT(SET,MAKE,TP_3,C,=0)
SPEDAT(SET,STORED,TP_3,C,=TP3[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_3,C,=TP_3)
SPEDAT(SET,MAKE,VP_3,C,=0)
SPEDAT(SET,STORED,VP_3,C,=VP3[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_3,C,=VP_3)
SPEDAT(SET,MAKE,TR_4,C,=0)
SPEDAT(SET,STORED,TR_4,C,=TR4[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,TR_4,C,=TR_4)
SPEDAT(SET,MAKE,VD_4,C,=0)
SPEDAT(SET,STORED,VD_4,C,=VD4[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,VD_4,C,=VD_4)
SPEDAT(SET,MAKE,TP_4,C,=0)
SPEDAT(SET,STORED,TP_4,C,=TP4[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_4,C,=TP_4)
SPEDAT(SET,MAKE,VP_4,C,=0)
SPEDAT(SET,STORED,VP_4,C,=VP4[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_4,C,=VP_4)
SPEDAT(SET,MAKE,TR_5,C,=0)
SPEDAT(SET,STORED,TR_5,C,=TR5[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,TR_5,C,=TR_5)
SPEDAT(SET,MAKE,VD_5,C,=0)
SPEDAT(SET,STORED,VD_5,C,=VD5[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,VD_5,C,=VD_5)
SPEDAT(SET,MAKE,TP_5,C,=0)
SPEDAT(SET,STORED,TP_5,C,=TP5[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_5,C,=TP_5)
SPEDAT(SET,MAKE,VP_5,C,=0)
SPEDAT(SET,STORED,VP_5,C,=VP5[60&15&1]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_5,C,=VP_5)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,OBJNAM,!OB1,C,INL)
SPEDAT(SET,OBJTYP,!OB1,C,INLET)
SPEDAT(SET,OBJNAM,!OB2,C,OUTL)
SPEDAT(SET,OBJTYP,!OB2,C,OPENING)
SPEDAT(SET,ARATIO,!OB2,R,1.)
SPEDAT(SET,OBJNAM,!OB3,C,HWAL)
SPEDAT(SET,OBJTYP,!OB3,C,PLATE)
SPEDAT(SET,OBJNAM,!PW3,C,HWAL)
SPEDAT(SET,OBJTYP,!PW3,C,USER_DEFINED)
SPEDAT(SET,OBJNAM,!OB4,C,SWAL)
SPEDAT(SET,OBJTYP,!OB4,C,PLATE)
SPEDAT(SET,OBJNAM,!PW4,C,SWAL)
SPEDAT(SET,OBJTYP,!PW4,C,USER_DEFINED)
SPEDAT(SET,OBJNAM,!OB5,C,NWAL)
SPEDAT(SET,OBJTYP,!OB5,C,PLATE)
SPEDAT(SET,OBJNAM,!PW5,C,NWAL)
SPEDAT(SET,OBJTYP,!PW5,C,USER_DEFINED)
SPEDAT(SET,OBJNAM,!OB6,C,LWAL)
SPEDAT(SET,OBJTYP,!OB6,C,PLATE)
SPEDAT(SET,OBJNAM,!PW6,C,LWAL)
SPEDAT(SET,OBJTYP,!PW6,C,USER_DEFINED)
SPEDAT(SET,FACETDAT,NUMOBJ,I,6)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=0.7165 ;EX(U1)=2.2
EX(V1)=1.884E-03 ;EX(W1)=1.884E-03
EX(KE)=0.02411 ;EX(EP)=0.3848
EX(C6)=0.9999 ;EX(C7)=0.9996
EX(C8)=0.999 ;EX(C9)=0.9978
EX(C10)=0.9584 ;EX(ENUN)=5.35E-03
EX(DI5)=3617. ;EX(DI4)=1.077E+04
EX(DI3)=1.038E+10 ;EX(DI2)=1.253E+11
EX(VP5)=3.239E-03 ;EX(TP5)=0.101
EX(VP4)=8.367E-04 ;EX(TP4)=0.03294
EX(VP3)=2.998E-04 ;EX(TP3)=0.01308
EX(VP2)=8.493E-05 ;EX(TP2)=3.894E-03
EX(VS5)=7.688E-03 ;EX(TR5)=7.837E-04
EX(VD5)=3.971E-04 ;EX(DEP5)=3.938E-04
EX(VS4)=2.506E-03 ;EX(TR4)=2.555E-04
EX(VD4)=1.013E-04 ;EX(DEP4)=1.177E-04
EX(VS3)=9.955E-04 ;EX(TR3)=1.015E-04
EX(VP1)=1.416E-05 ;EX(TP1)=6.379E-04
EX(DI1)=3.808E+21 ;EX(VS1)=4.816E-05
EX(TR1)=4.909E-06 ;EX(VD1)=1.7E-06
EX(VD3)=3.609E-05 ;EX(VSTR)=0.01587
EX(DEP1)=2.028E-06 ;EX(DEP3)=4.298E-05
EX(VS2)=2.963E-04 ;EX(ENUL)=1.526E-05
EX(TR2)=3.02E-05 ;EX(VD2)=1.019E-05
EX(DEP2)=1.219E-05 ;EX(YPLS)=2.619
EX(STRS)=1.98E-03 ;EX(EPKE)=9.658
EX(DEN1)=1.181 ;EX(EL1)=5.216E-03
EX(ENUT)=3.497E-04
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
OUTPUT(P1,N,N,Y,Y,Y,Y)
OUTPUT(U1,N,N,Y,Y,Y,Y)
OUTPUT(V1,N,N,Y,Y,Y,Y)
OUTPUT(W1,N,N,Y,Y,Y,Y)
OUTPUT(KE,N,N,Y,Y,Y,Y)
OUTPUT(EP,N,N,Y,Y,Y,Y)
OUTPUT(C6,N,N,Y,Y,Y,Y)
OUTPUT(C7,N,N,Y,Y,Y,Y)
OUTPUT(C8,N,N,Y,Y,Y,Y)
OUTPUT(C9,N,N,Y,Y,Y,Y)
OUTPUT(C10,N,N,Y,Y,Y,Y)
OUTPUT(ENUN,N,N,Y,N,N,N)
OUTPUT(DI5,N,N,Y,N,N,N)
OUTPUT(DI4,N,N,Y,N,N,N)
OUTPUT(DI3,N,N,Y,N,N,N)
OUTPUT(DI2,N,N,Y,N,N,N)
OUTPUT(VP5,N,N,Y,N,N,N)
OUTPUT(TP5,N,N,Y,N,N,N)
OUTPUT(VP4,N,N,Y,N,N,N)
OUTPUT(TP4,N,N,Y,N,N,N)
OUTPUT(VP3,N,N,Y,N,N,N)
OUTPUT(TP3,N,N,Y,N,N,N)
OUTPUT(VP2,N,N,Y,N,N,N)
OUTPUT(TP2,N,N,Y,N,N,N)
OUTPUT(VS5,N,N,Y,N,N,N)
OUTPUT(TR5,N,N,Y,N,N,N)
OUTPUT(VD5,N,N,Y,N,N,N)
OUTPUT(DEP5,N,N,Y,N,N,N)
OUTPUT(VS4,N,N,Y,N,N,N)
OUTPUT(TR4,N,N,Y,N,N,N)
OUTPUT(VD4,N,N,Y,N,N,N)
OUTPUT(DEP4,N,N,Y,N,N,N)
OUTPUT(VS3,N,N,Y,N,N,N)
OUTPUT(TR3,N,N,Y,N,N,N)
OUTPUT(VP1,N,N,Y,N,N,N)
OUTPUT(TP1,N,N,Y,N,N,N)
OUTPUT(DI1,N,N,Y,N,N,N)
OUTPUT(VS1,N,N,Y,N,N,N)
OUTPUT(TR1,N,N,Y,N,N,N)
OUTPUT(VD1,N,N,Y,N,N,N)
OUTPUT(VD3,N,N,Y,N,N,N)
OUTPUT(VSTR,N,N,Y,N,N,N)
OUTPUT(DEP1,N,N,Y,N,N,N)
OUTPUT(DEP3,N,N,Y,N,N,N)
OUTPUT(VS2,N,N,Y,N,N,N)
OUTPUT(ENUL,N,N,Y,N,N,N)
OUTPUT(TR2,N,N,Y,N,N,N)
OUTPUT(VD2,N,N,Y,N,N,N)
OUTPUT(DEP2,N,N,Y,N,N,N)
OUTPUT(YPLS,N,N,Y,N,Y,Y)
OUTPUT(STRS,N,N,Y,N,Y,Y)
OUTPUT(EPKE,N,N,Y,N,N,N)
OUTPUT(DEN1,N,N,Y,N,N,N)
OUTPUT(EL1,N,N,Y,N,N,N)
OUTPUT(ENUT,N,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 42 ;IYMON = 15 ;IZMON = 1
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
XZPR = T ;YZPR = F
IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP