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