TALK=T;RUN(1,1) ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(Flow Over Double-Wedge Airfoil. USP Test) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 1 XULAST =1. XFRAC(1)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 80 YVLAST =1. YFRAC(1)=0.0125 ;YFRAC(2)=0.025 YFRAC(3)=0.0375 ;YFRAC(4)=0.05 YFRAC(5)=0.0625 ;YFRAC(6)=0.075 YFRAC(7)=0.0875 ;YFRAC(8)=0.1 YFRAC(9)=0.1125 ;YFRAC(10)=0.125 YFRAC(11)=0.1375 ;YFRAC(12)=0.15 YFRAC(13)=0.1625 ;YFRAC(14)=0.175 YFRAC(15)=0.1875 ;YFRAC(16)=0.2 YFRAC(17)=0.2125 ;YFRAC(18)=0.225 YFRAC(19)=0.2375 ;YFRAC(20)=0.25 YFRAC(21)=0.2625 ;YFRAC(22)=0.275 YFRAC(23)=0.2875 ;YFRAC(24)=0.3 YFRAC(25)=0.3125 ;YFRAC(26)=0.325 YFRAC(27)=0.3375 ;YFRAC(28)=0.35 YFRAC(29)=0.3625 ;YFRAC(30)=0.375 YFRAC(31)=0.3875 ;YFRAC(32)=0.4 YFRAC(33)=0.4125 ;YFRAC(34)=0.425 YFRAC(35)=0.4375 ;YFRAC(36)=0.45 YFRAC(37)=0.4625 ;YFRAC(38)=0.475 YFRAC(39)=0.4875 ;YFRAC(40)=0.5 YFRAC(41)=0.5125 ;YFRAC(42)=0.525 YFRAC(43)=0.5375 ;YFRAC(44)=0.55 YFRAC(45)=0.5625 ;YFRAC(46)=0.575 YFRAC(47)=0.5875 ;YFRAC(48)=0.6 YFRAC(49)=0.6125 ;YFRAC(50)=0.625 YFRAC(51)=0.6375 ;YFRAC(52)=0.65 YFRAC(53)=0.6625 ;YFRAC(54)=0.675 YFRAC(55)=0.6875 ;YFRAC(56)=0.7 YFRAC(57)=0.7125 ;YFRAC(58)=0.725 YFRAC(59)=0.7375 ;YFRAC(60)=0.75 YFRAC(61)=0.7625 ;YFRAC(62)=0.775 YFRAC(63)=0.7875 ;YFRAC(64)=0.8 YFRAC(65)=0.8125 ;YFRAC(66)=0.825 YFRAC(67)=0.8375 ;YFRAC(68)=0.85 YFRAC(69)=0.8625 ;YFRAC(70)=0.875 YFRAC(71)=0.8875 ;YFRAC(72)=0.9 YFRAC(73)=0.9125 ;YFRAC(74)=0.925 YFRAC(75)=0.9375 ;YFRAC(76)=0.95 YFRAC(77)=0.9625 ;YFRAC(78)=0.975 YFRAC(79)=0.9875 ;YFRAC(80)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 160 ZWLAST =8. ZFRAC(1)=6.25E-03 ;ZFRAC(2)=0.0125 ZFRAC(3)=0.01875 ;ZFRAC(4)=0.025 ZFRAC(5)=0.03125 ;ZFRAC(6)=0.0375 ZFRAC(7)=0.04375 ;ZFRAC(8)=0.05 ZFRAC(9)=0.05625 ;ZFRAC(10)=0.0625 ZFRAC(11)=0.06875 ;ZFRAC(12)=0.075 ZFRAC(13)=0.08125 ;ZFRAC(14)=0.0875 ZFRAC(15)=0.09375 ;ZFRAC(16)=0.1 ZFRAC(17)=0.10625 ;ZFRAC(18)=0.1125 ZFRAC(19)=0.11875 ;ZFRAC(20)=0.125 ZFRAC(21)=0.13125 ;ZFRAC(22)=0.1375 ZFRAC(23)=0.14375 ;ZFRAC(24)=0.15 ZFRAC(25)=0.15625 ;ZFRAC(26)=0.1625 ZFRAC(27)=0.16875 ;ZFRAC(28)=0.175 ZFRAC(29)=0.18125 ;ZFRAC(30)=0.1875 ZFRAC(31)=0.19375 ;ZFRAC(32)=0.2 ZFRAC(33)=0.20625 ;ZFRAC(34)=0.2125 ZFRAC(35)=0.21875 ;ZFRAC(36)=0.225 ZFRAC(37)=0.23125 ;ZFRAC(38)=0.2375 ZFRAC(39)=0.24375 ;ZFRAC(40)=0.25 ZFRAC(41)=0.25625 ;ZFRAC(42)=0.2625 ZFRAC(43)=0.26875 ;ZFRAC(44)=0.275 ZFRAC(45)=0.28125 ;ZFRAC(46)=0.2875 ZFRAC(47)=0.29375 ;ZFRAC(48)=0.3 ZFRAC(49)=0.30625 ;ZFRAC(50)=0.3125 ZFRAC(51)=0.31875 ;ZFRAC(52)=0.325 ZFRAC(53)=0.33125 ;ZFRAC(54)=0.3375 ZFRAC(55)=0.34375 ;ZFRAC(56)=0.35 ZFRAC(57)=0.35625 ;ZFRAC(58)=0.3625 ZFRAC(59)=0.36875 ;ZFRAC(60)=0.375 ZFRAC(61)=0.38125 ;ZFRAC(62)=0.3875 ZFRAC(63)=0.39375 ;ZFRAC(64)=0.4 ZFRAC(65)=0.40625 ;ZFRAC(66)=0.4125 ZFRAC(67)=0.41875 ;ZFRAC(68)=0.425 ZFRAC(69)=0.43125 ;ZFRAC(70)=0.4375 ZFRAC(71)=0.44375 ;ZFRAC(72)=0.45 ZFRAC(73)=0.45625 ;ZFRAC(74)=0.4625 ZFRAC(75)=0.46875 ;ZFRAC(76)=0.475 ZFRAC(77)=0.48125 ;ZFRAC(78)=0.4875 ZFRAC(79)=0.49375 ;ZFRAC(80)=0.5 ZFRAC(81)=0.50625 ;ZFRAC(82)=0.5125 ZFRAC(83)=0.51875 ;ZFRAC(84)=0.525 ZFRAC(85)=0.53125 ;ZFRAC(86)=0.5375 ZFRAC(87)=0.54375 ;ZFRAC(88)=0.55 ZFRAC(89)=0.55625 ;ZFRAC(90)=0.5625 ZFRAC(91)=0.56875 ;ZFRAC(92)=0.575 ZFRAC(93)=0.58125 ;ZFRAC(94)=0.5875 ZFRAC(95)=0.59375 ;ZFRAC(96)=0.6 ZFRAC(97)=0.60625 ;ZFRAC(98)=0.6125 ZFRAC(99)=0.61875 ;ZFRAC(100)=0.625 ZFRAC(101)=0.63125 ;ZFRAC(102)=0.6375 ZFRAC(103)=0.64375 ;ZFRAC(104)=0.65 ZFRAC(105)=0.65625 ;ZFRAC(106)=0.6625 ZFRAC(107)=0.66875 ;ZFRAC(108)=0.675 ZFRAC(109)=0.68125 ;ZFRAC(110)=0.6875 ZFRAC(111)=0.69375 ;ZFRAC(112)=0.7 ZFRAC(113)=0.70625 ;ZFRAC(114)=0.7125 ZFRAC(115)=0.71875 ;ZFRAC(116)=0.725 ZFRAC(117)=0.73125 ;ZFRAC(118)=0.7375 ZFRAC(119)=0.74375 ;ZFRAC(120)=0.75 ZFRAC(121)=0.75625 ;ZFRAC(122)=0.7625 ZFRAC(123)=0.76875 ;ZFRAC(124)=0.775 ZFRAC(125)=0.78125 ;ZFRAC(126)=0.7875 ZFRAC(127)=0.79375 ;ZFRAC(128)=0.8 ZFRAC(129)=0.80625 ;ZFRAC(130)=0.8125 ZFRAC(131)=0.81875 ;ZFRAC(132)=0.825 ZFRAC(133)=0.83125 ;ZFRAC(134)=0.8375 ZFRAC(135)=0.84375 ;ZFRAC(136)=0.85 ZFRAC(137)=0.85625 ;ZFRAC(138)=0.8625 ZFRAC(139)=0.86875 ;ZFRAC(140)=0.875 ZFRAC(141)=0.88125 ;ZFRAC(142)=0.8875 ZFRAC(143)=0.89375 ;ZFRAC(144)=0.9 ZFRAC(145)=0.90625 ;ZFRAC(146)=0.9125 ZFRAC(147)=0.91875 ;ZFRAC(148)=0.925 ZFRAC(149)=0.93125 ;ZFRAC(150)=0.9375 ZFRAC(151)=0.94375 ;ZFRAC(152)=0.95 ZFRAC(153)=0.95625 ;ZFRAC(154)=0.9625 ZFRAC(155)=0.96875 ;ZFRAC(156)=0.975 ZFRAC(157)=0.98125 ;ZFRAC(158)=0.9875 ZFRAC(159)=0.99375 ;ZFRAC(160)=1. ************************************************************ Group 6. Body-Fitted Coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(5)=V1 ;NAME(7)=W1 NAME(149)=PRPS ;NAME(150)=POT * Y in SOLUTN argument list denotes: * 1-stored 2-solved 3-whole-field * 4-point-by-point 5-explicit 6-harmonic averaging SOLUTN(V1,Y,N,N,N,N,Y) SOLUTN(W1,Y,N,N,N,N,Y) SOLUTN(PRPS,Y,N,N,N,N,Y) SOLUTN(POT,Y,Y,N,N,N,Y) PRPS = 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(POT,N,Y,Y,Y,Y,Y) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = F ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1 USP = T ; USPDBG = F UTCPLT = T ; USPIMB = F UAUTO = F ; UGTGRD = F USPVTK = T ; USPTER = F USCANX = T ; USCANY = T USCANZ = T MXLEV = 2 ;MYLEV = 2 ;MZLEV = 2 DOMAT = -1 CELLST = 10 ; FACEST = 10 MINPRP = -1 ; MAXPRP = 250 UMXLEV = 1 ; UMXCLS = 10000 UNOLAY = 2 ; UIOLEV = 1 UWALEV = 1 ; UGRITR = 2 UMXCRS = 16 USPREL =0.7 ; UMNCEL =1. UDVREL =1. ************************************************************ Group 9. Properties used if PRPS is not stored, and where PRPS = -1.0 if it is! RHO1 =1. ;TMP1 =0. ;EL1 =0. TSURR =0. ;TEMP0 =0. ;PRESS0 =0. DVO1DT =0. ;DRH1DP =0. EMISS =0. ;SCATT =0. RADIA =0. ;RADIB =0. ENUL =1.0E-05 ;ENUT =0. PRNDTL(POT)=1. PRT(POT)=1. CP1 =1. ;CP2 =1. ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.0E-10 FIINIT(PRPS)=-1. ;FIINIT(POT)=1.0E-10 PATCH(HPOR1 ,INIVAL, 1, 1, 1, 1, 41, 119, 1, 1) INIT(HPOR1 ,PRPS,0. ,199. ) PATCH(HPOR2 ,INIVAL, 1, 1, 2, 2, 42, 118, 1, 1) INIT(HPOR2 ,PRPS,0. ,199. ) PATCH(HPOR3 ,INIVAL, 1, 1, 3, 3, 43, 117, 1, 1) INIT(HPOR3 ,PRPS,0. ,199. ) PATCH(HPOR4 ,INIVAL, 1, 1, 4, 4, 44, 116, 1, 1) INIT(HPOR4 ,PRPS,0. ,199. ) PATCH(HPOR5 ,INIVAL, 1, 1, 5, 5, 45, 115, 1, 1) INIT(HPOR5 ,PRPS,0. ,199. ) PATCH(HPOR6 ,INIVAL, 1, 1, 6, 6, 46, 114, 1, 1) INIT(HPOR6 ,PRPS,0. ,199. ) PATCH(HPOR7 ,INIVAL, 1, 1, 7, 7, 47, 113, 1, 1) INIT(HPOR7 ,PRPS,0. ,199. ) PATCH(HPOR8 ,INIVAL, 1, 1, 8, 8, 48, 112, 1, 1) INIT(HPOR8 ,PRPS,0. ,199. ) PATCH(HPOR9 ,INIVAL, 1, 1, 9, 9, 49, 111, 1, 1) INIT(HPOR9 ,PRPS,0. ,199. ) PATCH(HPOR10 ,INIVAL, 1, 1, 10, 10, 50, 110, 1, 1) INIT(HPOR10 ,PRPS,0. ,199. ) PATCH(HPOR11 ,INIVAL, 1, 1, 11, 11, 51, 109, 1, 1) INIT(HPOR11 ,PRPS,0. ,199. ) PATCH(HPOR12 ,INIVAL, 1, 1, 12, 12, 52, 108, 1, 1) INIT(HPOR12 ,PRPS,0. ,199. ) PATCH(HPOR13 ,INIVAL, 1, 1, 13, 13, 53, 107, 1, 1) INIT(HPOR13 ,PRPS,0. ,199. ) PATCH(HPOR14 ,INIVAL, 1, 1, 14, 14, 54, 106, 1, 1) INIT(HPOR14 ,PRPS,0. ,199. ) PATCH(HPOR15 ,INIVAL, 1, 1, 15, 15, 55, 105, 1, 1) INIT(HPOR15 ,PRPS,0. ,199. ) PATCH(HPOR16 ,INIVAL, 1, 1, 16, 16, 56, 104, 1, 1) INIT(HPOR16 ,PRPS,0. ,199. ) PATCH(HPOR17 ,INIVAL, 1, 1, 17, 17, 57, 103, 1, 1) INIT(HPOR17 ,PRPS,0. ,199. ) PATCH(HPOR18 ,INIVAL, 1, 1, 18, 18, 58, 102, 1, 1) INIT(HPOR18 ,PRPS,0. ,199. ) PATCH(HPOR19 ,INIVAL, 1, 1, 19, 19, 59, 101, 1, 1) INIT(HPOR19 ,PRPS,0. ,199. ) PATCH(HPOR20 ,INIVAL, 1, 1, 20, 20, 60, 100, 1, 1) INIT(HPOR20 ,PRPS,0. ,199. ) PATCH(HPOR21 ,INIVAL, 1, 1, 21, 21, 61, 99, 1, 1) INIT(HPOR21 ,PRPS,0. ,199. ) PATCH(HPOR22 ,INIVAL, 1, 1, 22, 22, 62, 98, 1, 1) INIT(HPOR22 ,PRPS,0. ,199. ) PATCH(HPOR23 ,INIVAL, 1, 1, 23, 23, 63, 97, 1, 1) INIT(HPOR23 ,PRPS,0. ,199. ) PATCH(HPOR24 ,INIVAL, 1, 1, 24, 24, 64, 96, 1, 1) INIT(HPOR24 ,PRPS,0. ,199. ) PATCH(HPOR25 ,INIVAL, 1, 1, 25, 25, 65, 95, 1, 1) INIT(HPOR25 ,PRPS,0. ,199. ) PATCH(HPOR26 ,INIVAL, 1, 1, 26, 26, 66, 94, 1, 1) INIT(HPOR26 ,PRPS,0. ,199. ) PATCH(HPOR27 ,INIVAL, 1, 1, 27, 27, 67, 93, 1, 1) INIT(HPOR27 ,PRPS,0. ,199. ) PATCH(HPOR28 ,INIVAL, 1, 1, 28, 28, 68, 92, 1, 1) INIT(HPOR28 ,PRPS,0. ,199. ) PATCH(HPOR29 ,INIVAL, 1, 1, 29, 29, 69, 91, 1, 1) INIT(HPOR29 ,PRPS,0. ,199. ) PATCH(HPOR30 ,INIVAL, 1, 1, 30, 30, 70, 90, 1, 1) INIT(HPOR30 ,PRPS,0. ,199. ) PATCH(HPOR31 ,INIVAL, 1, 1, 31, 31, 71, 89, 1, 1) INIT(HPOR31 ,PRPS,0. ,199. ) PATCH(HPOR32 ,INIVAL, 1, 1, 32, 32, 72, 88, 1, 1) INIT(HPOR32 ,PRPS,0. ,199. ) PATCH(HPOR33 ,INIVAL, 1, 1, 33, 33, 73, 87, 1, 1) INIT(HPOR33 ,PRPS,0. ,199. ) PATCH(HPOR34 ,INIVAL, 1, 1, 34, 34, 74, 86, 1, 1) INIT(HPOR34 ,PRPS,0. ,199. ) PATCH(HPOR35 ,INIVAL, 1, 1, 35, 35, 75, 85, 1, 1) INIT(HPOR35 ,PRPS,0. ,199. ) PATCH(HPOR36 ,INIVAL, 1, 1, 36, 36, 76, 84, 1, 1) INIT(HPOR36 ,PRPS,0. ,199. ) PATCH(HPOR37 ,INIVAL, 1, 1, 37, 37, 77, 83, 1, 1) INIT(HPOR37 ,PRPS,0. ,199. ) PATCH(HPOR38 ,INIVAL, 1, 1, 38, 38, 78, 82, 1, 1) INIT(HPOR38 ,PRPS,0. ,199. ) PATCH(HPOR39 ,INIVAL, 1, 1, 39, 39, 79, 81, 1, 1) INIT(HPOR39 ,PRPS,0. ,199. ) PATCH(HPOR40 ,INIVAL, 1, 1, 40, 40, 80, 80, 1, 1) INIT(HPOR40 ,PRPS,0. ,199. ) 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(UPSTRM ,LOW , 1, 1, 1, 80, 1, 1, 1, 1) COVAL(UPSTRM ,POT , FIXVAL ,4. ) PATCH(DWNSTRM ,HIGH , 1, 1, 1, 80, 160, 160, 1, 1) COVAL(DWNSTRM ,POT , FIXVAL ,-4. ) PATCH(^OB1 ,CELL , 1, 0, 0, 0, 0, 0, 1, 1) PATCH(^OB2 ,CELL , 2, 0, 0, 0, 0, 0, 1, 1) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 200 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SELREF = T RESFAC =1.0E-12 ************************************************************ Group 16. Terminate Iterations LITER(POT)=20 ENDIT(POT)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(V1,LINRLX,1.) RELAX(W1,LINRLX,1.) RELAX(PRPS,LINRLX,1.) RELAX(POT,FALSDT,1.0E+09) OVRRLX =0. EXPERT = F ;NNORSL = F ************************************************************ Group 18. Limits VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10 VARMAX(POT)=1.0E+10 ;VARMIN(POT)=-1.0E+10 ************************************************************ Group 19. Data transmitted to GROUND POTVEL = T PARSOL = F ISG62 = 1 SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,OBJNAM,^OB1,C,B1) SPEDAT(SET,OBJTYP,^OB1,C,USER_DEFINED) SPEDAT(SET,^OB1,LEVEL_X,I,0) SPEDAT(SET,^OB1,LEVEL_Y,I,0) SPEDAT(SET,^OB1,LEVEL_Z,I,0) SPEDAT(SET,B1,DATFILE,C,fine) SPEDAT(SET,^OB1,AXES,C,OYZ) SPEDAT(SET,OBJNAM,^OB2,C,B2) SPEDAT(SET,OBJTYP,^OB2,C,USER_DEFINED) SPEDAT(SET,^OB2,LEVEL_X,I,1) SPEDAT(SET,^OB2,LEVEL_Y,I,1) SPEDAT(SET,^OB2,LEVEL_Z,I,1) SPEDAT(SET,B2,DATFILE,C,fine) SPEDAT(SET,^OB2,AXES,C,OYZ) SPEDAT(SET,FACETDAT,NUMOBJ,I,2) SPEDAT(SET,MATERIAL,199,L,T) ************************************************************ Group 20. Preliminary Printout ************************************************************ 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(V1,Y,N,Y,N,N,N) OUTPUT(W1,Y,N,Y,N,N,N) OUTPUT(PRPS,Y,N,Y,N,N,N) OUTPUT(POT,Y,N,Y,Y,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 56 ;IZMON = 80 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1 UWATCH = T ;USTEER = T HIGHLO = F ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 ;NUMCLS = 5 NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 10000 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000 XZPR = F ;YZPR = F IPLTF = 2 ;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