TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2020, Date 16/06/20 CPVNAM=VDI; SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY Library case 638: Tutorial: Swirling flow through an orifice plate ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Swirling flow through an orifice plate ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,20,31,43) * 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 ) XCYIZ( 31, T ) ; XCYIZ( 32, T ) ; XCYIZ( 33, T ) XCYIZ( 34, T ) ; XCYIZ( 35, T ) ; XCYIZ( 36, T ) XCYIZ( 37, T ) ; XCYIZ( 38, T ) ; XCYIZ( 39, T ) XCYIZ( 40, T ) ; XCYIZ( 41, T ) ; XCYIZ( 42, T ) XCYIZ( 43, T ) ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(145)=HPOR ;NAME(146)=PRPS 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,PRPS,HPOR) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) TURMOD(KECHEN) ************************************************************ Group 8. Terms & Devices ************************************************************ Group 9. Properties PRESS0 =1.0E+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(U1)=20. ;FIINIT(W1)=1. FIINIT(HPOR)=1. ;FIINIT(PRPS)=-1. No PATCHes used for this Group INIADD = F IURINI = -1 ************************************************************ 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 XCYCLE = T EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 500 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 PARSOL = F CONWIZ = T IENUTA = 2 ISG62 = 1 SPEDAT(SET,FACETDAT,NUMOBJ,I,5) SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=10.94 ;EX(U1)=2.276 EX(V1)=0.2134 ;EX(W1)=0.9296 EX(KE)=0.4135 ;EX(EP)=7.717 EX(HPOR)=0.8192 ;EX(PRPS)=0.832 EX(EPKE)=13.8 ;EX(DEN1)=0.9892 EX(EL1)=6.73E-03 ;EX(ENUT)=2.74E-03 ************************************************************ Group 21. Print-out of Variables ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 23 ;IZMON = 9 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.814286,0.288411,0.503743) GVIEW(UP,0.162721,0.946451,-0.278843) GVIEW(VDIS,0.47571) GVIEW(CENTRE,0.,0.,0.5) > DOM, SIZE, 6.283185E+00, 2.000000E-01, 1.000000E+00 > DOM, MONIT, 1.570796E-01, 1.500000E-01, 1.627803E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > GRID, AUTO, F T T > GRID, RSET_X_1, 20, 1.000000E+00 > GRID, RSET_Y_1, 14,-1.100000E+00,G > GRID, RSET_Y_2, -17, 1.100000E+00,G > GRID, RSET_Z_1, -16, 1.100000E+00,G > GRID, RSET_Z_2, -12, 1.100000E+00,G > GRID, RSET_Z_3, -15, 1.100000E+00,G > OBJ, NAME, INLET > OBJ, POSITION, 0.000000E+00, 1.000000E-01, 0.000000E+00 > OBJ, SIZE, TO_END, 1.000000E-01, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, poldef > OBJ, TYPE, INLET > OBJ, PRESSURE, 0. > OBJ, VELOCITY, 20. ,0. ,1. > OBJ, TURB-INTENS, 5. > OBJ, NAME, OUTLET > OBJ, POSITION, 0.000000E+00, 0.000000E+00, AT_END > OBJ, SIZE, TO_END, 2.000000E-01, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, poldef > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0. > OBJ, COEFFICIENT, 1000. > OBJ, TURBULENCE, 1.0E-03 ,9.0E-03 > OBJ, NAME, CENTRE > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, 1.000000E-01, 3.000000E-01 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, poldef > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 198,Solid with smooth-wall friction > OBJ, NAME, ORIFICE > OBJ, POSITION, 0.000000E+00, 1.000000E-01, 6.000000E-01 > OBJ, SIZE, TO_END, TO_END, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, poldef > OBJ, TYPE, PLATE > OBJ, POROSITY, 0. > OBJ, SIDE, BOTH > OBJ, NAME, WALL > OBJ, POSITION, 0.000000E+00, AT_END, 0.000000E+00 > OBJ, SIZE, TO_END, 0.000000E+00, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, poldef > OBJ, VISIBLE, NO > OBJ, TYPE, PLATE STOP