GROUP 1. Run title and other preliminaries
TEXT(VAN DRIEST_2D-YX ELLIP CHANNL FLOW :T204
TITLE
  DISPLAY
  The problem considered is 2d y-x symmetric turbulent developing
  channel flow between two walls using Prandtl's mixing length
  plus Van-Driest damping with the Nikuradse's EL1=NIKURCH option.
  The main purpose of the calculation is to test the functionality
  of the van Driest option with the elliptic solution algorithm and
  also to test that PHOENICS returns a symmetric solution. About 100
  sweeps are required for convergence, but LSWEEP has been set equal
  to 20 so as to reduce computer time during PHOENICS testing
  procedures.
  ENDDIS
 
  AUTOPLOT USE
  file
  phi 5
 
  clear
  da 1 u1 x 12
  col3 1;plot 1
  @
   0.16021E+03 0.25859E+04 CR
  U1(m/s)@
  @
   0.92282E+03 0.16284E+03 CR
  Distance Between Two Walls of the Channel (m)@
  msg Press e to END
  ENDUSE
BOOLEAN(REYLOW,VARLAM); REYLOW=T;VARLAM=T
REAL(HEIGHT,REY,MIXL,DTF,WIN,MASIN);HEIGHT=0.1;WIN=1.0; REY=1.E5
    GROUP 4. Y-direction grid specification
ENULA=WIN*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
NX=20;GRDPWR(X,NX,5*HEIGHT,1.0)
    GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,U1);STORE(ENUT,LEN1);SOLUTN(P1,Y,Y,Y,N,N,N)
    GROUP 8. Terms (in differential equations) & devices
  ** use DIFCUT=0 to avoid error in EARTH with hybrid option
DIFCUT=0
    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(U1)=WIN
    GROUP 13. Boundary conditions and special sources
  ** upper wall
PATCH(WALLN,NWALL,1,NX,NY,NY,1,1,1,1);COVAL(WALLN,U1,LOGLAW,0.0)
  ** bottom wall
PATCH(WALLS,SWALL,1,NX,1,1,1,1,1,1);COVAL(WALLS,U1,LOGLAW,0.0)
  **  Inlet
INLET(INLET,WEST,1,1,1,NY,1,1,1,1)
VALUE(INLET,P1,WIN);VALUE(INLET,U1,WIN)
  **  Exit
PATCH(OUTLET,EAST,NX,NX,1,NY,1,1,1,1);COVAL(OUTLET,P1,1.E5,0.0)
    GROUP 15. Termination of sweeps
LSWEEP=20;TSTSWP=-1;LITHYD=1
    GROUP 16. Termination of iterations
MASIN=RHO1*WIN*HEIGHT; RESREF(P1)=1.E-12*WIN*HEIGHT
RESREF(U1)=1.E-12*MASIN*WIN; RESREF(V1)=RESREF(W1)
    GROUP 17. Under-relaxation devices
DTF=10.*ZWLAST/WIN
    GROUP 22. Spot-value print-out
IYMON=2;NPLT=1;NZPRIN=2;NYPRIN=2;IYPRF=1
NUMCLS=5;IZMON=5
    GROUP 24. Dumps for restarts
RELAX(U1,FALSDT,DTF); RELAX(V1,FALSDT,DTF);WALPRN=T