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 621: Tutorial: Simulation of a Labyrinth Flow ENDDIS VRV USE * Start of frame SLICE Y SLICE OUTLINE ON VARIABLE Velocity; VEC ON msg Velocity vectors on Y plane PAUSE * Start of frame * Setting display swiches VARIABLE Pressure VEC OFF CON ON CONTOUR BLANK ON CONTOUR AVERAGE ON msg Pressure contours on Y plane PAUSE * Start of frame VARIABLE Pressure;CON OFF * Generating streamlines STREAM CLEAR STREAM COLOUR Total STREAM TIME -1.000000E+10 1.000000E+10 STREAM DIRECTION Both STREAM ORIGIN Probe STREAM MODE Arrow STREAM 0.000000E+00 0.000000E+00 0.000000E+00 STREAM 1.428571E-01 5.000000E-01 7.142857E-02 STREAM 2.857143E-01 5.000000E-01 1.428571E-01 STREAM 4.285715E-01 5.000000E-01 2.142857E-01 STREAM 5.714286E-01 5.000000E-01 2.857143E-01 STREAM 7.142857E-01 5.000000E-01 3.571429E-01 STREAM 8.571429E-01 5.000000E-01 4.285714E-01 STREAM 1.000000E+00 5.000000E-01 5.000000E-01 STREAM 0.000000E+00 0.000000E+00 0.000000E+00 STREAM 0.000000E+00 0.000000E+00 0.000000E+00 STREAM 1.428572E+00 5.000000E-01 7.142858E-01 STREAM 1.571429E+00 5.000000E-01 7.857144E-01 STREAM 1.714286E+00 5.000000E-01 8.571430E-01 STREAM 1.857143E+00 5.000000E-01 9.285716E-01 STREAM 0.000000E+00 0.000000E+00 0.000000E+00 mgs Stream lines PAUSE ENDUSE ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Simulation of a Labyrinth Flow ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,39,1,60) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(146)=EPOR ;NAME(147)=PRPS NAME(148)=ENUT ;NAME(149)=EPKE NAME(150) =TEM1 * Solved variables list SOLVE(P1,U1,W1,TEM1) * Stored variables list STORE(EPKE,ENUT,PRPS,EPOR) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(TEM1,Y,Y,Y,N,N,Y) TURMOD(KECHEN) ************************************************************ Group 8. Terms & Devices TERMS (TEM1,Y,Y,Y,Y,Y,Y) ************************************************************ 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 PRNDTL(TEM1)=-0.0258 PRT(KE)=0.75 ;PRT(EP)=1.15 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1)=0. ;FIINIT(EPOR)=1. FIINIT(PRPS)=-1. ;FIINIT(TEM1)=20. 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 = 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 VARMAX(TEM1)=3000. ;VARMIN(TEM1)=-204.75 ************************************************************ Group 19. EARTH Calls To GROUND Station GENK = T PARSOL = F CONWIZ = T IENUTA = 2 ISG62 = 1 SPEDAT(SET,FACETDAT,NUMOBJ,I,9) SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=1.11 ;EX(U1)=0.283 EX(W1)=0.307 ;EX(KE)=0.01724 EX(EP)=0.03961 ;EX(EPOR)=0.97 EX(PRPS)=26.75 ;EX(ENUT)=1.027E-03 EX(EPKE)=1.3 ;EX(TEM1)=21.48 ************************************************************ Group 21. Print-out of Variables ************************************************************ Group 22. Monitor Print-Out IXMON = 34 ;IYMON = 1 ;IZMON = 18 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.175755) GVIEW(CENTRE,1.,0.5,0.5) > DOM, SIZE, 2.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, MONIT, 1.731310E+00, 5.000000E-01, 2.704330E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > GRID, AUTO, T T T > GRID, RSET_X_1, -11, 1.100000E+00,G > GRID, RSET_X_2, -8, 1.100000E+00,G > GRID, RSET_X_3, -8, 1.100000E+00,G > GRID, RSET_X_4, -12, 1.100000E+00,G > GRID, RSET_Y_1, 1, 1.000000E+00 > GRID, RSET_Z_1, 4,-1.100000E+00,G > GRID, RSET_Z_2, -16, 1.100000E+00,G > GRID, RSET_Z_3, -10, 1.100000E+00,G > GRID, RSET_Z_4, -10, 1.100000E+00,G > GRID, RSET_Z_5, -16, 1.100000E+00,G > GRID, RSET_Z_6, 4, 1.100000E+00,G > DOM, INI_AMB, YES > DOM, INI_BUOY, YES > OBJ, NAME, H-BLOCK > OBJ, POSITION, 0.000000E+00, 0.000000E+00, AT_END > OBJ, SIZE, TO_END, TO_END, 5.000000E-02 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 100,ALUMINIUM at 27 deg c > OBJ, NAME, L_BLOCK > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, TO_END, 5.000000E-02 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 100,ALUMINIUM at 27 deg c > OBJ, HEAT_FLUX, 0. ,100. > OBJ, NAME, WALL-E > OBJ, POSITION, 2.000000E+00, 0.000000E+00, 5.000000E-01 > OBJ, SIZE, 0.000000E+00, TO_END, 5.000000E-01 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, PLATE > OBJ, NAME, WALL-W > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, TO_END, 5.000000E-01 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, PLATE > OBJ, NAME, IN-PLATE > OBJ, POSITION, 5.000000E-01, 0.000000E+00, 3.000000E-01 > OBJ, SIZE, 0.000000E+00, 1.000000E+00, 6.500000E-01 > OBJ, GEOMETRY, box > OBJ, TYPE, PLATE > OBJ, POROSITY, 0. > OBJ, SIDE, BOTH > OBJ, NAME, IN-BLOCK > OBJ, POSITION, 1.000000E+00, 0.000000E+00, 5.000000E-02 > OBJ, SIZE, 4.000000E-01, TO_END, 6.500000E-01 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 103,COPPER at 27 deg C > OBJ, NAME, INLET > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 5.000000E-01 > OBJ, SIZE, 0.000000E+00, TO_END, 4.500000E-01 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, INLET > OBJ, PRESSURE, 0. > OBJ, VELOCITY, 0.5 ,0. ,0. > OBJ, TEMPERATURE, T_AMBIENT > OBJ, TURB-INTENS, 5. > OBJ, NAME, OUTLET > OBJ, POSITION, 2.000000E+00, 0.000000E+00, 5.000000E-02 > OBJ, SIZE, 0.000000E+00, TO_END, 4.500000E-01 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0. > OBJ, TEMPERATURE, T_AMBIENT > OBJ, COEFFICIENT, 1000. > OBJ, TURBULENCE, SAME , SAME STOP