#$r002
 
   PHOTON USE
   AUTOPLOT
   file
   phi 5
 
   cl
   msg Blue curve ---     Temperature at Top Plate
   msg Light Blue curve---Temperature at Lower Plate
   da 1 tem1 y 15;da 1 tem1 y 1
   col3 1;col5 2
   ENDUSE
   ************************************************************
   *  GROUP 1.  Run identifiers and other preliminaries.
LOAD(R102)
TEXT(TURB. FLOW BETWEEN PARALLEL PLATE: R103
TITLE
  DISPLAY
   The geometry and boundary conditions of this case are identical
   to those of case R102, but the flow is turbulent with a Reynolds
   number of 25,000.
  ENDDIS
 
REAL(KEIN,EPIN)
REYNO=25000.;GPRL=0.7;GCHNL=2.0;GCHNH=GCHNL/200.
GCOND=0.0263;gccp=1007.
GTIN=0.5*GCOND/(SIGMA*GCHNH);GTIN=GTIN**0.3333333
GQDOT=SIGMA*GTIN**4.
GENUL=GPRL*GCOND/GCCP;GUIN=REYNO*GENUL/(2.*GCHNH)
KEIN=0.02*GUIN**2;EPIN=KEIN**1.5*0.1643/(0.09*GCHNH)
STORE(ENUT)
TURMOD(KEMODL)
VALUE(IN,KE,KEIN);VALUE(IN,EP,EPIN)
COVAL(TOP,U1,LOGLAW,0.0)
COVAL(TOP,KE,LOGLAW,LOGLAW)
COVAL(TOP,EP,LOGLAW,LOGLAW)
COVAL(BOT,U1,LOGLAW,0.0)
COVAL(BOT,KE,LOGLAW,LOGLAW)
COVAL(BOT,EP,LOGLAW,LOGLAW)
LIBREF=103