PHOTON USE
  P



  gr ce ou x 1;set lhand
  vec x 1 sh
  msg                Velocity distribution
  msg  Press Enter to continue
  pause
  vec off;red
  con H1 x 1 fil;0.001
  msg                Temperature distribution
  msg  Press Enter to continue
  pause
  con off;red
  con KE x 1 fil;0.001
  msg                Turbulence energy distribution
  msg  Press Enter to continue
  pause
  con off;red
  con EL1 x 1 fil;0.001
  msg                Length scale distribution
  msg(Press Enter to continue
  pause
  con off;red
  con ENUT x 1 fil;0.001
  msg                Turbulence viscosity distribution
  ENDUSE
    GROUP 1. Run title and other preliminaries
TEXT(HEAT TRANSFER IN ABRUPT ENLARGEMENT AT PIPE ENTRANCE:140

  DISPLAY
  This case  represents  the  flow  in  a circular duct that has an
  abrupt enlargement at the  inlet.The  duct  wall  is  held  at  a
  constant temperature.

  Auxiliary relations  comprise  a turbulence model similar to that
  of Prandtl.  The length scale distributions are determined as  of
  the type proposed by Ranchal and Spalding in 1972 ( J. of Waerme-
  und Stoff Uebertragung, 5, 1972 )

  PLANT is   used   to  introduce  the  codings  for  length  scale
  distributions guided by the formulae provided in  Group  9 in the
  string marked PRPT01.
  ENDDIS

    GROUP 3. X-direction grid specification
  ** Activate polar coordinate system
CARTES=F
    GROUP 4. Y-direction grid specification
GRDPWR(Y,20,0.2,1.0)
    GROUP 5. Z-direction grid specification
GRDPWR(Z,20,2.0,1.0)
    GROUP 6. Body-fitted coordinates or grid distortion
    GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,W1,H1)
SOLUTN(P1,Y,Y,N,N,N,N)
TURMOD(KLMODL);STORE(ENUT,EL1)
    GROUP 8. Terms (in differential equations) & devices
TERMS(H1,N,Y,Y,Y,Y,Y)
    GROUP 9. Properties of the medium (or media)
REAL(RE);RE=5.e4
ENUL=2.*0.2/RE
PRNDTL(H1)=3.0;PRT(H1)=0.5;PRT(KE)=1.5
EL1=GRND

      PLANTBEGIN
   ** Length scale distribution
  LEN1=AMIN1(.0166*EXP(60.*(.2-RG2D)-1.),$
               .016*(ZGNZ/0.2)**1.09,0.12)
      PLANTEND

    GROUP 10. Inter-phase-transfer processes and properties
    GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=0.1
     ** Calculation of KE (where fric=0.03)...
REAL(TKEIN);TKEIN=2.*2.*0.03
    GROUP 12. Unused
    GROUP 13. Boundary conditions and special sources
PATCH(INLET,LOW,1,1,1,NY/2,1,1,1,1)
COVAL(INLET,P1,FIXFLU,2.);COVAL(INLET,W1,ONLYMS,2.)
COVAL(INLET,H1,ONLYMS,1.0);COVAL(INLET,KE,ONLYMS,TKEIN)

PATCH(OUTLET,HIGH,1,1,1,NY,NZ,NZ,1,1)
COVAL(OUTLET,P1,FIXVAL,0.0)

PATCH(LOWWAL,LWALL,1,1,NY/2+1,NY,1,1,1,1)
COVAL(LOWWAL,V1,LOGLAW,0.0)
COVAL(LOWWAL,KE,LOGLAW,LOGLAW)

PATCH(NORWAL,NWALL,1,1,NY,NY,1,NZ,1,1)
COVAL(NORWAL,W1,LOGLAW,0.0)
COVAL(NORWAL,H1,LOGLAW,0.0)
COVAL(NORWAL,KE,LOGLAW,LOGLAW)

    GROUP 15. Termination of sweeps
RESREF(W1)=1.E-6;RESREF(P1)=1.E-6;RESREF(V1)=1.E-6
RESREF(KE)=1.E-6
LSWEEP=50
RELAX(W1,FALSDT,1.0);RELAX(V1,FALSDT,1.0)
    GROUP 20. Preliminary print-out
    GROUP 22. Spot-value print-out
NPRMON=20;IYMON=5
    GROUP 22. Spot-value print-out
IYMON=19;IZMON=10;IPLTL=LSWEEP
    GROUP 23. Field print-out and plot control
NYPRIN=2;NZPRIN=5;NPLT=1
NAMSAT=MOSG

tstswp=-1
dmpstk=t
DISTIL=T
EX(P1)=1.492E-01; EX(V1)=2.481E-02; EX(W1)=6.473E-01
EX(H1)=9.695E-01; EX(EL1)=6.807E-02; EX(ENUT)=1.062E-02
 LIBREF=140
STOP