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
 ************************************************************
  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