TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2020, Date 16/06/20 CPVNAM=VDI; SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Inform: Inlet boundary layer profile ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,50,1,20) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(144)=DWDX ;NAME(145)=DUDZ NAME(146)=EPKE ;NAME(147)=PRPS NAME(148)=DEN1 ;NAME(149)=EL1 NAME(150) =ENUT * Solved variables list SOLVE(P1,U1,W1) * Stored variables list STORE(ENUT,EL1,DEN1,PRPS,EPKE,DUDZ,DWDX) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) TURMOD(KEMMK) ************************************************************ Group 8. Terms & Devices ************************************************************ 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 PRT(EP)=1.314 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(FOMG)=1. ;FIINIT(PRPS)=-1. 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 This patch is attached to object IN PATCH(INLET, WEST, -1, 0, 0, 0, 0, 0, 1, 1) EGWF = T ************************************************************ Echo save-block settings for Group 13 SAVE13BEGIN **DEFINE CONSTANT PARAMETER** REAL(GA,GD,GE) GA=5.0 GD=0.05 GE=20.0 **DEFINE VIN & MDOT** CHAR(VIN,MDOT) VIN=GA*(ZG)^(1.0/7.0) MDOT=:RHO1:*:VIN: **DEFINE SOURCE OF VARIABLES** (SOURCE OF P1 AT INLET IS :MDOT:) (SOURCE OF U1 AT INLET IS :MDOT:*:VIN:) (SOURCE OF KE AT INLET IS :MDOT:*(GD*:VIN:)^2.0) (SOURCE OF EP AT INLET IS :MDOT:*0.1643*KE^1.5/(GD*GE)) SAVE13END ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 500 SARAH =1. RESFAC =1.0E-03 ************************************************************ Group 16. Terminate Iterations LITER(P1)=50 ;LITER(KE)=10 LITER(EP)=10 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,1. ) RELAX(KE ,LINRLX,0.5 ) RELAX(EP ,LINRLX,0.5 ) KELIN = 3 ************************************************************ Group 18. Limits ************************************************************ Group 19. EARTH Calls To GROUND Station GENK = T PARSOL = F IENUTA = 12 ISG62 = 1 SPEDAT(SET,FACETDAT,NUMOBJ,I,6) SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=2.603 ;EX(U1)=5.488 EX(W1)=0.4156 ;EX(KE)=0.2738 EX(EP)=0.06319 ;EX(FOMG)=0.8053 EX(DWDX)=0.03762 ;EX(DUDZ)=0.4745 EX(EPKE)=0.1384 ;EX(PRPS)=0.965 EX(DEN1)=1.147 ;EX(EL1)=0.8312 EX(ENUT)=0.1644 ************************************************************ Group 21. Print-out of Variables OUTPUT(DWDX,Y,N,N,N,N,N) OUTPUT(DUDZ,Y,N,N,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 33 ;IYMON = 1 ;IZMON = 8 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,48.952942) GVIEW(CENTRE,50.,0.5,10.) > DOM, SIZE, 1.000000E+02, 1.000000E+00, 2.000000E+01 > DOM, MONIT, 6.500000E+01, 5.000000E-01, 5.072407E+00 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > GRID, RSET_X_1, 50, 1.000000E+00 > GRID, RSET_Y_1, 1, 1.000000E+00 > GRID, RSET_Z_1, 20, 1.400000E+00 > OBJ, NAME, IN > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 1.000000E+00, 2.000000E+01 > OBJ, GEOMETRY, box > OBJ, TYPE, USER_DEFINED > OBJ, PATCHES, INLET > OBJ, NAME, OUT > OBJ, POSITION, 1.000000E+02, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 1.000000E+00, 2.000000E+01 > OBJ, GEOMETRY, box > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0. > OBJ, COEFFICIENT, 1000. > OBJ, TURBULENCE, SAME , SAME > OBJ, NAME, SKY > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 2.000000E+01 > OBJ, SIZE, 1.000000E+02, 1.000000E+00, 0.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0. > OBJ, COEFFICIENT, 1000. > OBJ, VELOCITY, 9. ,0. ,0. > OBJ, TURBULENCE, SAME , SAME > OBJ, NAME, GROUND > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E+02, 1.000000E+00, 0.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, PLATE > OBJ, NAME, BLOCK > OBJ, POSITION, 4.000000E+01, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E+01, 1.000000E+00, 5.000000E+00 > OBJ, GEOMETRY, box > OBJ, GRID, N,N,N > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 198,Solid with smooth-wall friction STOP