****** TO LOAD CASE 107:TYPE L(N107) *****
    GROUP 1. Run title and other preliminaries
TEXT(2D TURBULENT BKWRD FACING STEP: N107
TITLE
mesg(PC486/50 time last reported as appx. 12 min
  DISPLAY
  This case is the same of T103 except that higher-order schemes
  are employed, and the calculation is performed in all six planes
  with the standard k-e model as a consistency test. Experiment
  indicates a reattachment length XR/H=7.2 downstream of the step,
  where H is the step height. The present calculations use the
  linear-upwind scheme (LUS) for the momentum equations, and the
  non-linear flux-limiter OSPRE scheme for KE and EP.
 
  The HYBRID scheme yields XR/H=5.12, whereas using LUS on the
  momentum equations and the upwind scheme on KE and EP predicts
  X/H=5.85. If LUS is used for momentum and OSPRE for KE and EP,
  then PHOENICS predicts XR/H=6.1. The higher-order schemes need
  700 sweeps for a converged solution, and HYBRID needs 250 sweeps.
  ENDDIS
  PHOTON USE
 
  p
 
 
   0.20443E+04 0.15633E+04 CR
  gr ou y 1;use patgeo
  stream 2d y 1
  -.99 0 10
  stream 2d y 1
  0. .0374 5
  msg Streamlines
  msg Press  to continue
  pause
  cl;use patgeo;vec y 1 sh
  msg Vectors
  msg Press  to continue
  pause
  ENDUSE
  AUTOPLOT USE
   FILE
   PHI 5
 
   D 1 W1 X 1;PLOT;DIV X .0381 1;SHIFT X -4 1;SCALE
   SCALE X 0 10;LEVEL Y 0
   msg Horizontal W1 velocity distribution along bottom wall
   msg Reattachment point where W1 passes through zero.
   msg Press  to continue
   pause
  ENDUSE
REAL(DTF,HEIGHT,WIDTH,CLEN,SLEN,REYNO,UIN,TKEIN,EPSIN,MIXL,FRIC)
REAL(DTKE);INTEGER(NYS,NXS,NXF,NYF)
CHAR(SCHM,DIRV,DIRH,VELV,VELH,PLANE)
     ** Calculation of domain specifications
HEIGHT=0.0381;WIDTH=3.*HEIGHT;SLEN=4.*HEIGHT;CLEN=20.*HEIGHT
REYNO=4.5E4;UIN=13.;NXS=5;NYS=10;NXF=20;NYF=10
MESG( Enter required vertical coordinate X,Y or Z
MESG( Default:  X
READVDU(DIRV,CHAR,X)
CASE :DIRV: OF
WHEN X,1
+ NREGX=2;VELV=U1
WHEN Y,1
+ NREGY=2;VELV=V1
WHEN Z,1
+ NREGZ=2;VELV=W1
ENDCASE
MESG( Enter required horizontal coordinate X,Y or Z
MESG( Default:  Z
READVDU(DIRH,CHAR,Z)
CASE :DIRH: OF
WHEN X,1
+ NREGX=2;VELH=U1
WHEN Y,1
+ NREGY=2;VELH=V1
WHEN Z,1
+ NREGZ=2;VELH=W1
ENDCASE
PLANE=:DIRV::DIRH:
TEXT(2D:PLANE: TURBULENT BKWRD FACING STEP: N107
    GROUP 3. X-direction grid specification
    GROUP 4. Y-direction grid specification
     ** channel length = 0.762 & channel width = 0.1143
IREG:DIRH:=1;GRDPWR(:DIRH:,NXS,SLEN,1.0)
IREG:DIRH:=2;GRDPWR(:DIRH:,NXF,CLEN-SLEN,1.1)
IREG:DIRV:=1;GRDPWR(:DIRV:,-NYS,HEIGHT,1.3)
IREG:DIRV:=2;GRDPWR(:DIRV:,-NYF,WIDTH-HEIGHT,1.4)
    GROUP 7. Variables stored, solved & named
SOLVE(P1,:VELH:,:VELV:);SOLUTN(P1,Y,Y,Y,N,N,N)
SOLUTN(:VELH:,P,P,P,P,P,N);SOLUTN(:VELV:,P,P,P,P,P,N)
TURMOD(KEMODL);STORE(ENUT)
    GROUP 8. Terms (in differential equations) & devices
MESG( Enter required convection scheme
MESG( Default: LUS - Linear upwind for momentum Scheme  OSPRE k-e
MESG( The alternative is:
MESG(  HYB  - Hybrid Differencing Scheme for all variables
READVDU(SCHM,CHAR,LUS)
CASE :SCHM: OF
WHEN HYB,3
+ MESG(Hybrid-differencing scheme
+ DIFCUT=0.5;DTF=CLEN/UIN;DTKE=DTF/5
+ LSWEEP=400
WHEN LUS,3
+ MESG(LUS for momentum  OSPRE for ke and ep
+ SCHEME(LUS,:VELH:,:VELV:);SCHEME(OSPRE,KE,EP)
+ DTF=0.2*CLEN/(UIN*N:DIRH:);DTKE=DTF
+ LSWEEP=700
ENDCASE
    GROUP 9. Properties of the medium (or media)
RHO1=1.0;ENUL=UIN*HEIGHT/REYNO
    GROUP 11. Initialization of variable or porosity fields
FRIC=0.018;TKEIN=0.25*UIN*UIN*FRIC;MIXL=0.09*HEIGHT
EPSIN=0.1643*TKEIN**1.5/MIXL;FIINIT(:VELH:)=UIN;FIINIT(P1)=1.3E-4
FIINIT(KE)=TKEIN;FIINIT(EP)=EPSIN
FIINIT(:VELV:)=0.1*UIN
     ** Initialization of variables in blocked region
CONPOR(STEP,0.0,CELL,#1,#1,#1,#1,#1,#1)
    GROUP 13. Boundary conditions and special sources
IF(:PLANE:.EQ.YX) THEN
+ INLET(INLET,WEST,#1,#1,#2,#NREGY,#1,#1,1,1)
+ PATCH(OUTLET,EAST,#NREGX,#NREGX,#1,#NREGY,#1,#1,1,1)
+ WALL (TOP,NORTH,#1,#NREGX,#NREGY,#NREGY,#1,#1,1,1)
+ WALL (BOT,SOUTH,#2,#NREGX,#1,#1,#1,#1,1,1)
+ PATCH(STEPL,WWALL,$NREGX,$NREGX,#1,#1,1,1,1,1)
+ PATCH(STEPT,SWALL,#1,#1,$NREGY,$NREGY,1,1,1,1)
ENDIF
IF(:PLANE:.EQ.XY) THEN
+ INLET(INLET,SOUTH,#2,#NREGX,#1,#1,#1,#1,1,1)
+ PATCH(OUTLET,NORTH,#1,#NREGX,#NREGY,#NREGY,#1,#1,1,1)
+ WALL (TOP,EAST,#NREGX,#NREGX,#1,#NREGY,#1,#1,1,1)
+ WALL (BOT,WEST,#1,#1,#2,#NREGY,#1,#1,1,1)
+ PATCH(STEPL,SWALL,#1,#1,$NREGY,$NREGY,1,1,1,1)
+ PATCH(STEPT,WWALL,$NREGX,$NREGX,#1,#1,1,1,1,1)
ENDIF
IF(:PLANE:.EQ.YZ) THEN
+ INLET(INLET,LOW,#1,#1,#2,#NREGY,#1,#1,1,1)
+ PATCH(OUTLET,HIGH,1,1,#1,#NREGY,#NREGZ,#NREGZ,1,1)
+ WALL (TOP,NORTH,1,1,#NREGY,#NREGY,#1,#NREGZ,1,1)
+ WALL (BOT,SOUTH,1,1,#1,#1,#NREGZ,#NREGZ,1,1)
+ PATCH(STEPL,LWALL,1,1,#1,#1,$NREGZ,$NREGZ,1,1)
+ PATCH(STEPT,SWALL,1,1,$NREGY,$NREGY,#1,#1,1,1)
ENDIF
IF(:PLANE:.EQ.ZY) THEN
+ INLET(INLET,SOUTH,#1,#1,#1,#1,#2,#NREGZ,1,1)
+ PATCH(OUTLET,NORTH,1,1,#NREGY,#NREGY,#1,#NREGZ,1,1)
+ WALL (TOP,HIGH,1,1,#1,#NREGY,#NREGZ,#NREGZ,1,1)
+ WALL (BOT,LOW,1,1,#2,#NREGY,#1,#1,1,1)
+ PATCH(STEPL,SWALL,1,1,$NREGY,$NREGY,#1,#1,1,1)
+ PATCH(STEPT,LWALL,1,1,#1,#1,$NREGZ,$NREGZ,1,1)
ENDIF
IF(:PLANE:.EQ.XZ) THEN
+ INLET(INLET,LOW,#2,#NREGX,#1,#1,#1,#1,1,1)
+ PATCH(OUTLET,HIGH,#1,#NREGX,1,1,#NREGZ,#NREGZ,1,1)
+ WALL (TOP,EAST,#NREGX,#NREGX,1,1,#1,#NREGZ,1,1)
+ WALL (BOT,WEST,#1,#1,1,1,#2,#NREGZ,1,1)
+ PATCH(STEPL,LWALL,#1,#1,1,1,$NREGZ,$NREGZ,1,1)
+ PATCH(STEPT,WWALL,$NREGX,$NREGX,1,1,#1,#1,1,1)
ENDIF
IF(:PLANE:.EQ.ZX) THEN
+ INLET(INLET,WEST,#1,#1,#1,#1,#2,#NREGZ,1,1)
+ PATCH(OUTLET,EAST,#NREGX,#NREGX,1,1,#1,#NREGZ,1,1)
+ WALL (TOP,HIGH,#1,#NREGX,1,1,#NREGZ,#NREGZ,1,1)
+ WALL (BOT,LOW,#2,#NREGX,1,1,#1,#1,1,1)
+ PATCH(STEPL,WWALL,$NREGX,$NREGX,1,1,#1,#1,1,1)
+ PATCH(STEPT,LWALL,#1,#1,1,1,$NREGZ,$NREGZ,1,1)
ENDIF
VALUE(INLET,P1,UIN);VALUE(INLET,:VELH:,UIN)
VALUE(INLET,KE,TKEIN);VALUE(INLET,EP,EPSIN)
 
COVAL(OUTLET,P1,1.E3,0.0)
COVAL(OUTLET,:VELH:,ONLYMS,0.0);COVAL(OUTLET,:VELV:,ONLYMS,0.0)
COVAL(OUTLET,KE,ONLYMS,0.0);COVAL(OUTLET,EP,ONLYMS,0.0)
 
COVAL(STEPT,:VELH:,LOGLAW,0.);COVAL(STEPT,KE,LOGLAW,LOGLAW)
COVAL(STEPT,EP,LOGLAW,LOGLAW)
 
COVAL(STEPL,:VELV:,LOGLAW,0.);COVAL(STEPL,KE,LOGLAW,LOGLAW)
COVAL(STEPL,EP,LOGLAW,LOGLAW)
    GROUP 15. Termination of sweeps
    GROUP 16. Termination of iterations
    GROUP 17. Under-relaxation devices
RELAX(:VELH:,FALSDT,DTF);RELAX(:VELV:,FALSDT,DTF)
KELIN=1;RELAX(KE,FALSDT,DTKE);RELAX(EP,FALSDT,DTKE)
I:DIRV:MON=NYS-2;I:DIRH:MON=NXS+2;NPRMON=100
    GROUP 23. Field print-out and plot control
ITABL=3;NPLT=10;IPLTL=LSWEEP;WALPRN=T;TSTSWP=-1
NSAVE=BS:PLANE: