TALK=T;RUN(1,1) DISPLAY The case considered is the normal impingement of a turbulent axisymmetric jet on a heated flat plate. The flow geometry and conditions approximate to those reported experimentally by D.Cooper et al ( Int.J.Heat Mass Transfer, Vol.36, No.10, p2675,[1993] ). The jet issues into ambient fluid at a Reynolds number of 2.3E4 from a pipe located 2 pipe diameters above the electrically heated plate. The jet injection temperature is ambient so that the purpose of the flow process is to effect cooling. For the case of orthogonal impaction, the flow is axisymmetric and is directed outward along the surface giving rise to a radial wall jet. The simulation may be performed with the standard k-e model or with the Reynolds stress transport model. For the latter, the model employs differential transport equations for determining the turbulent fluxes of heat. ENDDIS ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(RSTM_2D IMPINGING ROUND JET :T609) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 609 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = F NX = 1 XULAST =0.1 XFRAC(1)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 25 YVLAST =5. YFRAC(1)=1.0E-02 ;YFRAC(2)=0.02 YFRAC(3)=0.03 ;YFRAC(4)=0.04 YFRAC(5)=0.05 ;YFRAC(6)=0.06 YFRAC(7)=0.07 ;YFRAC(8)=0.08 YFRAC(9)=0.09 ;YFRAC(10)=0.1 YFRAC(11)=0.115492 ;YFRAC(12)=0.143818 YFRAC(13)=0.180498 ;YFRAC(14)=0.223935 YFRAC(15)=0.273205 ;YFRAC(16)=0.327684 YFRAC(17)=0.386915 ;YFRAC(18)=0.450542 YFRAC(19)=0.518282 ;YFRAC(20)=0.589898 YFRAC(21)=0.66519 ;YFRAC(22)=0.743988 YFRAC(23)=0.826141 ;YFRAC(24)=0.911517 YFRAC(25)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 30 ZWLAST =2. ZFRAC(1)=0.065812 ;ZFRAC(2)=0.114585 ZFRAC(3)=0.158489 ;ZFRAC(4)=0.199504 ZFRAC(5)=0.238495 ;ZFRAC(6)=0.275946 ZFRAC(7)=0.312163 ;ZFRAC(8)=0.347356 ZFRAC(9)=0.381678 ;ZFRAC(10)=0.415244 ZFRAC(11)=0.448144 ;ZFRAC(12)=0.48045 ZFRAC(13)=0.512221 ;ZFRAC(14)=0.543507 ZFRAC(15)=0.574349 ;ZFRAC(16)=0.604782 ZFRAC(17)=0.634837 ;ZFRAC(18)=0.66454 ZFRAC(19)=0.693914 ;ZFRAC(20)=0.722981 ZFRAC(21)=0.751759 ;ZFRAC(22)=0.780263 ZFRAC(23)=0.80851 ;ZFRAC(24)=0.836512 ZFRAC(25)=0.864281 ;ZFRAC(26)=0.891829 ZFRAC(27)=0.919166 ;ZFRAC(28)=0.946301 ZFRAC(29)=0.973243 ;ZFRAC(30)=1. ************************************************************ Group 6. Body-Fitted Coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(1)=P1 ;NAME(5)=V1 NAME(7)=W1 ;NAME(12)=KE NAME(13)=EP ;NAME(14)=H1 NAME(142)=TMP1 ;NAME(143)=WTRS NAME(144)=VTRS ;NAME(145)=U2RS NAME(146)=V2RS ;NAME(147)=W2RS NAME(148)=VWRS ;NAME(149)=LEN1 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,N) SOLUTN(V1,Y,Y,N,N,N,Y) SOLUTN(W1,Y,Y,N,N,N,Y) SOLUTN(KE,Y,N,N,N,N,N) SOLUTN(EP,Y,Y,N,N,N,N) SOLUTN(H1,Y,Y,Y,N,N,Y) SOLUTN(TMP1,Y,N,N,N,N,Y) SOLUTN(WTRS,Y,Y,N,N,N,N) SOLUTN(VTRS,Y,Y,N,N,N,N) SOLUTN(U2RS,Y,Y,N,N,N,N) SOLUTN(V2RS,Y,Y,N,N,N,N) SOLUTN(W2RS,Y,Y,N,N,N,N) SOLUTN(VWRS,Y,Y,N,N,N,N) SOLUTN(LEN1,Y,N,N,N,N,N) SOLUTN(ENUT,Y,N,N,N,N,N) VIST = 150 TEMP1 = 142 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(V1,Y,Y,Y,Y,Y,Y) TERMS(W1,Y,Y,Y,Y,Y,Y) TERMS(EP,Y,Y,Y,Y,N,Y) TERMS(H1,Y,Y,Y,Y,Y,Y) TERMS(WTRS,N,Y,Y,Y,N,Y) TERMS(VTRS,N,Y,Y,Y,Y,Y) TERMS(U2RS,N,Y,Y,Y,N,Y) TERMS(V2RS,N,Y,Y,Y,Y,Y) TERMS(W2RS,N,Y,Y,Y,N,Y) TERMS(VWRS,N,Y,Y,Y,Y,Y) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = F NEWENT = F UDIFF = T UCONNE = T USOURC = 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 =1. ;TMP1 = GRND2 ;EL1 = GRND TSURR =0. ;TEMP0 =0. ;PRESS0 =0. DVO1DT =0. ;DRH1DP =0. TMP1A =0. ;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 =1.408451E-05 ;ENUT = GRND ENUTA =0. ;ENUTB =0. ;ENUTC =0. IENUTA = 0 PRNDTL(V1)=1. ;PRNDTL(W1)=1. PRNDTL(EP)=1. ;PRNDTL(H1)=0.71 PRNDTL(WTRS)=1. ;PRNDTL(VTRS)=1. PRNDTL(U2RS)=1. ;PRNDTL(V2RS)=1. PRNDTL(W2RS)=1. ;PRNDTL(VWRS)=1. PRT(V1)=1. ;PRT(W1)=1. PRT(EP)=1.314 ;PRT(H1)=0.9 PRT(WTRS)=1. ;PRT(VTRS)=1. PRT(U2RS)=1. ;PRT(V2RS)=1. PRT(W2RS)=1. ;PRT(VWRS)=1. CP1 =1000. ;CP2 =1. ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(P1)=1.0E-10 ;FIINIT(V1)=0. FIINIT(W1)=1.0E-10 ;FIINIT(KE)=1.0E-02 FIINIT(EP)=0.0639 ;FIINIT(H1)=3.0E+05 FIINIT(TMP1)=1.0E-10 ;FIINIT(WTRS)=1.0E-10 FIINIT(VTRS)=1.0E-10 ;FIINIT(U2RS)=6.666666E-03 FIINIT(V2RS)=6.666666E-03 ;FIINIT(W2RS)=6.666666E-03 FIINIT(VWRS)=3.0E-03 ;FIINIT(LEN1)=1.0E-10 FIINIT(ENUT)=1.0E-10 PATCH(INWJET ,INIVAL, 1, 1, 1, 10, 1, 30, 1, 1) INIT(INWJET ,W1 ,0. ,1. ) 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(PLATE ,HWALL , 1, 1, 1, 25, 30, 30, 1, 1) COVAL(PLATE ,V1 , GRND2 ,0. ) COVAL(PLATE ,EP , GRND2 , GRND2 ) COVAL(PLATE ,H1 , FIXFLU ,1.0E+04 ) PATCH(WDMPH03 ,PHASEM, 1, 1, 1, 25, 1, 30, 1, 1) COVAL(WDMPH03 ,WTRS, FIXFLU , GRND1 ) COVAL(WDMPH03 ,VTRS, FIXFLU , GRND1 ) COVAL(WDMPH03 ,U2RS, FIXFLU , GRND1 ) COVAL(WDMPH03 ,V2RS, FIXFLU , GRND1 ) COVAL(WDMPH03 ,W2RS, FIXFLU , GRND1 ) COVAL(WDMPH03 ,VWRS, FIXFLU , GRND1 ) PATCH(JET1 ,LOW , 1, 1, 1, 10, 1, 1, 1, 1) COVAL(JET1 ,P1 , FIXFLU ,1. ) COVAL(JET1 ,V1 ,0. ,0. ) COVAL(JET1 ,W1 ,0. ,1. ) COVAL(JET1 ,EP ,0. ,1.825555E-03 ) COVAL(JET1 ,H1 ,0. ,3.0E+05 ) COVAL(JET1 ,WTRS,0. ,0. ) COVAL(JET1 ,VTRS,0. ,0. ) COVAL(JET1 ,U2RS,0. ,6.666666E-03 ) COVAL(JET1 ,V2RS,0. ,6.666666E-03 ) COVAL(JET1 ,W2RS,0. ,6.666666E-03 ) COVAL(JET1 ,VWRS,0. ,0. ) PATCH(TOP ,LOW , 1, 1, 11, 25, 1, 1, 1, 1) COVAL(TOP ,P1 ,1000. ,0. ) COVAL(TOP ,V1 ,0. ,0. ) COVAL(TOP ,W1 ,0. ,0. ) COVAL(TOP ,KE ,0. ,1.0E-10 ) COVAL(TOP ,EP ,0. ,1.0E-10 ) COVAL(TOP ,H1 ,0. ,3.0E+05 ) PATCH(SIDE ,NORTH , 1, 1, 25, 25, 1, 30, 1, 1) COVAL(SIDE ,P1 ,1000. ,0. ) COVAL(SIDE ,V1 ,0. ,0. ) COVAL(SIDE ,W1 ,0. ,0. ) COVAL(SIDE ,KE ,0. ,1.0E-10 ) COVAL(SIDE ,EP ,0. ,1.0E-10 ) COVAL(SIDE ,H1 ,0. ,3.0E+05 ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 1000 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SELREF = T RESFAC =1.0E-05 ************************************************************ Group 16. Terminate Iterations LITER(P1)=20 ;LITER(V1)=10 LITER(W1)=10 ;LITER(EP)=20 LITER(H1)=30 ;LITER(WTRS)=20 LITER(VTRS)=20 ;LITER(U2RS)=20 LITER(V2RS)=20 ;LITER(W2RS)=20 LITER(VWRS)=20 ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03 ENDIT(W1)=1.0E-03 ;ENDIT(EP)=1.0E-03 ENDIT(H1)=1.0E-03 ;ENDIT(WTRS)=1.0E-03 ENDIT(VTRS)=1.0E-03 ;ENDIT(U2RS)=1.0E-03 ENDIT(V2RS)=1.0E-03 ;ENDIT(W2RS)=1.0E-03 ENDIT(VWRS)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,1.) RELAX(V1,FALSDT,0.266667) RELAX(W1,FALSDT,0.266667) RELAX(KE,LINRLX,1.) RELAX(EP,FALSDT,0.066667) RELAX(H1,FALSDT,0.666667) RELAX(TMP1,LINRLX,1.) RELAX(WTRS,FALSDT,0.016667) RELAX(VTRS,FALSDT,0.016667) RELAX(U2RS,FALSDT,0.066667) RELAX(V2RS,FALSDT,0.066667) RELAX(W2RS,FALSDT,0.066667) RELAX(VWRS,FALSDT,0.066667) RELAX(LEN1,LINRLX,1.) RELAX(ENUT,LINRLX,1.) OVRRLX =0. EXPERT = F ;NNORSL = F ************************************************************ Group 18. Limits VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10 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(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=-1.0E+10 VARMAX(WTRS)=1.0E+10 ;VARMIN(WTRS)=-1.0E+10 VARMAX(VTRS)=1.0E+10 ;VARMIN(VTRS)=-1.0E+10 VARMAX(U2RS)=1.0E+10 ;VARMIN(U2RS)=-1.0E+10 VARMAX(V2RS)=1.0E+10 ;VARMIN(V2RS)=-1.0E+10 VARMAX(W2RS)=1.0E+10 ;VARMIN(W2RS)=-1.0E+10 VARMAX(VWRS)=1.0E+10 ;VARMIN(VWRS)=-1.0E+10 VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10 ************************************************************ Group 19. Data transmitted to GROUND GENK = T RSTM = T PARSOL = F IRSMSM = 2 ISG62 = 1 SPEDAT(SET,GXMONI,PLOTALL,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(P1,Y,N,Y,Y,Y,Y) OUTPUT(V1,Y,N,Y,Y,Y,Y) OUTPUT(W1,Y,N,Y,Y,Y,Y) OUTPUT(KE,Y,N,Y,N,N,N) OUTPUT(EP,Y,N,Y,Y,Y,Y) OUTPUT(H1,Y,N,Y,Y,Y,Y) OUTPUT(TMP1,Y,N,Y,N,N,N) OUTPUT(WTRS,Y,N,Y,Y,Y,Y) OUTPUT(VTRS,Y,N,Y,Y,Y,Y) OUTPUT(U2RS,Y,N,Y,Y,Y,Y) OUTPUT(V2RS,Y,N,Y,Y,Y,Y) OUTPUT(W2RS,Y,N,Y,Y,Y,Y) OUTPUT(VWRS,Y,N,Y,Y,Y,Y) OUTPUT(LEN1,Y,N,Y,N,N,N) OUTPUT(ENUT,Y,N,Y,N,N,N) WALPRN = T ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 12 ;IZMON = 28 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 = 25 NZPRIN = 1 ;IZPRF = 1 ;IZPRL = 10000 XZPR = F ;YZPR = F IPLTF = 1 ;IPLTL = -1 ;NPLT = 5 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