TALK=F;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2007, Date 16/01/07
CPVNAM=VDI;SPPNAM=Core
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(Test 21 Flow around of cylinder (fine grid).
This case is differed from test 20 by size of object.
This case is more correct because eddy after cylinder
lies into domain.
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,40,1,60)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(148) =DEN1
* Solved variables list
SOLVE(P1 ,U1 ,W1 )
* Stored variables list
STORE(DEN1)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
BOOLEAN(lAUTO)
mesg(Do you want to use automatic mode (y) or manual (n)?
readvdu(ans,char,n)
if(:ans:.eq.y)then
lAUTO=T
else
lAUTO=F
endif
INTEGER(LEVEL)
LEVEL = 0
if(lAUTO)then
mesg(Level of refinement is :LEVEL:
mesg(Do you want to change it? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
mesg(Enter level of refinement
readvdu(LEVEL,int,LEVEL)
mesg(New level of refinement is :LEVEL:
endif
endif
USP = T
UAUTO = lAUTO
USPDBG = F
UTCPLT = F
USPVTK = T
USPIMB = F
MXLEV = 0
MYLEV = 0
MZLEV = 0
DOMAT = -1
MINPRP = -1
MAXPRP = 100
CELLST = 10
FACEST = 10
UMXLEV = LEVEL ; UMXCLS= 100000
UNOLAY = 2
UIOLEV = 1
UWALEV = 3
UMXCRS = 16
if(lAUTO)then
mesg(Do you want to use PARSOL=T (y) or PARSOL=F (n)?
readvdu(ans,char,n)
if(:ans:.eq.y)then
PARSOL=T
else
PARSOL=F
endif
else
PARSOL = F
endif
************************************************************
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
SETPRPS(1, 0)
ENUT = 0.000000E+00
DVO1DT = 3.410000E-03
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
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 = 2000
RESFAC = 1.000000E-06
************************************************************
Group 16. Terminate Iterations
LITER (P1 ) = 200
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 7.E-01)
RELAX(U1 ,FALSDT, 1.E+05)
RELAX(W1 ,FALSDT, 1.E+05)
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, 0.3)
RELAX(U1 ,FALSDT, 1.E-01)
RELAX(W1 ,FALSDT, 1.E-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(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
ASAP = T
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 19 ;IYMON = 1 ;IZMON = 53
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
GVIEW(P,0.000000E+00,9.995500E-01,-2.999551E-02)
GVIEW(UP,1.000000E+00,0.000000E+00,0.000000E+00)
> DOM, SIZE, 2.000000E-01, 1.000000E-02, 3.000000E-01
> DOM, MONIT, 9.500000E-02, 5.000000E-03, 2.650000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> GRID, AUTO, F T F
> OBJ, NAME, CYL
> OBJ, POSITION, 6.000000E-02, 0.000000E+00, 5.000000E-02
> OBJ, SIZE, 8.000000E-02, 1.000000E-02, 8.000000E-02
> OBJ, GEOMETRY, cylinder
> OBJ, ROTATION24, 7
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, INLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 2.000000E-01, 9.999998E-03, 0.000000E+00
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 1.000000E-02
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, NAME, OUT
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 3.000000E-01
> OBJ, SIZE, 2.000000E-01, 9.999998E-03, 0.000000E+00
> OBJ, GEOMETRY, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, SAME
> OBJ, COEFFICIENT, 1.000000E+03
STOP