TALK=T;RUN(1,1)
  AUTOPLOT USE
  file
  phi 5


  da 1 u1; da 1 uliq
  col3 1;blb4 2;scale y 0 .015
  redr
  msg 1st-phase velocity
  msg blue PHOENICS   + analytical
  msg press  to continue
  pause
  cl
  screen
  da 1 r1; da 1 rliq
  col3 1;blb4 2;scale y 0 1.
  redr
  msg 1st-phase volume fraction
  msg blue PHOENICS   + analytical
  msg press  to continue
  pause
  cl
  msg press e to END
  enduse
  DISPLAY
  The problem considered is the one-dimensional vertical flow of
  spherical particles by water in a duct. The inlet void fraction
  is 0.4 and the ratio of the solids mass inflow rate to that of
  the gas is 0.5. The duct area is 1 m**2 and the height is 1m.
  The particles are 100 microns in diameter and the density ratio
  is 2.6. The inlet liquid velocity is taken as twice the particle
  terminal velocity. The task is to predict the volume fractions
  and phase velocities throughout the duct. The calculations employ
  the built-in particle-fluidisation drag model (CFIPD=7.), which
  uses the Ergun correlation for void fractions less than 0.8
  (dense fluidisation), and a modified spherical-particle drag
  correlation for void fraction greater than 0.8 (dilute
  fluidisation). The test case provides a test of both dense and
  dilute suspensions according to the setting of DENSE. In both
  cases the PHOENICS predictions show excellent agreement with the
  approximate analytical solutions.
  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(1D SOLIDS TRANSPORT IN WATER            )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 20
 XULAST =1.
 XFRAC(1)=0.05 ;XFRAC(2)=0.1
 XFRAC(3)=0.15 ;XFRAC(4)=0.2
 XFRAC(5)=0.25 ;XFRAC(6)=0.3
 XFRAC(7)=0.35 ;XFRAC(8)=0.4
 XFRAC(9)=0.45 ;XFRAC(10)=0.5
 XFRAC(11)=0.55 ;XFRAC(12)=0.6
 XFRAC(13)=0.65 ;XFRAC(14)=0.7
 XFRAC(15)=0.75 ;XFRAC(16)=0.8
 XFRAC(17)=0.85 ;XFRAC(18)=0.9
 XFRAC(19)=0.95 ;XFRAC(20)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 1
 YVLAST =1.
 YFRAC(1)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 1
 ZWLAST =1.
 ZFRAC(1)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = F
 EQUVEL = F
 NAME(1)=P1 ;NAME(3)=U1
 NAME(4)=U2 ;NAME(9)=R1
 NAME(10)=R2 ;NAME(144)=TMP1
 NAME(145)=RLIQ ;NAME(146)=ULIQ
 NAME(147)=SIZE ;NAME(148)=CFIP
 NAME(149)=VREL ;NAME(150)=REYN
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(P1,Y,Y,N,N,N,Y)
 SOLUTN(U1,Y,Y,N,N,N,Y)
 SOLUTN(U2,Y,Y,N,N,N,Y)
 SOLUTN(R1,Y,Y,N,N,N,Y)
 SOLUTN(R2,Y,Y,N,N,N,Y)
 SOLUTN(TMP1,Y,N,N,N,N,Y)
 SOLUTN(RLIQ,Y,N,N,N,N,Y)
 SOLUTN(ULIQ,Y,N,N,N,N,Y)
 SOLUTN(SIZE,Y,N,N,N,N,Y)
 SOLUTN(CFIP,Y,N,N,N,N,Y)
 SOLUTN(VREL,Y,N,N,N,N,Y)
 SOLUTN(REYN,Y,N,N,N,N,Y)
 TEMP1 = 144
 ************************************************************
  Group 8. Terms & Devices
    * Y in TERMS argument list denotes:
    * 1-built-in source 2-convection 3-diffusion 4-transient
    * 5-first phase variable 6-interphase transport         
 TERMS(P1,Y,Y,Y,N,Y,Y)
 TERMS(U1,Y,Y,Y,Y,Y,Y)
 TERMS(U2,Y,Y,Y,Y,N,Y)
 TERMS(R1,Y,Y,Y,Y,Y,Y)
 TERMS(R2,Y,Y,Y,Y,N,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWENL = T
 EQDVDP = F
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 =1000. ;TMP1 = GRND1
 EL1 =0.
 TSURR =0. ;TEMP0 =0.
 PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 TMP1A =0. ;TMP1B =0.
 TMP1C =0.
 TMP2A =0. ;TMP2B =0.
 TMP2C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 RHO2 =2600. ;TMP2 =0.
 EL2 =0.
 DVO2DT =0. ;DRH2DP =0.
 ENUL = GRND1 ;ENUT =0.
 ENULA =1.0E-06 ;ENULB =0.
 ENULC =0.
 PHINT(U1)=-2.022E+04 ;PHINT(U2)=-2.022E+04
 PHINT(R1)=-2.022E+04 ;PHINT(R2)=-2.022E+04
 PRNDTL(U1)=1. ;PRNDTL(U2)=1.
 PRNDTL(R1)=1. ;PRNDTL(R2)=1.
 PRT(U1)=1. ;PRT(U2)=1.
 PRT(R1)=1. ;PRT(R2)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 CFIPS = GRND7
 RLOLIM =0. ;CMDOT =0.
 CFIPA =1.0E-04 ;CFIPB =1.0E-04
 CFIPC =0. ;CFIPD =7.
 CINT(P1)=1. ;CINT(U1)=1.
 CINT(U2)=1. ;CINT(R1)=1.
 CINT(R2)=1.
 PHINT(P1)=-2.022E+04 ;PHINT(U1)=-2.022E+04
 PHINT(U2)=-2.022E+04 ;PHINT(R1)=-2.022E+04
 PHINT(R2)=-2.022E+04
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(U1)=0.01744
 FIINIT(U2)=2.235897E-03 ;FIINIT(R1)=0.6
 FIINIT(R2)=0.4 ;FIINIT(TMP1)=1.0E-10
 FIINIT(RLIQ)=0.779187 ;FIINIT(ULIQ)=8.952922E-03
 FIINIT(SIZE)=1.0E-10 ;FIINIT(CFIP)=1.0E-10
 FIINIT(VREL)=1.0E-10 ;FIINIT(REYN)=1.0E-10
   No PATCHes yet used for this Group
 INIADD = F
 FSWEEP = 1
 NAMFI =CHAM
 ************************************************************
  Group 12. Patchwise adjustment of terms
  Patches for this group are printed with those
  for Group 13.
  Their names begin either with GP12 or &
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH(IN ,CELL , 1, 1, 1, 1, 1, 1, 1, 1)
 COVAL(IN ,U1 ,0. ,0.01744 )
 COVAL(IN ,U2 ,0. ,2.235897E-03 )
 COVAL(IN ,R1 , FIXFLU ,6.976 )
 COVAL(IN ,R2 , FIXFLU ,3.488 )
 
 PATCH(OUT ,CELL , 20, 20, 1, 1, 1, 1, 1, 1)
 COVAL(OUT ,U1 ,0. ,0. )
 COVAL(OUT ,U2 ,0. ,0. )
 COVAL(OUT ,R1 ,895.292175 ,0. )
 COVAL(OUT ,R2 ,1579.616943 ,0. )
 
 PATCH(GRAVITY ,PHASEM, 1, 20, 1, 1, 1, 1, 1, 1)
 COVAL(GRAVITY ,U2 , FIXFLU ,-6.036923 )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 250 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 ISWR1 = 1 ;ISWR2 = 10000
 SELREF = T
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=10
 LITER(U2)=10 ;LITER(R1)=1
 LITER(R2)=1
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(U2)=1.0E-03 ;ENDIT(R1)=1.0E-03
 ENDIT(R2)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,0.286697)
 RELAX(U2,FALSDT,0.286697)
 RELAX(R1,LINRLX,0.3)
 RELAX(R2,LINRLX,0.3)
 RELAX(TMP1,LINRLX,1.)
 RELAX(RLIQ,LINRLX,1.)
 RELAX(ULIQ,LINRLX,1.)
 RELAX(SIZE,LINRLX,1.)
 RELAX(CFIP,LINRLX,0.3)
 RELAX(VREL,LINRLX,1.)
 RELAX(REYN,LINRLX,1.)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
 VARMAX(U2)=1.0E+06 ;VARMIN(U2)=-1.0E+06
 VARMAX(R1)=1. ;VARMIN(R1)=1.0E-06
 VARMAX(R2)=1. ;VARMIN(R2)=1.0E-10
 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=-1.0E+10
 VARMAX(RLIQ)=1.0E+10 ;VARMIN(RLIQ)=-1.0E+10
 VARMAX(ULIQ)=1.0E+10 ;VARMIN(ULIQ)=-1.0E+10
 VARMAX(SIZE)=1.0E+10 ;VARMIN(SIZE)=-1.0E+10
 VARMAX(CFIP)=1.0E+10 ;VARMIN(CFIP)=-1.0E+10
 VARMAX(VREL)=1.0E+10 ;VARMIN(VREL)=-1.0E+10
 VARMAX(REYN)=1.0E+10 ;VARMIN(REYN)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 PARSOL = F
 ISG62 = 1
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 INIFLD = F ;SUBWGR = F
    * Y in OUTPUT argument list denotes:
    * 1-field 2-correction-eq. monitor 3-selective dumping      
    * 4-whole-field residual 5-spot-value table 6-residual table
 OUTPUT(P1,Y,N,Y,Y,Y,Y)
 OUTPUT(U1,Y,N,Y,Y,Y,Y)
 OUTPUT(U2,Y,N,Y,Y,Y,Y)
 OUTPUT(R1,Y,N,Y,Y,Y,Y)
 OUTPUT(R2,Y,N,Y,Y,Y,Y)
 OUTPUT(TMP1,Y,N,Y,N,N,N)
 OUTPUT(RLIQ,Y,N,Y,N,N,N)
 OUTPUT(ULIQ,Y,N,Y,N,N,N)
 OUTPUT(SIZE,Y,N,Y,N,N,N)
 OUTPUT(CFIP,Y,N,Y,N,Y,N)
 OUTPUT(VREL,Y,N,Y,N,N,N)
 OUTPUT(REYN,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 10 ;IYMON = 1 ;IZMON = 1
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NXPRIN = 1 ;IXPRF = 1 ;IXPRL = 10000
 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.4
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
   No PATCHes yet used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP