GROUP 1. Run title and other preliminaries
TEXT(Viscous Heating In Couette Flow   
TITLE
  DISPLAY
         |             h   h
         |         h    u  | moving wall
  fixed  |      h    u     |
  wall   |   h    u        |
         | h   u           |
         |  u              |
         u -------> y      |ywall
 
  The flow is laminar, with uniform properties. The non-linearity
  of the temperature profile is due to viscous dissipation.
 
  The solution should be:  u/uwall = y/ywall;
     (h-hzero)/(hwall-hzero) = y/ywall + H*(y/ywall)*(1 - y/ywall)
  where:  H = 0.5 * Prandtl No * uwall**2 /(hwall -hzero)
 
  In this case, H = 1, which implies zero heat flux at the moving
  wall. The numerical calculations conform completely with the
  analytical solution.
 
  ENDDIS
REAL(HCONST)
HCONST=1.0
    GROUP 4. Y-direction grid specification
NY=40;GRDPWR(Y,NY,YVLAST,1.0)
    GROUP 7. Variables stored, solved & named
SOLVE(U1,H1)
 
    GROUP 8. Terms (in differential equations) & devices
TERMS(H1,Y,N,Y,N,Y,N);TERMS(U1,Y,N,Y,N,Y,N)
 
    GROUP 9. Properties of the medium (or media)
ENUL=1.0;PRNDTL(H1)=2.0*HCONST
 
    GROUP 13. Boundary conditions and special sources
WALL (INNER,SOUTH,1,1,1,1,1,1,1,1)
COVAL(INNER,U1,1.,0.0);COVAL(INNER,H1,1.0,0.0)
WALL (OUTER,NORTH,1,1,NY,NY,1,1,1,1)
COVAL(OUTER,U1,1.0,1.0);COVAL(OUTER,H1,1.0,1.0)
 
    GROUP 15. Termination of sweeps
LSWEEP=3
 
    GROUP 16. Termination of iterations
RESREF(U1)=1.E-10;RESREF(H1)=1.E-10
    GROUP 22. Spot-value print-out
;IYMON=20
    GROUP 23. Field print-out and plot control
PATCH(PROFUH,PROFIL,1,1,1,NY,1,1,1,1)
PLOT(PROFUH,U1,0.0,0.0);PLOT(PROFUH,H1,0.0,0.0)
ORSIZ=0.5;ITABL=2
 **END OF LIBRARY CASE 228