TALK=T;RUN( 1, 1)
DISPLAY
3-dimensional, polar, steady, elliptic simulation
This case calculates the accumulation of two groups of particles,
of different diameters, along the bottom of a circular pipe
through which liquid flows from low z to high z. The particle
slip velocities are normal to the continuous phase flow direction.
The larger particles settle to the bottom of the pipe, forcing
the smaller ones to rise with the displaced fluid.
ENDDIS
************************************************************
Q1 created by VDI menu, Version 2022, Date 24/08/22
CPVNAM=VDI; SPPNAM=Core
************************************************************
IRUNN = 1 ;LIBREF = 30
************************************************************
Group 1. Run Title
TEXT(3D flow in a circular pipe: P110 )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,20,20,30)
* Cylindrical-polar grid
CARTES=F
************************************************************
Group 6. Body-Fitted coordinates
* X-cyclic boundaries switched
XCYIZ( 1, T ) ; XCYIZ( 2, T ) ; XCYIZ( 3, T )
XCYIZ( 4, T ) ; XCYIZ( 5, T ) ; XCYIZ( 6, T )
XCYIZ( 7, T ) ; XCYIZ( 8, T ) ; XCYIZ( 9, T )
XCYIZ( 10, T ) ; XCYIZ( 11, T ) ; XCYIZ( 12, T )
XCYIZ( 13, T ) ; XCYIZ( 14, T ) ; XCYIZ( 15, T )
XCYIZ( 16, T ) ; XCYIZ( 17, T ) ; XCYIZ( 18, T )
XCYIZ( 19, T ) ; XCYIZ( 20, T ) ; XCYIZ( 21, T )
XCYIZ( 22, T ) ; XCYIZ( 23, T ) ; XCYIZ( 24, T )
XCYIZ( 25, T ) ; XCYIZ( 26, T ) ; XCYIZ( 27, T )
XCYIZ( 28, T ) ; XCYIZ( 29, T ) ; XCYIZ( 30, T )
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(16)=PT0 ;NAME(17)=PT1
NAME(18)=PT2 ;NAME(147)=IMB1
NAME(148)=VFOL ;NAME(149)=VISL
NAME(150) =DEN1
* Solved variables list
SOLVE(P1,U1,V1,W1,PT0,PT1,PT2,VFOL)
* Stored variables list
STORE(DEN1,VISL,IMB1)
************************************************************
Group 8. Terms & Devices
TERMS (VFOL,N,N,N,N,Y,Y)
GALA = T
NEWRH1 = T
NEWENL = T
************************************************************
Group 9. Properties
RHO1 = GRND ;RHO2 =1000.
ENUL = GRND
CP1 =1.
ENUT =0.
PHINT(PT0)=1000. ;PHINT(PT1)=2000.
PHINT(PT2)=2000.
PRNDTL(PT0)=1.0E-05 ;PRNDTL(PT1)=2.5E-06
PRNDTL(PT2)=2.5E-06
************************************************************
Group 10.Inter-Phase Transfer Processes
CINT(PT1)=1.0E-04 ;CINT(PT2)=2.0E-05
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1)=1.0E+04 ;FIINIT(U1)=0.
FIINIT(V1)=0. ;FIINIT(W1)=1.
FIINIT(PT0)=0.8 ;FIINIT(PT1)=0.1
FIINIT(PT2)=0.1 ;FIINIT(IMB1)=1.001E-10
FIINIT(VFOL)=1.001E-10 ;FIINIT(VISL)=1.001E-10
FIINIT(DEN1)=1.001E-10
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)
COVAL(BUOYANCY, W1, FIXFLU, GRND1)
BUOYA =0. ; BUOYB =-9.81 ; BUOYC =0.
XCYCLE = T
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 500
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(P1)=200 ;LITER(PT0)=10
LITER(PT1)=10 ;LITER(PT2)=10
LITER(VFOL)=10
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,0.5 )
RELAX(U1 ,LINRLX,0.5 )
RELAX(V1 ,LINRLX,0.5 )
RELAX(W1 ,LINRLX,0.5 )
RELAX(PT0 ,LINRLX,0.5 )
RELAX(PT1 ,LINRLX,0.5 )
RELAX(PT2 ,LINRLX,0.5 )
RELAX(VFOL,LINRLX,0.5 )
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
ASLP = T
PARSOL = F
CONWIZ = T
LASLPA = T ;LASLPB = T
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=1.2E+04 ;EX(U1)=4.59E-03
EX(V1)=3.302E-03 ;EX(W1)=1.003
EX(PT0)=0.8 ;EX(PT1)=0.09999
EX(PT2)=0.1 ;EX(IMB1)=5.243E-09
EX(VFOL)=1.001E-10 ;EX(VISL)=8.5E-06
EX(DEN1)=1200.
************************************************************
Group 21. Print-out of Variables
OUTPUT(P1 ,Y,Y,Y,Y,Y,Y)
OUTPUT(U1 ,Y,Y,Y,Y,Y,Y)
OUTPUT(V1 ,Y,Y,Y,Y,Y,Y)
OUTPUT(W1 ,Y,Y,Y,Y,Y,Y)
OUTPUT(VFOL,N,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 10 ;IYMON = 6 ;IZMON = 13
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NXPRIN = 1
NYPRIN = 1
NZPRIN = 1
IPLTF = 2 ;IPLTL = 35
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
IDISPA = 0 ;IDISPB = 0 ;IDISPC = 0
GVIEW(P,0.539702,0.564686,0.624381)
GVIEW(UP,-0.395026,0.824819,-0.404509)
GVIEW(VDIS,2.661729)
GVIEW(CENTRE,0.,0.,2.75)
> DOM, SIZE, 6.283185E+00, 2.500000E-01, 5.500000E+00
> DOM, MONIT, 2.984513E+00, 6.875000E-02, 2.291667E+00
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID, RSET_X_1, 20, 1.000000E+00
> GRID, RSET_Y_1, 20, 1.000000E+00
> GRID, RSET_Z_1, 30, 1.000000E+00
> DOM, T_AMBIENT, 0.000000E+00
> OBJ, NAME, ASMIN
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 6.283200E+00, 2.500000E-01, 0.000000E+00
> OBJ, GEOMETRY, polbox
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, VELOCITY, 0. ,0. ,1.
> OBJ, INLET_PT0, 0.8
> OBJ, INLET_PT1, 0.1
> OBJ, INLET_PT2, 0.1
> OBJ, NAME, ASMOU1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 5.500000E+00
> OBJ, SIZE, 6.283200E+00, 2.500000E-01, 0.000000E+00
> OBJ, GEOMETRY, polbox
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.
> OBJ, COEFFICIENT, 0.1
STOP