```
TALK=T;RUN( 1, 1)

************************************************************
Q1 created by VDI menu, Version 2009, Date 17/02/09
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
************************************************************
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) =VPOR
* Solved variables list
SOLVE(P1,U1,V1,W1)
* Stored variables list
STORE(VPOR)
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(VPOR)=1.
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 = 200
SARAH =1.0
RESFAC =1.0E-03
************************************************************
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
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
OUTPUT(VPOR,N,N,N,N,N,N)
************************************************************
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

PATCH(PROFILE, PROFIL, 7, 7, 6, 6, 1, 20, 1, 1)
PLOT(PROFILE, P1, 0., 0.)
PLOT(PROFILE, U1, 0., 0.)
PLOT(PROFILE, V1, 0., 0.)
PLOT(PROFILE, W1, 0., 0.)

PATCH(XEQ1, CONTUR, 1, 1, 1, 12, 1, 20, 1, 1)
PLOT(XEQ1, P1, 0., 20.)
PLOT(XEQ1, U1, 0., 20.)
PLOT(XEQ1, V1, 0., 20.)
PLOT(XEQ1, W1, 0., 20.)

PATCH(XEQ5, CONTUR, 5, 5, 1, 12, 1, 20, 1, 1)
PLOT(XEQ5, P1, 0., 20.)
PLOT(XEQ5, W1, 0., 20.)

PATCH(WAKE, CONTUR, 1, 8, 1, 12, 15, 15, 1, 1)
PLOT(WAKE, W1, 0., 20.)
************************************************************
Group 24. Dumps For Restarts

GVIEW(P,0.641234,0.496856,0.584767)
GVIEW(UP,-0.350212,0.867563,-0.353109)

> 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,    cube14
> 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,    cube14
> 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,    cube14
> 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,    cube14
> 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,    cube3t
> 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,    TEMPERATURE,  T_AMBIENT

> 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,    cube12t
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,     P_AMBIENT
> OBJ,    TEMPERATURE, SAME
> OBJ,    COEFFICIENT, 1.