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