TALK=T;RUN(1,1)
DISPLAY
The case considered is 2-phase solids-laden turbulent vertical
flow of air in a pipe, as studied experimentally by Tsuji et
al [J.Fluid Mech, Vol.139, p417, 1984].
The calculation is performed with the parabolic option.
For comparison with data the calculation should be carried out
until the flow is fully developed, i.e. 70 diameters downstream.
The calculation may be made with the standard k-e model, or
alternatively with the Chen-Wood or Mostafa-Mongia k-e variants
which allow for gas turbulence modulation due to the presence
of particles.
The pipe Reynolds number is 3E4, the density ratio
is 866, the particle mass-flow loading may be 1.0 or 2.1, and
the particle diameter is 200 microns. The task is to calculate
the fully-developed vertical velocity profiles for comparison
with the measured profiles.
For the results of a study made with a 1995 version of
PHOENICS, click
here.
ENDDIS
PHOTON USE
p
parphi
20 1
msg flow is from right to left; wall is at top, axis at bottom
msg radial dimension enlarged 20-fold
msg axial velocity contours for phase 1
con w1 x 1 fi ;.001
pause;con off; red
msg axial velocity contours for phase 2
con w2 x 1 fi;0.001
pause;con off;red
msg volume-fraction contours for phase 2
con r2 x 1 fi;0.001
pause;con off;red
msg interphase-friction factors
con cfip x 1 fi;0.001
pause; con off; red
msg effective (i.e. turbulent) viscosity
con enut x 1 fi;0.001
pause; con off; red
msg axial velocity contours for phase 1
vec x 1 sh
enduse
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(AIR-SOLIDS FLOW- No turb. Mods )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
NX = 1
XULAST =0.1
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 25
YVLAST =0.01525
YFRAC(1)=0.04 ;YFRAC(2)=0.08
YFRAC(3)=0.12 ;YFRAC(4)=0.16
YFRAC(5)=0.2 ;YFRAC(6)=0.24
YFRAC(7)=0.28 ;YFRAC(8)=0.32
YFRAC(9)=0.36 ;YFRAC(10)=0.4
YFRAC(11)=0.44 ;YFRAC(12)=0.48
YFRAC(13)=0.52 ;YFRAC(14)=0.56
YFRAC(15)=0.6 ;YFRAC(16)=0.64
YFRAC(17)=0.68 ;YFRAC(18)=0.72
YFRAC(19)=0.76 ;YFRAC(20)=0.8
YFRAC(21)=0.84 ;YFRAC(22)=0.88
YFRAC(23)=0.92 ;YFRAC(24)=0.96
YFRAC(25)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = T
NZ = 200
ZWLAST =8.54
ZFRAC(1)=5.0E-03 ;ZFRAC(41)=0.205
ZFRAC(81)=0.405 ;ZFRAC(121)=0.605
ZFRAC(161)=0.805
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = F
EQUVEL = F
NAME(1)=P1 ;NAME(5)=V1
NAME(6)=V2 ;NAME(7)=W1
NAME(8)=W2 ;NAME(9)=R1
NAME(10)=R2 ;NAME(12)=KE
NAME(13)=EP ;NAME(147)=REYN
NAME(148)=CD ;NAME(149)=ENUT
NAME(150)=CFIP
* 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(V1,Y,Y,N,N,N,Y)
SOLUTN(V2,Y,Y,N,N,N,Y)
SOLUTN(W1,Y,Y,N,N,N,Y)
SOLUTN(W2,Y,Y,N,N,N,Y)
SOLUTN(R1,Y,Y,N,N,N,Y)
SOLUTN(R2,Y,Y,N,N,N,Y)
SOLUTN(KE,Y,Y,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
SOLUTN(REYN,Y,N,N,N,N,Y)
SOLUTN(CD,Y,N,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
SOLUTN(CFIP,Y,N,N,N,N,Y)
VIST = 149
************************************************************
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(V1,Y,Y,Y,Y,Y,Y)
TERMS(V2,Y,Y,Y,Y,N,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(W2,Y,Y,Y,Y,N,Y)
TERMS(R1,Y,Y,Y,Y,Y,Y)
TERMS(R2,Y,Y,Y,Y,N,Y)
TERMS(KE,N,Y,Y,Y,Y,N)
TERMS(EP,N,Y,Y,Y,Y,N)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWENT = 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 =1.178 ;TMP1 =0. ;EL1 = GRND4
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0. ;EL1B =0. ;EL1C =0.
RHO2 =1020. ;TMP2 =0. ;EL2 =0.
DVO2DT =0. ;DRH2DP =0.
ENUL =1.585739E-05 ;ENUT = GRND3
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PHINT(V1)=-2.022E+04 ;PHINT(V2)=-2.022E+04
PHINT(W1)=-2.022E+04 ;PHINT(W2)=-2.022E+04
PHINT(R1)=-2.022E+04 ;PHINT(R2)=-2.022E+04
PRNDTL(V1)=1. ;PRNDTL(V2)=1.0E+10
PRNDTL(W1)=1. ;PRNDTL(W2)=1.0E+10
PRNDTL(R1)=1. ;PRNDTL(R2)=1.
PRNDTL(KE)=1. ;PRNDTL(EP)=1.
PRT(V1)=1. ;PRT(V2)=1.0E+10
PRT(W1)=1. ;PRT(W2)=1.0E+10
PRT(R1)=1. ;PRT(R2)=1.
PRT(KE)=1. ;PRT(EP)=1.314
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
CFIPS = GRND7
RLOLIM =1.0E-06 ;CMDOT =0.
CFIPA =1.0E-03 ;CFIPB =2.0E-04
CFIPC =0. ;CFIPD =0.
CINT(P1)=1. ;CINT(V1)=1.
CINT(V2)=1. ;CINT(W1)=1.
CINT(W2)=1. ;CINT(R1)=1.
CINT(R2)=1.
PHINT(P1)=-2.022E+04 ;PHINT(V1)=-2.022E+04
PHINT(V2)=-2.022E+04 ;PHINT(W1)=-2.022E+04
PHINT(W2)=-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(V1)=1.0E-10
FIINIT(V2)=1.0E-10 ;FIINIT(W1)=15.597429
FIINIT(W2)=15.597429 ;FIINIT(R1)=0.998845
FIINIT(R2)=1.154902E-03 ;FIINIT(KE)=1.0E-10
FIINIT(EP)=1.0E-10 ;FIINIT(REYN)=1.0E-10
FIINIT(CD)=1.0E-10 ;FIINIT(ENUT)=1.0E-10
FIINIT(CFIP)=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(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(KESOURCE,KE , GRND4 , GRND4 )
COVAL(KESOURCE,EP , GRND4 , GRND4 )
PATCH(IN ,LOW , 1, 1, 1, 25, 1, 1, 1, 1)
COVAL(IN ,V1 ,0. ,0. )
COVAL(IN ,V2 ,0. ,0. )
COVAL(IN ,W1 ,0. ,15.597429 )
COVAL(IN ,W2 ,0. ,15.597429 )
COVAL(IN ,R1 , FIXFLU ,18.352552 )
COVAL(IN ,R2 , FIXFLU ,18.373772 )
COVAL(IN ,KE ,0. ,0.608199 )
COVAL(IN ,EP ,0. ,51.101864 )
PATCH(GRAVITY ,PHASEM, 1, 1, 1, 25, 1, 200, 1, 1)
COVAL(GRAVITY ,W1 , FIXFLU ,-9.81 )
COVAL(GRAVITY ,W2 , FIXFLU ,-9.81 )
PATCH(NWALL ,NWALL , 1, 1, 25, 25, 1, 200, 1, 1)
COVAL(NWALL ,W1 , GRND2 ,0. )
COVAL(NWALL ,KE , GRND2 , GRND2 )
COVAL(NWALL ,EP , GRND2 , GRND2 )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
IPARAB = 0
AZPH =0. ;PBAR =0.
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1 ;ISWC1 = 1
LITHYD = 10 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
ISWR1 = 1 ;ISWR2 = 10000
SELREF = T
RESFAC =1.0E-02
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(V1)=10
LITER(V2)=10 ;LITER(W1)=10
LITER(W2)=10 ;LITER(R1)=1
LITER(R2)=1 ;LITER(KE)=20
LITER(EP)=20
ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
ENDIT(V2)=1.0E-03 ;ENDIT(W1)=1.0E-03
ENDIT(W2)=1.0E-03 ;ENDIT(R1)=1.0E-03
ENDIT(R2)=1.0E-03 ;ENDIT(KE)=1.0E-03
ENDIT(EP)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(V1,LINRLX,0.4)
RELAX(V2,LINRLX,0.4)
RELAX(W1,LINRLX,0.5)
RELAX(W2,LINRLX,0.5)
RELAX(R1,LINRLX,0.3)
RELAX(R2,LINRLX,0.3)
RELAX(KE,LINRLX,0.3)
RELAX(EP,LINRLX,0.3)
RELAX(REYN,LINRLX,1.)
RELAX(CD,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
RELAX(CFIP,LINRLX,1.)
KELIN = 0
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(V2)=1.0E+06 ;VARMIN(V2)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=1.0E-10
VARMAX(W2)=1.0E+06 ;VARMIN(W2)=1.0E-10
VARMAX(R1)=1. ;VARMIN(R1)=1.0E-10
VARMAX(R2)=1. ;VARMIN(R2)=1.0E-10
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
VARMAX(REYN)=1.0E+10 ;VARMIN(REYN)=-1.0E+10
VARMAX(CD)=1.0E+10 ;VARMIN(CD)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
VARMAX(CFIP)=1.0E+10 ;VARMIN(CFIP)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
GENK = T
PARSOL = F
ISG52 = 2
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
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(V1,Y,N,Y,Y,Y,Y)
OUTPUT(V2,Y,N,Y,Y,Y,Y)
OUTPUT(W1,Y,N,Y,Y,Y,Y)
OUTPUT(W2,Y,N,Y,Y,Y,Y)
OUTPUT(R1,Y,N,Y,Y,Y,Y)
OUTPUT(R2,Y,N,Y,Y,Y,Y)
OUTPUT(KE,Y,N,Y,Y,Y,Y)
OUTPUT(EP,Y,N,Y,Y,Y,Y)
OUTPUT(REYN,Y,N,Y,N,N,N)
OUTPUT(CD,Y,N,Y,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,N,N)
OUTPUT(CFIP,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 23 ;IZMON = 1
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -10
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = 40 ;IZPRF = 1 ;IZPRL = 10000
IPLTF = 1 ;IPLTL = -1 ;NPLT = 4
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
PATCH(ENDPROF ,PROFIL, 1, 1, 1, 25, 1, 200, 1, 1)
PLOT(ENDPROF ,W1 ,0. ,0. )
PLOT(ENDPROF ,W2 ,0. ,0. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
IDISPA = 10 ;IDISPB = 0 ;IDISPC = 0
STOP