GROUP 1. Run title and other preliminaries
TEXT(VAN DRIEST_1D PLANE COUETTE FLOW   :T200
TITLE
mesg(PC486/50 time last reported as 20.sec
  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.
 
  The turbulence is simulated by use of the Prandtl mixing length
  model with Nikuradse's prescription of the mixing length. If
  REYLOW=T in the Q1 file, the calculation integrates down to both
  walls by using van Driest's viscous-damping extension to the
  mixing length. The solution is performed one-dimensionally by use
  of the single-slab solver. A non-uniform grid is employed so as
  to concentrate cells very close to the walls.
  ENDDIS
 
  For both the low- and high-Reynolds-number turbulence closure
  models, PHOENICS predicts skin-friction coefficient of 2.77E-3
  which compares well with the experimental value of 3.07E-3
  suggested by the data of El Telbany and Reynolds [1982].
 
  AUTOPLOT USE
  file
  phi 5
 
  clear
  da 1 w1
  colour3
  plot 1
  @
   0.10254E+03 0.26250E+04 CR
  W1(m/s)@
  @
   0.18969E+04 0.16284E+03 CR
  Distance across the Channel(m)@
  @
   0.33644E+04 0.29051E+04 CR
    Moving Wall
 
            ^
            | W0
            |
            |
  @
  msg Press e to END
  ENDUSE
BOOLEAN(REYLOW,VARLAM); REYLOW=T;VARLAM=F
REAL(HEIGHT,WTOP,REY,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)
WSTAR
    GROUP 4. Y-direction grid specification
ENULA=WAV*HEIGHT/REY;YVLAST=HEIGHT
IF(REYLOW) THEN
+ NY=100;GRDPWR(Y,-NY,YVLAST,1.7)
ELSE
+ NY=30;GRDPWR(Y,NY,YVLAST,1.0)
ENDIF
    GROUP 7. Variables stored, solved & named
SOLVE(W1);STORE(ENUT,LEN1)
    GROUP 8. Terms (in differential equations) & devices
TERMS(W1,N,N,P,P,P,P)
    GROUP 9. Properties of the medium (or media)
ENUL=ENULA;ENUT=MIXLEN;EL1=NIKURCH;GENK=T
IF(REYLOW) THEN
+ IENUTA=5
ENDIF
   ** test for ground-set enul
IF(VARLAM) THEN
+ TMP1=CONST;TMP1A=0.0;ENUL=LINTEM
ENDIF
    GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=0.5*WTOP
    GROUP 13. Boundary conditions and special sources
  ** moving upper wall
WALL(WALLN,NORTH,1,1,NY,NY,1,NZ,1,1);COVAL(WALLN,W1,LOGLAW,WTOP)
  ** stationary bottom wall
WALL(WALLS,SOUTH,1,1,1,1,1,NZ,1,1)
    GROUP 15. Termination of sweeps
LSWEEP=30;TSTSWP=-1;LITHYD=6
    GROUP 16. Termination of iterations
MASIN=RHO1*WAV*HEIGHT; RESREF(W1)=1.E-12*MASIN*WAV
    GROUP 17. Under-relaxation devices
DTF=10.*ZWLAST/WAV;VARMIN(W1)=1.E-10
    GROUP 22. Spot-value print-out
IYMON=2;NPLT=2;NZPRIN=1;NYPRIN=4;IYPRF=1;TSTSWP=-1
    GROUP 24. Dumps for restarts
RELAX(W1,FALSDT,DTF);WALPRN=T