TALK=T;RUN( 1, 1)
DISPLAY
FLOW IN OPEN TURNAROUND CHANNEL
Steady Free Surface Flow of Water in an open channel with U-turn.
The task is to simulate the velocity distribution along with free
water surface elevation.
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 IN OPEN TURNAROUND CHANNEL )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,1,10,20)
* Set overall domain extent:
* xulast yvlast zwlast
name
XSI= 1.000000E+00; YSI= 9.000000E+01; ZSI= 1.450000E+02
RSET(D,CHAM )
************************************************************
Group 6. Body-Fitted coordinates
BFC=T
INTEGER(NI,NJ,NK,K1,K2);REAL(RR1,RR2,Y00)
NI=2;NJ=NY+1;NK=NZ+1;K1=5;K2=17;RR1=15.0
Y00=80.0+RR1;RR2=RR1+30.0
Inner radius of bend
DOMAIN(1,2,1,1,K1,K2);SETLIN(YC,Y00-RR1*COS(LNK*3.1416))
SETLIN(ZC,-RR1*SIN(LNK*3.1416))
Outer radius of bend
DOMAIN(1,2,NJ,NJ,K1,K2);SETLIN(YC,Y00-RR2*COS(LNK*3.1416))
SETLIN(ZC,-RR2*SIN(LNK*3.1416))
DOMAIN(1,2,1,NJ,1,K1);SETLIN(ZC,100.0*(1.0-LNK))
SETLIN(YC,80.0-LNJ*30.0)
DOMAIN(1,2,1,NJ,K2,NK);SETLIN(ZC,100.0*LNK)
SETLIN(YC,110.0+LNJ*30.0)
DOMAIN(1,2,1,NJ,K1,K2);MAGIC(T)
SETPT(1,1,1,0.0,80.0,100.0);SETPT(2,NJ,K1,1.0,50.0,0.0)
DOMAIN(1,2,1,NY,1,5);SETLIN(YC,YF+LNJ*(YL-YF))
SETLIN(ZC,ZF+LNK*(ZL-ZF))
DUMPC(grid)
READCO(grid)
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(147) =RHO1 ; NAME(148) =WCRT
NAME(149) =VCRT ; NAME(150) =UCRT
* Solved variables list
SOLVE(P1 ,V1 ,W1 )
* Stored variables list
STORE(UCRT,VCRT,WCRT,RHO1)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,N)
************************************************************
Group 8. Terms & Devices
NEWRH1 = T
************************************************************
Group 9. Properties
REAL(HIN,CON);ENUL=0.0;ENUT=0.0;RHO2=1.;HIN=1.
Compressibility is activated
RHO1=GRND3;DRH1DP=RHO1
CON=9.81/2.
RHO1A=CON**-0.5;PRESS0=9.81*HIN**2/2.
RHO1B=0.5;RHO1C=0.0
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(W1 ) = 1.000000E+00
No PATCHes used for this Group
INIADD = T
************************************************************
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 = 250
RESREF(P1 ) = 3.000000E-07 ;RESREF(V1 ) = 3.000000E-07
RESREF(W1 ) = 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(V1 ,FALSDT, 1.000000E+01)
RELAX(W1 ,FALSDT, 1.000000E+01)
************************************************************
Group 18. Limits
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
VARMAX(W1 ) = 5.000000E+00 ;VARMIN(W1 ) =-1.000000E+11
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
RSG1 = 7.500000E-01
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 5 ;IZMON = 10
NPRMON = 250
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 250
NPLT = 5
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 2
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
NOWIPE = T
GVIEW(P,-1.000000E+00,0.000000E+00,0.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
> DOM, SIZE, 1.000000E+00, 1.000000E+01, 2.000000E+01
> DOM, MONIT, 1.000000E+00, 5.000000E+00, 1.000000E+01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> DOM, RELAX, 3.000000E-01
> OBJ1, NAME, INLET
> OBJ1, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ1, SIZE, 1.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ1, CLIPART, INLET
> OBJ1, TYPE, INLET
> OBJ1, PRESSURE, 0.000000E+00
> OBJ1, VELOCITY, 0.000000E+00, 0.000000E+00,-1.000000E+00
> OBJ1, TEMPERATURE, 0.000000E+00
> OBJ2, NAME, OUTLET
> OBJ2, POSITION, 0.000000E+00, 0.000000E+00, 2.000000E+01
> OBJ2, SIZE, 1.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ2, CLIPART, OUTLET
> OBJ2, TYPE, OUTLET
> OBJ2, PRESSURE, 0.000000E+00
> OBJ2, TEMPERATURE, -1.026000E+04
> OBJ2, COEFFICIENT, 1.000000E+03
STOP
distil=t
EX(P1 )= 1.688E-01
EX(V1 )= 2.593E-02
EX(W1 )= 1.005E+00
EX(RHO1)= 1.000E+00
EX(WCRT)= 7.831E-01
EX(VCRT)= 4.137E-01