GROUP 1. Run title and other preliminaries
TEXT(Vortex Near A Moving Wall         
TITLE
  DISPLAY
  This unsteady calculation shows what happens in a cavity with
  a moving wall when an initial vortex is prescribed. The wall
  temperature differs from the initial temperature of the fluid.
  Interesting variations include: changes to the vortex intensity;
  changes to the temperature and velocity of the wall;and changes
  to the size of the vortex.
  ENDDIS
 
    GROUP 2. Transience; time-step specification
STEADY=F;GRDPWR(T,10,20.0,1.0)
 
    GROUP 3. X-direction grid specification
GRDPWR(X,10,1.0,1.0)
 
    GROUP 4. Y-direction grid specification
GRDPWR(Y,10,1.0,1.0)
 
    GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,H1)
 
    GROUP 8. Terms (in differential equations) & devices
TERMS(H1,N,Y,N,Y,Y,Y)
 
    GROUP 9. Properties of the medium (or media)
ENUL=1.E-2
 
    GROUP 11. Initialization of variable or porosity fields
FIINIT(U1)=0.0;FIINIT(V1)=0.0;FIINIT(P1)=0.0
FIINIT(H1)=0.0
  ** Vortex
PATCH(VORTEX,LINVLY,4,5,4,6,1,1,1,1)
INIT(VORTEX,U1,10.0,-1.0);INIT(VORTEX,H1,0.0,1.0)
 
    GROUP 13. Boundary conditions and special sources
XCYCLE=T
  ** Pressure relief
PATCH(RELIEF,CELL,1,1,1,1,1,1,1,LSTEP);COVAL(RELIEF,P1,FIXP,0.0)
COVAL(RELIEF,U1,ONLYMS,0.0);COVAL(RELIEF,V1,ONLYMS,0.0)
  ** Moving Wall
WALL (MOVING,SOUTH,1,NX,1,1,1,1,1,LSTEP);COVAL(MOVING,U1,1.0,.1)
 
    GROUP 15. Termination of sweeps
LSWEEP=20
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
    GROUP 16. Termination of iterations
 
    GROUP 21. Print-out of variables
OUTPUT(U1,Y,Y,Y,Y,Y,Y);OUTPUT(V1,Y,Y,Y,Y,Y,Y)
 
    GROUP 22. Spot-value print-out
IXMON=5;IYMON=5
 
    GROUP 23. Field print-out and plot control
NTPRIN=5;NXPRIN=2;NYPRIN=2
PATCH(Y15Z10,CONTUR,1,NX,1,NY,1,1,1,LSTEP)
PLOT(Y15Z10,U1,0.0,10.0);PLOT(Y15Z10,H1,0.0,10.0)
PATCH(TIMEPLOT,PROFIL,NX/2,NX/2,1,NY,1,1,1,LSTEP)
PLOT(TIMEPLOT,U1,-0.1,0.1)