TALK=T;RUN(1,1)
PHOTON USE
p;parphi
5 1
;;;;
up z
msg pressures
con p1 x 1 fi;0.001;pause;con off;red
msg phase-1 vol fract
con r1 x 1 fi;0.001;pause;con off;red
msg phase-1 radial vel.
con v1 x 1 fi;0.001;pause;con off;red
msg phase-2 radial vel.
con v2 x 1 fi;0.001;pause;con off;red
msg phase-1 axial vel.
con w1 x 1 fi;0.001;pause;con off;red
msg phase-2 axial vel.
con w2 x 1 fi;0.001;pause;con off;red
ENDUSE
DISPLAY
The case considered is 2-phase turbulent air-water flow in a
pipe, as studied experimentally by Seriwaza et al [1992] for
upward flow and by Lahey et al [1992] for both upward and
downward flow.
The input file is set up to run any one of these three cases.
Each calculation is performed with the parabolic option, and,
for testing purposes the calculation is terminated 5 diameters
downstream.
However, for comparison with data the calculation should be
continued to 35 diameters downstream, which corresponds to the
experimental measuring station.
The 2-phase model accounts for interfacial drag, lift, pressure
and virtual-mass forces.
The standard k-e model is employed with the option to select one
of two modifications to account for bubble-induced turbulence.
For upward flow, the predicted and measured void-fraction
profiles show that the gas is taken away from the centre and
towards the walls, while for downward flow the reverse is
observed.
For upward flow the predictions are in reasonable agreement with
the data, but much less so for downward flow. However, the
influences of mesh size, the interfacial-modelling coefficients
and the turbulence-modelling modifications need to be investiga$
ted.
For the results of a study made with a 1995 version of
PHOENICS, click
here.
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(BUBBLY AIR/WATER PIPE- Seriwaza upflow )
************************************************************
************************************************************
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 = 30
YVLAST =0.03
YFRAC(1)=0.033333 ;YFRAC(2)=0.066667
YFRAC(3)=0.1 ;YFRAC(4)=0.133333
YFRAC(5)=0.166667 ;YFRAC(6)=0.2
YFRAC(7)=0.233333 ;YFRAC(8)=0.266667
YFRAC(9)=0.3 ;YFRAC(10)=0.333333
YFRAC(11)=0.366667 ;YFRAC(12)=0.4
YFRAC(13)=0.433333 ;YFRAC(14)=0.466667
YFRAC(15)=0.5 ;YFRAC(16)=0.533333
YFRAC(17)=0.566667 ;YFRAC(18)=0.6
YFRAC(19)=0.633333 ;YFRAC(20)=0.666667
YFRAC(21)=0.7 ;YFRAC(22)=0.733333
YFRAC(23)=0.766667 ;YFRAC(24)=0.8
YFRAC(25)=0.833333 ;YFRAC(26)=0.866667
YFRAC(27)=0.9 ;YFRAC(28)=0.933333
YFRAC(29)=0.966667 ;YFRAC(30)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = T
NZ = 50
ZWLAST =0.3
ZFRAC(1)=0.02 ;ZFRAC(11)=0.22
ZFRAC(21)=0.42 ;ZFRAC(31)=0.62
ZFRAC(41)=0.82
************************************************************
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(144)=LISW
NAME(145)=LISV ;NAME(146)=VMSV
NAME(147)=VMSW ;NAME(148)=REYN
NAME(149)=CD ;NAME(150)=ENUT
* 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(LISW,Y,N,N,N,N,N)
SOLUTN(LISV,Y,N,N,N,N,N)
SOLUTN(VMSV,Y,N,N,N,N,Y)
SOLUTN(VMSW,Y,N,N,N,N,Y)
SOLUTN(REYN,Y,N,N,N,N,Y)
SOLUTN(CD,Y,N,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
VIST = 150
************************************************************
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 =1000. ;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 =1.34 ;TMP2 =0. ;EL2 =0.
DVO2DT =0. ;DRH2DP =0.
ENUL =1.0E-06 ;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.
PRNDTL(W1)=1. ;PRNDTL(W2)=1.
PRNDTL(R1)=1. ;PRNDTL(R2)=1.
PRNDTL(KE)=1. ;PRNDTL(EP)=1.
PRT(V1)=1. ;PRT(V2)=1.
PRT(W1)=1. ;PRT(W2)=1.
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-03 ;CMDOT =0.
CFIPA =1.0E-03 ;CFIPB =3.0E-03
CFIPC =0. ;CFIPD =5.
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)=1.437
FIINIT(W2)=1.437 ;FIINIT(R1)=0.946416
FIINIT(R2)=0.053584 ;FIINIT(KE)=4.827039E-03
FIINIT(EP)=0.018367 ;FIINIT(LISW)=1.0E-10
FIINIT(LISV)=1.0E-10 ;FIINIT(VMSV)=1.0E-10
FIINIT(VMSW)=1.0E-10 ;FIINIT(REYN)=1.0E-10
FIINIT(CD)=1.0E-10 ;FIINIT(ENUT)=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(LIFT ,CELL , 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(LIFT ,V1 ,1.0E-20 , GRND4 )
COVAL(LIFT ,V2 ,1.0E-20 , GRND4 )
COVAL(LIFT ,W1 ,1.0E-20 , GRND4 )
COVAL(LIFT ,W2 ,1.0E-20 , GRND4 )
PATCH(INTPL ,CELL , 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(INTPL ,V1 ,1.0E-20 , GRND4 )
COVAL(INTPL ,V2 ,1.0E-20 , GRND4 )
COVAL(INTPL ,W1 ,1.0E-20 , GRND4 )
COVAL(INTPL ,W2 ,1.0E-20 , GRND4 )
PATCH(IN ,LOW , 1, 1, 1, 30, 1, 1, 1, 1)
COVAL(IN ,V1 ,0. ,0. )
COVAL(IN ,V2 ,0. ,0. )
COVAL(IN ,W1 ,0. ,1.437 )
COVAL(IN ,W2 ,0. ,1.437 )
COVAL(IN ,R1 , FIXFLU ,1360. )
COVAL(IN ,R2 , FIXFLU ,0.10318 )
COVAL(IN ,KE ,0. ,4.827039E-03 )
COVAL(IN ,EP ,0. ,0.018367 )
PATCH(GRAVITY ,PHASEM, 1, 1, 1, 30, 1, 50, 1, 1)
COVAL(GRAVITY ,W2 , FIXFLU ,7311.085449 )
PATCH(NWALL ,NWALL , 1, 1, 30, 30, 1, 50, 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 = 40 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
ISWR1 = 1 ;ISWR2 = 10000
SELREF = F
RESREF(P1)=1.0E-08 ;RESREF(V1)=1.0E-08
RESREF(V2)=1.0E-08 ;RESREF(W1)=1.0E-08
RESREF(W2)=1.0E-08 ;RESREF(R1)=1.0E-08
RESREF(R2)=1.0E-08 ;RESREF(KE)=1.0E-08
RESREF(EP)=1.0E-08
************************************************************
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,FALSDT,2.087683E-06)
RELAX(V2,FALSDT,2.087683E-06)
RELAX(W1,FALSDT,2.087683E-04)
RELAX(W2,FALSDT,2.087683E-04)
RELAX(R1,LINRLX,0.3)
RELAX(R2,LINRLX,0.3)
RELAX(KE,LINRLX,0.3)
RELAX(EP,LINRLX,0.3)
RELAX(LISW,LINRLX,0.1)
RELAX(LISV,LINRLX,0.1)
RELAX(VMSV,LINRLX,1.)
RELAX(VMSW,LINRLX,1.)
RELAX(REYN,LINRLX,1.)
RELAX(CD,LINRLX,1.)
RELAX(ENUT,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)=10. ;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(LISW)=1.0E+10 ;VARMIN(LISW)=-1.0E+10
VARMAX(LISV)=1.0E+10 ;VARMIN(LISV)=-1.0E+10
VARMAX(VMSV)=1.0E+10 ;VARMIN(VMSV)=-1.0E+10
VARMAX(VMSW)=1.0E+10 ;VARMIN(VMSW)=-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
************************************************************
Group 19. Data transmitted to GROUND
GENK = T
PARSOL = F
ISG50 = 1
ISG51 = 2
ISG52 = 2
ISG62 = 1
CVM = GRND2
CVMA =0.5
CPIP = GRND2
CPIPA =0.25
CLIFT = GRND2
CLIFTA =0.075
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(LISW,Y,N,Y,Y,Y,Y)
OUTPUT(LISV,Y,N,Y,Y,Y,Y)
OUTPUT(VMSV,Y,N,Y,Y,Y,Y)
OUTPUT(VMSW,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)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 28 ;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
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = 10 ;IZPRF = 1 ;IZPRL = 10000
IPLTF = 1 ;IPLTL = -1 ;NPLT = 2
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
IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0
STOP