TALK=T;RUN(1,1) DISPLAY The case considered is 2D incompressible, turbulent flow past a thick flat rectangular plate with a sharp leading edge, as studied experimentally by Djilali,N. & Gartshore,I.S., "Turbulent flow around a bluff rectangular plate, Part I: Experimental Investigation", Journal of Fluid Mechanics, Vol.113, pp51-59, (1991). The flow separates at the leading edge of the plate, and then reattaches further downstream to form a long separation zone on top of the plate. The thickness of the plate H is taken as 0.1m, and the working fluid is air at standard temperature and pressure.The flow Reynolds number based on the free-stream velocity and plate thickness H is 50,000. The inlet and outlet planes are located 10H upstream and downstream of the leading edge of the plate. The height of the solution domain is taken as 10H, and for simplicity a zero flux boundary is assumed here. Symmetry is exploited so that only one half of the flow is simulated. A fixed-pressure boundary condition is applied at the outlet, and uniform flow profiles are specified at the inlet. Scalable wall functions are used at the walls of the plate. The default simulation is made with the realisable k-e model, but the case can also be run with the standard k-e and k-w models, and the following k-e variants: Chen-Kim, RNG, Kato-Launder, Murakami et al. The measured and predicted reattachment points are listed below: KE REAL_KE CHEN RNG MMK KL KO EXPT Xr/H = 1.08 4.85 4.7 3.9 3.2 3.2 0.73 4.7 These results are not grid independent. It can be seen that the standard k-e and k-w models seriously underestimate the separation length. This is because this model predicts excessive turbulence production at the stagnation zone on the front of the plate, and the high turbulence levels are then convected into the separation zone. The other models perform much better, with the Chen-Kim and Realisable models showing close agreement with the measurement. ENDDIS AUTOPLOT USE file phida 3 da 1 u1 y 11 divide x .1 1 shift x -10 1 colf 1 level y 0;level x 0 scale x 0 5 redraw msg Velocity (U1) profile msg Press e to END ENDUSE ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(REALISABLE_KE SCAL_WF BLUNT PLATE:T311 ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 14 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 115 XULAST =2. XFRAC(1)=0.021251 ;XFRAC(24)=0.336972 XFRAC(47)=0.465068 ;XFRAC(70)=0.531403 XFRAC(93)=0.673179 ************************************************************ Group 4. Y-Direction Grid Spacing NY = 75 YVLAST =1. YFRAC(1)=4.164548E-03 ;YFRAC(16)=0.071363 YFRAC(31)=0.152963 ;YFRAC(46)=0.29992 YFRAC(61)=0.564581 ************************************************************ 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(12)=KE NAME(13)=EP ;NAME(136)=PRPS NAME(138)=SKIN ;NAME(139)=STRS NAME(140)=YPLS ;NAME(141)=C1E NAME(142)=DVDY ;NAME(143)=DVDX NAME(144)=DUDY ;NAME(145)=DUDX NAME(146)=EPKE ;NAME(147)=CMU 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,N,N,N,Y) SOLUTN(U1,Y,Y,N,N,N,Y) SOLUTN(V1,Y,Y,N,N,N,Y) SOLUTN(KE,Y,Y,N,N,N,N) SOLUTN(EP,Y,Y,N,N,N,N) SOLUTN(PRPS,Y,N,N,N,N,Y) SOLUTN(SKIN,Y,N,N,N,N,Y) SOLUTN(STRS,Y,N,N,N,N,Y) SOLUTN(YPLS,Y,N,N,N,N,Y) SOLUTN(C1E,Y,N,N,N,N,Y) SOLUTN(DVDY,Y,N,N,N,N,N) SOLUTN(DVDX,Y,N,N,N,N,N) SOLUTN(DUDY,Y,N,N,N,N,N) SOLUTN(DUDX,Y,N,N,N,N,N) SOLUTN(EPKE,Y,N,N,N,N,N) SOLUTN(CMU,Y,N,N,N,N,N) 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 VIST = 150 LEN1 = 149 PRPS = 136 ************************************************************ 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(KE,N,Y,Y,Y,Y,N) TERMS(EP,N,Y,Y,Y,Y,N) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = T NEWENT = T ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1 DENPCO = T ************************************************************ Group 9. Properties used if PRPS is not stored, and where PRPS = -1.0 if it is! RHO1 =1.225 ;TMP1 =0. EL1 = GRND4 TSURR =0. ;TEMP0 =0. PRESS0 =0. DVO1DT =0. ;DRH1DP =0. EMISS =0. ;SCATT =0. RADIA =0. ;RADIB =0. EL1A =0. ;EL1B =0. EL1C =0. ENUL =1.46E-05 ;ENUT = GRND5 ENUTA =0. ;ENUTB =0. ENUTC =0. IENUTA = 14 PRNDTL(U1)=1. ;PRNDTL(V1)=1. PRNDTL(KE)=1. ;PRNDTL(EP)=1. PRT(U1)=1. ;PRT(V1)=1. PRT(KE)=1. ;PRT(EP)=1.2 CP1 =1. ;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)=7.300001E-03 ;FIINIT(KE)=5.329E-03 FIINIT(EP)=1.27831E-03 ;FIINIT(PRPS)=-1. FIINIT(SKIN)=1.0E-10 ;FIINIT(STRS)=1.0E-10 FIINIT(YPLS)=1.0E-10 ;FIINIT(C1E)=1.0E-10 FIINIT(DVDY)=1.0E-10 ;FIINIT(DVDX)=1.0E-10 FIINIT(DUDY)=1.0E-10 ;FIINIT(DUDX)=1.0E-10 FIINIT(EPKE)=1.0E-10 ;FIINIT(CMU)=0.09 FIINIT(DEN1)=1.0E-10 ;FIINIT(EL1)=1.0E-10 FIINIT(ENUT)=1.0E-10 PATCH(PLATE ,INIVAL, 61, 115, 1, 10, 1, 1, 1, 1) INIT(PLATE ,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(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(KESOURCE,KE , GRND4 , GRND4 ) PATCH(REKESO ,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(REKESO ,EP , GRND4 , GRND4 ) PATCH(INLET ,WEST , 1, 1, 1, 75, 1, 1, 1, 1) COVAL(INLET ,P1 , FIXFLU ,8.9425 ) COVAL(INLET ,U1 ,0. ,7.3 ) COVAL(INLET ,V1 ,0. ,0. ) COVAL(INLET ,KE ,0. ,5.329E-03 ) COVAL(INLET ,EP ,0. ,1.27831E-03 ) PATCH(OUTLET ,EAST , 115, 115, 11, 75, 1, 1, 1, 1) COVAL(OUTLET ,P1 ,1000. ,0. ) COVAL(OUTLET ,U1 ,0. ,0. ) COVAL(OUTLET ,V1 ,0. ,0. ) XCYCLE = F EGWF = T WALLCO = GRND2 SCALWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 1500 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SELREF = T RESFAC =1.0E-04 ************************************************************ Group 16. Terminate Iterations LITER(P1)=20 ;LITER(U1)=10 LITER(V1)=10 ;LITER(KE)=20 LITER(EP)=20 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03 ENDIT(V1)=1.0E-03 ;ENDIT(KE)=1.0E-03 ENDIT(EP)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,1.) RELAX(U1,FALSDT,1.) RELAX(V1,FALSDT,1.) RELAX(KE,LINRLX,0.5) RELAX(EP,LINRLX,0.5) RELAX(PRPS,LINRLX,1.) RELAX(SKIN,LINRLX,1.) RELAX(STRS,LINRLX,1.) RELAX(YPLS,LINRLX,1.) RELAX(C1E,LINRLX,1.) RELAX(DVDY,LINRLX,1.) RELAX(DVDX,LINRLX,1.) RELAX(DUDY,LINRLX,1.) RELAX(DUDX,LINRLX,1.) RELAX(EPKE,LINRLX,1.) RELAX(CMU,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(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10 VARMAX(SKIN)=1.0E+10 ;VARMIN(SKIN)=-1.0E+10 VARMAX(STRS)=1.0E+10 ;VARMIN(STRS)=-1.0E+10 VARMAX(YPLS)=1.0E+10 ;VARMIN(YPLS)=-1.0E+10 VARMAX(C1E)=1.0E+10 ;VARMIN(C1E)=-1.0E+10 VARMAX(DVDY)=1.0E+10 ;VARMIN(DVDY)=-1.0E+10 VARMAX(DVDX)=1.0E+10 ;VARMIN(DVDX)=-1.0E+10 VARMAX(DUDY)=1.0E+10 ;VARMIN(DUDY)=-1.0E+10 VARMAX(DUDX)=1.0E+10 ;VARMIN(DUDX)=-1.0E+10 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10 VARMAX(CMU)=1.0E+10 ;VARMIN(CMU)=-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 GENK = T PARSOL = F CONWIZ = T IENUTA = 14 ISG62 = 1 SPEDAT(SET,KECONST,C2E,R,1.9) SPEDAT(SET,MATERIAL,198,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=5.421 ;EX(U1)=6.629 EX(V1)=0.4262 ;EX(KE)=0.123 EX(EP)=4.972 ;EX(PRPS)=0.9362 EX(SKIN)=5.051E-05 ;EX(STRS)=2.497E-04 EX(YPLS)=0.1502 ;EX(C1E)=0.5868 EX(DVDY)=5.134 ;EX(DVDX)=3.958 EX(DUDY)=18.360001 ;EX(DUDX)=5.133 EX(EPKE)=5.302 ;EX(CMU)=0.07094 EX(DEN1)=1.147 ;EX(EL1)=0.04079 EX(ENUT)=1.602E-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,Y,N,Y,Y,Y,Y) OUTPUT(U1,Y,N,Y,Y,Y,Y) OUTPUT(V1,Y,N,Y,Y,Y,Y) OUTPUT(KE,Y,N,Y,Y,Y,Y) OUTPUT(EP,Y,N,Y,Y,Y,Y) OUTPUT(PRPS,Y,N,Y,N,N,N) OUTPUT(SKIN,Y,N,Y,N,N,N) OUTPUT(STRS,Y,N,Y,N,N,N) OUTPUT(YPLS,Y,N,Y,N,N,N) OUTPUT(C1E,Y,N,Y,N,Y,Y) OUTPUT(DVDY,Y,N,N,N,N,N) OUTPUT(DVDX,Y,N,N,N,N,N) OUTPUT(DUDY,Y,N,N,N,N,N) OUTPUT(DUDX,Y,N,N,N,N,N) OUTPUT(EPKE,Y,N,N,N,N,N) OUTPUT(CMU,Y,N,N,N,Y,Y) OUTPUT(DEN1,Y,N,Y,N,N,N) OUTPUT(EL1,Y,N,Y,N,N,N) OUTPUT(ENUT,Y,N,Y,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 73 ;IYMON = 12 ;IZMON = 1 NPRMON = 100 ;NPRMNT = 1 ;TSTSWP = -1 UWATCH = T ;USTEER = T HIGHLO = F ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 ;NUMCLS = 5 NXPRIN = 23 ;IXPRF = 1 ;IXPRL = 10000 NYPRIN = 15 ;IYPRF = 1 ;IYPRL = 10000 IPLTF = 1 ;IPLTL = 1500 ;NPLT = 10 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