TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.6, Date 02/11/04
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
PHOTON USE
p
gr ou x 1
MSG Velocity vectors
vec x 1 sh
msg
msg Press return to plot pressure contours
pause
cont p1 x 1 fil;.1
msg
msg Type e to End
ENDUSE
DISPLAY
This case concerns plane, two-dimensional, incompressible invscid
flow through a 180 degree turnaround duct with two straight legs.
Note the "parameterization" of the GSET commands, enabling the
grid to be changed easily.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 527
************************************************************
Group 1. Run Title
TEXT(Shallow water IN TURNAROUND DUCT: B)
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,1,20,72)
* Set overall domain extent:
* xulast yvlast zwlast
name
* Set overall domain extent:
* xulast yvlast zwlast
name
XSI= 1.000000E+00; YSI= 3.000000E+01; ZSI= 1.000000E+02
RSET(D,CHAM )
* Set objects: x0 y0 z0
* dx dy dz
name
XPO= 0.000000E+00; YPO= 0.000000E+00; ZPO= 0.000000E+00
XSI= 1.000000E+00; YSI= 3.000000E+01; ZSI= 0.000000E+00
RSET(B,IN )
XPO= 0.000000E+00; YPO= 0.000000E+00; ZPO= 1.000000E+02
XSI= 1.000000E+00; YSI= 3.000000E+01; ZSI= 0.000000E+00
RSET(B,OUTLET )
************************************************************
Group 6. Body-Fitted coordinates
BFC=T
* Copy/Transfer/Block grid planes
GSET(C,K73,F,K73,1,1,1,20,+,0,110,0)
GSET(C,K49,F,K73,1,1,1,20,+,0,0,-100,INC,1)
GSET(C,K25,F,K49,1,1,1,20,RX,-3.14159,95,0,INC,1)
GSET(C,K1,F,K25,1,1,1,20,+,0,0,100,INC,1)
**********
NONORT = T
UUP = F ;VUP = F ;WUP = T
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(147) =DEN1 ; NAME(148) =WCRT
NAME(149) =VCRT ; NAME(150) =UCRT
* Solved variables list
SOLVE(P1 ,V1 ,W1 )
* Stored variables list
STORE(UCRT,VCRT,WCRT,DEN1)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,N)
************************************************************
Group 8. Terms & Devices
NEWRH1 = T
************************************************************
Group 9. Properties
PRESS0 = 4.905000E+00
RHO1 = GRND3
RHO1A = 4.515000E-01 ;RHO1B = 5.000000E-01
RHO1C = 0.000000E+00
ENUL = 0.000000E+00
CP1 = 1.000000E+00
ENUT = 0.000000E+00
DRH1DP = GRND3
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(W1 ) = 1.000000E+00 ;FIINIT(WCRT) = 1.001000E-10
FIINIT(VCRT) = 1.001000E-10 ;FIINIT(UCRT) = 1.001000E-10
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
BFCA = 1.000000E+00
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 500
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
LITER (P1 ) = 10
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 5.000000E-01)
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 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
SYMBFC = T
PARSOL = T
SPEDAT(SET,FACETDAT,NUMOBJ,I,3)
************************************************************
Group 20. Preliminary Printout
ECHO = T
NULLPR = F
NDST = 0
DSTTOL = 1.000000E-02
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 10 ;IZMON = 36
NPRMON = 500
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 500
NPLT = 5
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 2
PATCH (YX ,CONTUR,1,1,1,20,1,72,1,1)
PLOT(YX ,P1 , 0.000000E+00, 2.000000E+01)
PATCH (INNER ,PROFIL,1,1,1,1,1,72,1,1)
PLOT(INNER ,P1 , 0.000000E+00, 0.000000E+00)
PLOT(INNER ,W1 , 0.000000E+00, 0.000000E+00)
PATCH (OUTER ,PROFIL,1,1,20,20,1,72,1,1)
PLOT(OUTER ,P1 , 0.000000E+00, 0.000000E+00)
PLOT(OUTER ,W1 , 0.000000E+00, 0.000000E+00)
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,-9.935490E-01,-5.311805E-02,1.001937E-01)
GVIEW(UP,-5.391878E-02,9.985313E-01,-5.298903E-03)
> DOM, SIZE, 1.000000E+00, 2.000000E+01, 7.200000E+01
> DOM, MONIT, 1.000000E+00, 1.000000E+01, 3.600000E+01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> OBJ, NAME, IN
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 2.000000E+01, 0.000000E+00
> OBJ, GEOMETRY, IN
> OBJ, TYPE, INLET
> OBJ, COLOR-MODE, DEFAULT2
> OBJ, COLOR-VAL, 0
> OBJ, OPAQUE, 50
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00,-1.000000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, NAME, OUTLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 7.200000E+01
> OBJ, SIZE, 1.000000E+00, 2.000000E+01, 0.000000E+00
> OBJ, GEOMETRY, OUTLET
> OBJ, TYPE, OUTLET
> OBJ, COLOR-MODE, DEFAULT2
> OBJ, COLOR-VAL, 0
> OBJ, OPAQUE, 50
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, SAME
> OBJ, COEFFICIENT, 1.000000E+04
STOP