TALK=T;RUN(1,1)
DISPLAY
The problem considered is the near-field of a sonic under-
expanded turbulent round jet issuing into stagnant surroundings.
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 turbulence
is represented by the k-e model, and the IPARAB=5 option of the
parabolic solver is used to perform the marching integration.
The calculations are started at the jet discharge plane, and
carried out until 6 diameters downstream. The calculations are
made with 35 radial grid cells and a forward step size of 1%
(DZW1) of the local grid width. The main input parameters are
the inlet and free-stream Mach numbers and the jet-to-ambient
static pressure ratio. This case has been tested for inlet
Mach numbers up to 2 and pressure ratios up to 2.5.
ENDDIS
PHOTON USE
P
PARPHI
VEC X 1 SH
pau;cl
con mach x 1 fi;.5
pau;cl
con p1 x 1 fi;.5
pau cl;
con tmp1 x 1 fi;.5
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(2D Sonic Underexpanded 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 = 35
YVLAST =1.5
YFRAC(1)=0.016667 ;YFRAC(3)=0.05
YFRAC(5)=0.083333 ;YFRAC(7)=0.116667
YFRAC(9)=0.15 ;YFRAC(11)=0.183333
YFRAC(13)=0.216667 ;YFRAC(15)=0.25
YFRAC(17)=0.283333 ;YFRAC(19)=0.316667
YFRAC(21)=0.336296 ;YFRAC(23)=0.36
YFRAC(25)=0.407407 ;YFRAC(27)=0.478519
YFRAC(29)=0.573333 ;YFRAC(31)=0.691852
YFRAC(33)=0.834074 ;YFRAC(35)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = T
NZ = 400
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(144)=EPKE ;NAME(145)=MACZ
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,Y,Y,N,Y)
SOLUTN(W1,Y,Y,Y,Y,N,Y)
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(MACZ,Y,N,N,N,N,Y)
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,N,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. ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWRH1 = T
NEWENT = T
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(W1)=1. ;PRNDTL(KE)=1.
PRNDTL(EP)=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.528282 ;FIINIT(V1)=1.0E-10
FIINIT(W1)=1.080123 ;FIINIT(KE)=2.916666E-03
FIINIT(EP)=5.176041E-04 ;FIINIT(H1)=3.5
FIINIT(EPKE)=1.0E-10 ;FIINIT(MACZ)=1.
FIINIT(TMP1)=0.833333 ;FIINIT(LEN1)=1.0E-10
FIINIT(ENUT)=1.0E-10 ;FIINIT(MACH)=1.0E-10
FIINIT(RHO1)=0.633938
PATCH(INITFS ,INIVAL, 1, 1, 21, 35, 1, 1, 1, 1)
INIT(INITFS ,P1 ,0. ,0.285558 )
INIT(INITFS ,W1 ,0. ,0. )
INIT(INITFS ,MACZ,0. ,1.0E-03 )
INIT(INITFS ,TMP1,0. ,1. )
INIT(INITFS ,RHO1,0. ,0.285558 )
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.684731 )
COVAL(IN ,V1 ,0. ,0. )
COVAL(IN ,W1 ,0. ,1.080123 )
COVAL(IN ,KE ,0. ,2.916666E-03 )
COVAL(IN ,EP ,0. ,5.176041E-04 )
PATCH(INFS ,LOW , 1, 1, 21, 35, 1, 1, 1, 1)
COVAL(INFS ,P1 , FIXFLU ,3.378764E-04 )
COVAL(INFS ,V1 ,0. ,0. )
COVAL(INFS ,W1 ,0. ,1.183216E-03 )
COVAL(INFS ,KE ,0. ,0. )
COVAL(INFS ,EP ,0. ,0. )
PATCH(NB ,NORTH , 1, 1, 35, 35, 1, 400, 1, 1)
COVAL(NB ,P1 ,1000. ,0.285558 )
COVAL(NB ,W1 ,0. ,1.183216E-03 )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
IPARAB = 5
AZPH =0. ;PBAR =0.
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1 ;ISWC1 = 1
LITHYD = 20 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = F
RESREF(P1)=8.559143E-15 ;RESREF(V1)=9.24493E-15
RESREF(W1)=9.24493E-15 ;RESREF(KE)=2.496416E-17
RESREF(EP)=4.430247E-18
************************************************************
Group 16. Terminate Iterations
LITER(P1)=100 ;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.7)
RELAX(V1,LINRLX,0.5)
RELAX(W1,LINRLX,0.5)
RELAX(KE,LINRLX,0.3)
RELAX(EP,LINRLX,0.3)
RELAX(H1,LINRLX,1.)
RELAX(EPKE,LINRLX,1.)
RELAX(MACZ,LINRLX,0.8)
RELAX(TMP1,LINRLX,1.)
RELAX(LEN1,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
RELAX(MACH,LINRLX,1.)
RELAX(RHO1,LINRLX,1.)
KELIN = 3
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=1.0E-04
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(MACZ)=1.0E+10 ;VARMIN(MACZ)=-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-03
************************************************************
Group 19. Data transmitted to GROUND
GENK = T
PARSOL = F
DZW1 =1.0E-02
ISG62 = 1
RMACHZ =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(MACZ,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,Y,N)
OUTPUT(RHO1,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 1 ;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 = 80 ;IZPRF = 1 ;IZPRL = 10000
IPLTF = 1 ;IPLTL = 20 ;NPLT = 2
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 = 400
STOP