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 627: Tutorial: Flow simulation in a pipe ENDDIS VRV USE * Start of frame VARIABLE Velocity; VEC ON msg Velocity vectors PAUSE * Start of frame VARIABLE Velocity; VEC OFF CON ON CONTOUR BLANK ON CONTOUR AVERAGE ON msg Velocity contours PAUSE * Start of frame VARIABLE Pressure; VEC OFF CON ON msg Pressure contours PAUSE ENDUSE ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Flow simulation in a pipe ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,40,10,55) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(145)=PRPS ;NAME(147)=WDIS NAME(148)=DEN1 ;NAME(149)=EL1 NAME(150) =ENUT * Solved variables list SOLVE(P1,U1,V1,W1) * Stored variables list STORE(ENUT,EL1,DEN1,WDIS,PRPS) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) TURMOD(LVEL) ************************************************************ Group 8. Terms & Devices ************************************************************ Group 9. Properties PRESS0 =1.0E+05 ;TEMP0 =273. * Domain material index is 5 signifying: * Carbon Dioxide at 27 deg C, 1 atm SETPRPS(1, 5) DVO1DT =3.33E-03 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(PRPS)=-1. ;FIINIT(WDIS)=0.03 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.0E-03 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,0.5 ) RELAX(U1 ,LINRLX,0.5 ) RELAX(V1 ,LINRLX,0.5 ) RELAX(W1 ,LINRLX,0.5 ) RELAX(LTLS,LINRLX,1. ) ************************************************************ Group 18. Limits ************************************************************ Group 19. EARTH Calls To GROUND Station PARSOL = F CONWIZ = T ISG62 = 1 SPEDAT(SET,CUTCELL,FRCMN,R,1.0E-02) SPEDAT(SET,CUTCELL,FRCMX,R,0.99) SPEDAT(SET,TOLERANCE,TOLNEWCT,R,1.0E-05) SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=0.1119 ;EX(U1)=0.1364 EX(V1)=6.016E-03 ;EX(W1)=0.1584 EX(PRPS)=0.2611 ;EX(LTLS)=9.305E-04 EX(WDIS)=0.01303 ;EX(DEN1)=0.4692 EX(EL1)=0.01303 ;EX(ENUT)=2.705E-04 ************************************************************ Group 21. Print-out of Variables OUTPUT(WDIS,Y,N,N,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 18 ;IYMON = 5 ;IZMON = 35 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.,-1.,0.) GVIEW(UP,0.,0.,1.) GVIEW(VDIS,1.208609) GVIEW(CENTRE,0.75,0.15,1.) > DOM, SIZE, 1.500000E+00, 3.000000E-01, 2.000000E+00 > DOM, MONIT, 6.750001E-01, 1.350000E-01, 1.385430E+00 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > GRID, RSET_X_1, 5, 1.000000E+00 > GRID, RSET_X_2, 10, 1.000000E+00 > GRID, RSET_X_3, 10, 1.000000E+00 > GRID, RSET_X_4, 10, 1.000000E+00 > GRID, RSET_X_5, 5, 1.000000E+00 > GRID, RSET_Y_1, 10, 1.000000E+00 > GRID, RSET_Z_1, 10, 1.000000E+00 > GRID, RSET_Z_2, 10, 1.000000E+00 > GRID, RSET_Z_3, -15, 1.300000E+00 > GRID, RSET_Z_4, 10, 1.000000E+00 > GRID, RSET_Z_5, 10, 1.000000E+00 > OBJ, NAME, SOLID > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.500000E+00, 3.000000E-01, 2.000000E+00 > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 198,Solid with smooth-wall friction > OBJ, NAME, PIPE-1 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 3.000000E-01, 3.000000E-01, 3.000000E-01 > OBJ, GEOMETRY, cylinder > OBJ, ROTATION24, 9 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, DOMAIN > OBJ, NAME, PIPE-2 > OBJ, POSITION, 1.200000E+00, 0.000000E+00, 1.700000E+00 > OBJ, SIZE, 3.000000E-01, 3.000000E-01, 3.000000E-01 > OBJ, GEOMETRY, cylinder > OBJ, ROTATION24, 9 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, DOMAIN > OBJ, NAME, PIPE-3 > OBJ, POSITION, 6.000000E-01, 0.000000E+00, 6.000000E-01 > OBJ, SIZE, 3.000000E-01, 3.000000E-01, 8.000000E-01 > OBJ, GEOMETRY, cylinder > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, DOMAIN > OBJ, NAME, BPIPE-1 > OBJ, POSITION, 3.000000E-01, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 6.000000E-01, 3.000000E-01, 6.000000E-01 > OBJ, GEOMETRY, pipebend > OBJ, ROTATION24, 9 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, DOMAIN > OBJ, NAME, BPIPE-2 > OBJ, POSITION, 6.000000E-01, 0.000000E+00, 1.400000E+00 > OBJ, SIZE, 6.000000E-01, 3.000000E-01, 6.000000E-01 > OBJ, GEOMETRY, pipebend > OBJ, ROTATION24, 3 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, DOMAIN > OBJ, NAME, INLET > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 3.000000E-01, 3.000000E-01 > OBJ, GEOMETRY, cylinder > OBJ, ROTATION24, 9 > OBJ, TYPE, INLET > OBJ, PRESSURE, 0. > OBJ, VELOCITY, 1. ,0. ,0. > OBJ, NAME, OUTLET > OBJ, POSITION, 1.500000E+00, 0.000000E+00, 1.700000E+00 > OBJ, SIZE, 0.000000E+00, 3.000000E-01, 3.000000E-01 > OBJ, GEOMETRY, cylinder > OBJ, ROTATION24, 9 > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0. > OBJ, COEFFICIENT, 1000. STOP