TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2025, Date 03/02/25
CPVNAM=VDI; SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
The case considered is 3D, steady, incompressible, turbulent flow
past a surface-mounted cube in a channel. The flow separates in
front of the cube to form a primary and secondary vortex, and
the main vortex wraps as a horse-shoe vortex around the cube into
the wake. The flow separates at the front corners of the cube on
the roof and the side walls; the reattaches on the side walls but
not on the roof. A large separation region develops behind the
cube which interacts with the horseshoe vortex. In the experiments
vortex shedding is observed from the side walls, and due to
momentum exchange with the wake, this will lead to a shorter
separation length than is reported here for a steady simulation.
This case is the equivalent of library case t308, for the high
Reynolds-number form of the Spalart-Allmaras turbulence model,
with scalable wall functions.
The height of the cube is 50% of that of the channel. The flow
Reynolds number based on channel bulk velocity and cube height
H is 40,000. The inlet plane is located 7H upstream of the cube,
and the outlet plane 10H downstream of the cube. Because of
symmetry conditions, only half of the width of the flow is
calculated. A fixed-pressure boundary condition is applied at the
outlet, and uniform flow profiles are specified at the inlet.
The case has been studied experimentally by Martinuzzi & Tropea
[J.Fluids Engng, 115, p85-92,1993] and numerically by Lakehal &
Rodi [J.Wind Eng. Ind.Aerodyn, 67 & 68, p65-78, 1997].
For this case, the main parameters that characterise separation
are the frontal stagnation point Ys/H, the primary upstream
separation point Zf/H, the roof reattachment point Zr/H and the
length of the separation zone behind the cube Zb/H. The
experimental and computed results for Zb are given below:
SA EXPT
Zb/H = 2.6 1.61
These results are not grid independent, and the mesh is not fine
enough to resolve the expected separation on the roof nor to
capture adequately the upstream and downstream separation regions.
For this rather coarse mesh, the model overpredicts Zb.
It is likely that more mesh and the inclusion of unsteady effects
are required for a much improved prediction of the separation
length behind the cube.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 308
************************************************************
Group 1. Run Title
TEXT(High Reynolds S-A Surface Cube Flow )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,38,36,84)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(134)=WDIS ;NAME(135)=PRPS
NAME(136)=VPOR ;NAME(137)=YPLS
NAME(139)=DWDZ ;NAME(140)=DWDY
NAME(141)=DWDX ;NAME(142)=DVDZ
NAME(143)=DVDY ;NAME(144)=DVDX
NAME(145)=DUDZ ;NAME(146)=DUDY
NAME(147)=DUDX ;NAME(148)=VOR1
NAME(150) =ENUT
* Solved variables list
SOLVE(P1,U1,V1,W1)
* Stored variables list
STORE(ENUT,VOR1,DUDX,DUDY,DUDZ,DVDX,DVDY,DVDZ)
STORE(DWDX,DWDY,DWDZ,YPLS,VPOR,PRPS,WDIS)
* Additional solver options
SOLUTN(P1,Y,Y,Y,N,N,N)
SOLUTN(U1,Y,Y,N,N,N,N)
SOLUTN(V1,Y,Y,N,N,N,N)
SOLUTN(W1,Y,Y,N,N,N,N)
TURMOD(SPALART-ALLMARAS)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
RHO1 =1.
ENUL =2.5E-05
CP1 =1.
PRT(ENTI)=0.666667
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1)=1.3E-04 ;FIINIT(V1)=1.0E-03
FIINIT(W1)=1. ;FIINIT(WDIS)=0.2
FIINIT(PRPS)=-1. ;FIINIT(VPOR)=1.
FIINIT(YPLS)=1.001E-10 ;FIINIT(DWDZ)=1.001E-10
FIINIT(DWDY)=1.001E-10 ;FIINIT(DWDX)=1.001E-10
FIINIT(DVDZ)=1.001E-10 ;FIINIT(DVDY)=1.001E-10
FIINIT(DVDX)=1.001E-10 ;FIINIT(DUDZ)=1.001E-10
FIINIT(DUDY)=1.001E-10 ;FIINIT(DUDX)=1.001E-10
FIINIT(ENTI)=6.6143E-03 ;FIINIT(ENUT)=1.001E-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
EGWF = T
SCALWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 2500
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(P1)=50
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
RELAX(U1 ,FALSDT,0.107143 )
RELAX(V1 ,FALSDT,0.107143 )
RELAX(W1 ,FALSDT,0.107143 )
RELAX(ENTI,LINRLX,0.5 )
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
PARSOL = F
ISG62 = 1
SPEDAT(SET,OUTPUT,NOFIELD,L,T)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=0.06523 ;EX(U1)=0.03132
EX(V1)=0.02635 ;EX(W1)=0.955
EX(WDIS)=0.5288 ;EX(PRPS)=0.9774
EX(VPOR)=0.9774 ;EX(YPLS)=4.995
EX(LTLS)=0.3321 ;EX(DWDZ)=0.07485
EX(DWDY)=0.3162 ;EX(DWDX)=0.1339
EX(DVDZ)=0.04259 ;EX(DVDY)=0.05339
EX(DVDX)=0.03444 ;EX(DUDZ)=0.03931
EX(DUDY)=0.03799 ;EX(DUDX)=0.05951
EX(VOR1)=0.3786 ;EX(ENTI)=3.512E-03
EX(ENUT)=3.51E-03
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 32 ;IZMON = 80
NPRMON = 100
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NYPRIN = 2
NZPRIN = 2
IPLTF = 1 ;IPLTL = 1200
NPLT = 10
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,-0.969549,0.075535,0.232957)
GVIEW(UP,0.082366,0.996407,0.019723)
GVIEW(WINDOW,974,1039)
GVIEW(DEPTH,1000.)
GVIEW(VDIS,8.9575)
GVIEW(CENTRE,2.25,1.,9.)
> DOM, SIZE, 4.500000E+00, 2.000000E+00, 1.800000E+01
> DOM, MONIT, 2.083333E-02, 1.645468E+00, 1.508752E+01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID, RSET_X_1, 12, 1.000000E+00
> GRID, RSET_X_2, 26, 1.080000E+00,G
> GRID, RSET_Y_1, 18,-1.060000E+00,G
> GRID, RSET_Y_2, 18, 1.060000E+00,G
> GRID, RSET_Z_1, 38,-1.050000E+00,G
> GRID, RSET_Z_2, 12, 1.000000E+00
> GRID, RSET_Z_3, 34, 1.070000E+00,G
> DOM, T_AMBIENT, 0.000000E+00
> OBJ, NAME, CUBE
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 7.000000E+00
> OBJ, SIZE, 5.000000E-01, 1.000000E+00, 1.000000E+00
> OBJ, GEOMETRY, box
> 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, 4.500000E+00, 2.000000E+00, 0.000000E+00
> OBJ, GEOMETRY, box
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, VELOCITY, 0. ,0. ,1.
> OBJ, ENTI_IN, 6.614300E-03
> OBJ, NAME, OUTL
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.800000E+01
> OBJ, SIZE, 4.500000E+00, 2.000000E+00, 0.000000E+00
> OBJ, GEOMETRY, box
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.
> OBJ, COEFFICIENT, 1000.
> OBJ, TURBULENCE, 0. ,0.
> OBJ, NAME, WALLN
> OBJ, POSITION, 0.000000E+00, 2.000000E+00, 0.000000E+00
> OBJ, SIZE, 4.500000E+00, 0.000000E+00, 1.800000E+01
> OBJ, GEOMETRY, box
> OBJ, TYPE, PLATE
> OBJ, NAME, WALLS
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 4.500000E+00, 0.000000E+00, 1.800000E+01
> OBJ, GEOMETRY, box
> OBJ, TYPE, PLATE
STOP