```
TALK=T;RUN(1,1)
DISPLAY
The problem considered is an axisymmetric sonic jet discharging
into stagnant surroundings from a nozzle at a pressure 3.56
times higher than the ambient pressure. 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 means of the k-e
turbulence model. The calculation may be made with the van-Leer
higher-order scheme or the hybrid differencing scheme. A
relatively coarse mesh of 30 radial by 90 axial cells is used in
the computations.
ENDDIS
PHOTON USE
P

0.20443E+04 0.15633E+04 CR
CON MACH X 1 FI;.5
msg Mach number contours
msg Press  to continue
pause
ENDUSE
AUTOPLOT USE
FILE
PHI 5

D 1 MACH Y 1;PLOT;REDR;LEVEL Y 1;level x 1.54
msg Mach number distribution along flow axis
msg Press  to continue
pause
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************

TEXT(2D TRANSONIC UNDEREXPANDED JET: N113    )

************************************************************
************************************************************

IRUNN = 1 ;LIBREF = 14
************************************************************
Group 2. Time dependence
************************************************************
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.
YFRAC(1)=0.025 ;YFRAC(2)=0.05
YFRAC(3)=0.075 ;YFRAC(4)=0.1
YFRAC(5)=0.125 ;YFRAC(6)=0.15
YFRAC(7)=0.175 ;YFRAC(8)=0.2
YFRAC(9)=0.225 ;YFRAC(10)=0.25
YFRAC(11)=0.275 ;YFRAC(12)=0.3
YFRAC(13)=0.325 ;YFRAC(14)=0.35
YFRAC(15)=0.375 ;YFRAC(16)=0.4
YFRAC(17)=0.425 ;YFRAC(18)=0.45
YFRAC(19)=0.475 ;YFRAC(20)=0.5
YFRAC(21)=0.525059 ;YFRAC(22)=0.561703
YFRAC(23)=0.604527 ;YFRAC(24)=0.651932
YFRAC(25)=0.703063 ;YFRAC(26)=0.757375
YFRAC(27)=0.814483 ;YFRAC(28)=0.874099
YFRAC(29)=0.935999 ;YFRAC(30)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 90
ZWLAST =5.
ZFRAC(1)=7.692308E-03 ;ZFRAC(2)=0.015385
ZFRAC(3)=0.023077 ;ZFRAC(4)=0.030769
ZFRAC(5)=0.038462 ;ZFRAC(6)=0.046154
ZFRAC(7)=0.053846 ;ZFRAC(8)=0.061538
ZFRAC(9)=0.069231 ;ZFRAC(10)=0.076923
ZFRAC(11)=0.084615 ;ZFRAC(12)=0.092308
ZFRAC(13)=0.1 ;ZFRAC(14)=0.107692
ZFRAC(15)=0.115385 ;ZFRAC(16)=0.123077
ZFRAC(17)=0.130769 ;ZFRAC(18)=0.138462
ZFRAC(19)=0.146154 ;ZFRAC(20)=0.153846
ZFRAC(21)=0.161538 ;ZFRAC(22)=0.169231
ZFRAC(23)=0.176923 ;ZFRAC(24)=0.184615
ZFRAC(25)=0.192308 ;ZFRAC(26)=0.2
ZFRAC(27)=0.207692 ;ZFRAC(28)=0.215385
ZFRAC(29)=0.223077 ;ZFRAC(30)=0.230769
ZFRAC(31)=0.238462 ;ZFRAC(32)=0.246154
ZFRAC(33)=0.253846 ;ZFRAC(34)=0.261538
ZFRAC(35)=0.269231 ;ZFRAC(36)=0.276923
ZFRAC(37)=0.284615 ;ZFRAC(38)=0.292308
ZFRAC(39)=0.3 ;ZFRAC(40)=0.307692
ZFRAC(41)=0.315385 ;ZFRAC(42)=0.323077
ZFRAC(43)=0.330769 ;ZFRAC(44)=0.338462
ZFRAC(45)=0.346154 ;ZFRAC(46)=0.353846
ZFRAC(47)=0.361538 ;ZFRAC(48)=0.369231
ZFRAC(49)=0.376923 ;ZFRAC(50)=0.384615
ZFRAC(51)=0.392308 ;ZFRAC(52)=0.4
ZFRAC(53)=0.407692 ;ZFRAC(54)=0.415385
ZFRAC(55)=0.423077 ;ZFRAC(56)=0.430769
ZFRAC(57)=0.438462 ;ZFRAC(58)=0.446154
ZFRAC(59)=0.453846 ;ZFRAC(60)=0.461538
ZFRAC(61)=0.469231 ;ZFRAC(62)=0.476923
ZFRAC(63)=0.484615 ;ZFRAC(64)=0.492308
ZFRAC(65)=0.5 ;ZFRAC(66)=0.511636
ZFRAC(67)=0.526733 ;ZFRAC(68)=0.543487
ZFRAC(69)=0.561416 ;ZFRAC(70)=0.580274
ZFRAC(71)=0.599906 ;ZFRAC(72)=0.620207
ZFRAC(73)=0.641098 ;ZFRAC(74)=0.662519
ZFRAC(75)=0.684422 ;ZFRAC(76)=0.706768
ZFRAC(77)=0.729525 ;ZFRAC(78)=0.752664
ZFRAC(79)=0.776163 ;ZFRAC(80)=0.8
ZFRAC(81)=0.812619 ;ZFRAC(82)=0.828991
ZFRAC(83)=0.84716 ;ZFRAC(84)=0.866604
ZFRAC(85)=0.887055 ;ZFRAC(86)=0.908346
ZFRAC(87)=0.930361 ;ZFRAC(88)=0.953016
ZFRAC(89)=0.976247 ;ZFRAC(90)=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(146)=TMP1 ;NAME(147)=VABS
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,Y,N,N,N)
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(TMP1,Y,N,N,N,N,Y)
SOLUTN(VABS,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
************************************************************
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. ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWRH1 = T
NEWENT = T
UCONV = 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.
EL1A =0. ;EL1B =0.
EL1C =0.
ENUL =2.0E-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.338275 ;FIINIT(V1)=0.
FIINIT(W1)=1.080123 ;FIINIT(KE)=2.916666E-03
FIINIT(EP)=5.176041E-04 ;FIINIT(H1)=3.5
FIINIT(TMP1)=1.0E-10 ;FIINIT(VABS)=1.0E-10
FIINIT(ENUT)=1.0E-10 ;FIINIT(MACH)=1.0E-10
FIINIT(RHO1)=0.316969
No PATCHes yet used for this Group
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(HOCS ,CELL , 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(HOCS ,V1 , FIXFLU , GRND2 )
COVAL(HOCS ,W1 , FIXFLU , GRND2 )
COVAL(HOCS ,KE , FIXFLU , GRND2 )
COVAL(HOCS ,EP , FIXFLU , GRND2 )

PATCH(IN ,LOW , 1, 1, 1, 15, 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(LB ,LOW , 1, 1, 16, 30, 1, 1, 1, 1)
COVAL(LB ,P1 ,1000. ,0.148269 )
COVAL(LB ,KE ,0. ,1.0E-10 )
COVAL(LB ,EP ,0. ,4.5E-15 )

PATCH(NB ,NORTH , 1, 1, 30, 30, 1, 90, 1, 1)
COVAL(NB ,P1 ,1000. ,0.148269 )
COVAL(NB ,KE ,0. ,1.0E-10 )
COVAL(NB ,EP ,0. ,4.5E-15 )

PATCH(OUT ,HIGH , 1, 1, 1, 30, 90, 90, 1, 1)
COVAL(OUT ,P1 ,1000. ,0.148269 )
COVAL(OUT ,KE ,0. , SAME )
COVAL(OUT ,EP ,0. , SAME )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1600 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
LITER(P1)=10 ;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,FALSDT,0.012859)
RELAX(W1,FALSDT,0.012859)
RELAX(KE,FALSDT,0.012859)
RELAX(EP,FALSDT,0.012859)
RELAX(H1,LINRLX,1.)
RELAX(TMP1,LINRLX,1.)
RELAX(VABS,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
RELAX(MACH,LINRLX,1.)
RELAX(RHO1,LINRLX,1.)
KELIN = 1
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=1.0E-10
VARMAX(V1)=100. ;VARMIN(V1)=-100.
VARMAX(W1)=100. ;VARMIN(W1)=-100.
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(TMP1)=1.0E+10 ;VARMIN(TMP1)=-1.0E+10
VARMAX(VABS)=1.0E+10 ;VARMIN(VABS)=-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)=1.0E+10 ;VARMIN(RHO1)=1.0E-10
************************************************************
Group 19. Data transmitted to GROUND
NAMGRD =CONV
GENK = T
PARSOL = F
ISG62 = 1
SPEDAT(SET,SCHEME,INLCS005,I,8)
SPEDAT(SET,SCHEME,INLCS007,I,8)
SPEDAT(SET,SCHEME,INLCS012,I,8)
SPEDAT(SET,SCHEME,INLCS013,I,8)
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(TMP1,Y,N,Y,N,N,N)
OUTPUT(VABS,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,Y,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 1 ;IZMON = 78
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 1600 ;NUMCLS = 5
NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 5
NZPRIN = 1 ;IZPRF = 43 ;IZPRL = 55
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = -1 ;NPLT = 20
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20

PATCH(PLOT4 ,PROFIL, 1, 1, 1, 1, 1, 90, 1, 1)
PLOT(PLOT4 ,W1 ,0. ,0. )

PATCH(PLOT5 ,PROFIL, 1, 1, 1, 1, 1, 90, 1, 1)
PLOT(PLOT5 ,MACH,0. ,0. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP
```