TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.6, Date 21/10/05
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
This unsteady case considers the prediction of vortex
shedding past a square cylinder at Re = 22000 by means
of a 2D Large Eddy Simulation. The cylinder dimension
is 0.4m and the free-stream water velocity is 0.535m/s.
The computations are initiated from solution fields
produced by a 2000-sweep, unconverged steady computation.
The basic Smagorinsky model is used with a coefficient
of 0.1 and no wall damping. The linear upwind scheme is
used for the momentum equations, and the near-wall region
is bridged by standard wall functions. The simulation
is advanced in time for 7 seconds using a uniform time
step of 1ms, corresponding to roughly 14 shedding cycles.
The INFORM facility is used to reduce gradually with
time the number of sweeps per time step.
The main results are listed below:
LES Data
Strouhal number 0.15 0.13
Mean Drag Coefficient 2.17 2.10
RMS Drag Coefficient 0.14 0.12
RMS Lift Coefficient 1.17 1.20
This run takes 2.25 hours to complete the 7000 time steps
on a 3GHZ PC with 2GB RAM.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(2D LES SQUARE CYLINDER- unsteady run )
************************************************************
Group 2. Transience
STEADY=F
* Set overall time and no. of steps
RSET(U,0.000000E+00,7.000000E+00,7000)
* Modify regions
RSET(T,1,7000,1.000000E+00)
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,120,85,1)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(137) =FMU ; NAME(138) =AREN
NAME(139) =AREE ; NAME(140) =GEN1
NAME(142) =SHRY ; NAME(143) =SHRX
NAME(144) =STRS ; NAME(145) =SKIN
NAME(146) =YPLS ; NAME(147) =WDIS
NAME(148) =DEN1 ; NAME(149) =EL1
NAME(150) =ENUT
* Solved variables list
SOLVE(P1 ,U1 ,V1 )
* Stored variables list
STORE(ENUT,EL1 ,DEN1,WDIS,YPLS,SKIN,STRS,SHRX)
STORE(SHRY,GEN1,AREE,AREN,FMU )
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
TURMOD(SGSMOD_NOWD)
************************************************************
Group 8. Terms & Devices
DIFCUT = 0.000000E+00
SCHEME(LUS ,U1 ,V1 )
************************************************************
Group 9. Properties
PRESS0 = 1.000000E+05 ;TEMP0 = 2.730000E+02
* Domain material index is 67 signifying:
* WATER at 20. deg C
SETPRPS(1, 67)
EL1A = 1.000000E-01 ;EL1B = 0.000000E+00
EL1C = 0.000000E+00
DVO1DT = 1.180000E-04
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
RESTRT(ALL)
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 = 1000
RESFAC = 1.000000E-03
************************************************************
Echo InForm settings for Group 15
inform15begin
(stored lsweep is 200 with IF(ISTEP.gt.3.and.ISTEP.le.6)!TSTSTR)
(stored lsweep is 100 with IF(ISTEP.GT.6.AND.ISTEP.LE.8)!TSTSTR)
(stored lsweep is 50 with IF(ISTEP.GT.8.AND.ISTEP.LE.10)!TSTSTR)
(stored lsweep is 40 with IF(ISTEP.GT.10.AND.ISTEP.LE.14)!TSTSTR)
(stored lsweep is 30 with IF(ISTEP.GT.14.AND.ISTEP.LE.20)!TSTSTR)
(stored lsweep is 20 with IF(ISTEP.GT.20)!TSTSTR)
inform15end
************************************************************
Group 16. Terminate Iterations
LITER (P1 ) = 300
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 1.000000E+00)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
************************************************************
Group 18. Limits
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
GENK = T
ASAP = T
PARSOL = F
CONWIZ = T
CALFOR = T
IDISPB = 1 ;IDISPC = 7000
CG( 1) ='DRAG'
************************************************************
Group 20. Preliminary Printout
ECHO = T
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL = 1.000000E-02
EX (P1 ) = 0.000000E+00 ;EX (U1 ) = 0.000000E+00
EX (V1 ) = 0.000000E+00 ;EX (FMU ) = 0.000000E+00
EX (AREN) = 0.000000E+00 ;EX (AREE) = 0.000000E+00
EX (GEN1) = 0.000000E+00 ;EX (LTLS) = 0.000000E+00
EX (SHRY) = 0.000000E+00 ;EX (SHRX) = 0.000000E+00
EX (STRS) = 0.000000E+00 ;EX (SKIN) = 0.000000E+00
EX (YPLS) = 0.000000E+00 ;EX (WDIS) = 0.000000E+00
EX (DEN1) = 0.000000E+00 ;EX (EL1 ) = 0.000000E+00
EX (ENUT) = 0.000000E+00
************************************************************
Group 21. Print-out of Variables
OUTPUT(GEN1,Y,N,N,N,N,N)
OUTPUT(WDIS,Y,N,N,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 68 ;IYMON = 43 ;IZMON = 1
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 10000
NXPRIN = 1
IXPRF = 20 ;IXPRL = 56
NYPRIN = 1
IYPRF = 25 ;IYPRL = 60
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
IDISPA = 500 ;IDISPB = 1 ;IDISPC = 7000
CSG1 ='m'
GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
> DOM, SIZE, 8.000000E-01, 5.600000E-01, 4.000000E-01
> DOM, MONIT, 2.607040E-01, 2.800000E-01, 2.000000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> GRID, RSET_X_1, 27, 9.000000E-01,G
> GRID, RSET_X_2, 22, 1.000000E+00
> GRID, RSET_X_3, 71, 1.040000E+00,G
> GRID, RSET_Y_1, 30, 9.000000E-01,G
> GRID, RSET_Y_2, 25, 1.000000E+00
> GRID, RSET_Y_3, 30, 1.100000E+00,G
> OBJ, NAME, SQCL
> OBJ, POSITION, 1.800000E-01, 2.600000E-01, 0.000000E+00
> OBJ, SIZE, 4.000000E-02, 4.000000E-02, 4.000000E-01
> OBJ, GEOMETRY, cube
> OBJ, ROTATION24, 1
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, NAME, INLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 5.600000E-01, 4.000000E-01
> OBJ, GEOMETRY, cube3t
> OBJ, ROTATION24, 1
> OBJ, TYPE, INLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, VELOCITY, 5.350000E-01, 0.000000E+00, 0.000000E+00
> OBJ, TEMPERATURE, 0.000000E+00
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
> OBJ, NAME, OUTLET
> OBJ, POSITION, 8.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 5.600000E-01, 4.000000E-01
> OBJ, GEOMETRY, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, SAME
> OBJ, COEFFICIENT, 1.000000E+03
> OBJ, TIME_LIMITS, ALWAYS_ACTIVE
STOP
** These parameters are used in GROUND and define the inlet veloc$
and cylinder dimension, respectively.
RG(1)=0.535;RG(2)=0.04