```
TALK=T;RUN( 1, 1)
DISPLAY

TWO-PHASE FLOW IN A BIFURCATION

The turbulent two-phase flow in bifurcating channel is
solved here on Cartesian grid by "cut-off" technique of PARSOL.

The case is aimed to demonstrate the PARSOL performance for
two-phase flows with a phase slip.

The demonstration case considers the turbulent flow of dispersed
air/water mixture in the branching channel designed as
two-dimensional asymmetrical bifurcation.

In this case the pressure distribution, phase volume fractions,
velocity fields of air and particulate water have to be
calculated.

ENDDIS

************************************************************
Q1 created by VDI menu, Version 3.3, Date 24/07/00
CPVNAM=VDI;SPPNAM=Core
************************************************************
IRUNN   =       1 ;LIBREF =     290
************************************************************
Group 1. Run Title
TEXT(Two-phase flow in a bifurcation)
************************************************************
Group 2. Transience
STEADY  =    T
************************************************************
Groups 3, 4, 5  Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,60,50,1,7.620000E-05)
* Set overall domain extent:
*        xulast  yvlast  zwlast
name
XSI= 1.200000E+00; YSI= 5.000000E-01; ZSI= 2.000000E-01
RSET(D,CHAM    )
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS  =    F
* Non-default variable names
NAME(147) =SMOK ; NAME(148) =EPKE
NAME(149) =LEN1 ; NAME(150) =ENUT
* Solved variables list
SOLVE(P1  ,U1  ,U2  ,V1  ,V2  ,R1  ,R2  ,SMOK)
* Stored variables list
STORE(ENUT,LEN1,EPKE)
* Additional solver options
SOLUTN(P1  ,Y,Y,Y,N,N,N)
TURMOD(KEMODL)

************************************************************
Group 8. Terms & Devices
TERMS (R1  ,Y,Y,N,N,Y,Y)
TERMS (R2  ,Y,Y,N,N,N,Y)
************************************************************
Group 9. Properties
RHO1    = 2.235000E+01
CP1     = 1.005000E+03
RHO2    = 1.279000E+03
CP2     = 4.181800E+03
ENUL    = 1.100667E-05
DVO1DT  = 3.330000E-03
PRT   (EP  ) =  1.314000E+00
************************************************************
Group 10.Inter-Phase Transfer Processes
CFIPS   = GRND7
RLOLIM  = 0.000000E+00 ;CMDOT  = 0.000000E+00
CFIPA   = 1.000000E-03 ;CFIPB  = 1.000000E-04
CMDTA   = 5.000000E-01 ;CMDTB  = 2.500000E-01
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1  ) =  1.300000E-04 ;FIINIT(U1  ) =  1.300000E+01
FIINIT(U2  ) =  1.300000E+01 ;FIINIT(R1  ) =  9.000000E-01
FIINIT(R2  ) =  1.000000E-01 ;FIINIT(KE  ) =  7.605000E-01
FIINIT(EP  ) =  3.177743E+01 ;FIINIT(SMOK) =  1.000000E+00
FIINIT(EPKE) =  1.000000E+00
No PATCHes used for this Group

INIADD  =    F
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources

PATCH (BUOYANCY,PHASEM,0,0,0,0,0,0,1,1)
COVAL (BUOYANCY,U1  , FIXFLU      , GRND1       )
COVAL (BUOYANCY,V1  , FIXFLU      , GRND1       )

BUOYA   = 0.000000E+00 ; BUOYB =-9.810000E+00
BUOYC   = 0.000000E+00
EGWF    =    T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP  =   10000
RESFAC  = 1.000000E-05
************************************************************
Group 16. Terminate Iterations
LITER (U1  ) =   20 ;LITER (V1  ) =   20
LITER (R1  ) =    1 ;LITER (R2  ) =    1
************************************************************
Group 17. Relaxation
RELAX(P1  ,LINRLX, 1.000000E+00)
RELAX(U1  ,FALSDT, 5.000000E-04)
RELAX(U2  ,FALSDT, 5.000000E-04)
RELAX(V1  ,FALSDT, 5.000000E-04)
RELAX(V2  ,FALSDT, 5.000000E-04)
RELAX(KE  ,LINRLX, 5.000000E-01)
RELAX(EP  ,LINRLX, 5.000000E-01)
RELAX(SMOK,LINRLX, 1.000000E+00)
KELIN   =       3
************************************************************
Group 18. Limits
VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
VARMAX(U2  ) = 1.000000E+06 ;VARMIN(U2  ) =-1.000000E+06
VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
VARMAX(V2  ) = 1.000000E+06 ;VARMIN(V2  ) =-1.000000E+06
VARMAX(R1  ) = 1.000000E+00 ;VARMIN(R1  ) = 1.000000E-06
VARMAX(R2  ) = 1.000000E+00 ;VARMIN(R2  ) = 1.000000E-06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD  =    T  ;USEGRX =    T
GENK    =    T
ASAP    =    T
LSG3    =    T
CSG3    ='FGEM'
************************************************************
Group 20. Preliminary Printout
ECHO    =    T
************************************************************
Group 21. Print-out of Variables
OUTPUT(U2  ,Y,Y,Y,N,Y,Y)
OUTPUT(V2  ,Y,Y,Y,N,Y,Y)
OUTPUT(SMOK,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON   =       4 ;IYMON  =      26 ;IZMON  =       1
NPRMON  =     100
NPRMNT  =       1
TSTSWP  =     -12
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT  =  100000
NXPRIN  =       1
NYPRIN  =       1
IPLTF   =       1 ;IPLTL  =     100
NPLT    =       1
ISWPRF  =       1 ;ISWPRL =  100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
NOWIPE  =    T

GVIEW(P,-3.899012E-02,0.000000E+00,9.992399E-01)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)

> DOM,    SIZE,        1.200000E+00, 5.000000E-01, 2.000000E-01
> DOM,    MONIT,       7.000000E-02, 2.550000E-01, 1.000000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-03
> GRID,   RSET_Z_1,      1, 1.000000E+00
> DOM,    VARIABLES,  SMOK
> DOM,    RELAX,       1.000000E+00

> OBJ1,   NAME,        OUTLET
> OBJ1,   POSITION,    1.200000E+00, 0.000000E+00, 0.000000E+00
> OBJ1,   SIZE,        0.000000E+00, 1.540000E-01, 2.000004E-01
> OBJ1,   CLIPART,     cube12t
> OBJ1,   ROTATION24,        1
> OBJ1,   GRID,        2.000000E+00
> OBJ1,   TYPE,        OUTLET
> OBJ1,   PRESSURE,      0.000000E+00
> OBJ1,   TEMPERATURE,  -1.026000E+04
> OBJ1,   COEFFICIENT,   2.235000E+01
> OBJ1,   TEMPERATURE-2, 0.000000E+00
> OBJ1,   COEFFICIENT-2, 1.279000E+03
> OBJ1,   TURBULENCE,   -1.026000E+04,-1.026000E+04
> OBJ1,   SMOKE,         0.000000E+00

> OBJ2,   NAME,        IN3
> OBJ2,   POSITION,    0.000000E+00, 2.000000E-01, 0.000000E+00
> OBJ2,   SIZE,        0.000000E+00, 1.000000E-01, 2.000004E-01
> OBJ2,   CLIPART,     cube3t
> OBJ2,   ROTATION24,        1
> OBJ2,   GRID,        2.000000E+00
> OBJ2,   TYPE,        INLET
> OBJ2,   DENSITY,       2.011500E+01
> OBJ2,   DENSITY-2,     1.279000E+02
> OBJ2,   VELOCITY,      5.000000E+01, 0.000000E+00, 0.000000E+00
> OBJ2,   VELOCITY-2,    5.000000E+01, 0.000000E+00, 0.000000E+00
> OBJ2,   TEMPERATURE,   0.000000E+00
> OBJ2,   TEMPERATURE-2, 0.000000E+00
> OBJ2,   SMOKE,         0.000000E+00
> OBJ2,   TURB-INTENS,   5.000000E+00

> OBJ3,   NAME,        B10
> OBJ3,   POSITION,    1.200000E+00, 4.300000E-01, 0.000000E+00
> OBJ3,   SIZE,        0.000000E+00, 7.000000E-02, 2.000004E-01
> OBJ3,   CLIPART,     cube12t
> OBJ3,   ROTATION24,        1
> OBJ3,   GRID,        2.000000E+00
> OBJ3,   TYPE,        OUTLET
> OBJ3,   PRESSURE,      0.000000E+00
> OBJ3,   TEMPERATURE,  -1.026000E+04
> OBJ3,   COEFFICIENT,   2.235000E+01
> OBJ3,   TEMPERATURE-2, 0.000000E+00
> OBJ3,   COEFFICIENT-2, 1.279000E+03
> OBJ3,   TURBULENCE,   -1.026000E+04,-1.026000E+04
> OBJ3,   SMOKE,         0.000000E+00

> OBJ4,   NAME,        B11
> OBJ4,   POSITION,    5.199996E-01, 2.200000E-01, 0.000000E+00
> OBJ4,   SIZE,        6.800004E-01, 2.100000E-01, 2.000004E-01
> OBJ4,   CLIPART,     wedge1
> OBJ4,   ROTATION24,       18
> OBJ4,   GRID,        2.000000E+00
> OBJ4,   TYPE,        BLOCKAGE
> OBJ4,   MATERIAL,      103
> OBJ4,   FIXED_TM1,     0.000000E+00, 0.000000E+00
> OBJ4,   ADIABATI2,     0.000000E+00, 0.000000E+00

> OBJ5,   NAME,        B17
> OBJ5,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ5,   SIZE,        1.200000E+00, 2.000000E-01, 2.000004E-01
> OBJ5,   CLIPART,     wedge1
> OBJ5,   ROTATION24,        9
> OBJ5,   GRID,        2.000000E+00
> OBJ5,   TYPE,        BLOCKAGE
> OBJ5,   MATERIAL,      103
> OBJ5,   FIXED_TM1,     0.000000E+00, 0.000000E+00
> OBJ5,   ADIABATI2,     0.000000E+00, 0.000000E+00

> OBJ6,   NAME,        B18
> OBJ6,   POSITION,    0.000000E+00, 3.000000E-01, 0.000000E+00
> OBJ6,   SIZE,        1.200000E+00, 2.000000E-01, 2.000004E-01
> OBJ6,   CLIPART,     wedge1
> OBJ6,   ROTATION24,       11
> OBJ6,   GRID,        2.000000E+00
> OBJ6,   TYPE,        BLOCKAGE
> OBJ6,   MATERIAL,      103
> OBJ6,   FIXED_TM1,     0.000000E+00, 0.000000E+00
> OBJ6,   ADIABATI2,     0.000000E+00, 0.000000E+00

> OBJ7,   NAME,        B7
> OBJ7,   POSITION,    5.199996E-01, 1.540000E-01, 0.000000E+00
> OBJ7,   SIZE,        6.800004E-01, 6.600000E-02, 2.000004E-01
> OBJ7,   CLIPART,     wedge1
> OBJ7,   ROTATION24,       23
> OBJ7,   GRID,        2.000000E+00
> OBJ7,   TYPE,        BLOCKAGE
> OBJ7,   MATERIAL,      103
> OBJ7,   FIXED_TM1,     0.000000E+00, 0.000000E+00
> OBJ7,   ADIABATI2,     0.000000E+00, 0.000000E+00
STOP
```