L(191)
 
  PHOTON USE
  p
  parphi
  1.0 1.0 0.2
 
 
 
 
  msg         FLAT-PLATE TURBULENT BOUNDARY LAYER ('k-e' model)
  msg
  set ref vec 75
  msg   Pressure contours & axial velocity vectors:
  gr ou x 1
  con p1 x 1 fi;.01;vec x 1
  msg
  msg Press  to continue
  pause
  con off;vec off;red
  msg        Axial velocity contours:
  con wcrt x 1 fi;.01
  msg
  msg Press  to continue
  pause
  con off;red
  msg        Temperature distribution:
  con temp x 1 fi;.01
  msg            -
  msg Press  to continue
  pause
  con off;red
  msg        Distribution of turbulent kinetic energy (k):
  con ke x 1 fi;.01
  msg
  msg Press  to continue
  pause
  con off;red
  msg        Eddy-viscosity distribution:
  con enut x 1 fi;.01
  msg
  msg Press e to END
  enduse
 
    GROUP 1. Run title and other preliminaries
TEXT(Boundary Layer K-E Turbulen Model 
TITLE
mesg(PC486/50 time last reported as 4.min
  The flow is as for cases 190 and 191; but now the two-equation
  k-e turbulence model..
    GROUP 7. Variables stored, solved & named
TURMOD(KEMODL)
 
    GROUP 11. Initialization of variable or porosity fields
FIINIT(EP)=FIINIT(KE)**1.5/YVLAST
 
    GROUP 13. Boundary conditions and special sources
  ** South Wall Boundary
COVAL(WFUN,KE,LOGLAW,LOGLAW);COVAL(WFUN,EP,LOGLAW,LOGLAW)
  ** North Free Boundary
COVAL(FREE,KE,ONLYMS,1.E-10);COVAL(FREE,EP,ONLYMS,1.E-10)
  ** Inlet Boundary
COVAL(PROF,KE,ONLYMS,GRND3);COVAL(PROF,EP,ONLYMS,GRND3)
 
    GROUP 21. Print-out of variables
OUTPUT(KE,Y,Y,Y,Y,Y,Y);OUTPUT(EP,Y,Y,Y,Y,Y,Y)
 
    GROUP 23. Field print-out and plot control
PLOT(IZEQNZ,KE,0.,0.0)
 
  ***actdem***
 
do ii=1,5
+ mesg(
enddo
mesg( Initial data that can be changed:
+ mesga(  Free stream velocity is set to 33. m/s
+ mesg(  Free stream temperature is set to 5. deg.
+ mesg(  Wall temperature is set to 10. deg.
mesga( Do you want to change settings (y/n)? (Default n)
readvdu(ans,char,n)
if(:ans:.eq.y) then
+ real(rt1)
+ do ii=1,5
+   mesg(
+ enddo
+ mesga(  Free stream velocity is 33. m/s. OK?
+ mesg(If not, insert new value.
+ readvdu(rt1,real,33.)
+ wfree=rt1
 
+ do ii=1,5
+   mesg(
+ enddo
+ mesga(  Free stream temperature is 5. deg. OK?
+ mesg(If not, insert new value.
+ readvdu(rt1,real,5.)
+ tfree=rt1
+ do ii=1,5
+   mesg(
+ enddo
+ mesga(  Wall temperature is 10. deg. OK? If not, insert new value.
+ readvdu(rt1,real,10.)
+ twall=rt1
 
+ FIINIT(W1)=WFREE;FIINIT(TEMP)=TFREE
+ COVAL(WFUN,TEMP,LOGLAW,TWALL)
 
+ COVAL(FREE,W1,ONLYMS,WFREE);COVAL(FREE,TEMP,ONLYMS,TFREE)
+ TEMP0=TWALL; EL1B=WFREE
+ FIINIT(KE)=WFREE**2*0.01
+ FIINIT(EP)=FIINIT(KE)**1.5/YVLAST
endif
 
IDISPA=1; TSTSWP=-1
selref=t; resfac=1.e-2