TALK=F;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.6, Date 31/01/06
CPVNAM=VDI;SPPNAM=Core
************************************************************
title
DISPLAY
This case solves a two-dimensional steady hydrodynamics
problem about flow in labirinth channel.
The analytical solution is absent. The solution can be compared
with solution of Structured PHOENICS.
Cases changes inlet velocity as Uin=0.1*caseno
Domain material is air.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(USP labirinth. )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
INTEGER(NCUBE,NCHAN)
NCUBE = 80
NCHAN = NCUBE/20
RSET(M,NCUBE,NCUBE,1)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(148) = PRPS
* Solved variables list
SOLVE(P1 ,U1 ,V1)
* Stored variables list
STORE(PRPS)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
Usp related variables
USP = T
UAUTO = F
USPDBG = F
UTCPLT = F
USPVTK = T
USPIMB = F
MXLEV = 0
MYLEV = 0
MZLEV = 0
DOMAT = -1
MINPRP = -1
MAXPRP = 250
CELLST = 100
FACEST = 100
PARSOL = F
************************************************************
Group 9. Properties
PRESS0 = 1.000000E+05 ;TEMP0 = 2.730000E+02
* Domain material index is 0 signifying:
* Air at 20 deg C, 1 atm, treated as incompressible
RHO1 = 1.
ENUL = 1.e-5
REAL(Uin)
Uin=0.1
REAL(REYNO)
REYNO = UIN*0.05/ENUL
mesg(REYNO = :REYNO:
DVO1DT = 3.410000E-03
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
No PATCHes used for this Group
PATCH(B1,INIVAL, 4*NX/5+NCHAN+1,NX,1,NY,1,NZ,1,1)
COVAL(B1, PRPS, 0.0, 198)
PATCH(B2,INIVAL, 3*NX/5+NCHAN+1,4*NX/5+NCHAN,NY/5+NCHAN+1,NY,1,NZ,1,1)
COVAL(B2, PRPS, 0.0, 198)
PATCH(B3,INIVAL, 2*NX/5+NCHAN+1,3*NX/5+NCHAN,2*NY/5+NCHAN+1,NY,1,NZ,1,1)
COVAL(B3, PRPS, 0.0, 198)
PATCH(B4,INIVAL, NX/5+NCHAN+1,2*NX/5+NCHAN,3*NY/5+NCHAN+1,NY,1,NZ,1,1)
COVAL(B4, PRPS, 0.0, 198)
PATCH(B5,INIVAL, 1,NX/5+NCHAN,4*NY/5+NCHAN+1,NY,1,NZ,1,1)
COVAL(B5, PRPS, 0.0, 198)
PATCH(B6,INIVAL, 1,4*NX/5,1,NY/5,1,NZ,1,1)
COVAL(B6, PRPS, 0.0, 198)
PATCH(B7,INIVAL, 1,3*NX/5,NY/5+1,2*NY/5,1,NZ,1,1)
COVAL(B7, PRPS, 0.0, 198)
PATCH(B8,INIVAL, 1,2*NX/5,2*NY/5+1,3*NY/5,1,NZ,1,1)
COVAL(B8, PRPS, 0.0, 198)
PATCH(B9,INIVAL, 1,NX/5,3*NY/5+1,4*NY/5,1,NZ,1,1)
COVAL(B9, PRPS, 0.0, 198)
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 = 1000
RESFAC = 1.000000E-06
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 8.000000E-01)
RELAX(U1 ,FALSDT, 1.000000E+00)
RELAX(V1 ,FALSDT, 1.000000E+00)
mesg(Do you want to use collocated arrangement (y) or staggered one (n)? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
SPEDAT(SET,USP,METHOD,I,1)
RELAX(P1 ,LINRLX, 5.000000E-01)
RELAX(U1 ,FALSDT, 1.000000E-01)
RELAX(V1 ,FALSDT, 1.000000E-01)
mesg(Do you want to use SIMPLEST (y) or SIMPLE (n)? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
SPEDAT(SET,USP,SIMPLEST,L,T)
endif
endif
************************************************************
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 = F
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 50 ;IYMON = 42 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
mesg(Do you want to use max&min(y) on monitor instead of default(n)? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
#maxmin
endif
mesg(Do you want to view results in the centres of cells? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
SPEDAT(SET,USPIO,VERTCENT,L,F)
endif
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
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, 1.000000E+00, 1.000000E+00
> DOM, MONIT, 4.950000E-01, 4.150000E-01, 5.000000E-02
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> OBJ, NAME, IN
> OBJ, POSITION, 0.000000E+00, 8.000000E-01, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 5.000000E-02, 1.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, Uin, 0.000000E+00, 0.000000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, NAME, OUT
> OBJ, POSITION, 8.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 5.000000E-02, 0.000000E+00, 1.000000E+00
> OBJ, GEOMETRY, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, SAME
> OBJ, COEFFICIENT, 1.000000E+07
STOP