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, steady, incompressible, turbulent flow
past a surface-mounted square rib in a channel. It is the
equivalent of library case t915, using the low Reynolds-number
form of the Spalart-Allmaras turbulence model.
This case has been studied experimentally by D.Crabb et al, Proc.
4th Brazilian Congress on Mech. Engng., Florianopolis, Brazil,
p415,(1997).
The height H of the rib is 8.5% of that of the channel. The flow
Reynolds number based on channel bulk velocity and rib height
H is 300,000. The inlet plane is located 6H upstream of the rib,
and the outlet plane 20H downstream of the rib. A fixed-pressure
boundary condition is applied at the outlet, and uniform flow
profiles are specified at the inlet.
For this case, the main parameter characterising separation is
the length of the separation zone behind the rib. The experimental
and computed results are:
SA EXPT
Lr/H = 14.7 12.3
where the separation length Lr is measured from the front of the
rib. These results have not been checked for grid independence.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(Low Reynolds S-A Square Rib Flow )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,1,271,442,1.0E-03,1.0E-05,1.0E-05)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(140)=PRPS ;NAME(141)=YPLS
NAME(142)=DWDY ;NAME(143)=DVDZ
NAME(145)=WDIS ;NAME(146)=VOR1
NAME(148)=DEN1 ;NAME(149)=EL1
NAME(150) =ENUT
* Solved variables list
SOLVE(P1,V1,W1)
* Stored variables list
STORE(ENUT,EL1,DEN1,VOR1,WDIS,DVDZ,DWDY,YPLS)
STORE(PRPS)
* Additional solver options
SOLUTN(P1,Y,Y,Y,N,N,Y)
SOLUTN(V1,Y,Y,Y,N,N,Y)
SOLUTN(W1,Y,Y,Y,N,N,Y)
TURMOD(SPALART-ALLMARAS-LOWRE)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
PRESS0 =1.01325E+05 ;TEMP0 =273.149994
* Domain material index is 0 signifying:
* Air at 20 deg C, 1 atm, treated as incompressible
SETPRPS(1, 0)
DVO1DT =3.41E-03
PRT(ENTI)=0.666667
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1)=0. ;FIINIT(V1)=1.0E-03
FIINIT(W1)=1. ;FIINIT(PRPS)=-1.
FIINIT(WDIS)=0.1 ;FIINIT(ENTI)=0.05
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 = 8000
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
LITER(P1)=200 ;LITER(V1)=60
LITER(W1)=60 ;LITER(ENTI)=500
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
RELAX(V1 ,FALSDT,0.06 )
RELAX(W1 ,FALSDT,0.06 )
RELAX(LTLS,LINRLX,1. )
RELAX(ENTI,LINRLX,0.3 )
************************************************************
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.2356 ;EX(V1)=0.1149
EX(W1)=0.5498 ;EX(PRPS)=0.8944
EX(YPLS)=3.509E-03 ;EX(DWDY)=9.813
EX(DVDZ)=3.342 ;EX(LTLS)=2.372
EX(WDIS)=0.6016 ;EX(VOR1)=12.78
EX(ENTI)=5.727E-03 ;EX(DEN1)=1.063
EX(EL1)=0. ;EX(ENUT)=5.72E-03
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 110 ;IZMON = 150
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.999975,4.999979E-03,4.999917E-03)
GVIEW(UP,4.999917E-03,0.999988,-2.499979E-05)
GVIEW(WINDOW,MAX)
GVIEW(DEPTH,100.)
GVIEW(VDIS,14.142192)
GVIEW(CENTRE,0.5,5.875,13.5)
> DOM, SIZE, 1.000000E+00, 1.175000E+01, 2.700000E+01
> DOM, MONIT, 5.000000E-01, 9.998668E-01, 5.999668E+00
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID, AUTO, T F F
> GRID, RSET_X_1, 1, 1.000000E+00
> GRID, RSET_Y_1, -110, 1.100000E+00,G
> GRID, RSET_Y_2, 1, 1.080000E+00,G
> GRID, RSET_Y_3, -160, 1.100000E+00,G
> GRID, RSET_Z_1, -150, 1.100000E+00,G
> GRID, RSET_Z_2, -1, 1.000000E+00,G
> GRID, RSET_Z_3, -115, 1.100000E+00,G
> GRID, RSET_Z_4, 1, 1.000000E+00,G
> GRID, RSET_Z_5, 175, 1.050000E+00,G
> DOM, INI_BUOY, YES
> OBJ, NAME, WALLS
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, TO_END, 0.000000E+00, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, box
> OBJ, TYPE, PLATE
> OBJ, NAME, WALLN
> OBJ, POSITION, 0.000000E+00, AT_END, 0.000000E+00
> OBJ, SIZE, TO_END, 0.000000E+00, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, box
> OBJ, TYPE, PLATE
> OBJ, NAME, INLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, TO_END, TO_END, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, box
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, VELOCITY, 0. ,0. ,1.
> OBJ, ENTI_IN, 1.338000E-02
> OBJ, NAME, OUT
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, AT_END
> OBJ, SIZE, TO_END, TO_END, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, box
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, COEFFICIENT, 1000.
> OBJ, TURBULENCE, 0. ,0.
> OBJ, NAME, RIB
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 6.000000E+00
> OBJ, SIZE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ, GEOMETRY, box
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, NULL1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 5.999785E+00
> OBJ, SIZE, 1.000000E+00, 1.000000E+00, 1.000400E+00
> OBJ, GEOMETRY, wirexyz
> OBJ, GRID, N,N,Y
> OBJ, TYPE, NULL
> OBJ, NAME, NULL2
> OBJ, POSITION, 0.000000E+00, 1.000000E+00, 6.000000E+00
> OBJ, SIZE, 1.000000E+00, 2.000000E-04, 1.000000E+00
> OBJ, GEOMETRY, wirexyz
> OBJ, GRID, N,Y,N
> OBJ, TYPE, NULL
STOP
SPEDAT(SET,SELREF,PRNSL2,L,T)