TALK=T;RUN(1,1) DISPLAY This case simulates a device for the flow turning. The turn device is described by a box objects with domain properties wich turned revolved under 55 degrees angle concerning the horizon. The turn is executed by means of the sources for U1 and V1 variables for U1 as coeff*(-V1), for V1 as coeff*U1. The In-Form is used for the setting of the source per unit volume as (SOURCE of U1 at TURNING is :COEF:*(-V1) with VOLU) (SOURCE of V1 at TURNING is :COEF:*U1 with VOLU) The Q1 contains PHOTON USE commands ENDDIS PHOTON USE p;;; ve z 1 sh use patgeo enduse ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(Device for turning flow ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 40 XULAST =1. XFRAC(1)=0.025 ;XFRAC(2)=0.05 XFRAC(3)=0.075 ;XFRAC(4)=0.1 XFRAC(5)=0.125 ;XFRAC(6)=0.15 XFRAC(7)=0.175 ;XFRAC(8)=0.2 XFRAC(9)=0.225 ;XFRAC(10)=0.25 XFRAC(11)=0.275 ;XFRAC(12)=0.3 XFRAC(13)=0.325 ;XFRAC(14)=0.35 XFRAC(15)=0.375 ;XFRAC(16)=0.4 XFRAC(17)=0.425 ;XFRAC(18)=0.45 XFRAC(19)=0.475 ;XFRAC(20)=0.5 XFRAC(21)=0.524346 ;XFRAC(22)=0.548691 XFRAC(23)=0.573037 ;XFRAC(24)=0.597382 XFRAC(25)=0.622546 ;XFRAC(26)=0.647709 XFRAC(27)=0.672873 ;XFRAC(28)=0.698037 XFRAC(29)=0.7232 ;XFRAC(30)=0.748364 XFRAC(31)=0.773528 ;XFRAC(32)=0.798691 XFRAC(33)=0.823855 ;XFRAC(34)=0.849018 XFRAC(35)=0.874182 ;XFRAC(36)=0.899346 XFRAC(37)=0.924509 ;XFRAC(38)=0.949673 XFRAC(39)=0.974836 ;XFRAC(40)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 40 YVLAST =1. YFRAC(1)=0.013353 ;YFRAC(2)=0.038966 YFRAC(3)=0.064579 ;YFRAC(4)=0.090192 YFRAC(5)=0.115805 ;YFRAC(6)=0.141418 YFRAC(7)=0.167031 ;YFRAC(8)=0.192644 YFRAC(9)=0.218257 ;YFRAC(10)=0.24387 YFRAC(11)=0.269483 ;YFRAC(12)=0.295096 YFRAC(13)=0.320709 ;YFRAC(14)=0.346322 YFRAC(15)=0.371935 ;YFRAC(16)=0.397548 YFRAC(17)=0.423161 ;YFRAC(18)=0.448774 YFRAC(19)=0.474387 ;YFRAC(20)=0.5 YFRAC(21)=0.521245 ;YFRAC(22)=0.546442 YFRAC(23)=0.57164 ;YFRAC(24)=0.596838 YFRAC(25)=0.622035 ;YFRAC(26)=0.647233 YFRAC(27)=0.672431 ;YFRAC(28)=0.697628 YFRAC(29)=0.722826 ;YFRAC(30)=0.748024 YFRAC(31)=0.773221 ;YFRAC(32)=0.798419 YFRAC(33)=0.823616 ;YFRAC(34)=0.848814 YFRAC(35)=0.874012 ;YFRAC(36)=0.899209 YFRAC(37)=0.924407 ;YFRAC(38)=0.949605 YFRAC(39)=0.974802 ;YFRAC(40)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 1 ZWLAST =1. ZFRAC(1)=1. ************************************************************ Group 6. Body-Fitted Coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(1)=P1 ;NAME(3)=U1 NAME(5)=V1 ;NAME(150)=PRPS * 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(PRPS,Y,N,N,N,N,Y) PRPS = 150 ************************************************************ 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) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = F ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1 ************************************************************ Group 9. Properties used if PRPS is not stored, and where PRPS = -1.0 if it is! RHO1 =1.189 ;TMP1 =0. ;EL1 =0. TSURR =0. ;TEMP0 =0. ;PRESS0 =0. DVO1DT =3.41E-03 ;DRH1DP =0. EMISS =0. ;SCATT =0. RADIA =0. ;RADIB =0. ENUL =1.544E-05 ;ENUT =0. PRNDTL(U1)=1. ;PRNDTL(V1)=1. PRT(U1)=1. ;PRT(V1)=1. CP1 =1005. ;CP2 =1. ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10 FIINIT(V1)=1.0E-10 ;FIINIT(PRPS)=-1. PATCH(^OB1 ,INIVAL, 1, 0, 0, 0, 0, 0, 1, 1) INIT(^OB1 ,PRPS,0. ,-1. ) PATCH(OB2 ,INIVAL, 1, 20, 21, 40, 1, 1, 1, 1) INIT(OB2 ,PRPS,0. ,198. ) 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(^OD1 ,CELL , 1, 18, 0, 0, 0, 0, 1, 1) COVAL(^OD1 ,U1 ,In-Form:source - see Grp 19) COVAL(^OD1 ,V1 ,In-Form:source - see Grp 19) PATCH(OB3 ,WEST , 1, 1, 1, 20, 1, 1, 1, 1) COVAL(OB3 ,P1 , FIXFLU ,1.189 ) COVAL(OB3 ,U1 ,0. ,1. ) COVAL(OB3 ,V1 ,0. ,0. ) PATCH(OB4 ,NORTH , 21, 40, 40, 40, 1, 1, 1, 1) COVAL(OB4 ,P1 ,1. ,0. ) COVAL(OB4 ,U1 ,0. ,0. ) COVAL(OB4 ,V1 ,0. ,0. ) PATCH(OB5 ,EWALL , 40, 40, 1, 40, 1, 1, 1, 1) COVAL(OB5 ,V1 ,1. ,0. ) PATCH(OB6 ,SWALL , 1, 40, 1, 1, 1, 1, 1, 1) COVAL(OB6 ,U1 ,1. ,0. ) 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-05 ************************************************************ Group 16. Terminate Iterations LITER(P1)=20 ;LITER(U1)=10 LITER(V1)=10 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03 ENDIT(V1)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,1.) RELAX(U1,FALSDT,10.) RELAX(V1,FALSDT,10.) RELAX(PRPS,LINRLX,1.) 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(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10 ************************************************************ Group 19. Data transmitted to GROUND PARSOL = F ISG62 = 1 SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,0) SPEDAT(SET,SOURCE,U1!TURNING,C,=10.*(-V1)!VOLU) SPEDAT(SET,SOURCE,V1!TURNING,C,=10.*U1!VOLU) SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,OBJNAM,^OB1,C,TURNING) SPEDAT(SET,OBJTYP,^OB1,C,BLOCKAGE) SPEDAT(SET,TURNING,MATERIAL,R,-1.) SPEDAT(SET,OBJNAM,^OD1,C,TURNING) SPEDAT(SET,OBJTYP,^OD1,C,BLOCKAGE) SPEDAT(SET,OBJNAM,!OB2,C,BLOCK) SPEDAT(SET,OBJTYP,!OB2,C,BLOCKAGE) SPEDAT(SET,BLOCK,MATERIAL,R,198.) SPEDAT(SET,OBJNAM,!OB3,C,INLET) SPEDAT(SET,OBJTYP,!OB3,C,INLET) SPEDAT(SET,OBJNAM,!OB4,C,OUTLET) SPEDAT(SET,OBJTYP,!OB4,C,OUTLET) SPEDAT(SET,ARATIO,!OB4,R,1.) SPEDAT(SET,OBJNAM,!OB5,C,E_WALL) SPEDAT(SET,OBJTYP,!OB5,C,PLATE) SPEDAT(SET,OBJNAM,!OB6,C,S_WALL) SPEDAT(SET,OBJTYP,!OB6,C,PLATE) SPEDAT(SET,FACETDAT,NUMOBJ,I,6) SPEDAT(SET,MATERIAL,198,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=1.023 ;EX(U1)=0.3851 EX(V1)=0.4004 ;EX(PRPS)=50.25 ************************************************************ 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,Y,N,Y,Y,Y,Y) OUTPUT(U1,Y,N,Y,Y,Y,Y) OUTPUT(V1,Y,N,Y,Y,Y,Y) OUTPUT(PRPS,Y,N,Y,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 26 ;IYMON = 14 ;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 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