TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2020, Date 16/06/20 CPVNAM=VDI; SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings PHOTON USE p msg pressure contours and velocity vectors msg (for case 276, pressure represents velocity potential) con p1 x 1 fi;0.05; con p1 y 1 fi;0.05 set ref vec 40.0; vec x 1; vec y 1 msg msg Press e to END ENDUSE DISPLAY This case is based on Library case 274, which deals with the flow pas a moving van. An impression of the geometry can be gained by examining the following sketch. | ______________ |______ |**************| |******| ^ __|**************| ^ |******| | |*****************| | |******| y| |*****************| y| |******| | ~**~~~~~~~~~~~~~**~ | |~~~~** ** ** | ** ------------------------------------------- z----> x---> Now, In-Form is used to make the inlet velocity a linear function of X and Y. VELin = 14*(xg/xulast+yg/yvlast) ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 274 ************************************************************ Group 1. Run Title TEXT(3D Nav.Sto.Flow Around A Moving Van ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,8,12,20) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd NAME(150) =PRPS * Solved variables list SOLVE(P1,U1,V1,W1) * Stored variables list STORE(PRPS) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,N) SOLUTN(U1,Y,Y,N,Y,N,Y) SOLUTN(V1,Y,Y,N,Y,N,Y) SOLUTN(W1,Y,Y,N,Y,N,Y) ************************************************************ Group 8. Terms & Devices ************************************************************ Group 9. Properties RHO1 =1. ENUL =1.0E-05 CP1 =1. ENUT =1.0E-03 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1)=0. ;FIINIT(W1)=14. 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 No PATCHes used for this Group EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 200 SARAH =1. RESFAC =1.0E-05 ************************************************************ Group 16. Terminate Iterations ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,1. ) RELAX(U1 ,FALSDT,1.0E-04 ) RELAX(V1 ,FALSDT,1.0E-04 ) RELAX(W1 ,FALSDT,1.0E-04 ) ************************************************************ Group 18. Limits ************************************************************ Group 19. EARTH Calls To GROUND Station PARSOL = F ISG62 = 1 SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=7.108 ;EX(U1)=0.2112 EX(V1)=0.1775 ;EX(W1)=13.41 EX(PRPS)=0.9646 ************************************************************ Group 21. Print-out of Variables ************************************************************ Group 22. Monitor Print-Out IXMON = 5 ;IYMON = 5 ;IZMON = 5 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 NXPRIN = 4 NYPRIN = 2 NZPRIN = 2 YZPR = T ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts GVIEW(P,0.641234,0.496856,0.584767) GVIEW(UP,-0.350212,0.867563,-0.353109) GVIEW(VDIS,10.306038) GVIEW(CENTRE,2.5,3.,10.) > DOM, SIZE, 5.000000E+00, 6.000000E+00, 2.000000E+01 > DOM, MONIT, 2.550000E+00, 2.250000E+00, 4.500000E+00 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02 > GRID, RSET_X_1, 2, 1.000000E+00 > GRID, RSET_X_2, 1, 1.000000E+00 > GRID, RSET_X_3, 5, 1.000000E+00 > GRID, RSET_Y_1, 1, 1.000000E+00 > GRID, RSET_Y_2, 2, 1.000000E+00 > GRID, RSET_Y_3, 3, 1.000000E+00 > GRID, RSET_Y_4, 6, 1.000000E+00 > GRID, RSET_Z_1, 5, 1.000000E+00 > GRID, RSET_Z_2, 1, 1.000000E+00 > GRID, RSET_Z_3, 3, 1.000000E+00 > GRID, RSET_Z_4, 1, 1.000000E+00 > GRID, RSET_Z_5, 10, 1.000000E+00 > DOM, T_AMBIENT, 0.000000E+00 > OBJ, NAME, BODY1 > OBJ, POSITION, 0.000000E+00, 5.000000E-01, 5.000000E+00 > OBJ, SIZE, 1.500000E+00, 1.000000E+00, 5.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 198,Solid with smooth-wall friction > OBJ, NAME, BODY2 > OBJ, POSITION, 0.000000E+00, 1.500000E+00, 6.000000E+00 > OBJ, SIZE, 1.500000E+00, 1.500000E+00, 4.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 198,Solid with smooth-wall friction > OBJ, NAME, FWHL > OBJ, POSITION, 1.000000E+00, 0.000000E+00, 5.000000E+00 > OBJ, SIZE, 5.000000E-01, 5.000000E-01, 1.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 199,Solid allowing fluid-slip at walls > OBJ, NAME, RWHL > OBJ, POSITION, 1.000000E+00, 0.000000E+00, 9.000000E+00 > OBJ, SIZE, 5.000000E-01, 5.000000E-01, 1.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 199,Solid allowing fluid-slip at walls > OBJ, NAME, UPSTR > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 5.000000E+00, 6.000000E+00, 0.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, INLET > OBJ, INFSRC_P1, den1*14*(xg/xulast+yg/yvlast) with area > OBJ, INFSRC_W1, 14*(xg/xulast+yg/yvlast) with onlyms > OBJ, PRESSURE, P_AMBIENT > OBJ, VELOCITY, 0. ,0. ,14. > OBJ, NAME, DWSTR > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 2.000000E+01 > OBJ, SIZE, 5.000000E+00, 6.000000E+00, 0.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, OUTLET > OBJ, PRESSURE, P_AMBIENT > OBJ, COEFFICIENT, 1. > OBJ, NAME, ROAD > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 5.000000E+00, 0.000000E+00, 2.000000E+01 > OBJ, GEOMETRY, box > OBJ, TYPE, PLATE > OBJ, VELOCITY, 0. ,0. ,14. STOP