TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2009, Date 26/05/09 CPVNAM=VDI;SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY "FLOW IN A$ SCREW PUMP" Shapemaker has been used here to create the helical coil object and the case presented makes use of it for the situation when the rotation of a helix attached to the shaft pumps the fluid. This case deals with the fluid flow driven by the tangential movement of screw rotor in a an open ended pipe creating the pumping effect. Note that the rotor and coils are treated as solid blockages, and the rotation is imparted by setting the slide velocity to a set angular velocity in rad/s. The chosen angular rotation results in a speed of 1m/s at the surface of the rotor. ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Flow In A Screw Pump ) ************************************************************ 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,20) * 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 ) ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd NAME(150) =PRPS * Solved variables list SOLVE(P1,U1,V1,W1) * Stored variables list STORE(PRPS) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) ************************************************************ Group 8. Terms & Devices ************************************************************ Group 9. Properties PRESS0 =1.0E+05 ;TEMP0 =273. RHO1 =1. ENUL =1.0E-02 CP1 =1005. ENUT =0. DVO1DT =3.41E-03 ************************************************************ 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 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)=50 ************************************************************ 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)=0.2881 ;EX(U1)=0.3682 EX(V1)=0.02747 ;EX(W1)=0.07981 EX(PRPS)=0.429 ************************************************************ Group 21. Print-out of Variables OUTPUT(P1 ,Y,Y,Y,N,Y,Y) OUTPUT(U1 ,Y,Y,Y,N,Y,Y) OUTPUT(V1 ,Y,Y,Y,N,Y,Y) OUTPUT(W1 ,Y,Y,Y,N,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 12 ;IZMON = 17 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.826719,0.01492,0.562418) GVIEW(UP,0.169339,0.946695,-0.274031) > DOM, SIZE, 6.283185E+00, 2.400000E-01, 6.800000E-01 > DOM, MONIT, 1.645621E-01, 1.380001E-01, 5.610000E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > GRID, AUTO, T T T > GRID, RSET_X_1, 8, 1.000000E+00 > GRID, RSET_X_2, 12, 1.000000E+00 > GRID, RSET_Y_1, 10, 1.000000E+00 > GRID, RSET_Y_2, 10, 1.000000E+00 > GRID, RSET_Z_1, 20, 1.000000E+00 > OBJ, NAME, INLET > OBJ, POSITION, 0.000000E+00, 1.200001E-01, 0.000000E+00 > OBJ, SIZE, 6.283000E+00, 1.199999E-01, 0.000000E+00 > OBJ, GEOMETRY, polcubet > OBJ, VISIBLE, NO > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0. > OBJ, TEMPERATURE, SAME > OBJ, COEFFICIENT, 1000. > OBJ, NAME, OUTLET > OBJ, POSITION, 0.000000E+00, 1.200000E-01, 6.800000E-01 > OBJ, SIZE, 6.283000E+00, 1.200000E-01, 0.000000E+00 > OBJ, GEOMETRY, polcubet > OBJ, VISIBLE, NO > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0. > OBJ, TEMPERATURE, SAME > OBJ, COEFFICIENT, 1000. > OBJ, NAME, B5 > OBJ, POSITION, 3.927000E+00, 1.697000E-01, 0.000000E+00 > OBJ, SIZE, 2.400000E-01, 2.400000E-01, 6.800000E-01 > OBJ, GEOMETRY, cylinder > OBJ, GRID, NO > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 198,Solid with smooth-wall friction > OBJ, VELOCITY, 8.33333 ,0. ,0. > OBJ, VEL_TYPE, ANG-VEL > OBJ, NAME, B7 > OBJ, POSITION, 3.927000E+00, 3.400000E-01, 0.000000E+00 > OBJ, SIZE, 4.800000E-01, 4.800000E-01, 4.000000E-01 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, coil > OBJ, GRID, NO > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 198,Solid with smooth-wall friction > OBJ, VELOCITY, 8.33333 ,0. ,0. > OBJ, VEL_TYPE, ANG-VEL > OBJ, NAME, B8 > OBJ, POSITION, 3.927000E+00, 3.400000E-01, 2.800000E-01 > OBJ, SIZE, 4.800000E-01, 4.800000E-01, 4.000000E-01 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, coil > OBJ, GRID, NO > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 198,Solid with smooth-wall friction > OBJ, VELOCITY, 8.33333 ,0. ,0. > OBJ, VEL_TYPE, ANG-VEL > OBJ, NAME, B6 > OBJ, POSITION, 0.000000E+00, 2.400000E-01, 0.000000E+00 > OBJ, SIZE, 2.632993E+00, 0.000000E+00, 6.800000E-01 > OBJ, GEOMETRY, polcu10 > OBJ, VISIBLE, NO > OBJ, TYPE, PLATE STOP