PHOTON USE
  p
  a2
 
 
  msg Circular pipe
  gr z 1
  msg Press  to continue
  pause
  cl
  msg Concentration contour of continuous phase
  msg after 2 seconds
  con pt0 z 1 fi;0.002
  msg Press  to continue
  pause
  cl
  msg Concentration contour of first dispersed phase
  msg after 2 seconds
  con pt1 z 1 fi;0.002
  msg Press  to continue
  pause
  cl
  msg Concentration contour of second dispersed phase
  msg after 2 seconds
  con pt2 z 1 fi;0.002
  msg Press  to continue
  pause
  cl
  p
  a4
 
 
  msg Concentration contour of continuous phase
  msg after 4 seconds
  con pt0 z 1 fi;0.002
  msg Press  to continue
  pause
  cl
  msg Concentration contour of first dispersed phase
  msg after 4 seconds
  con pt1 z 1 fi;0.002
  msg Press  to continue
  pause
  cl
  msg Concentration contour of second dispersed phase
  msg after 4 seconds
  con pt2 z 1 fi;0.002
  msg Press e to END
  enduse
 
    GROUP 1. Run title and other preliminaries
TEXT(SETTLING IN A CIRCULAR PIPE: P111
TITLE
  DISPLAY
  3-dimensional, polar, transient, elliptic simulation
 
  This case calculates the accumulation of two groups of particles,
  of different diameters, on the bottom of a closed circular pipe.
  The first dispersed phase is made up of larger droplets and so
  settles more quickly.
 
  enddis
 
steady=f
lstep=4; tlast=4.0
CARTES=F
nx=10; xulast=2.0*3.1416; ny=5; yvlast=0.25; nz=1; zwlast=0.5
#unigrid
 
char(setaslp); setaslp=$P001
#setaslp
 
    GROUP 7. Variables stored, solved & named
NAME(C3)=PT2
SOLVE(PT2)
SOLUTN(P1,Y,Y,N,N,N,N)
STORE(DEN1,VISL,IMB1)
    GROUP 9. Properties of the medium (or media)
RHO2=1000.0
REAL(dens); dens=2000.0
REAL(ENUP); ENUP=2.5E-6
REAL(VFRAC0,VFRAC1,VFRAC2)
VFRAC1=0.1; VFRAC2=0.1; VFRAC0=1.0-VFRAC1-VFRAC2
    ** Particle properties
   Densities use PHINT(), diameters CINT(), viscosities PRNDTL()
PRNDTL(PT0)=1.E-6
PRNDTL(PT1)=ENUP; PHINT(PT1)=dens; CINT(PT1)=1.E-3
PRNDTL(PT2)=ENUP; PHINT(PT2)=dens; CINT(PT2)=2.E-4
    GROUP 11. Initialisation of variable or porosity fields
FIINIT(P1)=1.E4; FIINIT(U1)=0.0; FIINIT(V1)=0.0
FIINIT(PT0)=VFRAC0; FIINIT(PT1)=VFRAC1; FIINIT(PT2)=VFRAC2
    GROUP 13. Boundary conditions and special sources
PATCH(REFP1,CELL,NX,NX,1,1,1,1,1,LSTEP)
COVAL(REFP1,P1,FIXP,0.0)
    ** Gravity
BUOYB=-9.81
PATCH(BUOY,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(BUOY,U1,FIXFLU,DENSITY)
COVAL(BUOY,V1,FIXFLU,DENSITY)
XCYCLE=T
    GROUP 15. Termination of sweeps
LITER(PT0)=2; LITER(PT1)=5
ENDIT(PT0)=1.E-8; ENDIT(PT1)=1.E-5
LSWEEP=50
RESREF(P1)=1.E-8;RESREF(U1)=1.e-8
RESREF(V1)=1.e-8;RESREF(PT0)=1.E-8
RESREF(PT1)=1.e-8;RESREF(PT2)=1.e-8
    GROUP 17. Under-relaxation devices
RELAX( P1,LINRLX,1.0);RELAX( U1,FALSDT,0.3)
RELAX( V1,FALSDT,0.3);RELAX(PT0,LINRLX,0.2)
RELAX(PT1,LINRLX,0.2)
    GROUP 19. Data communicated by satellite to GROUND
LASLPA=T; LASLPB=T
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
    GROUP 21. Print-out of variables
    GROUP 22. Spot-value print-out
TSTSWP=-1
    GROUP 23. Field print-out and plot control
IPLTF=2;IPLTL=LSWEEP;YZPR=F;IPROF=1;NTPRIN=2
IXMON=NX/2; IYMON=NY/2; NXPRIN=1; NYPRIN=1; NZPRIN=1
    GROUP 24. Dumps for restarts
IDISPA=2; CSG1=A