TALK=T;RUN( 1, 1) DISPLAY Air flows through a porous plate with an area of 2m^2, a porosity of 0.15, and a loss coefficient of 0.4. The expected pressure drop is: Dp = 0.5*K*rho*Vel^2 = 0.5*0.4*1.189*(18/0.15)^2 = 3424.32 The resistive force on the U1 should be dP*Area*porosity F = 3424.32*2*0.15 = 1027.29 These values are exactly those which are computed and printed in RESULT. ENDDIS ************************************************************ Q1 created by VDI menu, Version 3.3, Date 10/02/00 CPVNAM=VDI;SPPNAM=Core ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Pressure Drop Across Porous Plate ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,10,3,3) * Set overall domain extent: * xulast yvlast zwlast name XSI= 1.000000E+01; YSI= 2.000000E+00; ZSI= 1.000000E+00 RSET(D,CHAM ) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T * Solved variables list SOLVE(P1 ,U1 ,V1 ,W1 ) * Additional solver options SOLUTN(P1 ,Y,Y,Y,N,N,Y) ************************************************************ Group 8. Terms & Devices ************************************************************ Group 9. Properties SETPRPS(1, 0) RHO1 = 1.189000E+00 PRESS0 = 1.000000E+05 TEMP0 = 2.730000E+02 CP1 = 1.005000E+03 ENUL = 1.544000E-05 ;ENUT = 0.000000E+00 DVO1DT = 3.410000E-03 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields 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 = 100 RESFAC = 1.000000E-03 ************************************************************ Group 16. Terminate Iterations ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX, 1.000000E+00) RELAX(U1 ,FALSDT, 1.000000E+00) RELAX(V1 ,FALSDT, 1.000000E+00) RELAX(W1 ,FALSDT, 1.000000E+00) ************************************************************ Group 18. Limits VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06 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 SPEDAT(SET,SATLIT,TALK,L,T) ************************************************************ Group 20. Preliminary Printout ECHO = T ************************************************************ Group 21. Print-out of Variables OUTPUT(P1 ,Y,Y,Y,N,Y,Y) OUTPUT(U1 ,Y,Y,Y,N,Y,Y) OUTPUT(V1 ,Y,Y,Y,N,Y,Y) OUTPUT(W1 ,Y,Y,Y,N,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 2 ;IZMON = 2 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 NXPRIN = 1 ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts NOWIPE = T GVIEW(P,7.390102E-01,-6.736943E-01,0.000000E+00) GVIEW(UP,0.000000E+00,0.000000E+00,1.000000E+00) > DOM, SIZE, 1.000000E+01, 2.000000E+00, 1.000000E+00 > DOM, MONIT, 5.000000E-01, 1.000000E+00, 5.000000E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, SNAPSIZE, 1.000000E-02 > DOM, RELAX, 5.000000E-01 > OBJ1, NAME, B1 > OBJ1, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ1, SIZE, 0.000000E+00, 2.000000E+00, 1.000000E+00 > OBJ1, CLIPART, cube3t > OBJ1, ROTATION24, 1 > OBJ1, TYPE, INLET > OBJ1, PRESSURE, 0.000000E+00 > OBJ1, VELOCITY, 1.800000E+01, 0.000000E+00, 0.000000E+00 > OBJ1, TEMPERATURE, 0.000000E+00 > OBJ2, NAME, B2 > OBJ2, POSITION, 5.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ2, SIZE, 0.000000E+00, 2.000000E+00, 1.000000E+00 > OBJ2, CLIPART, cube11 > OBJ2, ROTATION24, 1 > OBJ2, TYPE, PLATE > OBJ2, POROSITY, 1.500000E-01 > OBJ2, RESISTANCE, 4.000000E-01 > OBJ2, SIDE, BOTH > OBJ3, NAME, B3 > OBJ3, POSITION, 1.000000E+01, 0.000000E+00, 0.000000E+00 > OBJ3, SIZE, 0.000000E+00, 2.000000E+00, 1.000000E+00 > OBJ3, CLIPART, cube12t > OBJ3, ROTATION24, 1 > OBJ3, TYPE, OUTLET > OBJ3, PRESSURE, 0.000000E+00 > OBJ3, TEMPERATURE, -1.026000E+04 > OBJ3, COEFFICIENT, 1.000000E+03 M STOP distil=t EX(P1 )= 2.277E+03 EX(U1 )= 2.933E+01 EX(V1 )= 2.097E-04 EX(W1 )= 1.461E-04 EX(EPOR)= 9.150E-01