TALK=T;RUN(1,1)
DISPLAY
Library Case I403: Deposition, horizontal duct @ 5.3m/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 highest air speed of 5.3m/s at a temperature of
23degC so as to replicate Runs 7 to 11 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(I403: 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.023171 ;YFRAC(2)=0.047501
YFRAC(3)=0.073047 ;YFRAC(4)=0.099871
YFRAC(5)=0.128035 ;YFRAC(6)=0.157608
YFRAC(7)=0.18866 ;YFRAC(8)=0.221264
YFRAC(9)=0.255498 ;YFRAC(10)=0.291444
YFRAC(11)=0.329188 ;YFRAC(12)=0.368818
YFRAC(13)=0.41043 ;YFRAC(14)=0.454123
YFRAC(15)=0.5 ;YFRAC(16)=0.545877
YFRAC(17)=0.58957 ;YFRAC(18)=0.631182
YFRAC(19)=0.670812 ;YFRAC(20)=0.708556
YFRAC(21)=0.744502 ;YFRAC(22)=0.778736
YFRAC(23)=0.81134 ;YFRAC(24)=0.842392
YFRAC(25)=0.871965 ;YFRAC(26)=0.90013
YFRAC(27)=0.926953 ;YFRAC(28)=0.952499
YFRAC(29)=0.976829 ;YFRAC(30)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 30
ZWLAST =0.152
ZFRAC(1)=0.023171 ;ZFRAC(2)=0.047501
ZFRAC(3)=0.073047 ;ZFRAC(4)=0.099871
ZFRAC(5)=0.128035 ;ZFRAC(6)=0.157608
ZFRAC(7)=0.18866 ;ZFRAC(8)=0.221264
ZFRAC(9)=0.255498 ;ZFRAC(10)=0.291444
ZFRAC(11)=0.329188 ;ZFRAC(12)=0.368818
ZFRAC(13)=0.41043 ;ZFRAC(14)=0.454123
ZFRAC(15)=0.5 ;ZFRAC(16)=0.545877
ZFRAC(17)=0.58957 ;ZFRAC(18)=0.631182
ZFRAC(19)=0.670812 ;ZFRAC(20)=0.708556
ZFRAC(21)=0.744502 ;ZFRAC(22)=0.778736
ZFRAC(23)=0.81134 ;ZFRAC(24)=0.842392
ZFRAC(25)=0.871965 ;ZFRAC(26)=0.90013
ZFRAC(27)=0.926953 ;ZFRAC(28)=0.952499
ZFRAC(29)=0.976829 ;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 =23. ;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)=5.3
FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.0E-10
FIINIT(KE)=0.129888 ;FIINIT(EP)=1.012128
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 ,6.321445 )
COVAL(OB1 ,U1 ,0. ,5.3 )
COVAL(OB1 ,V1 ,0. ,0. )
COVAL(OB1 ,W1 ,0. ,0. )
COVAL(OB1 ,KE ,0. ,0.129888 )
COVAL(OB1 ,EP ,0. ,1.012128 )
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,3.1E-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.8E-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)=3.43 ;EX(U1)=5.122
EX(V1)=4.131E-03 ;EX(W1)=4.132E-03
EX(KE)=0.1288 ;EX(EP)=5.773
EX(C6)=0.9999 ;EX(C7)=0.9996
EX(C8)=0.999 ;EX(C9)=0.9675
EX(C10)=0.9097 ;EX(ENUN)=0.01123
EX(DI5)=67.860001 ;EX(DI4)=1212.
EX(DI3)=5772. ;EX(DI2)=3.673E+09
EX(VP5)=6.433E-03 ;EX(TP5)=0.4688
EX(VP4)=1.509E-03 ;EX(TP4)=0.177
EX(VP3)=2.169E-04 ;EX(TP3)=0.06069
EX(VP2)=5.702E-05 ;EX(TP2)=0.02202
EX(VS5)=7.688E-03 ;EX(TR5)=7.837E-04
EX(VD5)=1.725E-03 ;EX(DEP5)=1.011E-03
EX(VS4)=2.903E-03 ;EX(TR4)=2.959E-04
EX(VD4)=4.076E-04 ;EX(DEP4)=3.803E-04
EX(VS3)=9.955E-04 ;EX(TR3)=1.015E-04
EX(VP1)=7.165E-06 ;EX(TP1)=2.959E-03
EX(DI1)=4.563E+10 ;EX(VS1)=4.816E-05
EX(TR1)=4.909E-06 ;EX(VD1)=1.849E-06
EX(VD3)=5.692E-05 ;EX(VSTR)=0.03374
EX(DEP1)=2.206E-06 ;EX(DEP3)=6.647E-05
EX(VS2)=3.584E-04 ;EX(ENUL)=1.563E-05
EX(TR2)=3.653E-05 ;EX(VD2)=1.477E-05
EX(DEP2)=1.758E-05 ;EX(YPLS)=3.917
EX(STRS)=9.183E-03 ;EX(EPKE)=28.889999
EX(DEN1)=1.181 ;EX(EL1)=4.675E-03
EX(ENUT)=7.338E-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