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

FLOW IMPINGEMENT ON A BLUNT BODY

The physical situation simulates the impingement of uniform
supercritical (Froude number is 2.8) water flow on a solid block.

One symmetric half of the flow is considered due to symmetric
flow pattern. In adjusting to the stagnation conditions on the
block face the   supercritical stream is expected to undergo an
oblique hydrailic jump.

A subcritical region is then created between the jump and the
body. The subcritical flow is deflected round the body and
accelerates to supercritical flow again.

The task is to compute the distributions of velocities and
depths.

ENDDIS

************************************************************
Q1 created by VDI menu, Version 3.3, Date 24/07/00
CPVNAM=VDI;SPPNAM=Core
************************************************************
IRUNN   =       1 ;LIBREF =       0
************************************************************
Group 1. Run Title
TEXT(Flow impingement on a blunt body        )
************************************************************
Group 2. Transience
************************************************************
Groups 3, 4, 5  Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,40,20,1)
* Set overall domain extent:
*        xulast  yvlast  zwlast
name
XSI= 2.000000E+00; YSI= 2.100000E+00; ZSI= 1.000000E+00
RSET(D,CHAM    )
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS  =    T
* Non-default variable names
NAME(149) =VPOR ; NAME(150) =RHO1
* Solved variables list
SOLVE(P1  ,U1  ,V1  )
* Stored variables list
STORE(RHO1,VPOR)
SOLUTN(P1  ,Y,Y,Y,N,N,N)

************************************************************
Group 8. Terms & Devices
NEWRH1  =    T
************************************************************
Group 9. Properties
REAL(UIN,HIN,RHOR,CON);RHOR=1.;ENUL=0.0;ENUT=0.0;RHO2=1.
HIN=0.03;UIN=1.5
RHO1=GRND3;DRH1DP=RHO1
CON=9.81/(2.*RHOR)
RHO1A=CON**-0.5;PRESS0=9.81*RHOR*HIN**2/2.
RHO1B=0.5;RHO1C=0.0
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(U1  ) =  1.250000E-01 ;FIINIT(VPOR) =  1.000000E+00
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

************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP  =     500
RESREF(P1  ) =  3.000000E-07 ;RESREF(U1  ) =  3.000000E-07
RESREF(V1  ) =  3.000000E-07
RESFAC  = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
LITER (P1  ) =   10
ENDIT (P1  ) =  1.000000E-04
************************************************************
Group 17. Relaxation
RELAX(P1  ,LINRLX, 1.000000E+00)
RELAX(U1  ,FALSDT, 1.000000E+00)
RELAX(V1  ,FALSDT, 1.000000E+00)
************************************************************
Group 18. Limits
VARMAX(P1  ) = 1.000000E+10 ;VARMIN(P1  ) =-3.905000E-03
VARMAX(U1  ) = 5.000000E+00 ;VARMIN(U1  ) =-1.000000E+11
VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD  =    T  ;USEGRX =    T
ASAP    =    T
************************************************************
Group 20. Preliminary Printout
ECHO    =    F
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON   =      10 ;IYMON  =       5 ;IZMON  =       1
NPRMON  =     250
NPRMNT  =       1
TSTSWP  =      -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT  =  100000
ISWPRF  =       1 ;ISWPRL =  100000
ITABL   =       0
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
NOWIPE  =    T

GVIEW(P,-5.936508E-02,0.000000E+00,9.982364E-01)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)

> DOM,    SIZE,        2.000000E+00, 2.100000E+00, 1.000000E+00
> DOM,    MONIT,       4.750000E-01, 4.531250E-01, 5.000000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   RSET_X_1,     24, 1.000000E+00
> GRID,   RSET_X_2,     16, 1.000000E+00
> GRID,   RSET_Y_1,      4, 1.000000E+00
> GRID,   RSET_Y_2,     16, 1.000000E+00
> DOM,    RELAX,       5.000000E-01

> OBJ1,   NAME,        CMP0
> OBJ1,   POSITION,    1.200000E+00, 0.000000E+00, 0.000000E+00
> OBJ1,   SIZE,        8.000000E-01, 4.000000E-01, 1.000000E+00
> OBJ1,   CLIPART,     cube14
> OBJ1,   ROTATION24,        1
> OBJ1,   TYPE,        BLOCKAGE
> OBJ1,   MATERIAL,      199

> OBJ2,   NAME,        INLET
> OBJ2,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ2,   SIZE,        0.000000E+00, 2.100000E+00, 1.000000E+00
> OBJ2,   CLIPART,     cube3t
> OBJ2,   ROTATION24,        1
> OBJ2,   TYPE,        INLET
> OBJ2,   PRESSURE,      0.000000E+00
> OBJ2,   VELOCITY,      1.500000E+00, 0.000000E+00, 0.000000E+00
> OBJ2,   TEMPERATURE,   0.000000E+00

> OBJ3,   NAME,        OUTLET
> OBJ3,   POSITION,    2.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ3,   SIZE,        0.000000E+00, 2.100000E+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+02

> OBJ4,   NAME,        TOP
> OBJ4,   POSITION,    0.000000E+00, 2.100000E+00, 0.000000E+00
> OBJ4,   SIZE,        2.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ4,   CLIPART,     cube12t
> OBJ4,   ROTATION24,        1
> OBJ4,   TYPE,        OUTLET
> OBJ4,   PRESSURE,      0.000000E+00
> OBJ4,   TEMPERATURE,  -1.026000E+04
> OBJ4,   COEFFICIENT,   1.000000E+02
STOP
distil=t
EX(P1  )=   6.091E-03
EX(U1  )=   1.212E+00
EX(V1  )=   1.045E-01
EX(PRPS)=   9.200E-01
EX(VPOR)=   9.200E-01
EX(RHO1)=   3.847E-02
```