TALK=T;RUN(1,1)
  DISPLAY
  The case considered is 2D, steady, incompressible, turbulent flow
  past a surface-mounted square rib in a channel. This case has
  been studied experimentally by D.Crabb et al, Proc. 4th Brazilian
  Congress on Mech. Engng., Florianopolis, Brazil, p415,(1997).

  The height H of the rib is 8.5% of that of the channel. The flow
  Reynolds number based on channel bulk velocity and rib height
  H is 300,000. The inlet plane is located 6H upstream of the rib,
  and the outlet plane 20H downstream of the rib. A fixed-pressure
  boundary condition is applied at the outlet, and uniform flow
  profiles are specified at the inlet.

  The case is set up to run one of six versions of the k-e model
  with scalable wall functions, namely the standard model, the MMK
  model, the Kato-Launder model, the Chen-Kim & RNG models, or
  the realisable k-e model. An option is provided to also run the
  standard k-w model.

  For this case, the main parameter characterising separation is
  the length of the separation zone behind the rib. The experimental
  and computed results are:

             K-E  K-O  KL    MMK   RKE   Chen  RNG   EXPT

     Lr/H =  5.2  6.3 11.4   11.7  12.8  14.1  13.0  12.3

  where the separation length Lr is measured from the front of the
  rib. These results are not grid independent, and the mesh is not
  fine enough to resolve the separation regions around the rib.
  The standard k-e & k-w models seriously underpredict the length
  of the separation behind the rib, and the realisable k-e model
  gives closest agreement with the data.
  ENDDIS
   AUTOPLOT USE
    AUTOPLOT
    FILE
    PHIDA 3

    D 1 W1 Y 1
    PLOT
    LEVEL Y 0
    SHIFT X -7 1
    REDR
    SCALE X 0 15
    msg Press e to END
  ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(RK K-E SQUARE RIB FLOW :T307            )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  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 = 64
 YVLAST =11.75
 YFRAC(1)=4.255319E-03 ;YFRAC(3)=0.012766
 YFRAC(5)=0.021277 ;YFRAC(7)=0.029787
 YFRAC(9)=0.038298 ;YFRAC(11)=0.046809
 YFRAC(13)=0.055319 ;YFRAC(15)=0.06383
 YFRAC(17)=0.07234 ;YFRAC(19)=0.080851
 YFRAC(21)=0.088544 ;YFRAC(23)=0.096159
 YFRAC(25)=0.104877 ;YFRAC(27)=0.114858
 YFRAC(29)=0.126285 ;YFRAC(31)=0.139369
 YFRAC(33)=0.154348 ;YFRAC(35)=0.171497
 YFRAC(37)=0.191131 ;YFRAC(39)=0.213611
 YFRAC(41)=0.239348 ;YFRAC(43)=0.268814
 YFRAC(45)=0.302549 ;YFRAC(47)=0.341173
 YFRAC(49)=0.385394 ;YFRAC(51)=0.436022
 YFRAC(53)=0.493986 ;YFRAC(55)=0.560349
 YFRAC(57)=0.636328 ;YFRAC(59)=0.723317
 YFRAC(61)=0.82291 ;YFRAC(63)=0.936934
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 110
 ZWLAST =27.
 ZFRAC(1)=0.01307 ;ZFRAC(3)=0.037372
 ZFRAC(5)=0.059415 ;ZFRAC(7)=0.079409
 ZFRAC(9)=0.097544 ;ZFRAC(11)=0.113992
 ZFRAC(13)=0.128912 ;ZFRAC(15)=0.142445
 ZFRAC(17)=0.154719 ;ZFRAC(19)=0.165852
 ZFRAC(21)=0.17595 ;ZFRAC(23)=0.18511
 ZFRAC(25)=0.193417 ;ZFRAC(27)=0.200953
 ZFRAC(29)=0.207787 ;ZFRAC(31)=0.213987
 ZFRAC(33)=0.21961 ;ZFRAC(35)=0.225309
 ZFRAC(37)=0.231481 ;ZFRAC(39)=0.237654
 ZFRAC(41)=0.243827 ;ZFRAC(43)=0.25
 ZFRAC(45)=0.256173 ;ZFRAC(47)=0.26188
 ZFRAC(49)=0.26744 ;ZFRAC(51)=0.273453
 ZFRAC(53)=0.279957 ;ZFRAC(55)=0.286992
 ZFRAC(57)=0.294602 ;ZFRAC(59)=0.302831
 ZFRAC(61)=0.311733 ;ZFRAC(63)=0.321361
 ZFRAC(65)=0.331774 ;ZFRAC(67)=0.343038
 ZFRAC(69)=0.35522 ;ZFRAC(71)=0.368396
 ZFRAC(73)=0.382648 ;ZFRAC(75)=0.398062
 ZFRAC(77)=0.414735 ;ZFRAC(79)=0.432767
 ZFRAC(81)=0.452272 ;ZFRAC(83)=0.473367
 ZFRAC(85)=0.496185 ;ZFRAC(87)=0.520864
 ZFRAC(89)=0.547557 ;ZFRAC(91)=0.576428
 ZFRAC(93)=0.607655 ;ZFRAC(95)=0.64143
 ZFRAC(97)=0.677961 ;ZFRAC(99)=0.717473
 ZFRAC(101)=0.760209 ;ZFRAC(103)=0.806433
 ZFRAC(105)=0.856428 ;ZFRAC(107)=0.910503
 ZFRAC(109)=0.96899
 ************************************************************
  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(140)=PRPS
 NAME(142)=YPLS ;NAME(143)=C1E
 NAME(144)=DWDZ ;NAME(145)=DWDY
 NAME(146)=DVDZ ;NAME(147)=DVDY
 NAME(148)=EPKE ;NAME(149)=CMU
 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,N)
 SOLUTN(W1,Y,Y,N,N,N,N)
 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(YPLS,Y,N,N,N,N,Y)
 SOLUTN(C1E,Y,N,N,N,N,Y)
 SOLUTN(DWDZ,Y,N,N,N,N,N)
 SOLUTN(DWDY,Y,N,N,N,N,N)
 SOLUTN(DVDZ,Y,N,N,N,N,N)
 SOLUTN(DVDY,Y,N,N,N,N,N)
 SOLUTN(EPKE,Y,N,N,N,N,N)
 SOLUTN(CMU,Y,N,N,N,N,N)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 VIST = 150
 PRPS = 140
 ************************************************************
  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(KE,N,Y,Y,Y,Y,N)
 TERMS(EP,N,Y,Y,Y,Y,N)
 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. ;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 =3.333333E-06 ;ENUT = GRND5
 ENUTA =0. ;ENUTB =0.
 ENUTC =0.
 IENUTA = 14
 PRNDTL(V1)=1. ;PRNDTL(W1)=1.
 PRNDTL(KE)=1. ;PRNDTL(EP)=1.
 PRT(V1)=1. ;PRT(W1)=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.3E-04 ;FIINIT(V1)=1.0E-03
 FIINIT(W1)=1. ;FIINIT(KE)=3.5E-03
 FIINIT(EP)=3.217061E-05 ;FIINIT(PRPS)=-1.
 FIINIT(YPLS)=1.0E-10 ;FIINIT(C1E)=1.0E-10
 FIINIT(DWDZ)=1.0E-10 ;FIINIT(DWDY)=1.0E-10
 FIINIT(DVDZ)=1.0E-10 ;FIINIT(DVDY)=1.0E-10
 FIINIT(EPKE)=1.0E-10 ;FIINIT(CMU)=0.09
 FIINIT(ENUT)=1.0E-10
 
 PATCH(RIB ,INIVAL, 1, 1, 1, 20, 35, 46, 1, 1)
 INIT(RIB ,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 ,LOW , 1, 1, 1, 64, 1, 1, 1, 1)
 COVAL(INLET ,P1 , FIXFLU ,1. )
 COVAL(INLET ,V1 ,0. ,0. )
 COVAL(INLET ,W1 ,0. ,1. )
 COVAL(INLET ,KE ,0. ,3.5E-03 )
 COVAL(INLET ,EP ,0. ,3.217061E-05 )
 
 PATCH(OUTL ,HIGH , 1, 1, 1, 64, 110, 110, 1, 1)
 COVAL(OUTL ,P1 ,1000. ,0. )
 COVAL(OUTL ,V1 ,0. ,0. )
 COVAL(OUTL ,W1 ,0. ,0. )
 COVAL(OUTL ,KE ,0. ,0. )
 COVAL(OUTL ,EP ,0. ,0. )
 
 PATCH(WALLN ,NWALL , 1, 1, 64, 64, 1, 110, 1, 1)
 COVAL(WALLN ,W1 , GRND2 ,0. )
 COVAL(WALLN ,KE , GRND2 , GRND2 )
 COVAL(WALLN ,EP , GRND2 , GRND2 )
 
 PATCH(WALLS ,SWALL , 1, 1, 1, 1, 1, 110, 1, 1)
 COVAL(WALLS ,W1 , GRND2 ,0. )
 COVAL(WALLS ,KE , GRND2 , GRND2 )
 COVAL(WALLS ,EP , GRND2 , GRND2 )
 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)=50 ;LITER(V1)=10
 LITER(W1)=10 ;LITER(KE)=5
 LITER(EP)=5
 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
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(V1,FALSDT,0.122727)
 RELAX(W1,FALSDT,0.122727)
 RELAX(KE,FALSDT,0.122727)
 RELAX(EP,FALSDT,0.122727)
 RELAX(PRPS,LINRLX,1.)
 RELAX(YPLS,LINRLX,1.)
 RELAX(C1E,LINRLX,1.)
 RELAX(DWDZ,LINRLX,1.)
 RELAX(DWDY,LINRLX,1.)
 RELAX(DVDZ,LINRLX,1.)
 RELAX(DVDY,LINRLX,1.)
 RELAX(EPKE,LINRLX,1.)
 RELAX(CMU,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(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(YPLS)=1.0E+10 ;VARMIN(YPLS)=-1.0E+10
 VARMAX(C1E)=1.0E+10 ;VARMIN(C1E)=-1.0E+10
 VARMAX(DWDZ)=1.0E+10 ;VARMIN(DWDZ)=-1.0E+10
 VARMAX(DWDY)=1.0E+10 ;VARMIN(DWDY)=-1.0E+10
 VARMAX(DVDZ)=1.0E+10 ;VARMIN(DVDZ)=-1.0E+10
 VARMAX(DVDY)=1.0E+10 ;VARMIN(DVDY)=-1.0E+10
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 VARMAX(CMU)=1.0E+10 ;VARMIN(CMU)=-1.0E+10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 GENK = T
 PARSOL = F
 IENUTA = 14
 ISG62 = 1
 SPEDAT(SET,KECONST,C2E,R,1.9)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 SPEDAT(SET,MATERIAL,198,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=0.1866 ;EX(V1)=0.1027
 EX(W1)=0.7593 ;EX(KE)=0.01526
 EX(EP)=2.468E-03 ;EX(PRPS)=0.9659
 EX(YPLS)=73.370003 ;EX(C1E)=0.5111
 EX(DWDZ)=0.08975 ;EX(DWDY)=0.28
 EX(DVDZ)=0.06251 ;EX(DVDY)=0.09097
 EX(EPKE)=0.09428 ;EX(CMU)=0.08361
 EX(ENUT)=0.02273
 ************************************************************
  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(YPLS,Y,N,Y,N,N,N)
 OUTPUT(C1E,Y,N,Y,N,N,N)
 OUTPUT(DWDZ,Y,N,N,N,N,N)
 OUTPUT(DWDY,Y,N,N,N,N,N)
 OUTPUT(DVDZ,Y,N,N,N,N,N)
 OUTPUT(DVDY,Y,N,N,N,N,N)
 OUTPUT(EPKE,Y,N,N,N,N,N)
 OUTPUT(CMU,Y,N,N,N,N,N)
 OUTPUT(ENUT,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 16 ;IZMON = 56
 NPRMON = 100 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NYPRIN = 2 ;IYPRF = 1 ;IYPRL = 10000
 NZPRIN = 2 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = F
 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