TALK=T;RUN(1,1)
DISPLAY
The case considered is the normal impingement of a turbulent
axisymmetric jet on a heated flat plate. The flow geometry
and conditions approximate to those reported experimentally
by D.Cooper et al ( Int.J.Heat Mass Transfer, Vol.36, No.10,
p2675,[1993] ). The jet issues into ambient fluid at a
Reynolds number of 2.3E4 from a pipe located 2 pipe diameters
above the electrically heated plate. The jet injection
temperature is ambient so that the purpose of the flow process
is to effect cooling. For the case of orthogonal impaction,
the flow is axisymmetric and is directed outward along the
surface giving rise to a radial wall jet. The simulation may
be performed with the standard k-e model or with the Reynolds
stress transport model. For the latter, the model employs
differential transport equations for determining the turbulent
fluxes of heat.
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(RSTM_2D IMPINGING ROUND JET :T609)
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 609
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
NX = 1
XULAST =0.1
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 25
YVLAST =5.
YFRAC(1)=1.0E-02 ;YFRAC(2)=0.02
YFRAC(3)=0.03 ;YFRAC(4)=0.04
YFRAC(5)=0.05 ;YFRAC(6)=0.06
YFRAC(7)=0.07 ;YFRAC(8)=0.08
YFRAC(9)=0.09 ;YFRAC(10)=0.1
YFRAC(11)=0.115492 ;YFRAC(12)=0.143818
YFRAC(13)=0.180498 ;YFRAC(14)=0.223935
YFRAC(15)=0.273205 ;YFRAC(16)=0.327684
YFRAC(17)=0.386915 ;YFRAC(18)=0.450542
YFRAC(19)=0.518282 ;YFRAC(20)=0.589898
YFRAC(21)=0.66519 ;YFRAC(22)=0.743988
YFRAC(23)=0.826141 ;YFRAC(24)=0.911517
YFRAC(25)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 30
ZWLAST =2.
ZFRAC(1)=0.065812 ;ZFRAC(2)=0.114585
ZFRAC(3)=0.158489 ;ZFRAC(4)=0.199504
ZFRAC(5)=0.238495 ;ZFRAC(6)=0.275946
ZFRAC(7)=0.312163 ;ZFRAC(8)=0.347356
ZFRAC(9)=0.381678 ;ZFRAC(10)=0.415244
ZFRAC(11)=0.448144 ;ZFRAC(12)=0.48045
ZFRAC(13)=0.512221 ;ZFRAC(14)=0.543507
ZFRAC(15)=0.574349 ;ZFRAC(16)=0.604782
ZFRAC(17)=0.634837 ;ZFRAC(18)=0.66454
ZFRAC(19)=0.693914 ;ZFRAC(20)=0.722981
ZFRAC(21)=0.751759 ;ZFRAC(22)=0.780263
ZFRAC(23)=0.80851 ;ZFRAC(24)=0.836512
ZFRAC(25)=0.864281 ;ZFRAC(26)=0.891829
ZFRAC(27)=0.919166 ;ZFRAC(28)=0.946301
ZFRAC(29)=0.973243 ;ZFRAC(30)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(5)=V1
NAME(7)=W1 ;NAME(12)=KE
NAME(13)=EP ;NAME(14)=H1
NAME(142)=TMP1 ;NAME(143)=WTRS
NAME(144)=VTRS ;NAME(145)=U2RS
NAME(146)=V2RS ;NAME(147)=W2RS
NAME(148)=VWRS ;NAME(149)=LEN1
NAME(150)=ENUT
* 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(V1,Y,Y,N,N,N,Y)
SOLUTN(W1,Y,Y,N,N,N,Y)
SOLUTN(KE,Y,N,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
SOLUTN(H1,Y,Y,Y,N,N,Y)
SOLUTN(TMP1,Y,N,N,N,N,Y)
SOLUTN(WTRS,Y,Y,N,N,N,N)
SOLUTN(VTRS,Y,Y,N,N,N,N)
SOLUTN(U2RS,Y,Y,N,N,N,N)
SOLUTN(V2RS,Y,Y,N,N,N,N)
SOLUTN(W2RS,Y,Y,N,N,N,N)
SOLUTN(VWRS,Y,Y,N,N,N,N)
SOLUTN(LEN1,Y,N,N,N,N,N)
SOLUTN(ENUT,Y,N,N,N,N,N)
VIST = 150
TEMP1 = 142
LEN1 = 149
************************************************************
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(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(EP,Y,Y,Y,Y,N,Y)
TERMS(H1,Y,Y,Y,Y,Y,Y)
TERMS(WTRS,N,Y,Y,Y,N,Y)
TERMS(VTRS,N,Y,Y,Y,Y,Y)
TERMS(U2RS,N,Y,Y,Y,N,Y)
TERMS(V2RS,N,Y,Y,Y,Y,Y)
TERMS(W2RS,N,Y,Y,Y,N,Y)
TERMS(VWRS,N,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWENT = F
UDIFF = T
UCONNE = T
USOURC = 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 = GRND2 ;EL1 = GRND
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
TMP1A =0. ;TMP1B =0. ;TMP1C =0.
TMP2A =0. ;TMP2B =0. ;TMP2C =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0. ;EL1B =0. ;EL1C =0.
ENUL =1.408451E-05 ;ENUT = GRND
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PRNDTL(V1)=1. ;PRNDTL(W1)=1.
PRNDTL(EP)=1. ;PRNDTL(H1)=0.71
PRNDTL(WTRS)=1. ;PRNDTL(VTRS)=1.
PRNDTL(U2RS)=1. ;PRNDTL(V2RS)=1.
PRNDTL(W2RS)=1. ;PRNDTL(VWRS)=1.
PRT(V1)=1. ;PRT(W1)=1.
PRT(EP)=1.314 ;PRT(H1)=0.9
PRT(WTRS)=1. ;PRT(VTRS)=1.
PRT(U2RS)=1. ;PRT(V2RS)=1.
PRT(W2RS)=1. ;PRT(VWRS)=1.
CP1 =1000. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(V1)=0.
FIINIT(W1)=1.0E-10 ;FIINIT(KE)=1.0E-02
FIINIT(EP)=0.0639 ;FIINIT(H1)=3.0E+05
FIINIT(TMP1)=1.0E-10 ;FIINIT(WTRS)=1.0E-10
FIINIT(VTRS)=1.0E-10 ;FIINIT(U2RS)=6.666666E-03
FIINIT(V2RS)=6.666666E-03 ;FIINIT(W2RS)=6.666666E-03
FIINIT(VWRS)=3.0E-03 ;FIINIT(LEN1)=1.0E-10
FIINIT(ENUT)=1.0E-10
PATCH(INWJET ,INIVAL, 1, 1, 1, 10, 1, 30, 1, 1)
INIT(INWJET ,W1 ,0. ,1. )
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(PLATE ,HWALL , 1, 1, 1, 25, 30, 30, 1, 1)
COVAL(PLATE ,V1 , GRND2 ,0. )
COVAL(PLATE ,EP , GRND2 , GRND2 )
COVAL(PLATE ,H1 , FIXFLU ,1.0E+04 )
PATCH(WDMPH03 ,PHASEM, 1, 1, 1, 25, 1, 30, 1, 1)
COVAL(WDMPH03 ,WTRS, FIXFLU , GRND1 )
COVAL(WDMPH03 ,VTRS, FIXFLU , GRND1 )
COVAL(WDMPH03 ,U2RS, FIXFLU , GRND1 )
COVAL(WDMPH03 ,V2RS, FIXFLU , GRND1 )
COVAL(WDMPH03 ,W2RS, FIXFLU , GRND1 )
COVAL(WDMPH03 ,VWRS, FIXFLU , GRND1 )
PATCH(JET1 ,LOW , 1, 1, 1, 10, 1, 1, 1, 1)
COVAL(JET1 ,P1 , FIXFLU ,1. )
COVAL(JET1 ,V1 ,0. ,0. )
COVAL(JET1 ,W1 ,0. ,1. )
COVAL(JET1 ,EP ,0. ,1.825555E-03 )
COVAL(JET1 ,H1 ,0. ,3.0E+05 )
COVAL(JET1 ,WTRS,0. ,0. )
COVAL(JET1 ,VTRS,0. ,0. )
COVAL(JET1 ,U2RS,0. ,6.666666E-03 )
COVAL(JET1 ,V2RS,0. ,6.666666E-03 )
COVAL(JET1 ,W2RS,0. ,6.666666E-03 )
COVAL(JET1 ,VWRS,0. ,0. )
PATCH(TOP ,LOW , 1, 1, 11, 25, 1, 1, 1, 1)
COVAL(TOP ,P1 ,1000. ,0. )
COVAL(TOP ,V1 ,0. ,0. )
COVAL(TOP ,W1 ,0. ,0. )
COVAL(TOP ,KE ,0. ,1.0E-10 )
COVAL(TOP ,EP ,0. ,1.0E-10 )
COVAL(TOP ,H1 ,0. ,3.0E+05 )
PATCH(SIDE ,NORTH , 1, 1, 25, 25, 1, 30, 1, 1)
COVAL(SIDE ,P1 ,1000. ,0. )
COVAL(SIDE ,V1 ,0. ,0. )
COVAL(SIDE ,W1 ,0. ,0. )
COVAL(SIDE ,KE ,0. ,1.0E-10 )
COVAL(SIDE ,EP ,0. ,1.0E-10 )
COVAL(SIDE ,H1 ,0. ,3.0E+05 )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1000 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(V1)=10
LITER(W1)=10 ;LITER(EP)=20
LITER(H1)=30 ;LITER(WTRS)=20
LITER(VTRS)=20 ;LITER(U2RS)=20
LITER(V2RS)=20 ;LITER(W2RS)=20
LITER(VWRS)=20
ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
ENDIT(W1)=1.0E-03 ;ENDIT(EP)=1.0E-03
ENDIT(H1)=1.0E-03 ;ENDIT(WTRS)=1.0E-03
ENDIT(VTRS)=1.0E-03 ;ENDIT(U2RS)=1.0E-03
ENDIT(V2RS)=1.0E-03 ;ENDIT(W2RS)=1.0E-03
ENDIT(VWRS)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(V1,FALSDT,0.266667)
RELAX(W1,FALSDT,0.266667)
RELAX(KE,LINRLX,1.)
RELAX(EP,FALSDT,0.066667)
RELAX(H1,FALSDT,0.666667)
RELAX(TMP1,LINRLX,1.)
RELAX(WTRS,FALSDT,0.016667)
RELAX(VTRS,FALSDT,0.016667)
RELAX(U2RS,FALSDT,0.066667)
RELAX(V2RS,FALSDT,0.066667)
RELAX(W2RS,FALSDT,0.066667)
RELAX(VWRS,FALSDT,0.066667)
RELAX(LEN1,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=-1.0E+10
VARMAX(WTRS)=1.0E+10 ;VARMIN(WTRS)=-1.0E+10
VARMAX(VTRS)=1.0E+10 ;VARMIN(VTRS)=-1.0E+10
VARMAX(U2RS)=1.0E+10 ;VARMIN(U2RS)=-1.0E+10
VARMAX(V2RS)=1.0E+10 ;VARMIN(V2RS)=-1.0E+10
VARMAX(W2RS)=1.0E+10 ;VARMIN(W2RS)=-1.0E+10
VARMAX(VWRS)=1.0E+10 ;VARMIN(VWRS)=-1.0E+10
VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
GENK = T
RSTM = T
PARSOL = F
IRSMSM = 2
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
************************************************************
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(V1,Y,N,Y,Y,Y,Y)
OUTPUT(W1,Y,N,Y,Y,Y,Y)
OUTPUT(KE,Y,N,Y,N,N,N)
OUTPUT(EP,Y,N,Y,Y,Y,Y)
OUTPUT(H1,Y,N,Y,Y,Y,Y)
OUTPUT(TMP1,Y,N,Y,N,N,N)
OUTPUT(WTRS,Y,N,Y,Y,Y,Y)
OUTPUT(VTRS,Y,N,Y,Y,Y,Y)
OUTPUT(U2RS,Y,N,Y,Y,Y,Y)
OUTPUT(V2RS,Y,N,Y,Y,Y,Y)
OUTPUT(W2RS,Y,N,Y,Y,Y,Y)
OUTPUT(VWRS,Y,N,Y,Y,Y,Y)
OUTPUT(LEN1,Y,N,Y,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,N,N)
WALPRN = T
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 12 ;IZMON = 28
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 25
NZPRIN = 1 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = -1 ;NPLT = 5
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