PHOTON USE
  p
  phi
 
 
 
  msg grid
  gr x 1
  msg Press  to continue
  pause
  cl
  msg 1st-phase (light phase) velocity vectors
  vec x 1 sh
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg 2nd-phase (dense phase) velocity vectors
  set vec comp
  - v2 w2
  vec x 1 sh
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of pressure
  con p1 x 1 fi;0.002
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of 1st-phase (light phase) volume fraction
  con r1 x 1 fi;0.002
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of 2nd-phase (dense phase) volume fraction
  con r2 x 1 fi;0.002
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of 1st-phase swirl velocity
  con u1/r x 1 fi;0.002
  gr ou x 1
  msg Press e to END
  enduse
    GROUP 1. Run title and other preliminaries
TEXT(EFFECT OF SWIRL :                   W575
TITLE
    GROUP 3. X-direction grid specification
    Activate polar coordinate system
CARTES=F;XULAST=.1
    GROUP 4. Y-direction grid specification
GRDPWR(Y,20,0.05,1.0);RINNER=0.1
    GROUP 5. Z-direction grid specification
GRDPWR(Z,20,0.40,2.0)
    GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,W1);SOLUTN(P1,Y,Y,Y,N,N,N)
ONEPHS=F;SOLVE(U2,V2,W2,R1,R2);STORE(CFIP)
    GROUP 9. Properties of the medium (or media)
ENUL=1.E-4;RHO2=20
    GROUP 10. Inter-phase-transfer processes and properties
CFIPS=1E2
    GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=0.1;FIINIT(W2)=0.1;FIINIT(R1)=0.5;FIINIT(R2)=0.5
    GROUP 13. Boundary conditions and special sources
IURVAL=-1
WALL(WALL,NORTH,1,1,NY,NY,1,NZ,1,1);VALUE(WALL,U1,6.0)
WALL(INSIDE,SOUTH,1,1,1,1,1,NZ,1,1);VALUE(INSIDE,U1,6.0)
INLET(IN,LOW,1,1,1,NY,1,1,1,1)
VALUE(IN,P1,0.5*RHO1*0.1);VALUE(IN,W1,0.1)
VALUE(IN,P2,0.5*RHO2*0.1);VALUE(IN,W2,0.1)
PATCH(OUTLET,HIGH,1,1,1,NY,NZ,NZ,1,1)
COVAL(OUTLET,P1,1E3,0.0)
COVAL(OUTLET,V1,ONLYMS,0.0);COVAL(OUTLET,W1,ONLYMS,0.0)
COVAL(OUTLET,U1,ONLYMS,0.0)
COVAL(OUTLET,P2,RHO2*1E3/(2*RHO1),0.0)
COVAL(OUTLET,V2,ONLYMS,0.0);COVAL(OUTLET,W2,ONLYMS,0.0)
COVAL(OUTLET,U2,ONLYMS,0.0)
    GROUP 15. Termination of sweeps
LSWEEP=75
    GROUP 16. Termination of iterations
RESREF(P1)=1E-10
RESREF(U1)=1E-10;RESREF(V1)=1E-10;RESREF(W1)=1E-10
RESREF(U2)=1E-10;RESREF(V2)=1E-10;RESREF(W2)=1E-10
    GROUP 17. Under-relaxation devices
REAL(DTF);DTF=ZWLAST/0.1/NZ
RELAX(W1,FALSDT,DTF);RELAX(V1,FALSDT,DTF)
RELAX(P1,LINRLX,0.5);RELAX(U1,FALSDT,DTF)
RELAX(W2,FALSDT,DTF);RELAX(V2,FALSDT,DTF)
RELAX(U2,FALSDT,DTF)
    GROUP 22. Spot-value print-out
IYMON=15;IZMON=18
    GROUP 23. Field print-out and plot control
IURPRN=-1;NPLT=1;tstswp=-1