TALK=T;RUN(1,1)
DISPLAY
The problem considered is turbulent couette flow in a plane
channel with one moving wall. 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. In the calculation made here
the Reynolds number is 1.E5 based on the channel height and
the average velocity.
ENDDIS
AUTOPLOT USE
file; phida 3
cl; da 1 w1; col9 1
msg Velocity (W1) profile; Press RETURN to continue
pause; cl; da 1 ke; col9 1; scale y 0. 2.e-3; redraw
msg KE profile. Press RETURN, then e to END
pause
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(CHEN-KIM K-E_1D PLANE COUETTE FLOW :T100)
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 100
************************************************************
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 = 30
YVLAST =0.1
YFRAC(1)=0.033333 ;YFRAC(2)=0.066667
YFRAC(3)=0.1 ;YFRAC(4)=0.133333
YFRAC(5)=0.166667 ;YFRAC(6)=0.2
YFRAC(7)=0.233333 ;YFRAC(8)=0.266667
YFRAC(9)=0.3 ;YFRAC(10)=0.333333
YFRAC(11)=0.366667 ;YFRAC(12)=0.4
YFRAC(13)=0.433333 ;YFRAC(14)=0.466667
YFRAC(15)=0.5 ;YFRAC(16)=0.533333
YFRAC(17)=0.566667 ;YFRAC(18)=0.6
YFRAC(19)=0.633333 ;YFRAC(20)=0.666667
YFRAC(21)=0.7 ;YFRAC(22)=0.733333
YFRAC(23)=0.766667 ;YFRAC(24)=0.8
YFRAC(25)=0.833333 ;YFRAC(26)=0.866667
YFRAC(27)=0.9 ;YFRAC(28)=0.933333
YFRAC(29)=0.966667 ;YFRAC(30)=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(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(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)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
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 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =5.0E-07 ;ENUT =0.
PRNDTL(W1)=1.
PRT(W1)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(W1)=1.0E-10 ;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, 30, 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(WALLN ,NWALL , 1, 1, 30, 30, 1, 1, 1, 1)
COVAL(WALLN ,W1 , GRND2 ,1. )
PATCH(WALLS ,SWALL , 1, 1, 1, 1, 1, 1, 1, 1)
COVAL(WALLS ,W1 ,1. ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 10 ;ISWC1 = 1
LITHYD = 20 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(W1)=10
ENDIT(W1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(W1,FALSDT,40.)
RELAX(CF,LINRLX,1.)
RELAX(STRS,LINRLX,1.)
RELAX(LEN1,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=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
PARSOL = F
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. ;EX(CF)=0.
EX(STRS)=0. ;EX(LEN1)=0.
EX(ENUT)=0.
************************************************************
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(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 = 28 ;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 = 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