DISPLAY Library case 647: Tutorial: Turbulence Model KE-CHEM ENDDIS TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 3.6, Date 02/11/04 CPVNAM=VDI;SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings PHOTON USE p;;;; use patgeo gr ou z 1 msg contours of velocity, and velocity vectors con u1 z 1 fi;0.01;vec z 1 msg press RETURN for length scale pause;con off;vec off;red msg contours of turbulence length scale con len1 z 1 fi;0.01;pause;con off;red msg contours of effective viscosity con enut z 1 fi;0.01 enduse DISPLAY This case involves steady, turbulent flow over a backward-facing step, of height h, in a two-dimensional channel of width 3*h. The calculations are started at a distance 4h upstream of the step and terminated at a distance 16h downstream. //////////////////////// wall /////////////////////// ----------------------------------------------------- Pressure Inlet -------> ---------> fixed at zero /________________ | ////////////////| Exit | wall /| Recirculation | /| <---- -----> \ constant /|____________________________________ prescribed mass ////////////////////////////////////// inflow rate wall ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 290 ************************************************************ Group 1. Run Title TEXT(KECHEN model; Re = 45000.00; EGWF= T ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,36,15,1,7.620000E-05) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T * Non-default variable names NAME(147) =VPOR ; NAME(148) =EPKE NAME(149) =LEN1 ; NAME(150) =ENUT * Solved variables list SOLVE(P1 ,U1 ,V1 ) * Stored variables list STORE(ENUT,LEN1,EPKE,VPOR) * Additional solver options SOLUTN(P1 ,Y,Y,Y,N,N,N) TURMOD(KECHEN) ************************************************************ Group 8. Terms & Devices ************************************************************ Group 9. Properties RHO1 = 1.000000E+00 ENUL = 1.100667E-05 CP1 = 1.000000E+00 PRT (KE ) = 7.500000E-01 ;PRT (EP ) = 1.150000E+00 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1 ) = 1.300000E-04 ;FIINIT(U1 ) = 1.300000E+01 FIINIT(KE ) = 7.605000E-01 ;FIINIT(EP ) = 3.177743E+01 FIINIT(VPOR) = 1.000000E+00 ;FIINIT(EPKE) = 1.001000E-10 FIINIT(LEN1) = 1.001000E-10 ;FIINIT(ENUT) = 1.001000E-10 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 = 300 RESFAC = 1.000000E-02 ************************************************************ Group 16. Terminate Iterations LITER (U1 ) = 20 ;LITER (V1 ) = 20 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX, 1.000000E+00) RELAX(U1 ,FALSDT, 1.172308E-02) RELAX(V1 ,FALSDT, 1.172308E-02) RELAX(KE ,LINRLX, 5.000000E-01) RELAX(EP ,LINRLX, 5.000000E-01) KELIN = 3 ************************************************************ Group 18. Limits VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06 VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06 ************************************************************ Group 19. EARTH Calls To GROUND Station USEGRD = T ;USEGRX = T GENK = T ASAP = T PARSOL = T CONWIZ = T IENUTA = 2 SPEDAT(SET,FACETDAT,NUMOBJ,I,7) ************************************************************ Group 20. Preliminary Printout ECHO = T NULLPR = F NDST = 0 DSTTOL = 1.000000E-02 ************************************************************ Group 21. Print-out of Variables ************************************************************ Group 22. Monitor Print-Out IXMON = 9 ;IYMON = 8 ;IZMON = 1 NPRMON = 100 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 IPLTF = 1 ;IPLTL = 100 NPLT = 1 ISWPRF = 1 ;ISWPRL = 100000 ITABL = 1 PATCH (MAP ,CONTUR,1,20,1,15,1,1,1,1) PLOT(MAP ,P1 , 0.000000E+00, 1.000000E+01) PLOT(MAP ,U1 , 0.000000E+00, 1.000000E+01) PLOT(MAP ,V1 , 0.000000E+00, 1.000000E+01) ************************************************************ Group 24. Dumps For Restarts GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00) GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00) > DOM, SIZE, 7.620000E-01, 1.143000E-01, 1.000000E+00 > DOM, MONIT, 1.817660E-01, 2.857500E-02, 5.000000E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, SNAPSIZE, 1.000000E-02 > GRID, RSET_X_1, 4, 1.000000E+00 > GRID, RSET_X_2, 16, 1.300000E+00 > GRID, RSET_X_3, 16, 1.300000E+00 > GRID, RSET_Y_1, 10, 1.000000E+00 > GRID, RSET_Y_2, 5, 1.000000E+00 > OBJ, NAME, X2 > OBJ, POSITION, 3.048000E-01, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 1.143000E-01, 1.000000E+00 > OBJ, GEOMETRY, wirexyz > OBJ, ROTATION24, 1 > OBJ, TYPE, NULL > OBJ, NAME, STEP > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.524000E-01, 3.810000E-02, 1.000000E+00 > OBJ, GEOMETRY, cube14 > OBJ, ROTATION24, 1 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 198,Solid with smooth-wall friction > OBJ, NAME, INLET > OBJ, POSITION, 0.000000E+00, 3.810000E-02, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 7.619999E-02, 1.000000E+00 > OBJ, GEOMETRY, cube3t > OBJ, ROTATION24, 1 > OBJ, TYPE, INLET > OBJ, PRESSURE, 1.000000E+00 > OBJ, VELOCITY, 1.300000E+01, 0.000000E+00, 0.000000E+00 > OBJ, TEMPERATURE, 0.000000E+00 > OBJ, KE_IN, 7.605000E-01 > OBJ, EP_IN, 3.177743E+01 > OBJ, NAME, OUTLET > OBJ, POSITION, 7.619999E-01, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 1.143000E-01, 1.000000E+00 > OBJ, GEOMETRY, cube12t > OBJ, ROTATION24, 1 > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0.000000E+00 > OBJ, TEMPERATURE, SAME > OBJ, COEFFICIENT, 1.000000E+05 > OBJ, TURBULENCE, 0.000000E+00, 0.000000E+00 > OBJ, NAME, WFNN > OBJ, POSITION, 0.000000E+00, 1.143000E-01, 0.000000E+00 > OBJ, SIZE, 7.619999E-01, 0.000000E+00, 1.000000E+00 > OBJ, GEOMETRY, cube11 > OBJ, ROTATION24, 1 > OBJ, TYPE, PLATE > OBJ, NAME, WFNS > OBJ, POSITION, 1.524000E-01, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 6.095999E-01, 0.000000E+00, 1.000000E+00 > OBJ, GEOMETRY, cube11 > OBJ, ROTATION24, 1 > OBJ, TYPE, PLATE STOP