TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 3.4, Date 10/08/01 CPVNAM=VDI;SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY Library case 30 loaded, to activate EXPERT xprtbegin ... start of the inputs to expert via Q1 expsol f .... expert is not used in the solver expdtf t .... expert is used to adjust false time step indtf 2 ... adjust the false time step of all velocities inres 3 ... monitor the residuals of variable 3 ifrsts 20 ... first adjustment sweep ifrequ 10 ... adjustments will be made every ifrequ sweeps facdec 0.5 .. decreases will be by a factor of facdec facinc 2.0 .. increases will be by a factor of facinc exprin t .... print adjusted false time steps in RESULT file expend ..... end of the inputs to expert via Q1 xprtend ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Case13<-09 Sutherland viscosity m=1/2 ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,1,100,200) * Set overall domain extent: * xulast yvlast zwlast name XSI= 1.000000E-01; YSI= 5.000000E-03; ZSI= 1.200000E-02 RSET(D,CHAM ) * Cylindrical-polar grid CARTES=F ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T * Non-default variable names NAME(147) =MACH ; NAME(148) =DEN1 NAME(149) =TEM1 ; NAME(150) =PRPS * Solved variables list SOLVE(P1 ,V1 ,W1 ,TEM1) * Stored variables list STORE(PRPS,DEN1,MACH) * Additional solver options SOLUTN(P1 ,Y,Y,Y,N,N,Y) SOLUTN(V1 ,Y,Y,Y,N,N,Y) SOLUTN(W1 ,Y,Y,Y,N,N,Y) SOLUTN(TEM1,Y,Y,Y,N,N,Y) ************************************************************ Group 8. Terms & Devices TERMS (TEM1,Y,Y,Y,Y,N,Y) NEWRH1 = T NEWENL = T DENPCO = T ************************************************************ Group 9. Properties SETPRPS(1, 19) RHO1 = GRND5 RHO1A = 0.000000E+00 ;RHO1B = 3.369000E-03 RHO1C = 7.143000E-01 CP1 = 1.039000E+03 ENUL = GRND6 ;ENUT = 0.000000E+00 ENULA = 1.390000E-06 ;ENULB = 1.020000E+02 ENULC = 0.000000E+00 DRH1DP = GRND5 DVO1DT = 3.330000E-03 PRNDTL(TEM1) = -2.610000E-02 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1 ) = 3.000000E+03 ;FIINIT(TEM1) = 2.000000E+02 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 = 2000 RESFAC = 1.000000E-03 ************************************************************ Group 16. Terminate Iterations ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX, 1.000000E+00) RELAX(V1 ,FALSDT, 1.000000E-09) RELAX(W1 ,FALSDT, 1.000000E-09) RELAX(TEM1,FALSDT, 1.000000E-04) EXPERT = T ************************************************************ Group 18. Limits VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06 VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06 ************************************************************ Group 19. EARTH Calls To GROUND Station USEGRD = T ;USEGRX = T ASAP = T RADIA = 0.000000E+00 ;RADIB = 0.000000E+00 EMISS = 0.000000E+00 ************************************************************ Group 20. Preliminary Printout ECHO = T ************************************************************ Group 21. Print-out of Variables ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 1 ;IZMON = 168 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 NOWIPE = T GVIEW(P,-1.000000e+000,0.000000e+000,0.000000e+000) GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00) > DOM, SIZE, 1.000000E-01, 5.000000E-03, 1.200000E-02 > DOM, MONIT, 5.000000E-02, 3.216200E-06, 1.004220E-02 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, SNAPSIZE, 1.000000E-02 > GRID, RSET_Y_1, 50, 1.300000E+00 > GRID, RSET_Y_2, 50, 1.300000E+00 > DOM, RELAX, 5.000000E-01 > OBJ1, NAME, NOZZLE > OBJ1, POSITION, 0.000000E+00, 0.000000E+00, 2.000000E-03 > OBJ1, SIZE, 1.000000E-01, 5.000000E-03, 1.000000E-02 > OBJ1, CLIPART, nozzleshape > OBJ1, ROTATION24, 1 > OBJ1, TYPE, BLOCKAGE > OBJ1, MATERIAL, 198 > OBJ1, ADIABATIC_W, 0.000000E+00, 0.000000E+00 > OBJ1, ADIABATIC_E, 0.000000E+00, 0.000000E+00 > OBJ1, ADIABATIC_S, 0.000000E+00, 0.000000E+00 > OBJ1, ADIABATIC_N, 0.000000E+00, 0.000000E+00 > OBJ1, ADIABATIC_L, 0.000000E+00, 0.000000E+00 > OBJ1, ADIABATIC_H, 0.000000E+00, 0.000000E+00 > OBJ2, NAME, IN > OBJ2, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ2, SIZE, 1.000000E-01, 4.900000E-03, 0.000000E+00 > OBJ2, CLIPART, polcu5t > OBJ2, ROTATION24, 1 > OBJ2, TYPE, INLET > OBJ2, DENSITY, 6.000000E-02 > OBJ2, VOLUFLOW, 2.400000E-05 > OBJ2, TEMPERATURE, 3.000000E+02 > OBJ3, NAME, OUT > OBJ3, POSITION, 0.000000E+00, 0.000000E+00, 1.200000E-02 > OBJ3, SIZE, 1.000000E-01, 1.040000E-03, 0.000000E+00 > OBJ3, CLIPART, polcubet > OBJ3, ROTATION24, 1 > OBJ3, TYPE, OUTLET > OBJ3, PRESSURE, 5.000000E+02 > OBJ3, TEMPERATURE, -1.026000E+04 > OBJ3, COEFFICIENT, 1.000000E+00 > OBJ3, VELOCITY, -1.026000E+04,-1.026000E+04,-1.026000E+04 > OBJ4, NAME, PIPE > OBJ4, POSITION, 0.000000E+00, 4.900000E-03, 0.000000E+00 > OBJ4, SIZE, 1.000000E-01, 1.000000E-04, 2.000000E-03 > OBJ4, CLIPART, polcu8 > OBJ4, ROTATION24, 1 > OBJ4, TYPE, BLOCKAGE > OBJ4, MATERIAL, 198 > OBJ4, ADIABATIC_W, 0.000000E+00, 0.000000E+00 > OBJ4, ADIABATIC_E, 0.000000E+00, 0.000000E+00 > OBJ4, ADIABATIC_S, 0.000000E+00, 0.000000E+00 > OBJ4, ADIABATIC_N, 0.000000E+00, 0.000000E+00 > OBJ4, ADIABATIC_L, 0.000000E+00, 0.000000E+00 > OBJ4, ADIABATIC_H, 0.000000E+00, 0.000000E+00 STOP LG(1)=T