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
STEADY=F
* 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
INIADD = F
************************************************************
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