GROUP 1. Run title and other preliminaries
TEXT(2S K-E MODEL_1D PLANE COUETTE FLOW :T400
TITLE
mesg(PC486/50 time last reported as 26.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.
 
  Calculations are made with the high-Re forms of the standard
  k-e model and 2-scale k-e model. The latter is selected by
  setting TSKE=T. For this case, the standard and 2-scale forms of
  the k-e model produce very similar results.
  ENDDIS
 
   AUTOPLOT USE
   file
   phi 5
 
   clear
   da 1 w1
   col9 1
   msg Velocity (W1) profile
   msg Press RETURN to continue
   pause
   clear
   da 1 ke
   col9 1
   scale y 0. 2.e-3
   redraw
   msg KE profile
   msg Press e to END
   ENDUSE
 
BOOLEAN(TSKE);TSKE=T
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=30;YVLAST=HEIGHT;GRDPWR(Y,NY,YVLAST,1.0)
    GROUP 7. Variables stored, solved & named
SOLVE(W1);STORE(ENUT,LEN1)
IF(TSKE) THEN
+ TURMOD(TSKEMO)
ELSE
+ TURMOD(KEMODL)
ENDIF
    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)
IF(TSKE) THEN
+ TERMS(KP,N,N,P,P,P,P);TERMS(KT,Y,N,P,P,P,P)
+ TERMS(ET,Y,N,P,P,P,P)
ELSE
+ TERMS(KE,N,N,P,P,P,P)
ENDIF
    GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=0.5*WTOP;FIINIT(EP)=EPSIN
IF(TSKE) THEN
+ REAL(KTDKP);KTDKP=0.25
+ FIINIT(ET)=EPSIN;FIINIT(KP)=TKEIN/(1.+KTDKP)
+ FIINIT(KT)=KTDKP*FIINIT(KP)
ENDIF
FIINIT(KE)=TKEIN
    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=20;TSTSWP=-1;LITHYD=10
    GROUP 16. Termination of iterations
MASIN=RHO1*WAV*HEIGHT
RESREF(W1)=1.E-12*MASIN*WAV
RESREF(EP)=RESREF(W1)*EPSIN
IF(TSKE) THEN
+ RESREF(KP)=RESREF(W1)*TKEIN; RESREF(ET)=RESREF(EP)
+ RESREF(KT)=RESREF(KP)
ELSE
+ RESREF(KE)=RESREF(W1)*TKEIN
ENDIF
    GROUP 17. Under-relaxation devices
DTF=20.*ZWLAST/WAV
    GROUP 18. Limits on variables or increments to them
VARMIN(W1)=1.E-10
    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,DTF); RELAX(EP,FALSDT,DTF)
IF(TSKE) THEN
+ RELAX(KP,FALSDT,DTF); RELAX(KT,FALSDT,DTF)
+ RELAX(ET,FALSDT,DTF)
ELSE
+ RELAX(KE,FALSDT,DTF)
ENDIF
WALPRN=T