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

WATER POURING INTO A BUND

FREE SURFACE EXAMPLE - Water rise in a channel with
obstruction: 3D, elliptic, Scalar-Equation-Model.

The case represents a water flow within an open channel.
The channel has an immersed solid obstruction. The model
predicts the elevation of the water surface over the
blockage and the gravity-driven sideway 'spills'.

ENDDIS

************************************************************
Q1 created by VDI menu, Version 3.3, Date 24/07/00
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
PHOTON USE
p
n4

up z
gr x 1
use patgeo
msg contours of den1 at time = 0.08s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n8

up z
gr x 1
use patgeo
msg contours of den1 at time = 0.16s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n12

up z
gr x 1
use patgeo
msg contours of den1 at time = 0.24s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n16

up z
gr x 1
use patgeo
msg contours of den1 at time = 0.32s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n20

up z
gr x 1
use patgeo
msg contours of den1 at time = 0.40s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n24

up z
gr x 1
use patgeo
msg contours of den1 at time = 0.48s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n28

up z
gr x 1
use patgeo
msg contours of den1 at time = 0.56s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n32

up z
gr x 1
use patgeo
msg contours of den1 at time = 0.64s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n36

up z
gr x 1
use patgeo
msg contours of den1 at time = 0.72s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n40

up z
gr x 1
use patgeo
msg contours of den1 at time = 0.8s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n44

up z
gr x 1
use patgeo
msg contours of den1 at time = 0.88s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n48

up z
gr x 1
use patgeo
msg contours of den1 at time = 0.96s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n52

up z
gr x 1
use patgeo
msg contours of den1 at time = 1.04s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n56

up z
gr x 1
use patgeo
msg contours of den1 at time = 1.12s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n60

up z
gr x 1
use patgeo
msg contours of den1 at time = 1.2s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n64

up z
gr x 1
use patgeo
msg contours of den1 at time = 1.28s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n68

up z
gr x 1
use patgeo
msg contours of den1 at time = 1.36s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n72

up z
gr x 1
use patgeo
msg contours of den1 at time = 1.44s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n76

up z
gr x 1
use patgeo
msg contours of den1 at time = 1.52s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n80

up z
gr x 1
use patgeo
msg contours of den1 at time = 1.6s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n84

up z
gr x 1
use patgeo
msg contours of den1 at time = 1.68s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n88

up z
gr x 1
use patgeo
msg contours of den1 at time = 1.76s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n92

up z
gr x 1
use patgeo
msg contours of den1 at time = 1.84s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n96

up z
gr x 1
use patgeo
msg contours of den1 at time = 1.92s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
pause
p
n100

up z
gr x 1
use patgeo
msg contours of den1 at time = 2.00s
con den1 z 1 fi;.001
con den1 x 1 fi;.001
con den1 x 4 fi;.001
con den1 x 8 fi;.001
msg velocity vectors
vec x 1;vec z 1
msg            -
msg Press e to END
enduse

************************************************************
IRUNN   =       1 ;LIBREF =     106
************************************************************
Group 1. Run Title
TEXT(Water pouring into a bund)
************************************************************
Group 2. Transience
* Set overall time and no. of steps
RSET(U,0.000000E+00,2.000000E+00,100)
* Modify regions
************************************************************
Groups 3, 4, 5  Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,10,10,10)
* Set overall domain extent:
*        xulast  yvlast  zwlast
name
XSI= 1.000000E+00; YSI= 1.000000E+00; ZSI= 3.125000E-01
RSET(D,CHAM    )
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS  =    T
* Non-default variable names
NAME(147) =SURN ; NAME(148) =VFOL
NAME(149) =PRPS ; NAME(150) =DEN1
* Solved variables list
SOLVE(P1  ,U1  ,V1  ,W1  ,SURN,VFOL)
* Stored variables list
STORE(DEN1,PRPS)
* Additional solver options
SOLUTN(P1  ,Y,Y,Y,N,N,N)
SOLUTN(V1  ,Y,Y,N,N,N,N)
SOLUTN(W1  ,Y,Y,N,N,N,N)

************************************************************
Group 8. Terms & Devices
TERMS (SURN,N,N,N,N,N,Y)
TERMS (VFOL,N,N,N,N,Y,Y)
GALA    =    T
NEWRH1  =    T
NEWENL  =    T
RUPLIM  = 6.000000E-01 ;RLOLIM = 4.000000E-01
************************************************************
Group 9. Properties
DVO1DT  = 3.330000E-03
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1  ) =  0.000000E+00 ;FIINIT(V1  ) =  0.000000E+00
FIINIT(W1  ) =  0.000000E+00 ;FIINIT(SURN) =  0.000000E+00
FIINIT(PRPS) =  0.000000E+00 ;FIINIT(DEN1) =  1.189000E+00
No PATCHes used for this Group

************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources

PATCH (BUOYANCY,PHASEM,0,0,0,0,0,0,1,100)
COVAL (BUOYANCY,U1  , FIXFLU      , GRND1       )
COVAL (BUOYANCY,V1  , FIXFLU      , GRND1       )
COVAL (BUOYANCY,W1  , FIXFLU      , GRND1       )

BUOYA   = 0.000000E+00 ; BUOYB = 0.000000E+00
BUOYC   =-9.810000E+00
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP  =      40
RESREF(P1  ) =  1.000000E-06 ;RESREF(V1  ) =  1.000000E-06
RESREF(W1  ) =  1.000000E-06
RESFAC  = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
LITER (SURN) =    1
************************************************************
Group 17. Relaxation
RELAX(P1  ,LINRLX, 5.000000E-01)
RELAX(U1  ,FALSDT, 3.000000E-02)
RELAX(V1  ,FALSDT, 3.000000E-02)
RELAX(W1  ,FALSDT, 3.000000E-02)
RELAX(SURN,FALSDT, 1.000000E+01)
RELAX(VFOL,FALSDT, 1.000000E+00)
RELAX(DEN1,LINRLX, 4.000000E-01)
************************************************************
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
VARMAX(SURN) = 1.000000E+00 ;VARMIN(SURN) = 0.000000E+00
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD  =    T  ;USEGRX =    T
SURF    =    T
ASAP    =    T
IPRPSA  =      67
CSG1    ='N'
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
************************************************************
Group 20. Preliminary Printout
ECHO    =    T
distil=t
EX(P1  )=   4.275E+02
EX(U1  )=   3.480E-01
EX(V1  )=   7.489E-01
EX(W1  )=   3.694E-01
EX(SURN)=   4.914E-01
EX(VFOL)=   4.885E-01
EX(PRPS)=   3.839E+01
EX(DEN1)=   5.117E+02
************************************************************
Group 21. Print-out of Variables
OUTPUT(P1  ,Y,Y,Y,Y,Y,Y)
OUTPUT(U1  ,Y,Y,Y,N,Y,Y)
OUTPUT(V1  ,Y,Y,Y,Y,Y,Y)
OUTPUT(W1  ,Y,Y,Y,Y,Y,Y)
OUTPUT(SURN,Y,Y,Y,Y,Y,Y)
OUTPUT(VFOL,Y,N,Y,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
NTPRIN  =      10 ;ISTPRF =       1 ;ISTPRL =   10000
ISWPRF  =       1 ;ISWPRL =  100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
NOWIPE  =    T
IDISPA  =       2 ;IDISPB =       0 ;IDISPC =       0
CSG1    ='N'

GVIEW(P,-9.841644E-01,0.000000E+00,1.772584E-01)
GVIEW(UP,1.772584E-01,0.000000E+00,9.841644E-01)

> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 3.125000E-01
> DOM,    MONIT,       4.500000E-01, 4.500000E-01, 1.406250E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> DOM,    RELAX,       1.000000E+00

> OBJ1,   NAME,        SOUTHOUT
> OBJ1,   POSITION,    0.000000E+00, 0.000000E+00, 9.375000E-02
> OBJ1,   SIZE,        1.000000E+00, 0.000000E+00, 2.187500E-01
> OBJ1,   CLIPART,     cube13
> OBJ1,   ROTATION24,        1
> OBJ1,   GRID,        2.000000E+00
> OBJ1,   TYPE,        PLATE
> OBJ1,   SURF_TEMP,     0.000000E+00,-1.026000E+04

> OBJ2,   NAME,        HIGHOUT
> OBJ2,   POSITION,    0.000000E+00, 0.000000E+00, 3.125000E-01
> OBJ2,   SIZE,        1.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ2,   CLIPART,     cube12t
> OBJ2,   ROTATION24,        1
> OBJ2,   TYPE,        OUTLET
> OBJ2,   PRESSURE,      0.000000E+00
> OBJ2,   TEMPERATURE,  -1.026000E+04
> OBJ2,   COEFFICIENT,   1.000000E+00
> OBJ2,   TIME_LIMITS,   ALWAYS_ACTIVE

> OBJ3,   NAME,        NORTHOUT
> OBJ3,   POSITION,    8.000000E-01, 1.000000E+00, 0.000000E+00
> OBJ3,   SIZE,        2.000000E-01, 0.000000E+00, 3.125000E-01
> OBJ3,   CLIPART,     cube13
> OBJ3,   ROTATION24,        1
> OBJ3,   TYPE,        PLATE
> OBJ3,   SURF_TEMP,     0.000000E+00,-1.026000E+04

> OBJ4,   NAME,        HEAVYINL
> OBJ4,   POSITION,    5.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ4,   SIZE,        5.000000E-01, 0.000000E+00, 9.375000E-02
> OBJ4,   CLIPART,     cube13
> OBJ4,   ROTATION24,        1
> OBJ4,   GRID,        2.000000E+00
> OBJ4,   TYPE,        PLATE
> OBJ4,   SURF_TEMP,     0.000000E+00, 0.000000E+00
> OBJ4,   VELOCITY,      0.000000E+00, 2.000000E+00, 0.000000E+00

> OBJ5,   NAME,        BLOCK
> OBJ5,   POSITION,    0.000000E+00, 5.000000E-01, 0.000000E+00
> OBJ5,   SIZE,        1.000000E+00, 9.999999E-02, 9.450000E-02
> OBJ5,   CLIPART,     cube4
> OBJ5,   ROTATION24,        1
> OBJ5,   GRID,        2.000000E+00
> OBJ5,   TYPE,        BLOCKAGE
> OBJ5,   MATERIAL,      100
> OBJ5,   TIME_LIMITS,   ALWAYS_ACTIVE
> OBJ5,   FIXED_TMP,     0.000000E+00, 0.000000E+00

> OBJ6,   NAME,        EASTOUT
> OBJ6,   POSITION,    1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ6,   SIZE,        0.000000E+00, 1.000000E+00, 3.125000E-01
> OBJ6,   CLIPART,     cube13
> OBJ6,   ROTATION24,        1
> OBJ6,   GRID,        2.000000E+00
> OBJ6,   TYPE,        PLATE
> OBJ6,   FIXED_TMP,     0.000000E+00, 0.000000E+00

> OBJ7,   NAME,        B7
> OBJ7,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ7,   SIZE,        5.000000E-01, 0.000000E+00, 9.375000E-02
> OBJ7,   CLIPART,     cube3t
> OBJ7,   ROTATION24,        1
> OBJ7,   GRID,        2.000000E+00
> OBJ7,   TYPE,        INLET
> OBJ7,   FLUID,         1.000000E+00
> OBJ7,   VELOCITY,      0.000000E+00, 2.000000E+00, 0.000000E+00
> OBJ7,   TEMPERATURE,   0.000000E+00
> OBJ7,   TIME_LIMITS,   ALWAYS_ACTIVE

> OBJ8,   NAME,        B8
> OBJ8,   POSITION,    2.000000E-01, 1.000000E+00, 0.000000E+00
> OBJ8,   SIZE,        6.000000E-01, 0.000000E+00, 6.000000E-02
> OBJ8,   CLIPART,     cube12t
> OBJ8,   ROTATION24,        1
> OBJ8,   GRID,        2.000000E+00
> OBJ8,   TYPE,        OUTLET
> OBJ8,   PRESSURE,      0.000000E+00
> OBJ8,   TEMPERATURE,  -1.026000E+04
> OBJ8,   COEFFICIENT,   1.000000E+00
> OBJ8,   TIME_LIMITS,   ALWAYS_ACTIVE

> OBJ9,   NAME,        B10
> OBJ9,   POSITION,    0.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ9,   SIZE,        2.000000E-01, 0.000000E+00, 3.125000E-01
> OBJ9,   CLIPART,     cube13
> OBJ9,   ROTATION24,        1
> OBJ9,   GRID,        2.000000E+00
> OBJ9,   TYPE,        PLATE
> OBJ9,   SURF_TEMP,     0.000000E+00,-1.026000E+04

> OBJ10,  NAME,        B12
> OBJ10,  POSITION,    2.000000E-01, 1.000000E+00, 6.000000E-02
> OBJ10,  SIZE,        6.000000E-01, 0.000000E+00, 2.525000E-01
> OBJ10,  CLIPART,     cube13
> OBJ10,  ROTATION24,        1
> OBJ10,  GRID,        2.000000E+00
> OBJ10,  TYPE,        PLATE
> OBJ10,  SURF_TEMP,     0.000000E+00,-1.026000E+04
STOP
```