TALK=T;RUN(1,1) DISPLAY The problem concerns 2D incompressible, laminar/turbulent flow over a back-facing step. The case is similar to that described for library case T103. The main objectives of it are: 1. Provide comparative test against calculations by stag- gered algorithm for the case set whether as Cartesian (CARTES= T, or cylindrical geometry (CARTES= F). 2. Provide test for the use of turbulence model with CCM- method (LTURB=T). The calculation is performed with either the Lam-Bremhorst low-Re K-E model (LTWOL=F) or the 2-layer low-Re K-E model (LTWOL=T), rather than with the standard high-Re form plus wall functions. User can switch from the default colocated computational algorithm (CCM) to the staggered one (STAG) by setting LCCM = F. While using CCM-method, it is possible to use higher order convective schemes (see settings at the end of Q1-file). A fully developed parabolic laminar flow profile can be prescribed at the inflow by setting LUNIF = F. Otherwise inflow is uniform. ---------------------------------------------------------- ENDDIS PHOTON USE p ; ; ; ; ; msg Computational Domain: gr i 1 msg Press Any Key to Continue... pause cl set vec av off msg Velocity Vectors: vec i 1 sh msg Press Any Key to Continue... pause cl msg Contours of Pressure: con p1 i 1 fi;0.1 pause cl msg Contours of W1-velocity: con w1 i 1 fi;0.1 pause cl msg Contours of V1-velocity: con v1 i 1 fi;0.1 msg Press Eto exit PHOTON ... ENDUSE SCHMBEGIN VARNAM VC1 SCHEME MINMOD VARNAM WC1 SCHEME MINMOD VARNAM KE SCHEME MINMOD VARNAM EP SCHEME MINMOD SCHMEND ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(CCM: Back-facing step (2 layer K-E). ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 150 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 1 XULAST =0.1 XFRAC(1)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 25 YVLAST =1.5 YFRAC(1)=0.01751 ;YFRAC(2)=0.04621 YFRAC(3)=0.081519 ;YFRAC(4)=0.121948 YFRAC(5)=0.166667 ;YFRAC(6)=0.211385 YFRAC(7)=0.251814 ;YFRAC(8)=0.287124 YFRAC(9)=0.315823 ;YFRAC(10)=0.333333 YFRAC(11)=0.353185 ;YFRAC(12)=0.385722 YFRAC(13)=0.425753 ;YFRAC(14)=0.471587 YFRAC(15)=0.522285 ;YFRAC(16)=0.577229 YFRAC(17)=0.635976 ;YFRAC(18)=0.697357 YFRAC(19)=0.756104 ;YFRAC(20)=0.811048 YFRAC(21)=0.861746 ;YFRAC(22)=0.907581 YFRAC(23)=0.947612 ;YFRAC(24)=0.980148 YFRAC(25)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 45 ZWLAST =22. ZFRAC(1)=0.02647 ;ZFRAC(2)=0.051226 ZFRAC(3)=0.07407 ;ZFRAC(4)=0.094722 ZFRAC(5)=0.112759 ;ZFRAC(6)=0.127419 ZFRAC(7)=0.136364 ;ZFRAC(8)=0.143995 ZFRAC(9)=0.155155 ;ZFRAC(10)=0.168196 ZFRAC(11)=0.182633 ;ZFRAC(12)=0.198204 ZFRAC(13)=0.214744 ;ZFRAC(14)=0.232136 ZFRAC(15)=0.250291 ;ZFRAC(16)=0.269142 ZFRAC(17)=0.288633 ;ZFRAC(18)=0.308718 ZFRAC(19)=0.329359 ;ZFRAC(20)=0.350524 ZFRAC(21)=0.372183 ;ZFRAC(22)=0.394311 ZFRAC(23)=0.416887 ;ZFRAC(24)=0.43989 ZFRAC(25)=0.463302 ;ZFRAC(26)=0.487109 ZFRAC(27)=0.511295 ;ZFRAC(28)=0.535846 ZFRAC(29)=0.56075 ;ZFRAC(30)=0.585997 ZFRAC(31)=0.611575 ;ZFRAC(32)=0.637475 ZFRAC(33)=0.663688 ;ZFRAC(34)=0.690205 ZFRAC(35)=0.717018 ;ZFRAC(36)=0.74412 ZFRAC(37)=0.771504 ;ZFRAC(38)=0.799164 ZFRAC(39)=0.827092 ;ZFRAC(40)=0.855283 ZFRAC(41)=0.883732 ;ZFRAC(42)=0.912433 ZFRAC(43)=0.941381 ;ZFRAC(44)=0.970572 ZFRAC(45)=1. ************************************************************ Group 6. Body-Fitted Coordinates CCM = T ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(1)=P1 ;NAME(5)=V1 NAME(7)=W1 ;NAME(12)=KE NAME(13)=EP ;NAME(143)=PRPS NAME(145)=WC1 ;NAME(146)=VC1 NAME(147)=LTLS ;NAME(148)=WDIS 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,Y) SOLUTN(V1,Y,Y,N,N,N,Y) SOLUTN(W1,Y,Y,N,N,N,Y) SOLUTN(KE,Y,Y,Y,N,N,N) SOLUTN(EP,Y,Y,Y,N,N,N) SOLUTN(PRPS,Y,N,N,N,N,N) SOLUTN(WC1,Y,Y,Y,N,N,Y) SOLUTN(VC1,Y,Y,Y,N,N,Y) SOLUTN(LTLS,Y,Y,Y,N,N,Y) SOLUTN(WDIS,Y,N,N,N,N,N) SOLUTN(LEN1,Y,N,N,N,N,Y) SOLUTN(ENUT,Y,N,N,N,N,Y) VIST = 150 LEN1 = 149 PRPS = 143 ************************************************************ 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,N,N,N,N,N,N) TERMS(W1,N,N,N,N,N,N) TERMS(KE,N,Y,Y,Y,Y,N) TERMS(EP,N,Y,Y,Y,Y,N) TERMS(WC1,N,Y,Y,Y,N,Y) TERMS(VC1,N,Y,Y,Y,Y,Y) TERMS(LTLS,N,N,Y,N,Y,Y) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = F 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 =1.189 ;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.444444E-04 ;ENUT = GRND3 ENUTA =0. ;ENUTB =0. ENUTC =0. IENUTA = 8 PRNDTL(KE)=1. ;PRNDTL(EP)=1. PRNDTL(WC1)=1. ;PRNDTL(VC1)=1. PRNDTL(LTLS)=1. PRT(KE)=1. ;PRT(EP)=1.314 PRT(WC1)=1. ;PRT(VC1)=1. PRT(LTLS)=1. CP1 =1. ;CP2 =1. ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(P1)=1.3E-04 ;FIINIT(V1)=1.0E-10 FIINIT(W1)=1.0E-10 ;FIINIT(KE)=1.8E-04 FIINIT(EP)=8.817263E-06 ;FIINIT(PRPS)=-1. FIINIT(WC1)=1.0E-05 ;FIINIT(VC1)=1.0E-05 FIINIT(LTLS)=1.0E-10 ;FIINIT(WDIS)=1.0E-02 FIINIT(LEN1)=1.0E-10 ;FIINIT(ENUT)=1.0E-10 PATCH(STEP ,INIVAL, 1, 1, 1, 10, 1, 7, 1, 1) INIT(STEP ,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(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(KESOURCE,KE , GRND4 , GRND4 ) COVAL(KESOURCE,EP , GRND4 , GRND4 ) PATCH(INL1 ,LOW , 1, 1, 11, 11, 1, 1, 1, 1) COVAL(INL1 ,P1 , FIXFLU ,15.457 ) COVAL(INL1 ,V1 ,0. ,0. ) COVAL(INL1 ,W1 ,0. ,0. ) COVAL(INL1 ,KE ,0. ,1.8E-04 ) COVAL(INL1 ,EP ,0. ,8.817263E-06 ) COVAL(INL1 ,WC1 ,0. ,13. ) COVAL(INL1 ,VC1 ,0. ,0. ) PATCH(INL2 ,LOW , 1, 1, 12, 12, 1, 1, 1, 1) COVAL(INL2 ,P1 , FIXFLU ,15.457 ) COVAL(INL2 ,V1 ,0. ,0. ) COVAL(INL2 ,W1 ,0. ,0. ) COVAL(INL2 ,KE ,0. ,1.8E-04 ) COVAL(INL2 ,EP ,0. ,8.817263E-06 ) COVAL(INL2 ,WC1 ,0. ,13. ) COVAL(INL2 ,VC1 ,0. ,0. ) PATCH(INL3 ,LOW , 1, 1, 13, 13, 1, 1, 1, 1) COVAL(INL3 ,P1 , FIXFLU ,15.457 ) COVAL(INL3 ,V1 ,0. ,0. ) COVAL(INL3 ,W1 ,0. ,0. ) COVAL(INL3 ,KE ,0. ,1.8E-04 ) COVAL(INL3 ,EP ,0. ,8.817263E-06 ) COVAL(INL3 ,WC1 ,0. ,13. ) COVAL(INL3 ,VC1 ,0. ,0. ) PATCH(INL4 ,LOW , 1, 1, 14, 14, 1, 1, 1, 1) COVAL(INL4 ,P1 , FIXFLU ,15.457 ) COVAL(INL4 ,V1 ,0. ,0. ) COVAL(INL4 ,W1 ,0. ,0. ) COVAL(INL4 ,KE ,0. ,1.8E-04 ) COVAL(INL4 ,EP ,0. ,8.817263E-06 ) COVAL(INL4 ,WC1 ,0. ,13. ) COVAL(INL4 ,VC1 ,0. ,0. ) PATCH(INL5 ,LOW , 1, 1, 15, 15, 1, 1, 1, 1) COVAL(INL5 ,P1 , FIXFLU ,15.457 ) COVAL(INL5 ,V1 ,0. ,0. ) COVAL(INL5 ,W1 ,0. ,0. ) COVAL(INL5 ,KE ,0. ,1.8E-04 ) COVAL(INL5 ,EP ,0. ,8.817263E-06 ) COVAL(INL5 ,WC1 ,0. ,13. ) COVAL(INL5 ,VC1 ,0. ,0. ) PATCH(INL6 ,LOW , 1, 1, 16, 16, 1, 1, 1, 1) COVAL(INL6 ,P1 , FIXFLU ,15.457 ) COVAL(INL6 ,V1 ,0. ,0. ) COVAL(INL6 ,W1 ,0. ,0. ) COVAL(INL6 ,KE ,0. ,1.8E-04 ) COVAL(INL6 ,EP ,0. ,8.817263E-06 ) COVAL(INL6 ,WC1 ,0. ,13. ) COVAL(INL6 ,VC1 ,0. ,0. ) PATCH(INL7 ,LOW , 1, 1, 17, 17, 1, 1, 1, 1) COVAL(INL7 ,P1 , FIXFLU ,15.457 ) COVAL(INL7 ,V1 ,0. ,0. ) COVAL(INL7 ,W1 ,0. ,0. ) COVAL(INL7 ,KE ,0. ,1.8E-04 ) COVAL(INL7 ,EP ,0. ,8.817263E-06 ) COVAL(INL7 ,WC1 ,0. ,13. ) COVAL(INL7 ,VC1 ,0. ,0. ) PATCH(INL8 ,LOW , 1, 1, 18, 18, 1, 1, 1, 1) COVAL(INL8 ,P1 , FIXFLU ,15.457 ) COVAL(INL8 ,V1 ,0. ,0. ) COVAL(INL8 ,W1 ,0. ,0. ) COVAL(INL8 ,KE ,0. ,1.8E-04 ) COVAL(INL8 ,EP ,0. ,8.817263E-06 ) COVAL(INL8 ,WC1 ,0. ,13. ) COVAL(INL8 ,VC1 ,0. ,0. ) PATCH(INL9 ,LOW , 1, 1, 19, 19, 1, 1, 1, 1) COVAL(INL9 ,P1 , FIXFLU ,15.457 ) COVAL(INL9 ,V1 ,0. ,0. ) COVAL(INL9 ,W1 ,0. ,0. ) COVAL(INL9 ,KE ,0. ,1.8E-04 ) COVAL(INL9 ,EP ,0. ,8.817263E-06 ) COVAL(INL9 ,WC1 ,0. ,13. ) COVAL(INL9 ,VC1 ,0. ,0. ) PATCH(INL10 ,LOW , 1, 1, 20, 20, 1, 1, 1, 1) COVAL(INL10 ,P1 , FIXFLU ,15.457 ) COVAL(INL10 ,V1 ,0. ,0. ) COVAL(INL10 ,W1 ,0. ,0. ) COVAL(INL10 ,KE ,0. ,1.8E-04 ) COVAL(INL10 ,EP ,0. ,8.817263E-06 ) COVAL(INL10 ,WC1 ,0. ,13. ) COVAL(INL10 ,VC1 ,0. ,0. ) PATCH(INL11 ,LOW , 1, 1, 21, 21, 1, 1, 1, 1) COVAL(INL11 ,P1 , FIXFLU ,15.457 ) COVAL(INL11 ,V1 ,0. ,0. ) COVAL(INL11 ,W1 ,0. ,0. ) COVAL(INL11 ,KE ,0. ,1.8E-04 ) COVAL(INL11 ,EP ,0. ,8.817263E-06 ) COVAL(INL11 ,WC1 ,0. ,13. ) COVAL(INL11 ,VC1 ,0. ,0. ) PATCH(INL12 ,LOW , 1, 1, 22, 22, 1, 1, 1, 1) COVAL(INL12 ,P1 , FIXFLU ,15.457 ) COVAL(INL12 ,V1 ,0. ,0. ) COVAL(INL12 ,W1 ,0. ,0. ) COVAL(INL12 ,KE ,0. ,1.8E-04 ) COVAL(INL12 ,EP ,0. ,8.817263E-06 ) COVAL(INL12 ,WC1 ,0. ,13. ) COVAL(INL12 ,VC1 ,0. ,0. ) PATCH(INL13 ,LOW , 1, 1, 23, 23, 1, 1, 1, 1) COVAL(INL13 ,P1 , FIXFLU ,15.457 ) COVAL(INL13 ,V1 ,0. ,0. ) COVAL(INL13 ,W1 ,0. ,0. ) COVAL(INL13 ,KE ,0. ,1.8E-04 ) COVAL(INL13 ,EP ,0. ,8.817263E-06 ) COVAL(INL13 ,WC1 ,0. ,13. ) COVAL(INL13 ,VC1 ,0. ,0. ) PATCH(INL14 ,LOW , 1, 1, 24, 24, 1, 1, 1, 1) COVAL(INL14 ,P1 , FIXFLU ,15.457 ) COVAL(INL14 ,V1 ,0. ,0. ) COVAL(INL14 ,W1 ,0. ,0. ) COVAL(INL14 ,KE ,0. ,1.8E-04 ) COVAL(INL14 ,EP ,0. ,8.817263E-06 ) COVAL(INL14 ,WC1 ,0. ,13. ) COVAL(INL14 ,VC1 ,0. ,0. ) PATCH(INL15 ,LOW , 1, 1, 25, 25, 1, 1, 1, 1) COVAL(INL15 ,P1 , FIXFLU ,15.457 ) COVAL(INL15 ,V1 ,0. ,0. ) COVAL(INL15 ,W1 ,0. ,0. ) COVAL(INL15 ,KE ,0. ,1.8E-04 ) COVAL(INL15 ,EP ,0. ,8.817263E-06 ) COVAL(INL15 ,WC1 ,0. ,13. ) COVAL(INL15 ,VC1 ,0. ,0. ) PATCH(WS1 ,SWALL , 1, 1, 11, 11, 1, 7, 1, 1) COVAL(WS1 ,KE ,1. ,0. ) COVAL(WS1 ,WC1 , GRND2 ,0. ) COVAL(WS1 ,VC1 , GRND2 ,0. ) COVAL(WS1 ,LTLS,1. ,0. ) PATCH(WN ,NWALL , 1, 1, 25, 25, 1, 45, 1, 1) COVAL(WN ,KE ,1. ,0. ) COVAL(WN ,WC1 , GRND2 ,0. ) COVAL(WN ,VC1 , GRND2 ,0. ) COVAL(WN ,LTLS,1. ,0. ) PATCH(WL ,LWALL , 1, 1, 1, 10, 8, 8, 1, 1) COVAL(WL ,KE ,1. ,0. ) COVAL(WL ,WC1 , GRND2 ,0. ) COVAL(WL ,VC1 , GRND2 ,0. ) COVAL(WL ,LTLS,1. ,0. ) PATCH(WS2 ,SWALL , 1, 1, 1, 1, 8, 45, 1, 1) COVAL(WS2 ,KE ,1. ,0. ) COVAL(WS2 ,WC1 , GRND2 ,0. ) COVAL(WS2 ,VC1 , GRND2 ,0. ) COVAL(WS2 ,LTLS,1. ,0. ) PATCH(OUT1 ,HIGH , 1, 1, 1, 25, 45, 45, 1, 1) COVAL(OUT1 ,P1 ,1000. ,0. ) COVAL(OUT1 ,KE ,0. ,0. ) COVAL(OUT1 ,EP ,0. ,0. ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 250 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SELREF = T RESFAC =1.0E-03 ************************************************************ Group 16. Terminate Iterations LITER(P1)=20 ;LITER(V1)=10 LITER(W1)=10 ;LITER(KE)=20 LITER(EP)=20 ;LITER(WC1)=20 LITER(VC1)=20 ;LITER(LTLS)=20 ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03 ENDIT(W1)=1.0E-03 ;ENDIT(KE)=1.0E-03 ENDIT(EP)=1.0E-03 ;ENDIT(WC1)=1.0E-03 ENDIT(VC1)=1.0E-03 ;ENDIT(LTLS)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,0.5) RELAX(V1,FALSDT,1.) RELAX(W1,FALSDT,1.) RELAX(KE,FALSDT,0.037607) RELAX(EP,FALSDT,0.037607) RELAX(PRPS,LINRLX,1.) RELAX(WC1,FALSDT,0.376068) RELAX(VC1,FALSDT,0.376068) RELAX(LTLS,FALSDT,1.0E+09) RELAX(WDIS,LINRLX,1.) RELAX(LEN1,LINRLX,1.) RELAX(ENUT,LINRLX,1.) KELIN = 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(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10 VARMAX(WC1)=1.0E+10 ;VARMIN(WC1)=-1.0E+10 VARMAX(VC1)=1.0E+10 ;VARMIN(VC1)=-1.0E+10 VARMAX(LTLS)=1.0E+10 ;VARMIN(LTLS)=-1.0E+10 VARMAX(WDIS)=1.0E+10 ;VARMIN(WDIS)=-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 PARSOL = F CONWIZ = T IENUTA = 8 ISG62 = 1 CSG3 ='LCRU' 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(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,Y,Y,Y) OUTPUT(EP,Y,N,Y,Y,Y,Y) OUTPUT(PRPS,Y,N,Y,N,N,N) OUTPUT(WC1,Y,N,Y,Y,Y,Y) OUTPUT(VC1,Y,N,Y,Y,Y,Y) OUTPUT(LTLS,Y,N,Y,Y,Y,Y) OUTPUT(WDIS,Y,N,N,N,N,N) OUTPUT(LEN1,Y,N,Y,N,N,N) OUTPUT(ENUT,Y,N,Y,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 13 ;IZMON = 23 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 = 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