TALK=T;RUN(1,1)
DISPLAY
Library Case I402: Deposition, horizontal duct @ 9m/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 9m/s at a temperature of
26degC so as to replicate Runs 12 to 16 of the experimental serie$
s,
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(I402: Deposition, horizontal duct @9m/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.018415 ;YFRAC(2)=0.038303
YFRAC(3)=0.059782 ;YFRAC(4)=0.082979
YFRAC(5)=0.108032 ;YFRAC(6)=0.135089
YFRAC(7)=0.164311 ;YFRAC(8)=0.195871
YFRAC(9)=0.229956 ;YFRAC(10)=0.266767
YFRAC(11)=0.306523 ;YFRAC(12)=0.349459
YFRAC(13)=0.395831 ;YFRAC(14)=0.445912
YFRAC(15)=0.5 ;YFRAC(16)=0.554088
YFRAC(17)=0.604169 ;YFRAC(18)=0.650541
YFRAC(19)=0.693477 ;YFRAC(20)=0.733233
YFRAC(21)=0.770044 ;YFRAC(22)=0.804129
YFRAC(23)=0.835689 ;YFRAC(24)=0.86491
YFRAC(25)=0.891968 ;YFRAC(26)=0.917021
YFRAC(27)=0.940218 ;YFRAC(28)=0.961697
YFRAC(29)=0.981585 ;YFRAC(30)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 30
ZWLAST =0.152
ZFRAC(1)=0.018415 ;ZFRAC(2)=0.038303
ZFRAC(3)=0.059782 ;ZFRAC(4)=0.082979
ZFRAC(5)=0.108032 ;ZFRAC(6)=0.135089
ZFRAC(7)=0.164311 ;ZFRAC(8)=0.195871
ZFRAC(9)=0.229956 ;ZFRAC(10)=0.266767
ZFRAC(11)=0.306523 ;ZFRAC(12)=0.349459
ZFRAC(13)=0.395831 ;ZFRAC(14)=0.445912
ZFRAC(15)=0.5 ;ZFRAC(16)=0.554088
ZFRAC(17)=0.604169 ;ZFRAC(18)=0.650541
ZFRAC(19)=0.693477 ;ZFRAC(20)=0.733233
ZFRAC(21)=0.770044 ;ZFRAC(22)=0.804129
ZFRAC(23)=0.835689 ;ZFRAC(24)=0.86491
ZFRAC(25)=0.891968 ;ZFRAC(26)=0.917021
ZFRAC(27)=0.940218 ;ZFRAC(28)=0.961697
ZFRAC(29)=0.981585 ;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 =26. ;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)=9.
FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.0E-10
FIINIT(KE)=0.374544 ;FIINIT(EP)=4.956048
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 ,10.626824 )
COVAL(OB1 ,U1 ,0. ,9. )
COVAL(OB1 ,V1 ,0. ,0. )
COVAL(OB1 ,W1 ,0. ,0. )
COVAL(OB1 ,KE ,0. ,0.374544 )
COVAL(OB1 ,EP ,0. ,4.956048 )
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.4E-06)
SPEDAT(SET,DFLUX,DENP4,R,1000.)
SPEDAT(SET,DFLUX,DIAP4,R,8.7E-06)
SPEDAT(SET,DFLUX,DENP5,R,1000.)
SPEDAT(SET,DFLUX,DIAP5,R,1.5E-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)=8.866 ;EX(U1)=8.527
EX(V1)=6.682E-03 ;EX(W1)=6.575E-03
EX(KE)=0.3554 ;EX(EP)=30.26
EX(C6)=0.9999 ;EX(C7)=0.9996
EX(C8)=0.999 ;EX(C9)=0.9462
EX(C10)=0.8038 ;EX(ENUN)=0.01738
EX(DI5)=2.807 ;EX(DI4)=202.899994
EX(DI3)=1530. ;EX(DI2)=5125.
EX(VP5)=0.01561 ;EX(TP5)=1.025
EX(VP4)=2.728E-03 ;EX(TP4)=0.3475
EX(VP3)=3.741E-04 ;EX(TP3)=0.1625
EX(VP2)=4.82E-05 ;EX(TP2)=0.05474
EX(VS5)=6.71E-03 ;EX(TR5)=6.84E-04
EX(VD5)=6.773E-03 ;EX(DEP5)=3.133E-03
EX(VS4)=2.275E-03 ;EX(TR4)=2.319E-04
EX(VD4)=1.202E-03 ;EX(DEP4)=8.564E-04
EX(VS3)=1.064E-03 ;EX(TR3)=1.085E-04
EX(VP1)=4.942E-06 ;EX(TP1)=7.356E-03
EX(DI1)=2.151E+07 ;EX(VS1)=4.816E-05
EX(TR1)=4.909E-06 ;EX(VD1)=2.038E-06
EX(VD3)=1.623E-04 ;EX(VSTR)=0.0545
EX(DEP1)=2.407E-06 ;EX(DEP3)=1.781E-04
EX(VS2)=3.584E-04 ;EX(ENUL)=1.563E-05
EX(TR2)=3.653E-05 ;EX(VD2)=2.014E-05
EX(DEP2)=2.368E-05 ;EX(YPLS)=4.88
EX(STRS)=0.02342 ;EX(EPKE)=56.93
EX(DEN1)=1.181 ;EX(EL1)=4.326E-03
EX(ENUT)=1.136E-03
************************************************************
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