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 2D incompressible, turbulent flow
past a thick flat rectangular plate with a sharp leading
edge, as studied experimentally by Djilali,N. & Gartshore,I.S.,
"Turbulent flow around a bluff rectangular plate, Part I:
Experimental Investigation", Journal of Fluid Mechanics,
Vol.113, pp51-59, (1991).
This case is the equivalent of library case t311, for the high
Reynolds-number form of the Spalart-Allmaras turbulence model.
The flow separates at the leading edge of the plate, and then
reattaches further downstream to form a long separation zone on
top of the plate. The thickness of the plate H is taken as 0.1m,
and the working fluid is air at standard temperature and
pressure.The flow Reynolds number based on the free-stream
velocity and plate thickness H is 50,000.
The inlet and outlet planes are located 10H upstream and
downstream of the leading edge of the plate. The height of the
solution domain is taken as 10H, and for simplicity a zero
flux boundary is assumed here. Symmetry is exploited so that
only one half of the flow is simulated. A fixed-pressure boundary
condition is applied at the outlet, and uniform flow profiles are
specified at the inlet. Scalable wall functions are used at
the walls of the plate.
The measured and predicted reattachment points are as follows:
SA EXPT
Xr/H= 4.5 4.7
These results are not grid independent. However, the
Spalart-Allmaras model shows good agreement with the
experimental result.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 311
************************************************************
Group 1. Run Title
TEXT(High Reynolds S-A Blunt Flat Plate )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,115,75,1)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(136)=WDIS ;NAME(137)=PRPS
NAME(138)=SKIN ;NAME(139)=STRS
NAME(140)=YPLS ;NAME(142)=DVDY
NAME(143)=DVDX ;NAME(144)=DUDY
NAME(145)=DUDX ;NAME(146)=VOR1
NAME(148)=DEN1 ;NAME(150)=ENUT
* Solved variables list
SOLVE(P1,U1,V1)
* Stored variables list
STORE(ENUT,DEN1,VOR1,DUDX,DUDY,DVDX,DVDY,YPLS)
STORE(STRS,SKIN,PRPS,WDIS)
TURMOD(SPALART-ALLMARAS)
************************************************************
Group 8. Terms & Devices
DENPCO = T
************************************************************
Group 9. Properties
RHO1 =1.225
ENUL =1.46E-05
CP1 =1.
PRT(ENTI)=0.666667
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(U1)=7.3 ;FIINIT(V1)=7.300001E-03
FIINIT(WDIS)=0.1 ;FIINIT(PRPS)=-1.
FIINIT(SKIN)=1.001E-10 ;FIINIT(STRS)=1.001E-10
FIINIT(YPLS)=1.001E-10 ;FIINIT(DVDY)=1.001E-10
FIINIT(DVDX)=1.001E-10 ;FIINIT(DUDY)=1.001E-10
FIINIT(DUDX)=1.001E-10 ;FIINIT(DEN1)=1.001E-10
FIINIT(ENTI)=1.99939E-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 = 1200
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
RELAX(ENTI,LINRLX,0.5 )
************************************************************
Group 18. Limits
VARMAX(ENUT)=0.1 ;VARMIN(ENUT)=-1.0E+10
************************************************************
Group 19. EARTH Calls To GROUND Station
PARSOL = F
CONWIZ = T
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=5.333 ;EX(U1)=6.654
EX(V1)=0.42 ;EX(WDIS)=0.4522
EX(PRPS)=0.9362 ;EX(SKIN)=8.363E-05
EX(STRS)=2.857E-04 ;EX(YPLS)=0.1538
EX(DVDY)=5.153 ;EX(DVDX)=4.
EX(DUDY)=18.98 ;EX(DUDX)=5.137
EX(VOR1)=17.77 ;EX(LTLS)=0.4368
EX(DEN1)=1.147 ;EX(ENTI)=1.976E-03
EX(ENUT)=1.975E-03
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 73 ;IYMON = 12 ;IZMON = 1
NPRMON = 100
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NXPRIN = 23
NYPRIN = 15
IPLTF = 1 ;IPLTL = 1500
NPLT = 10
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,0.,0.,1.)
GVIEW(UP,0.,1.,0.)
GVIEW(WINDOW,974,1039)
GVIEW(DEPTH,1000.)
GVIEW(VDIS,1.211028)
GVIEW(CENTRE,1.,0.5,0.5)
> DOM, SIZE, 2.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, MONIT, 1.082791E+00, 5.489546E-02, 5.000000E-01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID, RSET_X_1, 60,-1.040000E+00,G
> GRID, RSET_X_2, 55, 1.040000E+00,G
> GRID, RSET_Y_1, 10, 1.040000E+00,G
> GRID, RSET_Y_2, 65, 1.040000E+00,G
> GRID, RSET_Z_1, 1, 1.000000E+00
> DOM, T_AMBIENT, 0.000000E+00
> OBJ, NAME, PLATE
> OBJ, POSITION, 9.999999E-01, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 5.000001E-02, 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, 0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ, GEOMETRY, box
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, VELOCITY, 7.3 ,0. ,0.
> OBJ, TURB-INTENS, 1.
> OBJ, NAME, OUTLET
> OBJ, POSITION, 2.000000E+00, 5.000001E-02, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 9.500000E-01, 1.000000E+00
> OBJ, GEOMETRY, box
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.
> OBJ, COEFFICIENT, 1000.
> OBJ, TURBULENCE, 0. ,0.
STOP