PHOTON USE
  p
  phi
  10 1
 
 
  msg grid (scale enlarged 10 times in y direction)
  gr x 1
  msg Press  to continue
  pause
  cl
  msg 1st-phase (air) velocity vectors
  vec x 1 sh
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg 2nd-phase (liquid) 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 (air) volume fraction
  con r1 x 1 fi;0.002
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of 1st-phase (air) enthalpy
  con h1 x 1 fi;0.002
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of 2nd-phase (liquid) enthalpy
  con h2 x 1 fi;0.002
  gr ou x 1
  msg Press e to END
  enduse
    GROUP 1. Run title and other preliminaries
TEXT(AREA RESTRICTION AT QUARTER LENGTH: W371
TITLE
  DISPLAY
  This run examines a two-phase, fluid-solid flow in a pipe with an
  orifice at one quarter of the pipe length from the entrance. The
  pipe wall is held at constant temperature.
  Interesting variants include: changes to the mass-inflow rate;
  changes to the density of both phases; changes to the interphase
  friction coefficient; and addition of gravitational effects.
  ENDDIS
    GROUP 3. X-direction grid specification
CARTES=F;XULAST=0.01
    GROUP 4. Y-direction grid specification
GRDPWR(Y,20,0.01,1.0)
    GROUP 5. Z-direction grid specification
GRDPWR(Z,20,0.20,1.0)
    GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,W1,H1);ONEPHS=F;SOLVE(R1,R2,V2,W2,H2)
TERMS(H1,N,Y,Y,Y,Y,Y)
TERMS(H2,N,Y,Y,Y,N,Y)
    GROUP 9. Properties of the medium (or media)
ENUL=1.E-5;PRNDTL(H1)=0.7;RHO2=1000.0
    GROUP 10. Inter-phase-transfer processes and properties
  ** Set constant inter-phase friction coefficient.
CFIPS=1.E2
    GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=0.1;FIINIT(W2)=0.1
FIINIT(R1)=0.999;FIINIT(R2)=0.001
    GROUP 13. Boundary conditions and special sources
  ** Wall
WALL (WALL,NORTH,1,1,NY,NY,1,NZ,1,1)
COVAL(WALL,W1,1.0,0.0);COVAL(WALL,H1,1./PRNDTL(H1),1.0)
  ** Inlet
INLET(INLET,LOW,1,1,1,NY,1,1,1,1)
VALUE(INLET,P1,0.1);VALUE(INLET,W1,0.1)
VALUE(INLET,H1,0.0);VALUE(INLET,P2,0.1)
VALUE(INLET,W2,0.1)
  ** Outlet
PATCH(OUTLET,HIGH,1,1,1,NY,NZ,NZ,1,1)
COVAL(OUTLET,P1,FIXP,0.0);COVAL(OUTLET,P2,1.E3,0.0)
COVAL(OUTLET,W1,ONLYMS,0.0);COVAL(OUTLET,W2,ONLYMS,0.0)
  ** Orifice
PATCH(ORIFICE,HIGH,1,1,11,NY,NZ/4,NZ/4,1,1)
COVAL(ORIFICE,W1,FIXVAL,0.0);COVAL(ORIFICE,W2,FIXVAL,0.0)
    GROUP 15. Termination of sweeps
LSWEEP=20;RESREF(P1)=1.E-6;RESREF(V1)=1.E-6;RESREF(H1)=1.E-6
RESREF(R1)=1.E-6;RESREF(V2)=1.E-6;RESREF(H2)=1.E-6
RESREF(W1)=1.E-6;RESREF(W2)=1.E-6;RESREF(R2)=1.E-6
    GROUP 17. Under-relaxation devices
RELAX(W1,FALSDT,1.0);RELAX(V1,FALSDT,1.0)
RELAX(W2,FALSDT,1.0);RELAX(V2,FALSDT,1.0)
RELAX(R1,LINRLX,0.75);RELAX(R2,LINRLX,0.75)
    GROUP 22. Spot-value print-out
NPRMON=20;IYMON=5
    GROUP 22. Spot-value print-out
IYMON=19;IZMON=10;IPLTL=30
    GROUP 23. Field print-out and plot control
NPRINT=10
IYPRF=9;IZPRF=4;IZPRL=8;TSTSWP=-1
PATCH(IZEQ19,PROFIL,1,1,1,NY,NZ-1,NZ-1,1,1)
PLOT(IZEQ19,W1,0.0,0.4);PLOT(IZEQ19,W2,0.0,0.4)
PATCH(CONT,CONTUR,1,1,1,NY,1,NZ,1,1)
PLOT(CONT,W1,0.0,10);PLOT(CONT,W2,0.0,10)
PLOT(CONT,H1,0.0,10);PLOT(CONT,R2,0.0,10)
PLOT(CONT,H2,0.0,10)
PATCH(IZ1TONZ,PROFIL,1,1,1,1,1,NZ,1,1)
PLOT(IZ1TONZ,W1,0.0,0.0);PLOT(IZ1TONZ,W2,0.0,0.0)
PLOT(IZ1TONZ,R2,0.0,0.0)