TALK=T;RUN(1,1)
DISPLAY
The problem concerns 2d incompressible, turbulent flow and heat
transfer over a backward-facing step at a Reynolds number of
45,000 based on step height Hs and inlet velocity. The expansion
ratio Hc/Hs = 3.0, where Hc is the downstream channel height 3.0.
The case is similar to that described for library case T103,
excepting that heat transfer is considered and the calculation is
performed with any one of the following low-Re models: the Lam-
Bremhorst k-e model, the Chen-Kim k-e model, the 2-layer k-e mod$
el,
and the Wilcox or Menter k-w models, or the k-w-SST model. The
thermal field is calculated for a constant heat flux through the
south wall downstream of the step and an adiabatic condition on
all other walls. The laminar Prandtl number is taken as 1.0 and
the turbulent Prandtl as 0.86.
ENDDIS
AUTOPLOT USE
file
phida 3
da 1 u1 y 1
divide x .0381 1
shift x -4 1
col9 1
level y 0;level x 0
scale x -4 15
msg Velocity (U1) profile
msg Press e to END
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Low-Re LB k-e BF Step Y-X :T208 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 14
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 80
XULAST =0.762
XFRAC(1)=0.027551 ;XFRAC(2)=0.053061
XFRAC(3)=0.076531 ;XFRAC(4)=0.097959
XFRAC(5)=0.117347 ;XFRAC(6)=0.134694
XFRAC(7)=0.15 ;XFRAC(8)=0.163265
XFRAC(9)=0.17449 ;XFRAC(10)=0.183673
XFRAC(11)=0.190816 ;XFRAC(12)=0.195918
XFRAC(13)=0.19898 ;XFRAC(14)=0.2
XFRAC(15)=0.200124 ;XFRAC(16)=0.200498
XFRAC(17)=0.20112 ;XFRAC(18)=0.201992
XFRAC(19)=0.203112 ;XFRAC(20)=0.204482
XFRAC(21)=0.2061 ;XFRAC(22)=0.207967
XFRAC(23)=0.210084 ;XFRAC(24)=0.212449
XFRAC(25)=0.215063 ;XFRAC(26)=0.217927
XFRAC(27)=0.221039 ;XFRAC(28)=0.2244
XFRAC(29)=0.22801 ;XFRAC(30)=0.23187
XFRAC(31)=0.235978 ;XFRAC(32)=0.240335
XFRAC(33)=0.244941 ;XFRAC(34)=0.249796
XFRAC(35)=0.2549 ;XFRAC(36)=0.260254
XFRAC(37)=0.265856 ;XFRAC(38)=0.271707
XFRAC(39)=0.277807 ;XFRAC(40)=0.284156
XFRAC(41)=0.290754 ;XFRAC(42)=0.297601
XFRAC(43)=0.304697 ;XFRAC(44)=0.312042
XFRAC(45)=0.319636 ;XFRAC(46)=0.327479
XFRAC(47)=0.33557 ;XFRAC(48)=0.343911
XFRAC(49)=0.352501 ;XFRAC(50)=0.36134
XFRAC(51)=0.370428 ;XFRAC(52)=0.379765
XFRAC(53)=0.38935 ;XFRAC(54)=0.399185
XFRAC(55)=0.409269 ;XFRAC(56)=0.419602
XFRAC(57)=0.430183 ;XFRAC(58)=0.441014
XFRAC(59)=0.452094 ;XFRAC(60)=0.463422
XFRAC(61)=0.475 ;XFRAC(62)=0.490334
XFRAC(63)=0.510228 ;XFRAC(64)=0.532306
XFRAC(65)=0.555933 ;XFRAC(66)=0.580784
XFRAC(67)=0.606655 ;XFRAC(68)=0.633406
XFRAC(69)=0.660936 ;XFRAC(70)=0.689164
XFRAC(71)=0.718027 ;XFRAC(72)=0.747475
XFRAC(73)=0.777463 ;XFRAC(74)=0.807956
XFRAC(75)=0.838922 ;XFRAC(76)=0.870334
XFRAC(77)=0.902168 ;XFRAC(78)=0.934403
XFRAC(79)=0.967019 ;XFRAC(80)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 80
YVLAST =0.1143
YFRAC(1)=1.023523E-03 ;YFRAC(2)=3.325437E-03
YFRAC(3)=6.625275E-03 ;YFRAC(4)=0.010804
YFRAC(5)=0.015789 ;YFRAC(6)=0.021526
YFRAC(7)=0.027975 ;YFRAC(8)=0.035104
YFRAC(9)=0.042885 ;YFRAC(10)=0.051298
YFRAC(11)=0.060321 ;YFRAC(12)=0.069937
YFRAC(13)=0.080131 ;YFRAC(14)=0.09089
YFRAC(15)=0.1022 ;YFRAC(16)=0.114052
YFRAC(17)=0.126433 ;YFRAC(18)=0.139335
YFRAC(19)=0.152749 ;YFRAC(20)=0.166667
YFRAC(21)=0.180584 ;YFRAC(22)=0.193998
YFRAC(23)=0.2069 ;YFRAC(24)=0.219282
YFRAC(25)=0.231133 ;YFRAC(26)=0.242444
YFRAC(27)=0.253202 ;YFRAC(28)=0.263397
YFRAC(29)=0.273013 ;YFRAC(30)=0.282036
YFRAC(31)=0.290448 ;YFRAC(32)=0.29823
YFRAC(33)=0.305359 ;YFRAC(34)=0.311808
YFRAC(35)=0.317545 ;YFRAC(36)=0.322529
YFRAC(37)=0.326708 ;YFRAC(38)=0.330008
YFRAC(39)=0.33231 ;YFRAC(40)=0.333333
YFRAC(41)=0.334167 ;YFRAC(42)=0.336667
YFRAC(43)=0.340833 ;YFRAC(44)=0.346667
YFRAC(45)=0.354167 ;YFRAC(46)=0.363333
YFRAC(47)=0.374167 ;YFRAC(48)=0.386667
YFRAC(49)=0.400833 ;YFRAC(50)=0.416667
YFRAC(51)=0.434167 ;YFRAC(52)=0.453333
YFRAC(53)=0.474167 ;YFRAC(54)=0.496667
YFRAC(55)=0.520833 ;YFRAC(56)=0.546667
YFRAC(57)=0.574167 ;YFRAC(58)=0.603333
YFRAC(59)=0.634167 ;YFRAC(60)=0.666667
YFRAC(61)=0.699167 ;YFRAC(62)=0.73
YFRAC(63)=0.759167 ;YFRAC(64)=0.786667
YFRAC(65)=0.8125 ;YFRAC(66)=0.836667
YFRAC(67)=0.859167 ;YFRAC(68)=0.88
YFRAC(69)=0.899167 ;YFRAC(70)=0.916667
YFRAC(71)=0.9325 ;YFRAC(72)=0.946667
YFRAC(73)=0.959167 ;YFRAC(74)=0.97
YFRAC(75)=0.979167 ;YFRAC(76)=0.986667
YFRAC(77)=0.9925 ;YFRAC(78)=0.996667
YFRAC(79)=0.999167 ;YFRAC(80)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST =1.
ZFRAC(1)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(12)=KE
NAME(13)=EP ;NAME(140)=EPKE
NAME(141)=PRPS ;NAME(142)=YPLS
NAME(143)=FTWO ;NAME(144)=FMU
NAME(145)=REYN ;NAME(146)=LTLS
NAME(147)=WDIS ;NAME(148)=LEN1
NAME(149)=ENUT ;NAME(150)=TEM1
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
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(KE,Y,Y,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
SOLUTN(EPKE,Y,N,N,N,N,N)
SOLUTN(PRPS,Y,N,N,N,N,Y)
SOLUTN(YPLS,Y,N,N,N,N,Y)
SOLUTN(FTWO,Y,N,N,N,N,Y)
SOLUTN(FMU,Y,N,N,N,N,Y)
SOLUTN(REYN,Y,N,N,N,N,Y)
SOLUTN(LTLS,Y,Y,Y,N,N,Y)
SOLUTN(WDIS,Y,N,N,N,N,N)
SOLUTN(LEN1,Y,N,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
SOLUTN(TEM1,Y,Y,N,N,N,Y)
VIST = 149
LEN1 = 148
PRPS = 141
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
TERMS(P1,Y,Y,Y,N,Y,Y)
TERMS(U1,Y,Y,Y,Y,Y,Y)
TERMS(V1,Y,Y,Y,Y,Y,Y)
TERMS(KE,N,Y,Y,Y,Y,N)
TERMS(EP,N,Y,Y,Y,Y,N)
TERMS(LTLS,N,N,Y,N,Y,Y)
TERMS(TEM1,N,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWENT = T
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 =1. ;TMP1 =0. ;EL1 = GRND4
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0. ;EL1B =0. ;EL1C =0.
ENUL =1.100667E-05 ;ENUT = GRND3
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 3
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(KE)=1. ;PRNDTL(EP)=1.
PRNDTL(LTLS)=1. ;PRNDTL(TEM1)=1.
PRT(U1)=1. ;PRT(V1)=1.
PRT(KE)=1. ;PRT(EP)=1.314
PRT(LTLS)=1. ;PRT(TEM1)=0.86
CP1 =1000. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.3E-04 ;FIINIT(U1)=1.3
FIINIT(V1)=0. ;FIINIT(KE)=6.76
FIINIT(EP)=31.777433 ;FIINIT(EPKE)=1.0E-10
FIINIT(PRPS)=-1. ;FIINIT(YPLS)=1.0E-10
FIINIT(FTWO)=1.0E-10 ;FIINIT(FMU)=1.0E-10
FIINIT(REYN)=1.0E-10 ;FIINIT(LTLS)=1.0E-10
FIINIT(WDIS)=0.01143 ;FIINIT(LEN1)=1.0E-10
FIINIT(ENUT)=1.0E-10 ;FIINIT(TEM1)=0.5
PATCH(BLOCK ,INIVAL, 1, 14, 1, 40, 1, 1, 1, 1)
INIT(BLOCK ,PRPS,0. ,198. )
INIADD = F
FSWEEP = 1
NAMFI =CHAM
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources
PATCH(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(KESOURCE,KE , GRND4 , GRND4 )
COVAL(KESOURCE,EP , GRND4 , GRND4 )
PATCH(INLET ,WEST , 1, 1, 41, 80, 1, 1, 1, 1)
COVAL(INLET ,P1 , FIXFLU ,13. )
COVAL(INLET ,U1 ,0. ,13. )
COVAL(INLET ,V1 ,0. ,0. )
COVAL(INLET ,KE ,0. ,0.7605 )
COVAL(INLET ,EP ,0. ,31.777433 )
COVAL(INLET ,TEM1,0. ,0. )
PATCH(OUTLET ,EAST , 80, 80, 1, 80, 1, 1, 1, 1)
COVAL(OUTLET ,P1 ,1.0E+05 ,0. )
COVAL(OUTLET ,TEM1,0. , SAME )
PATCH(WALLN ,NWALL , 1, 80, 80, 80, 1, 1, 1, 1)
COVAL(WALLN ,U1 , GRND2 ,0. )
COVAL(WALLN ,KE ,1. ,0. )
PATCH(WALLS ,SWALL , 15, 80, 1, 1, 1, 1, 1, 1)
COVAL(WALLS ,U1 , GRND2 ,0. )
COVAL(WALLS ,KE ,1. ,0. )
PATCH(HEATIN ,SOUTH , 15, 80, 1, 1, 1, 1, 1, 1)
COVAL(HEATIN ,TEM1, FIXFLU ,1625. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1500 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(U1)=2
LITER(V1)=2 ;LITER(KE)=5
LITER(EP)=5 ;LITER(LTLS)=20
LITER(TEM1)=10
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(KE)=1.0E-03
ENDIT(EP)=1.0E-03 ;ENDIT(LTLS)=1.0E-03
ENDIT(TEM1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(U1,FALSDT,7.326923E-03)
RELAX(V1,FALSDT,7.326923E-03)
RELAX(KE,LINRLX,0.3)
RELAX(EP,LINRLX,0.3)
RELAX(EPKE,LINRLX,1.)
RELAX(PRPS,LINRLX,1.)
RELAX(YPLS,LINRLX,1.)
RELAX(FTWO,LINRLX,1.)
RELAX(FMU,LINRLX,1.)
RELAX(REYN,LINRLX,1.)
RELAX(LTLS,FALSDT,1.0E+09)
RELAX(WDIS,LINRLX,1.)
RELAX(LEN1,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
RELAX(TEM1,FALSDT,0.732692)
KELIN = 3
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
VARMAX(YPLS)=1.0E+10 ;VARMIN(YPLS)=-1.0E+10
VARMAX(FTWO)=1.0E+10 ;VARMIN(FTWO)=-1.0E+10
VARMAX(FMU)=1.0E+10 ;VARMIN(FMU)=-1.0E+10
VARMAX(REYN)=1.0E+10 ;VARMIN(REYN)=-1.0E+10
VARMAX(LTLS)=1.0E+10 ;VARMIN(LTLS)=-1.0E+10
VARMAX(WDIS)=1.0E+10 ;VARMIN(WDIS)=-1.0E+10
VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
VARMAX(TEM1)=1.0E+10 ;VARMIN(TEM1)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
GENK = T
PARSOL = F
IENUTA = 3
ISG62 = 1
SPEDAT(SET,OUTPUT,NOFIELD,L,T)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,MATERIAL,198,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=19.76 ;EX(U1)=6.848
EX(V1)=0.2304 ;EX(KE)=1.992
EX(EP)=401.799988 ;EX(EPKE)=214.699997
EX(PRPS)=0.9125 ;EX(YPLS)=0.07523
EX(FTWO)=0.9083 ;EX(FMU)=0.825
EX(REYN)=2129. ;EX(LTLS)=6.673E-04
EX(WDIS)=0.01796 ;EX(LEN1)=3.382E-03
EX(ENUT)=2.477E-03 ;EX(TEM1)=4.132
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
OUTPUT(P1,Y,N,Y,Y,Y,Y)
OUTPUT(U1,Y,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(KE,Y,N,Y,Y,Y,Y)
OUTPUT(EP,Y,N,Y,Y,Y,Y)
OUTPUT(EPKE,Y,N,Y,N,N,N)
OUTPUT(PRPS,Y,N,Y,N,N,N)
OUTPUT(YPLS,Y,N,Y,N,N,N)
OUTPUT(FTWO,Y,N,Y,N,N,N)
OUTPUT(FMU,Y,N,Y,N,N,N)
OUTPUT(REYN,Y,N,Y,N,N,N)
OUTPUT(LTLS,Y,N,Y,Y,Y,Y)
OUTPUT(WDIS,Y,N,Y,N,N,N)
OUTPUT(LEN1,Y,N,Y,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,N,N)
OUTPUT(TEM1,Y,N,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 16 ;IYMON = 38 ;IZMON = 1
NPRMON = 10000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
IPLTF = 1 ;IPLTL = 2000 ;NPLT = -1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP