```
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
************************************************************
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  )
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

************************************************************
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
```