TALK=T;RUN(1,1)
DISPLAY
The problem considered is the near-field of a subsonic turbulent
round jet issuing into stagnant surroundings at Mach 0.95. The
stagnation enthalpy of the nozzle fluid is equal to that of the
free stream, so that with the assumption of unit Prandtl numbers
the energy equation need not be solved. The k-e model is used
and the calculations are started at the jet discharge plane. The
parabolic marching integration is carried out until about 10
diameters downstream. The calculations are made with 30 radial
grid cells and a forward step size of 5% (DZW1) of the local grid
width. The radial width of the grid increases linearly with
downstream distance so as to accommodate the jet spread.
ENDDIS
PHOTON USE
P
PARPHI
vec x 1 sh
pau;cl
con mach x 1 fi;.1
pau;cl
con enut x 1 fi;.1
pau cl;
con tmp1 x 1 fi;.1
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(2D Subsonic Turbulent Round Jet )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 14
************************************************************
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 = 30
YVLAST =1.
AZYV =1. ;AZRI =0. ;AZAL =0.
YFRAC(1)=0.025 ;YFRAC(3)=0.075
YFRAC(5)=0.125 ;YFRAC(7)=0.175
YFRAC(9)=0.225 ;YFRAC(11)=0.275
YFRAC(13)=0.325 ;YFRAC(15)=0.375
YFRAC(17)=0.425 ;YFRAC(19)=0.475
YFRAC(21)=0.525059 ;YFRAC(23)=0.604527
YFRAC(25)=0.703063 ;YFRAC(27)=0.814483
YFRAC(29)=0.935999
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = T
NZ = 140
ZWADD =7.843137
ZWLAST =1.
AZDZ = GRND2
ZFRAC(1)=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(145)=EPKE ;NAME(146)=TMP1
NAME(147)=LEN1 ;NAME(148)=ENUT
NAME(149)=MACH ;NAME(150)=RHO1
* 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,N,N,N,Y)
SOLUTN(V1,Y,Y,N,N,N,N)
SOLUTN(W1,Y,Y,N,N,N,N)
SOLUTN(KE,Y,Y,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
SOLUTN(H1,Y,N,N,N,N,Y)
SOLUTN(EPKE,Y,N,N,N,N,N)
SOLUTN(TMP1,Y,N,N,N,N,Y)
SOLUTN(LEN1,Y,N,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
SOLUTN(MACH,Y,N,N,N,N,Y)
SOLUTN(RHO1,Y,N,N,N,N,Y)
DEN1 = 150
VIST = 148
TEMP1 = 146
LEN1 = 147
************************************************************
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(KE,N,Y,Y,Y,Y,N)
TERMS(EP,N,Y,Y,Y,Y,N)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWRH1 = T
NEWENT = T
U1AD =0. ;V1AD = GRND1 ;W1AD =0.
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
DENPCO = T
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 = GRND5 ;TMP1 = GRND6 ;EL1 = GRND4
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP = GRND5
RHO1A =0. ;RHO1B =1. ;RHO1C =0.714286
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 =3.042555E-05 ;ENUT = GRND3
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PRNDTL(V1)=1. ;PRNDTL(W1)=1.
PRNDTL(KE)=1. ;PRNDTL(EP)=1.
PRT(V1)=1. ;PRT(W1)=1.
PRT(KE)=1. ;PRT(EP)=1.314
CP1 =3.5 ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=0.55946 ;FIINIT(V1)=1.0E-10
FIINIT(W1)=1.034557 ;FIINIT(KE)=0.010703
FIINIT(EP)=3.638576E-03 ;FIINIT(H1)=3.5
FIINIT(EPKE)=1.0E-10 ;FIINIT(TMP1)=1.0E-10
FIINIT(LEN1)=1.0E-10 ;FIINIT(ENUT)=1.0E-10
FIINIT(MACH)=1.0E-10 ;FIINIT(RHO1)=0.660442
PATCH(INITFS ,INIVAL, 1, 1, 21, 30, 1, 1, 1, 1)
INIT(INITFS ,W1 ,0. ,0. )
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(IN ,LOW , 1, 1, 1, 20, 1, 1, 1, 1)
COVAL(IN ,P1 , FIXFLU ,0.683266 )
COVAL(IN ,V1 ,0. ,0. )
COVAL(IN ,W1 ,0. ,1.034557 )
COVAL(IN ,KE ,0. ,0.010703 )
COVAL(IN ,EP ,0. ,3.638576E-03 )
PATCH(NB ,NORTH , 1, 1, 30, 30, 1, 140, 1, 1)
COVAL(NB ,P1 ,1000. ,0.55946 )
COVAL(NB ,W1 ,0. ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
IPARAB = 1
AZPH =0. ;PBAR =0.
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1 ;ISWC1 = 1
LITHYD = 60 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = F
RESREF(P1)=8.540821E-15 ;RESREF(V1)=8.835971E-15
RESREF(W1)=8.835971E-15 ;RESREF(KE)=9.141319E-17
RESREF(EP)=3.107642E-17
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(V1)=10
LITER(W1)=10 ;LITER(KE)=20
LITER(EP)=20
ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
ENDIT(W1)=1.0E-03 ;ENDIT(KE)=1.0E-03
ENDIT(EP)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,0.15)
RELAX(V1,FALSDT,0.483298)
RELAX(W1,FALSDT,0.483298)
RELAX(KE,LINRLX,0.2)
RELAX(EP,LINRLX,0.2)
RELAX(H1,LINRLX,1.)
RELAX(EPKE,LINRLX,1.)
RELAX(TMP1,LINRLX,1.)
RELAX(LEN1,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
RELAX(MACH,LINRLX,1.)
RELAX(RHO1,LINRLX,0.05)
KELIN = 3
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=5. ;VARMIN(P1)=-1.0E+10
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=1.0E-10
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
VARMAX(H1)=3.5 ;VARMIN(H1)=1.0E-10
VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=-1.0E+10
VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
VARMAX(MACH)=1.0E+10 ;VARMIN(MACH)=-1.0E+10
VARMAX(RHO1)=5. ;VARMIN(RHO1)=1.0E-04
************************************************************
Group 19. Data transmitted to GROUND
GENK = T
PARSOL = F
DZW1 =0.05
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,Y,Y,Y)
OUTPUT(EP,Y,N,Y,Y,Y,Y)
OUTPUT(H1,Y,N,Y,N,N,N)
OUTPUT(EPKE,Y,N,Y,N,N,N)
OUTPUT(TMP1,Y,N,Y,N,N,N)
OUTPUT(LEN1,Y,N,Y,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,N,N)
OUTPUT(MACH,Y,N,Y,N,N,N)
OUTPUT(RHO1,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 20 ;IZMON = 1
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NYPRIN = 2 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = 28 ;IZPRF = 1 ;IZPRL = 10000
IPLTF = 1 ;IPLTL = 60 ;NPLT = 5
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 2 ;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
IDISPA = 2 ;IDISPB = 1 ;IDISPC = 140
STOP