TALK=T;RUN( 1, 1)
DISPLAY
"WATER CURRENTS IN A RIVER COURSE "
VRGEOM shapemaker utility has been used here to create the
geometry object representing the topography of a river course
by converting BFC grid XYZ file.
The case presented makes use of the river geometry object to
simulate the water flow distribution on Cartesian grid.
This case deals with the 3D turbulent water flow in the river
course containing two islands and having nonlinear bed topography.
The free surface elevations are neglected assumed to be small
for essentially subcritical flow.
ENDDIS
************************************************************
Q1 created by VDI menu, Version 3.3, Date 27/01/00
CPVNAM=VDI;SPPNAM=Core
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(Water Currents In A River Course )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,45,45,10)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(149) =WDIS ; NAME(150) =PRPS
* Solved variables list
SOLVE(P1 ,U1 ,V1 ,W1 )
* Stored variables list
STORE(PRPS,WDIS)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
TURMOD(LVEL)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
PRESS0 = 1.000000E+05 ;TEMP0 = 2.730000E+02
SETPRPS(1, 67)
DVO1DT = 1.180000E-04
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(WDIS) = 1.000000E-03 ;FIINIT(PRPS) = 6.700000E+01
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
No PATCHes used for this Group
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 400
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 5.000000E-01)
RELAX(U1 ,FALSDT, 1.000000E-01)
RELAX(V1 ,FALSDT, 1.000000E-01)
RELAX(W1 ,FALSDT, 1.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
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
PARSOL = 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)
OUTPUT(WDIS,Y,N,N,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 20 ;IYMON = 25 ;IZMON = 10
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
NOWIPE = T
GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
> DOM, SIZE, 1.000000E+00, 4.600000E-01, 1.000000E-02
> DOM, MONIT, 4.333330E-01, 2.504440E-01, 9.500000E-03
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+01
> DOM, SNAPSIZE, 5.000000E-03
> OBJ, NAME, B2
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 4.600000E-01, 1.000000E-02
> OBJ, CLIPART, box
> OBJ, ROTATION24, 1
> OBJ, VISIBLE, -1
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198
> OBJ, NAME, B3
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 4.600000E-01, 1.000000E-02
> OBJ, CLIPART, vrg4
> OBJ, ROTATION24, 1
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 67
> OBJ, SCAL_FIXF, 0.000000E+00
> OBJ, NAME, B4
> OBJ, POSITION, 8.660001E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.339999E-01, 1.400000E-01, 1.000000E-02
> OBJ, CLIPART, wedge1
> OBJ, ROTATION24, 18
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 67
> OBJ, SCAL_FIXF, 0.000000E+00
> OBJ, NAME, OUTLETS
> OBJ, POSITION, 8.590000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.410000E-01, 0.000000E+00, 1.000000E-02
> OBJ, CLIPART, cube12t
> OBJ, ROTATION24, 1
> OBJ, VISIBLE, -1
> OBJ, GRID, 2
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, SAME
> OBJ, COEFFICIENT, 1.000000E+03
> OBJ, NAME, OUTLETE
> OBJ, POSITION, 1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 1.400000E-01, 1.000000E-02
> OBJ, CLIPART, cube12t
> OBJ, ROTATION24, 1
> OBJ, VISIBLE, -1
> OBJ, GRID, 2
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, SAME
> OBJ, COEFFICIENT, 1.000000E+03
> OBJ, NAME, B6
> OBJ, POSITION, 0.000000E+00, 2.990000E-01, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 1.610000E-01, 1.000000E-02
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, VISIBLE, -1
> OBJ, GRID, 2
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, 1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, TEMPERATURE, 0.000000E+00
STOP
distil=t
EX(P1 )= 3.051E+01
EX(U1 )= 1.733E-01
EX(V1 )= 6.934E-02
EX(W1 )= 2.604E-03
EX(PRPS)= 1.603E+02
EX(LTLS)= 1.163E-05
EX(WDIS)= 5.603E-04