PHOTON USE
  p
 
 
 
  msg             Numerical solution  velocity contours
  con U Z 1 fi;0.001
  msg Press Enter to continue
  PAUSE
  CLEAR
  con UEXC Z 1 fi;0.001
  msg               Exact solution  velocity contours
  msg Press e to END
  ENDUSE
 
    GROUP 1. Run title and other preliminaries

  DISPLAY 
    The non-isothermal Coette flow with exponentially temperature
    dependent viscosity has an exact solution published in Numerical
    Heat Transfer, vol. 13, p. 272, 1988. PLANT is used to perform
    the codings for out-put preparations and to introduce the linear
    temperature distribution.  The alternative is  to initialise the
    latter inserting PLANT instructions in Group 11.
  ENDDIS
 
TEXT( NON-ISOTHERMAL COUETTE FLOW:113
NX=10;NY=2
    GROUP 3. X-direction grid specification
GRDPWR(X,NX,1.,1.0)
    GROUP 4. Y-direction grid specification
GRDPWR(Y,NY,2.,1.0)
    GROUP 7. Variables stored, solved & named
SOLVE(U,TEMP);STORE(UEXC,TEMP)
    GROUP 8. Terms (in differential equations) & devices
TERMS(U,N,N,Y,P,P,P)
    GROUP 9. Properties of the medium (or media)
RHO1=1.0
ENUL=GRND

  PLANTBEGIN    
   VISL=EXP(-1.*TEMP)
   UEXC=(EXP(1.)+1.)/(EXP(1.)-1.)*(EXP(XG2D)-1.)-$
    XG2D*EXP(XG2D)
   TEMP=XG2D
  PLANTEND

    GROUP 13. Boundary conditions and special sources
PATCH(DIFWES,WWALL,1,1,1,NY,1,NZ,1,LSTEP)
COVAL(DIFWES,U,1./PRNDTL(U),0.0)
 
PATCH(DIFEAS,EWALL,NX,NX,1,NY,1,NZ,1,LSTEP)
COVAL(DIFEAS,U,1./PRNDTL(U),1.0)
 
PATCH(DPDX,NORTH,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(DPDX,U,FIXFLU,1.)
    GROUP 15. Termination of sweeps
LSWEEP=20
    GROUP 19.
NAMSAT=MOSG
    GROUP 20. Preliminary print-out
    GROUP 23. Field print-out and plot control
tstswp=-1
dmpstk=t
DISTIL=T
 EX(UEXC)=5.546E-01; EX(TEMP)=5.000E-01; EX(U   )=5.580E-01
 LIBREF=113
STOP