TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2013, Date 10/01/13
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
This case is 45 degree duct channel with cut-cell forming
its friction walls. The pressure drop can be calculated
analytically as:
dP = 12*Mu*L*U/(H^2)
where Mu is viscosity, 0.01 [kg/m.s]
L is duct length, 13.4 [m]
H is duct gap, 0.707 [m]
U is mean velocity, 1.41e-2 [m/s]
dP = 4.53e-2
Click here to see P1 from simulation
The velocity solution appeared in very good agreement with
analytical solution.
Click here to see the comparision
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 142
************************************************************
Group 1. Run Title
TEXT(Inclined Channel; grid 152,152,1
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,152,152,1)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
NAME(150) =PRPS
* Solved variables list
SOLVE(P1,U1,V1)
* Stored variables list
STORE(PRPS)
* Additional solver options
SOLUTN(P1,Y,Y,Y,N,N,Y)
SOLUTN(U1,Y,Y,Y,N,N,Y)
SOLUTN(V1,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
RHO1 =1.
ENUL =1.0E-02
CP1 =1.
ENUT =0.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(PRPS)=-1.
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
No PATCHes used for this Group
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 500
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
CONWIZ = T
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=2.53E-03 ;EX(U1)=1.019E-03
EX(V1)=1.019E-03 ;EX(PRPS)=1.077E-01
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 74 ;IYMON = 73 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,0.,0.,1.)
GVIEW(UP,0.,1.,0.)
GVIEW(DEPTH,46.221096)
GVIEW(NEARPLANE,0.02)
GVIEW(VDIS,6.805175)
GVIEW(CENTRE,5.,5.,0.)
> DOM, SIZE, 1.000000E+01, 1.000000E+01, 1.000000E+00
> DOM, MONIT, 4.843750E+00, 4.781250E+00, 5.000000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID, RSET_X_1, 4, 1.000000E+00
> GRID, RSET_X_2, 144, 1.000000E+00
> GRID, RSET_X_3, 4, 1.000000E+00
> GRID, RSET_Y_1, 4, 1.000000E+00
> GRID, RSET_Y_2, 144, 1.000000E+00
> GRID, RSET_Y_3, 4, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> OBJ, NAME, B3
> OBJ, POSITION, 9.500000E+00, 1.000000E+01, 0.000000E+00
> OBJ, SIZE, 5.000000E-01, 0.000000E+00, 1.000000E+00
> OBJ, GEOMETRY, cube12t
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.
> OBJ, COEFFICIENT, 1000.
> OBJ, NAME, B5
> OBJ, POSITION, 0.000000E+00, 5.000000E-01, 0.000000E+00
> OBJ, SIZE, 9.500000E+00, 9.500000E+00, 1.000000E+00
> OBJ, GEOMETRY, wedge1
> OBJ, ROTATION24, 6
> OBJ, GRID, N,N,N
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, B6
> OBJ, POSITION, 5.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 9.500000E+00, 9.500000E+00, 1.000000E+00
> OBJ, GEOMETRY, wedge1
> OBJ, ROTATION24, 18
> OBJ, GRID, N,N,N
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, B7
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 5.000000E-01, 0.000000E+00, 1.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 11
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.
> OBJ, VELOCITY, 1.0E-02 ,1.0E-02 ,0.
> OBJ, NAME, B8
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 5.000000E-01, 1.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.
> OBJ, VELOCITY, 1.0E-02 ,1.0E-02 ,0.
> OBJ, NAME, B9
> OBJ, POSITION, 1.000000E+01, 9.500000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 5.000000E-01, 1.000000E+00
> OBJ, GEOMETRY, cube12t
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.
> OBJ, COEFFICIENT, 1000.
STOP