TEXT( HORIZENTAL DUCT
DISPLAY
This is the case of horizental duct with plate at the bottom
and blockage at the top, the grid was made to cut the object
of blockage, so the frition wall of of blcokage is located
in cut-cell. The preesure drop could be calculated 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 velocity contour
Click here to see the comparision
ENDDIS
TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.5, Date 16/12/02
CPVNAM=VDI;SPPNAM=Core
************************************************************
IRUNN = 1 ;LIBREF = 144
************************************************************
Group 1. Run Title
TEXT(No title has been set for this run. )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,100,11,1)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Solved variables list
SOLVE(P1 ,U1 ,V1 )
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
RHO1 = 1.000000E+00
ENUL = 1.000000E-02
CP1 = 1.000000E+00
ENUT = 0.000000E+00
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
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 = 200
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000000E+00)
************************************************************
Group 18. Limits
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
PARSOL = T
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 92 ;IYMON = 6 ;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.000000E+00,0.000000E+00,1.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
> DOM, SIZE, 1.343300E+01, 8.777000E-01, 1.000000E+00
> DOM, MONIT, 1.227015E+01, 4.000000E-01, 5.000000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> GRID, RSET_X_1, 100, 1.000000E+00
> GRID, RSET_Y_1, 9, 1.000000E+00
> GRID, RSET_Y_2, 1, 1.000000E+00
> GRID, RSET_Y_3, 1, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> OBJ, NAME, B1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.343300E+01, 0.000000E+00, 1.000000E+00
> OBJ, CLIPART, cube11
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, NAME, B3
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 8.777000E-01, 1.000000E+00
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, 1.414000E-02, 0.000000E+00, 0.000000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, NAME, B4
> OBJ, POSITION, 1.343300E+01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 8.777000E-01, 1.000000E+00
> OBJ, CLIPART, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, SAME
> OBJ, COEFFICIENT, 1.000000E+03
> OBJ, NAME, B5
> OBJ, POSITION, 0.000000E+00, 6.363000E-01, 0.000000E+00
> OBJ, SIZE, 1.343300E+01, 1.414000E-01, 1.000000E+00
> OBJ, CLIPART, wirexyz
> OBJ, ROTATION24, 1
> OBJ, TYPE, NULL
> OBJ, NAME, B6
> OBJ, POSITION, 0.000000E+00, 7.070000E-01, 0.000000E+00
> OBJ, SIZE, 1.343300E+01, 1.707000E-01, 1.000000E+00
> OBJ, CLIPART, cube14
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198
STOP
distil=t
EX(P1 )= 2.014E-02
EX(U1 )= 1.286E-02
EX(V1 )= 4.530E-05
EX(PRPS)= 1.891E+01