PHOTON USE
  p
 
 
 
 
  gr ou x 1
  MSG Velocity vectors
  vec x 1 sh
  msg
  msg Press return to plot pressure contours
  pause
  cont p1 x 1 fil;.01
  msg
  msg Type e to End
  ENDUSE
 
#$B517
TEXT(As for 517 with Navier-Stokes Soltn:B518
TITLE
l(pause
  DISPLAY
  Navier Stokes equations over ellipse.
    This case considers creeping flow over the ellipse: the
  Reynolds number is 1.0. The is done by de-activating the
  Darcy formulation used for the potential-flow simulations in
  the preceding cases, and restoring solution for v1 and w1
  ENDDIS
DARCY=F; COVAL(DARCY,V1,0.0,0.0); COVAL(DARCY,W1,0.0,0.0)
SOLUTN(V1,Y,Y,Y,P,P,P);SOLUTN(W1,Y,Y,Y,P,P,P)
REAL(WIN)
    GROUP 6. Body-fitted coordinates or grid distortion
  ** Reset Algebraic Grid
GSET(M,F1,+J+K,1,1,1,TRANS)
GSET(C,I2,F,I1,+,1,0,0)
NONORT=T; SYMBFC=T
    GROUP 8. Terms (in differential equations) & devices
TERMS(V1,Y,Y,Y,N,N,N); TERMS(W1,Y,Y,Y,N,N,N)
    GROUP 9. Properties of the medium (or media)
ENUL=1.0
    GROUP 11. Initialization of variable or porosity fields
WIN=1.25;FIINIT(W1)=WIN
    GROUP 13. Boundary conditions and special sources
COVAL(INLET,P1,FIXFLU,WIN);COVAL(INLET,W1,ONLYMS,WIN)
COVAL(INLET,V1,ONLYMS,0)
   ** Outlet
COVAL(OUTLET,P1,FIXP,0.0)
COVAL(OUTLET,V1,ONLYMS,0.0);COVAL(OUTLET,W1,ONLYMS,0.0)
   ** Cylinder surface friction
WALL (WALL,SOUTH,1,NX,1,1,IZ1,IZ2-1,1,1)
COVAL(WALL,W1,1.0,0.)
    GROUP 16. Termination of iterations
LSWEEP=200
LITER(V1)=2;LITER(W1)=2
    GROUP 17. Under-relaxation devices
relax(p1,linrlx,0.5)
real(dtf)
dtf=1.e-3
RELAX(V1,FALSDT,dtf); RELAX(W1,FALSDT,dtf)
    GROUP 22. Spot-value print-out
IYMON=3;IZMON=10