AUTOPLOT USE
  FI
  PHI 5
 
  cl
  DA 1 W1;PL
  msg Profile of velocity, W1
  msg            -
  msg Press  to continue
  pause
  cl
  DA 1 U2RS;DA 1 V2RS;DA 1 VWRS
  scale
  col4 1
  col8 2
  cole 3
  msg Profiles of U2RS - blue
  msg             V2RS - green
  msg             VWRS - red
  msg            -
  msg Press e to end
  ENDUSE
 
    GROUP 1. Run title and other preliminaries
TEXT(RSTM_1DY PLANE COUETTE FLOW        :T606
TITLE
mesg(PC486/50 time last reported as 1.min
 
  DISPLAY
  The problem considered is plane turbulent couette flow in a
  channel at a Reynolds number of 1.E5, as described in detail
  for library case T100.
  ENDDIS
 
  The calculations are performed with the RSTM plus wall functions,
  and like the k-e model they show that the turbulent kinetic
  energy k is essentially uniform across the flow, whereas the
  experimental data reported in the literature show that k is about
  70% larger near the walls than it is at the centre plane.
 
  The RSTM predicts a skin friction coefficient of 3.13E-3,
  which compares well with the experimental value of 3.07E-3.
  The standard k-e model predicts a value of 3.38E-3 and
  here the skin friction coefficient is based on the average
  velocity.
 
REAL(HEIGHT,WTOP,REY,TKEIN,EPSIN,MIXL,DTF,WAV,WSTAR,MASIN)
HEIGHT=0.1;WTOP=1.0; REY=1.E5;WAV=0.5*WTOP
  ** wstar from data of El Telbany & Reynolds [1982]
WSTAR=WAV*0.196/LOG10(REY);TKEIN=WSTAR*WSTAR/.3
MIXL=0.045*HEIGHT;EPSIN=TKEIN**1.5/MIXL*0.1643
    GROUP 4. Y-direction grid specification
ENUL=WAV*HEIGHT/REY
NY=60;YVLAST=HEIGHT;GRDPWR(Y,NY,YVLAST,1.0)
    GROUP 7. Variables stored, solved & named
SOLVE(W1);STORE(ENUT,LEN1)
  ** moving upper wall
PATCH(WAL1,NWALL,1,1,NY,NY,1,NZ,1,1);COVAL(WAL1,W1,1.0,WTOP)
  ** stationary bottom wall
PATCH(WAL2,SWALL,1,1,1,1,1,NZ,1,1);COVAL(WAL2,W1,1.0,0.0)
STORE(V1,KE,DWDY,PVW,PW2,PV2,PU2,DVW)
STORE(PK,EPDK,FWAL,VWDK,U2DK,V2DK,W2DK)
DTF=0.05;TURMOD(REYSTRS,DTF,WAL1,WAL2)
    GROUP 8. Terms (in differential equations) & devices
  ** deactivate convection
TERMS(W1,N,N,P,P,P,P);TERMS(EP,Y,N,P,P,P,P)
TERMS(U2RS,N,N,P,P,P,P);TERMS(V2RS,Y,N,P,P,P,P)
TERMS(W2RS,N,N,P,P,P,P);TERMS(VWRS,Y,N,P,P,P,P)
    GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=0.5*WTOP;FIINIT(EP)=EPSIN;FIINIT(KE)=TKEIN
FIINIT(W2RS)=2.*TKEIN/3.;FIINIT(V2RS)=FIINIT(W2RS)
FIINIT(U2RS)=FIINIT(W2RS);FIINIT(VWRS)=-WSTAR*WSTAR
    GROUP 15. Termination of sweeps
LSWEEP=50;TSTSWP=-1;LITHYD=20
    GROUP 16. Termination of iterations
MASIN=RHO1*WAV*HEIGHT; RESREF(W1)=1.E-12*MASIN*WAV
RESREF(KE)=RESREF(W1)*TKEIN; RESREF(EP)=RESREF(W1)*EPSIN
    GROUP 22. Spot-value print-out
IYMON=NY-2;ITABL=3;NPLT=1;NZPRIN=1;NYPRIN=1;IYPRF=1
    GROUP 24. Dumps for restarts
RELAX(W1,FALSDT,0.3);WALPRN=T;VARMIN(W1)=1.E-10