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, and four different versions of the
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:
KE MMK KL CK RKE RNG KW KWR MKW SST EXPT
Lr/H = 8.3 11.4 11.2 15.1 12.9 13.7 8.0 11.7 8.3 12.3 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, whilst the k-w-SST, the Wilcox
2008 k-w and realisable k-e models give 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 = 90
YVLAST =11.75
YFRAC(1)=3.039514E-03 ;YFRAC(5)=0.015198
YFRAC(9)=0.027356 ;YFRAC(13)=0.039514
YFRAC(17)=0.051672 ;YFRAC(21)=0.06383
YFRAC(25)=0.075988 ;YFRAC(29)=0.088815
YFRAC(33)=0.106864 ;YFRAC(37)=0.131418
YFRAC(41)=0.164825 ;YFRAC(45)=0.210274
YFRAC(49)=0.272107 ;YFRAC(53)=0.356231
YFRAC(57)=0.47068 ;YFRAC(61)=0.614427
YFRAC(65)=0.728876 ;YFRAC(69)=0.812999
YFRAC(73)=0.874832 ;YFRAC(77)=0.920281
YFRAC(81)=0.953688 ;YFRAC(85)=0.978243
YFRAC(89)=0.996291
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 110
ZWLAST =27.
ZFRAC(1)=0.01307 ;ZFRAC(5)=0.059415
ZFRAC(9)=0.097544 ;ZFRAC(13)=0.128912
ZFRAC(17)=0.154719 ;ZFRAC(21)=0.17595
ZFRAC(25)=0.193417 ;ZFRAC(29)=0.207787
ZFRAC(33)=0.21961 ;ZFRAC(37)=0.231481
ZFRAC(41)=0.243827 ;ZFRAC(45)=0.256173
ZFRAC(49)=0.26744 ;ZFRAC(53)=0.279957
ZFRAC(57)=0.294602 ;ZFRAC(61)=0.311733
ZFRAC(65)=0.331774 ;ZFRAC(69)=0.35522
ZFRAC(73)=0.382648 ;ZFRAC(77)=0.414735
ZFRAC(81)=0.452272 ;ZFRAC(85)=0.496185
ZFRAC(89)=0.547557 ;ZFRAC(93)=0.607655
ZFRAC(97)=0.677961 ;ZFRAC(101)=0.760209
ZFRAC(105)=0.856428 ;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(141)=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,Y,N,N,N)
SOLUTN(W1,Y,Y,Y,N,N,N)
SOLUTN(KE,Y,Y,Y,N,N,N)
SOLUTN(EP,Y,Y,Y,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 = 141
************************************************************
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)=2.133969E-03
FIINIT(EP)=3.063158E-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, 28, 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, 90, 1, 1, 1, 1)
COVAL(INLET ,P1 , FIXFLU ,1. )
COVAL(INLET ,V1 ,0. ,0. )
COVAL(INLET ,W1 ,0. ,1. )
COVAL(INLET ,KE ,0. ,2.133969E-03 )
COVAL(INLET ,EP ,0. ,3.063158E-05 )
PATCH(OUTL ,HIGH , 1, 1, 1, 90, 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, 90, 90, 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 = 1800 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
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,LINRLX,0.3)
RELAX(EP,LINRLX,0.3)
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,0.5)
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,OUTPUT,NOFIELD,L,T)
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.1512 ;EX(V1)=0.07132
EX(W1)=0.7815 ;EX(KE)=0.01017
EX(EP)=1.776E-03 ;EX(PRPS)=0.9661
EX(YPLS)=4.931 ;EX(C1E)=0.4664
EX(DWDZ)=0.07146 ;EX(DWDY)=0.2955
EX(DVDZ)=0.04759 ;EX(DVDY)=0.07292
EX(EPKE)=0.1117 ;EX(CMU)=0.1096
EX(ENUT)=0.01484
************************************************************
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,Y,N,N,N)
OUTPUT(DWDY,Y,N,Y,N,N,N)
OUTPUT(DVDZ,Y,N,Y,N,N,N)
OUTPUT(DVDY,Y,N,Y,N,N,N)
OUTPUT(EPKE,Y,N,Y,N,N,N)
OUTPUT(CMU,Y,N,Y,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 67 ;IYMON = 30 ;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 = 4 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = 4 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = 1800 ;NPLT = 25
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