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*Vel^2 = 0.5*0.4*(18/0.15)^2 = 2880



  The resistive force on the U1 should be dP*Area*porosity



        F = 2880*2*0.15 = 864



  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