AUTOPLOT USE
  file
  phi 5
 
  da 1
  u1
  screen
  msg 1st-phase velocity
  pl 1
  msg press  to continue
  pause
  cl
  da 1
  u2
  screen
  msg 2nd-phase velocity
  pl 1
  msg press  to continue
  pause
  cl
  da 1
  r1
  screen
  msg 1st-phase volume fraction
  pl 1
  msg press  to continue
  pause
  cl
  da 1
  c1
  screen
  msg 1st-group of particles
  pl 1
  msg press  to continue
  pause
  cl
  da 1
  p1
  screen
  msg pressure
  pl 1
  msg press e to END
  enduse
TEXT(1D PARTICLE FLOW (cint=grnd7 test): W571
TITLE
    GROUP 3. X-direction grid specification
GRDPWR(X,10,1,1)
    GROUP 7. Variables stored, solved & named
ONEPHS=F;SOLVE(P1,U1,U2,R1,R2,C1,C2)
    Activate the following line if storage for the properties is
    required
    STORE(NUSS,REYN,VREL,CD,APRJ,CFIP)
    GROUP 9. Properties of the medium (or media)
RHO1=2;RHO2=1500
    Note that ENUL=LINTEM (GRND1) selects ENUL=ENULA+ENULB*TMP1
ENUL=LINTEM;ENULA=1E-3;ENULB=0.0;STORE(TMP1);TMP1=CONST
PRNDTL(C1)=-GRND1;STORE(PRL);FIINIT(PRL)=1.E-3/0.7
    GROUP 10. Inter-phase-transfer processes and properties
CFIPS=GRND7;CFIPA=0.01;CFIPB=1E-3;STORE(SIZE)
CINT(C1)=GRND7;CINT(C2)=1E15
    GROUP 11. Initialization of variable or porosity fields
FIINIT(R1)=0.5;FIINIT(R2)=0.5
    GROUP 13. Boundary conditions and special sources
REAL(UIN1,UIN2,R1IN,R2IN)
UIN1=20;UIN2=10;R1IN=0.25;R2IN=1-R1IN
INLET(IN,WEST,1,1,1,1,1,1,1,1)
VALUE(IN,P1,RHO1*UIN1*R1IN)
VALUE(IN,U1,UIN1);VALUE(IN,C1,1.0)
VALUE(IN,P2,RHO2*UIN2*R2IN)
VALUE(IN,U2,UIN2);VALUE(IN,C2,0.0)
OUTLET(OUT,EAST,NX,NX,1,1,1,1,1,1)
COVAL(OUT,P1,1E3,0);COVAL(OUT,P2,RHO2/(2*RHO1)*1E3,0)
    GROUP 15. Termination of sweeps
LSWEEP=40
    GROUP 17. Under-relaxation devices
RELAX(R1,LINRLX,1);RELAX(R2,LINRLX,1)
    GROUP 21. Print-out of variables
NXPRIN=1
    GROUP 22. Spot-value print-out
IXMON=7;TSTSWP=-1
    GROUP 24. Dumps for restarts