TALK=T;RUN(1,1)
DISPLAY
The problem considered is plane turbulent couette flow in a
channel with one moving wall, as simulated for library case T100
using various high-Re turbulence models with wall functions. Here,
the turbulence is simulated by use of various low-Re k-e and
k-w models. The calculation integrates down to the wall and the
solution is performed by use of the single-slab solver. A non-
uniform grid is employed so as to concentrate cells very close to
the walls. The Reynolds number is 1.E5 based on the channel height
and the average velocity.
For this flow, the shear stress is uniform across the flow and
hence the axial pressure-gradient is zero. The velocity profile
is S-shaped and symmetrical about the central plane, and so the
average velocity is one half of the velocity of the moving wall.
ENDDIS
AUTOPLOT USE
file
phi 5
da 1 w1;col9 1
msg Velocity (W1) profile
msg Press RETURN to continue
pause
clear
da 1 ke;col9 1
msg KE profile
msg Press RETURN to continue
pause
clear
da 1 enut;col9 1
msg ENUT profile
msg Press e to END
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(LAM-BRE KE_1D PLANE COUETTE FLOW :T205 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 205
************************************************************
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 = 132
YVLAST =0.1
YFRAC(1)=2.55102E-04 ;YFRAC(2)=5.306122E-04
YFRAC(3)=8.281632E-04 ;YFRAC(4)=1.149518E-03
YFRAC(5)=1.496582E-03 ;YFRAC(6)=1.87141E-03
YFRAC(7)=2.276225E-03 ;YFRAC(8)=2.713425E-03
YFRAC(9)=3.185601E-03 ;YFRAC(10)=3.695552E-03
YFRAC(11)=4.246298E-03 ;YFRAC(12)=4.841104E-03
YFRAC(13)=5.483495E-03 ;YFRAC(14)=6.177276E-03
YFRAC(15)=6.926561E-03 ;YFRAC(16)=7.735788E-03
YFRAC(17)=8.609753E-03 ;YFRAC(18)=9.553635E-03
YFRAC(19)=0.010573 ;YFRAC(20)=0.011674
YFRAC(21)=0.012863 ;YFRAC(22)=0.014147
YFRAC(23)=0.015534 ;YFRAC(24)=0.017032
YFRAC(25)=0.018649 ;YFRAC(26)=0.020397
YFRAC(27)=0.022283 ;YFRAC(28)=0.024321
YFRAC(29)=0.026522 ;YFRAC(30)=0.028899
YFRAC(31)=0.031466 ;YFRAC(32)=0.034238
YFRAC(33)=0.037232 ;YFRAC(34)=0.040466
YFRAC(35)=0.043958 ;YFRAC(36)=0.04773
YFRAC(37)=0.051804 ;YFRAC(38)=0.056203
YFRAC(39)=0.060954 ;YFRAC(40)=0.066086
YFRAC(41)=0.071628 ;YFRAC(42)=0.077613
YFRAC(43)=0.084077 ;YFRAC(44)=0.091059
YFRAC(45)=0.098599 ;YFRAC(46)=0.106741
YFRAC(47)=0.115536 ;YFRAC(48)=0.125034
YFRAC(49)=0.135292 ;YFRAC(50)=0.14637
YFRAC(51)=0.158335 ;YFRAC(52)=0.171257
YFRAC(53)=0.185212 ;YFRAC(54)=0.200285
YFRAC(55)=0.216562 ;YFRAC(56)=0.234143
YFRAC(57)=0.253129 ;YFRAC(58)=0.273634
YFRAC(59)=0.29578 ;YFRAC(60)=0.319698
YFRAC(61)=0.345529 ;YFRAC(62)=0.373426
YFRAC(63)=0.403556 ;YFRAC(64)=0.436095
YFRAC(65)=0.471238 ;YFRAC(66)=0.5
YFRAC(67)=0.528762 ;YFRAC(68)=0.563905
YFRAC(69)=0.596444 ;YFRAC(70)=0.626574
YFRAC(71)=0.654471 ;YFRAC(72)=0.680302
YFRAC(73)=0.70422 ;YFRAC(74)=0.726366
YFRAC(75)=0.746871 ;YFRAC(76)=0.765857
YFRAC(77)=0.783438 ;YFRAC(78)=0.799715
YFRAC(79)=0.814788 ;YFRAC(80)=0.828743
YFRAC(81)=0.841665 ;YFRAC(82)=0.85363
YFRAC(83)=0.864708 ;YFRAC(84)=0.874966
YFRAC(85)=0.884464 ;YFRAC(86)=0.893259
YFRAC(87)=0.901402 ;YFRAC(88)=0.908941
YFRAC(89)=0.915923 ;YFRAC(90)=0.922387
YFRAC(91)=0.928372 ;YFRAC(92)=0.933914
YFRAC(93)=0.939046 ;YFRAC(94)=0.943797
YFRAC(95)=0.948196 ;YFRAC(96)=0.95227
YFRAC(97)=0.956042 ;YFRAC(98)=0.959534
YFRAC(99)=0.962768 ;YFRAC(100)=0.965762
YFRAC(101)=0.968534 ;YFRAC(102)=0.971101
YFRAC(103)=0.973478 ;YFRAC(104)=0.975679
YFRAC(105)=0.977717 ;YFRAC(106)=0.979603
YFRAC(107)=0.981351 ;YFRAC(108)=0.982968
YFRAC(109)=0.984466 ;YFRAC(110)=0.985853
YFRAC(111)=0.987137 ;YFRAC(112)=0.988326
YFRAC(113)=0.989427 ;YFRAC(114)=0.990446
YFRAC(115)=0.99139 ;YFRAC(116)=0.992264
YFRAC(117)=0.993073 ;YFRAC(118)=0.993823
YFRAC(119)=0.994516 ;YFRAC(120)=0.995159
YFRAC(121)=0.995754 ;YFRAC(122)=0.996304
YFRAC(123)=0.996814 ;YFRAC(124)=0.997287
YFRAC(125)=0.997724 ;YFRAC(126)=0.998129
YFRAC(127)=0.998503 ;YFRAC(128)=0.99885
YFRAC(129)=0.999172 ;YFRAC(130)=0.999469
YFRAC(131)=0.999745 ;YFRAC(132)=1.
************************************************************
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(7)=W1 ;NAME(12)=KE
NAME(13)=EP ;NAME(140)=FTWO
NAME(141)=FONE ;NAME(142)=REYN
NAME(143)=REYT ;NAME(144)=FMU
NAME(145)=LTLS ;NAME(146)=WDIS
NAME(147)=CF ;NAME(148)=STRS
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(W1,Y,Y,N,N,N,N)
SOLUTN(KE,Y,Y,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
SOLUTN(FTWO,Y,N,N,N,N,Y)
SOLUTN(FONE,Y,N,N,N,N,Y)
SOLUTN(REYN,Y,N,N,N,N,Y)
SOLUTN(REYT,Y,N,N,N,N,Y)
SOLUTN(FMU,Y,N,N,N,N,Y)
SOLUTN(LTLS,Y,Y,Y,N,N,Y)
SOLUTN(WDIS,Y,N,N,N,N,N)
SOLUTN(CF,Y,N,N,N,N,N)
SOLUTN(STRS,Y,N,N,N,N,Y)
SOLUTN(LEN1,Y,N,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
VIST = 150
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(W1,N,N,Y,Y,Y,Y)
TERMS(KE,Y,N,Y,Y,Y,N)
TERMS(EP,Y,N,Y,Y,Y,N)
TERMS(LTLS,N,N,Y,N,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWENL = T
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 = GRND1 ;EL1 = GRND4
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 = GRND1 ;ENUT = GRND3
ENULA =5.0E-07 ;ENULB =0. ;ENULC =0. ;ENULD =0.
ENULE =0. ;ENULF =0. ;ENULG =0.
IENULA = 0 ;IENULB = 0
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 3
PRNDTL(W1)=1. ;PRNDTL(KE)=1.
PRNDTL(EP)=1. ;PRNDTL(LTLS)=1.
PRT(W1)=1. ;PRT(KE)=1.
PRT(EP)=1.314 ;PRT(LTLS)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(W1)=1.0E-10 ;FIINIT(KE)=1.280533E-03
FIINIT(EP)=1.673059E-03 ;FIINIT(FTWO)=1.0E-10
FIINIT(FONE)=1.0E-10 ;FIINIT(REYN)=1.0E-10
FIINIT(REYT)=1.0E-10 ;FIINIT(FMU)=1.0E-10
FIINIT(LTLS)=1.0E-10 ;FIINIT(WDIS)=1.0E-02
FIINIT(CF)=1.0E-10 ;FIINIT(STRS)=1.0E-10
FIINIT(LEN1)=1.0E-10 ;FIINIT(ENUT)=1.0E-10
PATCH(ICOUF ,LINVLY, 1, 1, 1, 132, 1, 1, 1, 1)
INIT(ICOUF ,W1 ,10. ,0. )
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(WALLN ,NWALL , 1, 1, 132, 132, 1, 1, 1, 1)
COVAL(WALLN ,W1 , GRND2 ,1. )
COVAL(WALLN ,KE ,1. ,0. )
PATCH(WALLS ,SWALL , 1, 1, 1, 1, 1, 1, 1, 1)
COVAL(WALLS ,W1 , GRND2 ,0. )
COVAL(WALLS ,KE ,1. ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 80 ;ISWC1 = 1
LITHYD = 6 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(W1)=10 ;LITER(KE)=20
LITER(EP)=20 ;LITER(LTLS)=20
ENDIT(W1)=1.0E-03 ;ENDIT(KE)=1.0E-03
ENDIT(EP)=1.0E-03 ;ENDIT(LTLS)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(W1,FALSDT,0.2)
RELAX(KE,FALSDT,0.05)
RELAX(EP,FALSDT,0.05)
RELAX(FTWO,LINRLX,1.)
RELAX(FONE,LINRLX,1.)
RELAX(REYN,LINRLX,1.)
RELAX(REYT,LINRLX,1.)
RELAX(FMU,LINRLX,1.)
RELAX(LTLS,FALSDT,1.0E+09)
RELAX(WDIS,LINRLX,1.)
RELAX(CF,LINRLX,1.)
RELAX(STRS,LINRLX,1.)
RELAX(LEN1,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
KELIN = 1
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=1.0E-10
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
VARMAX(FTWO)=1.0E+10 ;VARMIN(FTWO)=-1.0E+10
VARMAX(FONE)=1.0E+10 ;VARMIN(FONE)=-1.0E+10
VARMAX(REYN)=1.0E+10 ;VARMIN(REYN)=-1.0E+10
VARMAX(REYT)=1.0E+10 ;VARMIN(REYT)=-1.0E+10
VARMAX(FMU)=1.0E+10 ;VARMIN(FMU)=-1.0E+10
VARMAX(LTLS)=1.0E+10 ;VARMIN(LTLS)=-1.0E+10
VARMAX(WDIS)=1.0E+10 ;VARMIN(WDIS)=-1.0E+10
VARMAX(CF)=1.0E+10 ;VARMIN(CF)=-1.0E+10
VARMAX(STRS)=1.0E+10 ;VARMIN(STRS)=-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
IENUTA = 3
ISG62 = 1
SPEDAT(SET,STORED,CF,C,=2.*STRS/(RHO1*0.5*0.5)!ZSLFIN)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(W1)=0.5 ;EX(KE)=1.589E-03
EX(EP)=0.0169 ;EX(FTWO)=0.9653
EX(FONE)=2.24 ;EX(REYN)=773.700012
EX(REYT)=1833. ;EX(FMU)=0.7624
EX(LTLS)=3.596E-04 ;EX(WDIS)=9.707E-03
EX(CF)=5.737E-05 ;EX(STRS)=7.172E-06
EX(LEN1)=3.769E-03 ;EX(ENUT)=8.077E-05
************************************************************
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(W1,Y,N,Y,Y,Y,Y)
OUTPUT(KE,Y,N,Y,Y,Y,Y)
OUTPUT(EP,Y,N,Y,Y,Y,Y)
OUTPUT(FTWO,Y,N,Y,N,N,N)
OUTPUT(FONE,Y,N,Y,N,N,N)
OUTPUT(REYN,Y,N,Y,N,N,N)
OUTPUT(REYT,Y,N,Y,N,N,N)
OUTPUT(FMU,Y,N,Y,N,N,N)
OUTPUT(LTLS,Y,N,Y,Y,Y,Y)
OUTPUT(WDIS,Y,N,Y,N,N,N)
OUTPUT(CF,Y,N,Y,N,N,N)
OUTPUT(STRS,Y,N,Y,N,N,N)
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 = 2 ;IZMON = 1
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
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