TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2009, Date 26/05/09
CPVNAM=VDI; SPPNAM = Core
display
This is one of six test cases which check the operation
of the INLET object in polar coordinates. They test all
combinations of flow direction (in or out), and IURVAL (-1,+1,0)
In each case there are eight inlet objects, four at the
low-X end and 4 at the high-X end. Each uses one of the
four inflow options (grid vel, volume flow, mass flow and
cartesian components). The IURVAL setting should only
influence the inlets using grid velocity.
Area of inlet = 0.5
Rinner = 2. Yg = 0.125, 0.375. Radii = 2.125, 2.375
density = 1.189
Inlet grid-directed velocity U=10
IURVAL= 0 - U is velocity
mass source = 0.5*1.189*10 = 5.945
U1 source = (0.5*1.189*10)*10 = 59.45
IURVAL=-1 - U is angular velocity
mass source = 0.25*1.189*10*(2.125 + 2.375) = 13.376
U1 source = 0.25*1.189*10*10*(2.125*2.125 + 2.375*2.375) = 301.89
IURVAL=+1 - U is angular momentum
mass source = 0.25*1.189*10*(1/2.125 + 1/2.375) = 2.65
U1 source = 0.25*1.189*10*10*(1/(2.125*2.125) + 1/(2.375*2.375)) = 11.85
Inlet volume flow = 10 - unaffected by IURVAL
Mass source = 1.189*10 = 23.78
U1 source = (1.189*10)*(10/0.5) = 237.8
Inlet mass flow = 10 - unaffected by IURVAL
mass source = 10
U1 source = 10*(10/1.189/0.5) = 168.2
Cartesian inlet velocities (10,0,0) - unaffected by IURVAL
mass source = 0.5*1.189*10 = 5.945
U1 source = (0.5*1.189*10)*10 = 59.45
Expected results are included.
enddis
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(Angular Momentum IURVAL=+1 Inflow )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,10,2,4)
* Cylindrical-polar grid
CARTES=F;RINNER=2.000000E+00
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(145)=PRPS ;NAME(146)=HPOR
NAME(147)=EPKE ;NAME(148)=DEN1
NAME(149)=EL1 ;NAME(150)=ENUT
* Solved variables list
SOLVE(P1,U1,V1,W1)
* Stored variables list
STORE(ENUT,EL1,DEN1,EPKE,HPOR,PRPS)
* Additional solver options
SOLUTN(P1,Y,Y,Y,N,N,Y)
TURMOD(KECHEN)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
PRESS0 =1.01325E+05 ;TEMP0 =273.
* Domain material index is 0 signifying:
* Air at 20 deg C, 1 atm, treated as incompressible
SETPRPS(1, 0)
DVO1DT =3.41E-03
PRT(KE)=0.75 ;PRT(EP)=1.15
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1)=0. ;FIINIT(PRPS)=-1.
FIINIT(HPOR)=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
IURVAL = 1
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 100
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
LITER(P1)=200
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
RELAX(KE ,LINRLX,0.5 )
RELAX(EP ,LINRLX,0.5 )
KELIN = 3
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
GENK = T
CONWIZ = T
IENUTA = 2
SPEDAT(SET,KECONST,C1E,R,1.44)
SPEDAT(SET,KECONST,C2E,R,1.92)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=168.199997 ;EX(U1)=9.959
EX(V1)=1.7 ;EX(W1)=0.
EX(KE)=0.7129 ;EX(EP)=5.505
EX(PRPS)=0.8 ;EX(HPOR)=0.2
EX(EPKE)=4.874 ;EX(DEN1)=0.9512
EX(EL1)=0.02463 ;EX(ENUT)=0.01213
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 4 ;IYMON = 1 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NXPRIN = 1
NYPRIN = 1
NZPRIN = 1
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,0.416368,0.461451,0.783391)
GVIEW(UP,0.190372,0.798283,-0.571405)
> DOM, SIZE, 1.570800E+00, 5.000000E-01, 4.000000E+00
> DOM, MONIT, 5.100001E-01, 1.250000E-01, 5.000000E-01
> 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, 1, 1.000000E+00
> GRID, RSET_X_2, 3, 1.000000E+00
> GRID, RSET_X_3, 1, 1.000000E+00
> GRID, RSET_X_4, 4, 1.000000E+00
> GRID, RSET_X_5, 1, 1.000000E+00
> GRID, RSET_Y_1, 2, 1.000000E+00
> GRID, RSET_Z_1, 1, 1.000000E+00
> GRID, RSET_Z_2, 1, 1.000000E+00
> GRID, RSET_Z_3, 1, 1.000000E+00
> GRID, RSET_Z_4, 1, 1.000000E+00
> DOM, INI_AMB, YES
> DOM, INI_BUOY, YES
> OBJ, NAME, BLOCK1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E-02, TO_END, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu8
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, BLOCK2
> OBJ, POSITION, AT_END, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E-02, TO_END, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu8
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, RAD/S
> OBJ, POSITION, 1.000000E-02, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, TO_END, 1.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu5t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, VELOCITY, 10. ,0. ,0.
> OBJ, TURB-INTENS, 5.
> OBJ, OBJECT-SIDE, HIGH
> OBJ, NAME, VOLFLO
> OBJ, POSITION, 1.000000E-02, 0.000000E+00, 1.000000E+00
> OBJ, SIZE, 0.000000E+00, TO_END, 1.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu5t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, VOLUFLOW, 10.
> OBJ, TURB-INTENS, 5.
> OBJ, OBJECT-SIDE, HIGH
> OBJ, NAME, MASFLO
> OBJ, POSITION, 1.000000E-02, 0.000000E+00, 2.000000E+00
> OBJ, SIZE, 0.000000E+00, TO_END, 1.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu5t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, MASSFLOW, 10.
> OBJ, TURB-INTENS, 5.
> OBJ, OBJECT-SIDE, HIGH
> OBJ, NAME, CARTVEL
> OBJ, POSITION, 1.000000E-02, 0.000000E+00, 3.000000E+00
> OBJ, SIZE, 0.000000E+00, TO_END, 1.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu5t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, CART-VEL, 10. ,0. ,0.
> OBJ, TURB-INTENS, 5.
> OBJ, OBJECT-SIDE, HIGH
> OBJ, NAME, PL1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ, SIZE, TO_END, TO_END, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu10
> OBJ, TYPE, PLATE
> OBJ, POROSITY, 0.
> OBJ, SIDE, BOTH
> OBJ, WALL_CO_L, FRICTIONLESS
> OBJ, WALL_CO_H, FRICTIONLESS
> OBJ, NAME, PL2
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 2.000000E+00
> OBJ, SIZE, TO_END, TO_END, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu10
> OBJ, TYPE, PLATE
> OBJ, POROSITY, 0.
> OBJ, SIDE, BOTH
> OBJ, WALL_CO_L, FRICTIONLESS
> OBJ, WALL_CO_H, FRICTIONLESS
> OBJ, NAME, PL3
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 3.000000E+00
> OBJ, SIZE, TO_END, TO_END, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu10
> OBJ, TYPE, PLATE
> OBJ, POROSITY, 0.
> OBJ, SIDE, BOTH
> OBJ, WALL_CO_L, FRICTIONLESS
> OBJ, WALL_CO_H, FRICTIONLESS
> OBJ, NAME, RAD/S-2
> OBJ, POSITION, 1.560800E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, TO_END, 1.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu5t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, VELOCITY, -10. ,0. ,0.
> OBJ, TURB-INTENS, 5.
> OBJ, NAME, VOLFLO-2
> OBJ, POSITION, 1.560800E+00, 0.000000E+00, 1.000000E+00
> OBJ, SIZE, 0.000000E+00, TO_END, 1.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu5t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, VOLUFLOW, 10.
> OBJ, TURB-INTENS, 5.
> OBJ, NAME, MASFLO-2
> OBJ, POSITION, 1.560800E+00, 0.000000E+00, 2.000000E+00
> OBJ, SIZE, 0.000000E+00, TO_END, 1.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu5t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, MASSFLOW, 10.
> OBJ, TURB-INTENS, 5.
> OBJ, NAME, CARTVEL2
> OBJ, POSITION, 1.560800E+00, 0.000000E+00, 3.000000E+00
> OBJ, SIZE, 0.000000E+00, TO_END, 1.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu5t
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, CART-VEL, 0. ,10. ,0.
> OBJ, TURB-INTENS, 5.
> OBJ, NAME, OUT
> OBJ, POSITION, 6.100000E-01, AT_END, 0.000000E+00
> OBJ, SIZE, 1.570800E-01, 0.000000E+00, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcubet
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, COEFFICIENT, 1000.
> OBJ, TURBULENCE, SAME , SAME
STOP