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
This case concerns plane, two-dimensional, incompressible
flow through a 180 degree turnaround duct, as studied
experimentally by Monson & Seegmiller (1988). The duct has a
width W=0.0381m and an aspect ratio of 10. The inlet and
outlet planes are located 3.5W upstream and downstream of the
bend, respectively. The inner radius of the bend is 0.01905m,
which corresponds to a curvature ratio of 0.5. The Reynolds
number based on duct width is 1.E5, and the Dean number is
7.071E4. The flow exhibits large streamline curvature with
flow relaminarisation along the inner convex wall together
with the formation of a separation zone near the bend exit
on the inner (convex) surface of the duct. The turnaround duct
is representative of many flows of engineering interest, such
as flow in the turnaround duct of the Space-Shuttle-Main-Engine
powerhead.
Calculations are made with the high-Re form of the
Spalart-Allmaras model, and the calculation employs a relatively
coarse non-uniform mesh of NY=25 and NZ=80. This is
the equivalent of Library Case t106 which uses alternative
turbulence models.
The model predicts the occurence of separation along the
inner wall at the exit of the U-bend. The current solution
remains sensitive to mesh numbers.
Low-Reynolds-number modelling, as exemplified by Library Case
t919, is recommended to predict the important features of
this flow.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 106
************************************************************
Group 1. Run Title
TEXT(High Reynolds S-A Turnaround Duct 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,25,80)
* Set overall domain extent:
* xulast yvlast zwlast
name
* Set overall domain extent:
* xulast yvlast zwlast
name
XSI= 1.000000E+00; YSI= 3.810000E-02; ZSI= 1.333500E-01
RSET(D,CHAM )
* Set objects: x0 y0 z0
* dx dy dz
name
XPO= 0.000000E+00; YPO= 0.000000E+00; ZPO= 0.000000E+00
XSI= 1.000000E+00; YSI= 3.810000E-02; ZSI= 0.000000E+00
RSET(B,IN )
XPO= 0.000000E+00; YPO= 0.000000E+00; ZPO= 1.333500E-01
XSI= 1.000000E+00; YSI= 3.810000E-02; ZSI= 0.000000E+00
RSET(B,OUTLET )
XPO= 0.000000E+00; YPO= 3.810000E-02; ZPO= 0.000000E+00
XSI= 1.000000E+00; YSI= 0.000000E+00; ZSI= 1.333500E-01
RSET(B,WFNN )
XPO= 0.000000E+00; YPO= 0.000000E+00; ZPO= 0.000000E+00
XSI= 1.000000E+00; YSI= 0.000000E+00; ZSI= 1.333500E-01
RSET(B,WFNS )
************************************************************
Group 6. Body-Fitted coordinates
BFC=T
* Set points
XPO= 0.000000E+00;YPO= 0.000000E+00;ZPO= 1.333500E-01;GSET(P,A )
XPO= 0.000000E+00;YPO= 3.810000E-02;ZPO= 1.333500E-01;GSET(P,B )
XPO= 1.000000E+00;YPO= 0.000000E+00;ZPO= 1.333500E-01;GSET(P,C )
XPO= 1.000000E+00;YPO= 3.810000E-02;ZPO= 1.333500E-01;GSET(P,D )
* Set lines/arcs
GSET(L,LAB,A,B,25,S1.7)
GSET(L,LBD,B,D,1,1.0)
GSET(L,LCD,C,D,25,S1.7)
GSET(L,LCA,C,A,1,1.0)
* Set frames
GSET(F,FABC,A,-,B,-,D,-,C,-)
* Match a grid mesh
GSET(M,FABC,+J+I,1,1,1,TRANS)
* Copy/Transfer/Block grid planes
GSET(C,K81,F,K1,1,1,1,25,+,0,0,0)
GSET(C,K81,F,K81,1,1,1,25,+,0,.0762,0)
GSET(C,K51,F,K81,1,1,1,25,+,0,0,-.13335,INC,.8)
GSET(C,K21,F,K51,1,1,1,25,RX,-3.14159,.05715,0,INC,1)
GSET(C,K1,F,K21,1,1,1,25,+,0,0,.13335,INC,1.2)
**********
NONORT = T
UUP = F ;VUP = F ;WUP = T
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(133)=DUDZ ;NAME(134)=DUDY
NAME(135)=DWDY ;NAME(136)=DVDZ
NAME(138)=WDIS ;NAME(139)=VOR1
NAME(141)=W1NR ;NAME(142)=TINT
NAME(143)=CP ;NAME(144)=STRS
NAME(145)=YPLS ;NAME(146)=LEN1
NAME(147)=ENUT ;NAME(148)=WCRT
NAME(149)=VCRT ;NAME(150)=UCRT
* Solved variables list
SOLVE(P1,V1,W1)
* Stored variables list
STORE(UCRT,VCRT,WCRT,ENUT,LEN1,YPLS,STRS,CP)
STORE(TINT,W1NR,VOR1,WDIS,DVDZ,DWDY,DUDY,DUDZ)
* Additional solver options
SOLUTN(P1,Y,Y,Y,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 =1.000125E-05
CP1 =1.
PRT(ENTI)=0.666667
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(W1)=26.25 ;FIINIT(WDIS)=0.1
FIINIT(ENTI)=2.459E-03 ;FIINIT(W1NR)=1.001E-10
FIINIT(TINT)=1.001E-10 ;FIINIT(STRS)=1.001E-10
FIINIT(YPLS)=1.001E-10 ;FIINIT(LEN1)=1.001E-10
FIINIT(ENUT)=1.001E-10 ;FIINIT(WCRT)=1.001E-10
FIINIT(VCRT)=1.001E-10 ;FIINIT(UCRT)=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
BFCA =1.
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1000
RESREF(P1)=1.000125E-12 ;RESREF(V1)=2.625328E-11
RESREF(W1)=2.625328E-11
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
RELAX(V1 ,FALSDT,6.35E-05 )
RELAX(W1 ,FALSDT,6.35E-05 )
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)=122. ;EX(V1)=0.6114
EX(W1)=25.08 ;EX(DUDZ)=0.
EX(DUDY)=0. ;EX(DWDY)=1735.
EX(DVDZ)=589.5 ;EX(LTLS)=9.508E-05
EX(WDIS)=7.157E-03 ;EX(VOR1)=2085.
EX(ENTI)=1.004E-03 ;EX(W1NR)=1.001E-10
EX(TINT)=1.001E-10 ;EX(CP)=1.0E-10
EX(STRS)=0.1668 ;EX(YPLS)=1.378
EX(LEN1)=1.001E-10 ;EX(ENUT)=1.001E-03
EX(WCRT)=21.450001 ;EX(VCRT)=6.764
EX(UCRT)=0.
************************************************************
Group 21. Print-out of Variables
OUTPUT(ENUT,Y,N,Y,N,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 2 ;IZMON = 56
NPRMON = 1000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 1000
NYPRIN = 1
NZPRIN = 1
IZPRF = 18 ;IZPRL = 24
NPLT = 2
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,-0.999937,9.99952E-03,-5.000544E-03)
GVIEW(UP,9.999771E-03,0.99995,-2.499602E-05)
GVIEW(WINDOW,MAX)
GVIEW(DEPTH,1000.)
GVIEW(VDIS,0.076185)
GVIEW(CENTRE,0.499654,0.0409,0.03015)
> DOM, SIZE, 1.000000E+00, 2.500000E+01, 8.000000E+01
> DOM, MONIT, 1.000000E+00, 2.000000E+00, 5.600000E+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, 1, 1.000000E+00
> GRID, RSET_Y_1, 25, 1.000000E+00
> GRID, RSET_Z_1, 20, 1.000000E+00
> GRID, RSET_Z_2, 30, 1.000000E+00
> GRID, RSET_Z_3, 30, 1.000000E+00
> DOM, T_AMBIENT, 0.000000E+00
> OBJ, NAME, IN
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 2.500000E+01, 0.000000E+00
> OBJ, GEOMETRY, IN
> OBJ, TYPE, INLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, VELOCITY, 0. ,0. ,-26.25
> OBJ, ENTI_IN, 2.459000E-03
> OBJ, NAME, OUTLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 8.000000E+01
> OBJ, SIZE, 1.000000E+00, 2.500000E+01, 0.000000E+00
> OBJ, GEOMETRY, OUTLET
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.
> OBJ, COEFFICIENT, 1.0E+04
> OBJ, TURBULENCE, 0. ,0.
> OBJ, NAME, WFNN
> OBJ, POSITION, 0.000000E+00, 2.500000E+01, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 0.000000E+00, 8.000000E+01
> OBJ, GEOMETRY, WFNN
> OBJ, TYPE, PLATE
> OBJ, NAME, WFNS
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, 1.000000E+00, 0.000000E+00, 8.000000E+01
> OBJ, GEOMETRY, WFNS
> OBJ, TYPE, PLATE
STOP