TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2023, Date 05/06/24 CPVNAM=VDI; SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY This run illustrates the activation of a power-law laminar viscosity formula via the PROPS file ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 576 ************************************************************ Group 1. Run Title TEXT(2D Channel Flow;Non Newtonian Fluid ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,1,10,15) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(147)=GEN1 ;NAME(148)=ENUL NAME(149)=PRPS ;NAME(150)=TEM1 * Solved variables list SOLVE(P1,V1,W1,TEM1) * Stored variables list STORE(PRPS,ENUL,GEN1) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(TEM1,Y,Y,Y,N,N,Y) ************************************************************ Group 8. Terms & Devices NEWENL = T ************************************************************ Group 9. Properties * Domain material index is 90 signifying: * NOTSET SETPRPS(1, 90,Q1) ENUT =0. PRNDTL(TEM1)=-0.1428 * List of user-defined materials to be read by EARTH MATFLG=T;IMAT=1 * Name *Ind. Dens. Viscos. Spec.heat Conduct. Expans. Compr. *90 1000.0 GRND4 1000 .1428 0 * constants for GRND option no 1 0.001 2 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(GEN1)=1.001E-10 ;FIINIT(ENUL)=1.001E-10 FIINIT(PRPS)=-1. ;FIINIT(TEM1)=1.001E-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 = 500 RESFAC =1.0E-05 ************************************************************ Group 16. Terminate Iterations ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,0.7 ) RELAX(V1 ,FALSDT,0.010667 ) RELAX(W1 ,FALSDT,0.010667 ) RELAX(TEM1,LINRLX,0.5 ) ************************************************************ Group 18. Limits ************************************************************ Group 19. EARTH Calls To GROUND Station GENK = T PARSOL = F ISG62 = 1 CSG10 ='q1' SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=1.568E+04 ;EX(V1)=0.5726 EX(W1)=4.987 ;EX(GEN1)=451.5 EX(ENUL)=0.06668 ;EX(PRPS)=17.67 EX(TEM1)=8.995 ************************************************************ Group 21. Print-out of Variables ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 5 ;IZMON = 13 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,-1.,0.,0.) GVIEW(UP,0.,1.,0.) GVIEW(VDIS,2.097562) GVIEW(CENTRE,0.5,0.5,2.) > DOM, SIZE, 1.000000E+00, 1.000000E+00, 4.000000E+00 > DOM, MONIT, 5.000000E-01, 4.683772E-01, 2.520000E+00 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > GRID, RSET_X_1, 1, 1.000000E+00 > GRID, RSET_Y_1, -5, 1.500000E+00 > GRID, RSET_Y_2, -5, 1.500000E+00 > GRID, RSET_Z_1, 5,-1.200000E+00 > GRID, RSET_Z_2, -5, 1.500000E+00 > GRID, RSET_Z_3, 5, 2.000000E+00 > DOM, T_AMBIENT, 0.000000E+00 > OBJ, NAME, BLK > OBJ, POSITION, 0.000000E+00, 5.000000E-01, 1.000000E+00 > OBJ, SIZE, 1.000000E+00, 5.000000E-01, 1.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 101,BAKELITE at 20 deg c > OBJ, NAME, IN > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E+00, 1.000000E+00, 0.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, INLET > OBJ, DENSITY, 1000. > OBJ, VELOCITY, 0. ,0. ,5. > OBJ, TEMPERATURE, 9. > OBJ, NAME, OUTLET > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 4.000000E+00 > OBJ, SIZE, 1.000000E+00, 1.000000E+00, 0.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0. > OBJ, TEMPERATURE, SAME > OBJ, COEFFICIENT, 1000. > OBJ, NAME, NORTH > OBJ, POSITION, 0.000000E+00, 1.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E+00, 0.000000E+00, 1.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, PLATE > OBJ, SURF_TEMP, 0. ,0. > OBJ, NAME, NORTH1 > OBJ, POSITION, 0.000000E+00, 1.000000E+00, 2.000000E+00 > OBJ, SIZE, 1.000000E+00, 0.000000E+00, 2.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, PLATE > OBJ, SURF_TEMP, 0. ,0. ************************************************************ STOP